SEARCH · Search NASA

Results for “upwind schemes”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for CFD. This is especially important for solution of complex three-dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for both their capability of resolving discontinuities and their sound theoretical basis in characteristic theory for hyperbolic systems. With this in mind, two new flux splitting techniques are presented for upwind differencing.

Liou, MENG-S.↗

NASA NTRS1992

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

NASA NTRS1993

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

NASA NTRS1986

Three-dimensional, conservative, Euler computations using patched grid systems and explicit methods

The method of 'zonal approach' (in which the flow field is partitioned into regions with independent grids) for computation of flow over complex geometries, such as aircraft configuration, requires application of a grid-interfacing procedure. A three-dimensional conservative, boundary scheme for patched grids, applicable in generalized coordinates for arbitrary point distributions on a planar zonal surface is presented. The computation technique is derived within the framework of the Osher (1983) upwind scheme, using Euler equations. The three-dimensional interfacing method is applied to the computation of flow about a wing-canard combination using a two-zone, patched grid.

Hessenius, K. A.↗

NASA NTRS1988

Investigation of supersonic chemically reacting and radiating channel flow

The 2-D time-dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multistage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multidimensional radiative transfer equations for a nongray model are provided for a general configuration and then reduced for a planar geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. The supersonic inviscid and viscous, nonreacting flows are solved by employing the finite volume technique of Jameson and the unsplit finite difference scheme of MacCormack. The specified problem considered is of the flow in a channel with a 10 deg compression-expansion ramp. The calculated results are compared with those of an upwind scheme. The problem of chemically reacting and radiating flows are solved for the flow of premixed hydrogen-air through a channel with parallel boundaries, and a channel with a compression corner. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the entire flow field.

Mani, Mortaza↗

NASA NTRS1992

Parabolized Navier-Stokes methods for hypersonic flows

The development and numerical solution of the Parabolized Navier-Stokes (PNS) equations is discussed. Special emphasis is placed on numerical integration algorithms for PNS solvers, most notably recent applications of upwind schemes. Application of this technology to high speed, high Reynolds number flows was the focus, and modeling techniques for real gas effects and turbulence were discussed. Experiences in validating the UPS upwind PNS flow solver for relatively simple body shapes are presented. The test cases range in complexity from two dimensional perfect gas cases to three dimensional finite rate chemistry cases.

Lawrence, S. L.↗

NASA NTRS1993

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗

NASA NTRS1995

Numerical simulation of supersonic and hypersonic inlet flow fields

This report summarizes the research performed by North Carolina State University and NASA Ames Research Center under Cooperative Agreement NCA2-719, 'Numerical Simulation of Supersonic and Hypersonic Inlet Flow Fields". Four distinct rotated upwind schemes were developed and investigated to determine accuracy and practicality. The scheme found to have the best combination of attributes, including reduction to grid alignment with no rotation, was the cell centered non-orthogonal (CCNO) scheme. In 2D, the CCNO scheme improved rotation when flux interpolation was extended to second order. In 3D, improvements were less dramatic in all cases, with second order flux interpolation showing the least improvement over grid aligned upwinding. The reduction in improvement is attributed to uncertainty in determining optimum rotation angle and difficulty in performing accurate and efficient interpolation of the angle in 3D. The CCNO rotational technique will prove very useful for increasing accuracy when second order interpolation is not appropriate and will materially improve inlet flow solutions.

Mcrae, D. Scott↗

NASA NTRS1985

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities.

Davis, S. F.↗

NASA NTRS1986

Kinetic theory based new upwind methods for inviscid compressible flows

Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.

Deshpande, S. M.↗

NASA NTRS1986

Implicit hybrid schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Implicit hybrid algorithms employing symmetric planar Gauss-Seidel (SPGS) relaxation and either block-tridiagonally structured coefficient matrices (AF-SPGS) or block-triangular coefficient matrices (LU-SPGS) are derived to solve the flux-difference-split Navier-Stokes equations for three-dimensional incompressible flow in an upwind scheme. The physical basis of the approach is discussed, and results for problems involving vortex flow around a thin delta wing at Reynolds numbers 900,000 and 10,000 are presented graphically. It is found that AF-SPGS converges faster on vector computers which depend on long vector lengths to achieve optimum performance, whereas LU-SPGS is preferable on sequentially operating machines and vector computers using shorter vector lengths.

Hartwich, P.-M.↗

NASA NTRS1987

Vectorizable implicit algorithms for the flux-difference split, three-dimensional Navier-Stokes equations

The computational efficiency of four vectorizable implicit algorithms is assessed when applied to calculate steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in general coordinates. Two of these algorithms are characterized as hybrid schemes; that is, they combine some approximate factorization in two coordinate directions with relaxation in the remaining spatial direction. The other two algorithms utilize an approximate factorization approach which yields two-factor algorithms for three-dimensional systems. All four algorithms are implemented in identical high-resolution upwind schemes for the flux-difference split Navier-Stokes equations. These highly nonlinear schemes are obtained by extending an implicit Total Variation Diminishing (TVD) scheme recently developed for linear one-dimensional systems of hyperbolic conservation laws to the three-dimensional Navier-Stokes equations. The computation of vortical flow over a sharp-edged, thin delta wing has been chosen as a common numerical test case. The convergence of the algorithms is discussed and the accuracy of the computed flow-field results is assessed. The validity of the present results are demonstrated by a comparison with experimental data.

Hartwich, P. M.↗

NASA NTRS1989

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗

NASA NTRS1989

Numerical simulation of flow over a hypersonic aircraft using an explicit upwind PNS solver

A hypersonic flow field over a generic airplane configuration is simulated by solving the Parabolized Navier-Stokes (PNS) equations. The finite difference solution of the PNS equations is calculated using a noniterative space marching, explicit, upwind scheme recently developed by the authors. Special gridding techniques are used which allowed the sharp changes in surface geometry of the airplane configuration to be modelled without smoothing of corners. Comparisons of the PNS results to a solution of the Navier-Stokes equations demonstrates a good agreement of the numerical results in approximately 1/6 of the cpu time. This paper demonstrates that the explicit upwind algorithm for solving the PNS equations is an efficient method for simulating hypersonic flow fields about complete airplane configurations and should be considered as an alternative to solving the Navier-Stokes equations for flow fields where the PNS equations are valid.

Korte, John J.↗

NASA NTRS1990

An interactive adaptive remeshing algorithm for the two-dimensional Euler equations

An interactive adaptive remeshing algorithm utilizing a frontal grid generator and a variety of time integration schemes for the two-dimensional Euler equations on unstructured meshes is presented. Several device dependent interactive graphics interfaces have been developed along with a device independent DI-3000 interface which can be employed on any computer that has the supporting software including the Cray-2 supercomputers Voyager and Navier. The time integration methods available include: an explicit four stage Runge-Kutta and a fully implicit LU decomposition. A cell-centered finite volume upwind scheme utilizing Roe's approximate Riemann solver is developed. To obtain higher order accurate results a monotone linear reconstruction procedure proposed by Barth is utilized. Results for flow over a transonic circular arc and flow through a supersonic nozzle are examined.

Slack, David C.↗

NASA NTRS1990

Efficient real gas Navier-Stokes computations of high speed flows using an LU scheme

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three dimensional Navier-Stokes calculations was demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature and transport properties are directly related to the flow variables. Since the pressure is explicitly known as a general function of the flow variables no assumptions are made regarding the pressure derivatives in the construction of the flux Jacobians. The method is efficient since no sub-iterations are required to deduce the pressure and temperature from the flux variables and allows different equations of state to be easily supplied to the code. The flexibility of the EOS approach is demonstrated by implementing a high order TVD upwinding scheme based upon flux differencing and Van Leer's flux vector splitting. The EOS approach is demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates and through a simplified model of a scramjet module gap-seal configuration.

Coirier, William J.↗

NASA NTRS1990

Efficient real gas Navier-Stokes computations of high speed flows using an LU scheme

Coirier, William J.↗

DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗