Engineering PapersSearch

SEARCH · Engineering Papers

Results for “GLOBAL COARSE GRID”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

A truncation error injection approach to viscous-inviscid interaction

An approach to viscous-inviscid interaction which is based on truncation error injection is presented in the context of solving flow over an airfoil. A two-dimensional interpolation scheme is used to restrict the fine grid solutions to the global coarse grid. Details on the current implementation of the approach are given, and the boundary conditions being used are discussed. Inviscid results from a NACA0012 airfoil test case and the viscous results are presented.

Goble, B. D.

Three-dimensional elliptic grid generation for an F-16

A case history depicting the effort to generate a computational grid for the simulation of transonic flow about an F-16 aircraft at realistic flight conditions is presented. The flow solver for which this grid is designed is a zonal one, using the Reynolds averaged Navier-Stokes equations near the surface of the aircraft, and the Euler equations in regions removed from the aircraft. A body conforming global grid, suitable for the Euler equation, is first generated using 3-D Poisson equations having inhomogeneous terms modeled after the 2-D GRAPE code. Regions of the global grid are then designated for zonal refinement as appropriate to accurately model the flow physics. Grid spacing suitable for solution of the Navier-Stokes equations is generated in the refinement zones by simple subdivision of the given coarse grid intervals. That grid generation project is described, with particular emphasis on the global coarse grid.

Reese L. Sorenson

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on a multiple-instruction/multiple-data stream parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed: one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All of the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the Intel iPSC/860.

TRT-THEORETICAL

Three-dimensional elliptic grid generation about fighter aircraft for zonal finite-difference computations

An elliptic grid-generation method for finite-difference computations about complex aerodynamic configurations is developed. A zonal approach is used, which involves first making a coarse global grid filling the entire physical domain and then subdividing regions of that grid to make the individual zone grids. The details of the grid-generation method are presented along with results of the present application, a wing-body configuration based on the F-16 fighter aircraft.

Sorenson, R. L.

Viscous-inviscid interaction and local grid refinement via truncation error injection

A methodology is presented which makes it possible to decouple a complex problem having multiple disparate length scales into problems of single length scale so that they can be solved more efficiently on a computer. The method is applied to a viscous transonic flow over an airfoil. It is found that accurate prediction of the flow over an airfoil can be obtained by solving the Euler equations on a relatively coarse global grid with viscous effects computed separately on a boundary-layer type grid and injected into the global grid solution as a combination of vorticity and trucation error.

Goble, Brian D.

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on an MIMD (multiple instruction/multiple data stream) parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed, one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the iPSC/860.

Venkatakrishnan, V.

Advances in the computation of transonic separated flows over finite wings

The transonic flow around a low-aspect-ratio infinite wing in a wind tunnel is investigated by means of numerical simulations, with a focus on shock-induced separated flows. A coarse global grid with far-field boundaries matching those of the test section is subdivided into zones: the flow in the clustered zones near the wing is analyzed by solving the Reynolds-averaged Navier-Stokes equations, while that farther from the wing is modeled with the Euler equations. The results are presented graphically, and it is shown that a mushroomlike separated flow with two counterrotating vortices can be simulated when the correct shock strength is imposed (by careful selection of the artificial dissipation, the boundary conditions, the grid refinement, the algebraic turbulence model, and the geometry representation).

Kaynak, Unver

3-D multilateration for measurement of earth crustal deformation and network densification

Discussion of how range and range-difference data types can make possible precise three-dimensional measurement of ground station positions and the position of an artificial signal source, without explicit dependence upon the signal source trajectory. An effective strategy for such measurement is to combine the multilateration approach with a VLBI system using natural radio sources. The VLBI methods would provide a coarse grid of three-dimensional benchmark locations on a regional and global scale. Multilateration stations would then occupy these coarse grid locations and provide a means for highly portable, relatively low-cost units to then densify networks on a regional and local scale. Because a multilateration approach can make use of strong artificial radio sources, it makes possible the use of relatively low-cost, highly mobile stations. Such mobile stations are virtually essential for three dimensional surveying in heavily urbanized areas or in rugged terrain.

Ong, K. M.

Development of a Navier-Stokes algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.

Development of an aerodynamics algorithm for parallel-processing supercomputers

An explicit flow solver, applicable to the hierarchy of model equations ranging from Euler to full Navier-Stokes, is combined with several techniques designed to reduce computational expense. The computational domain consists of local grid refinements embedded in a global coarse mesh, where the locations of these refinements are defined by the physics of the flow. Flow characteristics are also used to determine which set of model equations is appropriate for solution in each region, thereby reducing not only the number of grid points at which the solution must be obtained, but also the computational effort required to get that solution. Acceleration to steady-state is achieved by applying multigrid on each of the subgrids, regardless of the particular model equations being solved. Since each of these components is explicit, advantage can readily be taken of the vector- and parallel-processing capabilities of machines such as the Cray X-MP and Cray-2.

Swisshelm, Julie M.

Numerical simulation of shock wave diffraction by TVD schemes

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

Young, Victor Y. C.

Zonal multigrid solution of compressible flow problems on unstructured and adaptive meshes

The simultaneous use of adaptive meshing techniques with a multigrid strategy for solving the 2-D Euler equations in the context of unstructured meshes is studied. To obtain optimal efficiency, methods capable of computing locally improved solutions without recourse to global recalculations are pursued. A method for locally refining an existing unstructured mesh, without regenerating a new global mesh is employed, and the domain is automatically partitioned into refined and unrefined regions. Two multigrid strategies are developed. In the first, time-stepping is performed on a global fine mesh covering the entire domain, and convergence acceleration is achieved through the use of zonal coarse grid accelerator meshes, which lie under the adaptively refined regions of the global fine mesh. Both schemes are shown to produce similar convergence rates to each other, and also with respect to a previously developed global multigrid algorithm, which performs time-stepping throughout the entire domain, on each mesh level. However, the present schemes exhibit higher computational efficiency due to the smaller number of operations on each level.

Mavriplis, Dimitri J.

High-Resolution Gridded Level 3 Aerosol Optical Depth Data from MODIS

The state-of-art satellite observations of atmospheric aerosols over the last two decades from NASA's MODIS instruments have been extensively utilized in climate change and air quality research and applications. The operational algorithms now produce level 2 aerosol data at varying spatial resolutions (1, 3, and 10 km) and level 3 data at 1 degree. The local and global applications have been benefited from the coarse resolution gridded data sets (i.e., level 3, 1 degree), as it is easier to use since data volume is low and, several online and offline tools are readily available to access and analyze the data with minimal computing resources. At the same time, researchers who require data at much finer spatial scales have to go through a challenging process of obtaining, processing, and analyzing larger volumes of data sets that require high-end computing resources and coding skills. Therefore, we have created a high spatial resolution (HRG, 0.1x0.1 degree) daily and monthly aerosol optical depth (AOD) product by combining two MODIS operational algorithms, namely Deep Blue (DB) and Dark Target (DT). The new HRG AODs meets the accuracy requirements of level 2 AOD data and provide either the same or more spatial coverage on daily and monthly scales. The data sets are provided in daily and monthly files through open Ftp server with python scripts to read and map the data. The reduced data volume with an easy to use format and tools to access the data will encourage more users to utilize the data for research and applications.

aerosol

On the effect of high latitude filtering in global grid point models

The application of various Fourier filtering techniques in a coarse resolution grid point model and the poor simulation of the nonlinear effects responsible for change in the mean zonal flow is discussed. The poor simulation results in spurious energy and momentum transfers. An operational filter, a second order filter on a staggered grid, a fourth order filter on a nonstaggered grid, and a biharmonic diffusion filter were considered.

Takacs, L. L.

Global Soil Water Estimates as Landslide Predictor: The Effectiveness of SMOS, SMAP, and GRACE Observations, Land Surface Simulations, and Data Assimilation

This global feasibility study assesses the potential of coarse-scale, gridded soil water estimates for the probabilistic modeling of hydrologically-triggered landslides, using Soil Moisture Ocean Salinity (SMOS), Soil Moisture Active Passive (SMAP) and Gravity Recovery and Climate Experiment (GRACE) remote sensing data, Catchment Land Surface Model (CLSM) simulations and six data products based on the assimilation of SMOS, SMAP, and/or GRACE observations into CLSM. SMOS or SMAP observations (~40-km resolution) are only available for less than 20% of the globally reported landslide events, because they are intermittent and uncertain in regions with complex terrain. GRACE terrestrial water storage estimates include 75% of the reported landslides but have coarse spatial and temporal resolutions (monthly, ~300-km). CLSM soil water simulations have the added advantage of complete spatial and temporal coverage, and are found to be able to distinguish between “stable slope” (no landslide) conditions and landslide-inducing conditions in a probabilistic way. Assimilating SMOS and/or GRACE data increases the landslide probability estimates based on soil water percentiles for the reported landslides, relative to model-only estimates at 36-km resolution for the period 2011-2016, unless the CLSM model-only soil water content is already high (≥ 50th percentile). The SMAP Level 4 data assimilation product (at 9-km resolution, period 2015-2019) more generally updates the soil water conditions towards higher landslide probabilities for the reported landslides, but is similar to model-only estimates for the majority of landslides where SMAP data cannot easily be converted to soil moisture owing to complex terrain.

Anne Felsberg

Computation of Flow Over a Drag Prediction Workshop Wing/Body Transport Configuration Using CFL3D

A Drag Prediction Workshop was held in conjunction with the 19th AIAA Applied Aerodynamics Conference in June 2001. The purpose of the workshop was to assess the prediction of drag by computational methods for a wing/body configuration (DLR-F4) representative of subsonic transport aircraft. This report details computed results submitted to this workshop using the Reynolds-averaged Navier-Stokes code CFL3D. Two supplied grids were used: a point-matched 1-to-1 multi-block grid, and an overset multi-block grid. The 1-to-1 grid, generally of much poorer quality and with less streamwise resolution than the overset grid, is found to be too coarse to adequately resolve the surface pressures. However, the global forces and moments are nonetheless similar to those computed using the overset grid. The effect of three different turbulence models is assessed using the 1-to-1 grid. Surface pressures are very similar overall, and the drag variation due to turbulence model is 18 drag counts. Most of this drag variation is in the friction component, and is attributed in part to insufficient grid resolution of the 1-to-1 grid. The misnomer of 'fully turbulent' computations is discussed; comparisons are made using different transition locations and their effects on the global forces and moments are quantified. Finally, the effect of two different versions of a widely used one-equation turbulence model is explored.

Rumsey, Christopher L.

Improving Estimation of Ground Casualty Risk From Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the Earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, Chris L.

Improving Estimation of Ground Casualty Risk from Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination, and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, C.