SEARCH · Search NASA

Results for “Variable Grid Method”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Accuracy of schemes with nonuniform meshes for compressible fluid flows

The accuracy of the space discretization for time-dependent problems when a nonuniform mesh is used is considered. Many schemes reduce to first-order accuracy while a popular finite volume scheme is even inconsistent for general grids. This accuracy is based on physical variables. However, when accuracy is measured in computational variables then second-order accuracy can be obtained. This is meaningful only if the mesh accurately reflects the properties of the solution. In addition, the stability properties of some improved accurate schemes are analyzed and it can be shown that they also allow for larger time steps when Runge-Kutta type methods are used to advance in time.

Turkel, E.↗

NASA NTRS1986

Accuracy of schemes with nonuniform meshes for compressible fluid flows

Turkel, Eli↗

NASA NTRS1982

Developing lake sampling strategies with existing Landsat data

On the basis of the assumptions that Landsat image densities relate to water quality, that lake densities can be modeled, and that homogeneous water quality areas can therefore be predicted, a method is developed for the selection of optimum water quality sampling locations on Oneida Lake in central New York. Patterns of Oneida Lake density on 24 Landsat band 5 images were traced from the screen of a density slicer, and the lake was gridded into 161 1.25 X 1.25-km cells, with image density levels encoded for each cell. Cell density levels were then employed as the dependent variable in correlations and multiple regression equations with degree days, resultant wind speed, and resultant direction. Cell density levels were averaged over two-week periods, and the averages were correlated with water quality measurements.

Philipson, W. R.↗

NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

NASA NTRS1994

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

NASA NTRS1994

Approaches of comparison for clear-sky radiative fluxes from general circulation models with Earth Radiation Budget Experiment data

In order to compare the clear-sky greenhouse effect and cloud-radiative forcing from general circulation models with Earth Radiation Budget Experiments (ERBE) data, it is necessary to calculate the general circulation model (GCM) clear-sky radiative fluxes in a way consistent with ERBE. This study discusses problems associated with the available methods for clear-sky radiative flux computations in GCMs and proposes a new approach, which uses a statistical relationship between the grid cloud cover and the availability of ERBE clear-sky measurement, established from ERBE pixel data, to sample the model radiative fluxes. Calculations with version 2 of the National Center for Atmospheric Research (NCAR) Community Climate Model using observed sea surface temperature (SST) show good agreement of clear-sky sampling from the proposed method with ERBE sampling. It is also shown that large improvements are achieved in the spatial variability of the model clear-sky radiative fluxes over ocean, with reference to ERBE, by using the new clear-sky sampling method.

Zhang, M. H.↗

NASA NTRS1994

Visualizing Flow Over Parametric Surfaces Using Line Integral Convolution

Line Integral Convolution (LIC) is a powerful technique for imaging and animating vector fields. We extend the LIC paradigm in three ways: (1) The existing technique is limited to vector fields over a regular Cartesian grid. We extend it to vector fields over parametric surfaces, such as those found in curvilinear grids, used in computational fluid dynamics simulations; (2) Periodic motion filters can be used to animate the flow visualization. When the flow lies on a parametric surface, however, the motion appears misleading. We explain why this problem arises and show how to adjust the LIC algorithm to handle it; (3) We introduce a technique to visualize vector magnitudes as well as vector direction. Cabral and Leedom have suggested a method for variable-speed animation, which is based on varying the frequency of the filter function. We develop a different technique based on kernel phase shifts which we have found to show substantially better results. Our implementation of these algorithms utilizes texture-mapping hardware to run in real time, which allows them to be included in interactive applications.

Forssell, Lisa↗

NASA NTRS1987

Statistics and dynamics of persistent anomalies

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Mo, Kingtse C.↗

NASA NTRS1994

Multiblock, Multigrid Solution Of Euler Equations

Method of numerical solution of Euler equations of three-dimensional flow of compressible fluid involves combination of multiblock and multigrid strategies. In multiblock strategy, flow field divided, into multiple smaller, more computationally-convenient zones and computational grid fitted to applicable flow boundaries generated in each block. In multigrid strategy used here, different quantities computed, variously, on finer or coarser grids. Minimizing cost of computation by using fewest grid points yielding acceptably accurate values of affected variable. Multigrid strategy found effective in accelerating convergence to steady state, while multiblock strategy provides geometric flexibility.

Melson, N. Duane↗

NASA NTRS1995

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗

NASA NTRS1980

Effect of grid system on finite element calculation

Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

Lee, K. D.↗

NASA NTRS2000

Accuracy Improvement in Magnetic Field Modeling for an Axisymmetric Electromagnet

This paper examines the accuracy and calculation speed for the magnetic field computation in an axisymmetric electromagnet. Different numerical techniques, based on an adaptive nonuniform grid, high order finite difference approximations, and semi-analitical calculation of boundary conditions are considered. These techniques are being applied to the modeling of the Variable Specific Impulse Magnetoplasma Rocket. For high-accuracy calculations, a fourth-order scheme offers dramatic advantages over a second order scheme. For complex physical configurations of interest in plasma propulsion, a second-order scheme with nonuniform mesh gives the best results. Also, the relative advantages of various methods are described when the speed of computation is an important consideration.

Ilin, Andrew V.↗

NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

NASA NTRS1991

A multiple-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations

A multiple block multigrid method for the solution of the three dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell vertex method which employs central difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques and variable coefficient numerical dissipation results in an efficient and robust scheme is discussed. The multiblock strategy places the block loop within the Runge-Kutta Loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one and two block calculations in which the two block grid is generated by splitting the one block grid. Results are presented for both Euler and Navier-Stokes computations of wing/fuselage combinations.

Atkins, Harold↗

NASA NTRS1976

Diffusion analysis for two-phase metal-matrix composite

Diffusion controlled filament matrix interaction in a metal matrix composite, where the filaments and matrix comprise a two phase binary alloy system, was mathematically modeled. The problem of a diffusion controlled, two phase moving interface by means of a one dimensional, variable grid, finite difference technique was analyzed. Concentration dependent diffusion coefficients and equilibrium solubility limits were used, and the change in filament diameter and compositional changes in the matrix were calculated as a function of exposure time at elevated temperatures. With the tungsten nickel (W-Ni) system as a model composite system, unidirectional composites containing from 0.06 to 0.44 initial filament volume fraction were modeled. Compositional changes in the matrix were calculated by superposition of the contributions from neighboring filaments. Alternate methods for determining compositional changes between first and second nearest neighbor filaments were also considered. The results show the relative importance of filament volume fraction, filament diameter, exposure temperature, and exposure time as they affect the rate and extent of filament matrix interaction.

Tenney, D. R.↗

NASA NTRS2001

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗

NASA NTRS1986

Improving the convergence rate to steady state of parabolic ADI methods

The present, residuals' L(2)-norms analysis of the rate of convergence to steady state for parabolic ADI solvers allows the prediction of the number of iterations required for convergence, as a function of the Courant number alpha. A modification of current ADI codes is presented which significantly improves the convergence rate and is insensitive to the Courant number over a large range of alpha. This corrected algorithm is tested for the cases of Dirichlet problems for uniform grids of many mesh sizes, mixed Dirichlet-Neumann problems, and problems defined on stretched grids and/or problems with variable coefficients.

Abarbanel, Saul S.↗

NASA NTRS1975

Sasaki's variational optimization analysis of temperature and moisture advection in a severe storm environment

Temperature and moisture advection patterns of dry air intrusion were analyzed by a variational optimization method to relate mesoscale advection to a severe thunderstorm environment. The analysis results presented here were obtained from a case study, May 26, 1973, when a tornado squall line developed in Oklahoma in conjunction with a clearly delineated dry air intrusion. Observations from the synoptic upper air network were used for inputs of horizontal wind components, potential temperature, and mixing ratio at each 50 mb increment from 950 to 100 mb. An objective analysis was performed to assign values of the variables to grid points of a 15 x 14 horizontal grid with approximate 125 km grid spacing at each pressure level. After the analysis the data fields were smoothed using a variational formalism.

Mcfarland, M. J.↗