Search NASA⌕ Search

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.

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

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

White, J. A.↗
NASA NTRS1986

Higher-order flux difference splitting schemes for the Euler equations using upstream interpolations

A class of explicit two time-level, 2p + 1 space-point, (2p 1)-th order, upwind-biased flux difference splitting schemes are proposed for the numerical advection based on Lagrange's interpolation, and the method is an accord with the physical domain of dependence. A normalized Jacobian coefficient matrix is introduced to convert the schemes to hyperbolic systems of conservation laws, and approaches to make the higher-order schemes total variation stable are discussed. Accuracy and stability of the present schemes are examined, and implicit total variation diminishing schemes are developed for steady-state calculations.Application to gasdynamic problems for both steady and unsteady flows covering a wide range of Mach numbers is considered, and results for a blast wave passing a cylinder, and head-on collision of two blast waves over a circular arc, are presented. The flow patterns were found to be symmetric, and good resolution of flow structures was obtained.

Yang, J. Y.↗
NASA NTRS1991

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗
NASA NTRS1992

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗
DOE OSTI2021

A high-order WENO-limited finite-volume algorithm for atmospheric flow using the ADER-differential transform time discretization

A high-order-accurate weighted essentially non-oscillatory (WENO) limited upwind finite-volume scheme is detailed for the compressible, nonhydrostatic, inviscid Euler equations using an arbitrary derivatives (ADER) time-stepping scheme based on differential transforms (DTs). A second-order-accurate alternating Strang dimensional splitting is compared against multidimensional simulation with 2D transport using solid body rotation of various data. The two were found to give nearly identical accuracy in orthogonal, Cartesian coordinates. Orders of convergence are demonstrated at up to ninth-order accuracy with 2D transport. 1D transport is used to confirm that error decreases monotonically with increasing order of accuracy with WENO limiting even for discontinuous data. Further, WENO limiting always decreased the error compared with simulation without limiting in the L 1 norm. A series of standard 2D compressible nonhydrostatic Euler equation test cases were validated against previous results from literature. Finally, it was demonstrated that increasing the order of accuracy led to better resolved features and increased power for kinetic energy at small wavelengths.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Second-order closure analysis of turbulent reacting flows

A reacting shear layer program for a complete second-order closure of turbulent reacting flows is developed. The program uses a typical eddy box model and many of the fluid mechanical turbulence models developed in earlier studies. The computer program actually solves a total of 30 equations, where the numerical integration of the equations is performed by a forward-time-centered space quasi-implicit upwind finite-difference scheme. Results for nonreacting and simple isothermal reacting flows are presented along with preliminary results for propane-air flames. It is shown that the program predictions for simple nonreacting flows are in good agreement with available data. The calculations for isothermal reacting flows demonstrate the importance of including a mixedness correlation in finite rate chemistry calculations, where significant errors in the flowfield properties can result from neglect of the interaction between turbulence and the chemistry. Results for propane-air diffusion flames further support this conclusion.

Varma, A. K.↗
NASA NTRS1983

A numerical study of the steady scalar convective diffusion equation for small viscosity

A time-independent convection diffusion equation is studied by means of a compact finite difference scheme and numerical solutions are compared to the analytic inviscid solutions. The correct internal and external boundary layer behavior is observed, due to an inherent feature of the scheme which automatically produces upwind differencing in inviscid regions and the correct viscous behavior in viscous regions.

Giles, M. B.↗
NASA NTRS1984

A numerical study of the steady scalar convective diffusion equation for small viscosity

A time-independent convection diffusin equation is studied by means of a compact finite difference scheme and numerical solutions are compared to the analytic inviscid solutions. The correct internal and external boundary layer behavior is observed, due to an inherent feature of the scheme which automatically produces upwind differencing in inviscid regions and the correct viscous behavior in viscous regions.

Giles, M. B.↗
NASA NTRS1987

Computation of viscous blast wave solutions with an upwind finite volume method

A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.

Molvik, Gregory A.↗
NASA NTRS1988

A longitudinally-patched grid approach with applications to high speed flows

The combined use of an implicit, upwind finite-volume scheme with an 'equivalent' gamma formulation (for real gas effects) and with a three-dimensional conservative patched grid formulation is discussed. Results are presented over a wide Mach number range on both single and patched grids.

Walters, R. W.↗
NASA NTRS1988

Improving the accuracy of central difference schemes

General difference approximations to the fluid dynamic equations require an artificial viscosity in order to converge to a steady state. This artificial viscosity serves two purposes. One is to suppress high frequency noise which is not damped by the central differences. The second purpose is to introduce an entropy-like condition so that shocks can be captured. These viscosities need a coefficient to measure the amount of viscosity to be added. In the standard scheme, a scalar coefficient is used based on the spectral radius of the Jacobian of the convective flux. However, this can add too much viscosity to the slower waves. Hence, it is suggested that a matrix viscosity be used. This gives an appropriate viscosity for each wave component. With this matrix valued coefficient, the central difference scheme becomes closer to upwind biased methods.

Turkel, Eli↗
NASA NTRS1989

Unsteady Euler airfoil solutions using unstructured dynamic meshes

Two algorithms for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis of oscillating airfoils. Both algorithms were developed for use on an unstructured grid made up of triangles. The first flow solver involves a Runge-Kutta time-stepping scheme with a finite-volume spatial discretization that reduces to central differencing on a rectangular mesh. The second flow solver involves a modified Euler time-integration scheme with an upwind-biased spatial discretization based on the flux-vector splitting of Van Leer. The paper presents descriptions of the Euler solvers and dynamic mesh algorithm along with results which assess the capability.

Batina, John T.↗
NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗
NASA NTRS1988

Navier-Stokes calculations of scramjet-afterbody flowfields

Scramjet simulant-gas exhaust flows are computed using two-dimensional Navier-Stokes (NS) equations. The NS equations are solved using an implicit, upwind, finite-volume scheme, and the Reynolds stresses are modeled algebraically. The solutions are compared with experimental data where the freestream and the nozzle exhaust fluids are air. The NS equations for the flow, where the scramjet exhaust simulant gas is a mixture of Freon-12 and argon, are solved by an explicit, finite difference scheme. The results are compared with experimental data, and certain discrepancies are explained.

Baysal, O.↗
NASA NTRS1989

Studies of vortex flow aerodynamics using CFD flow visualizations

Navier-Stokes computations of three-dimensional vortical flows over a round-edged double-delta wing and a tangent-ogive forebody are performed using an implicit upwind finite-difference scheme. Computed particle traces are compared with experimental oil-flow streaks.

Hsu, C.-H.↗
NASA NTRS1990

Upwind Algorithm For Parabolized Navier-Stokes Equations

Supersonic flow about cone calculated accurately. Report presents theoretical basis of computer code solving parabolized Navier-Stokes equations of supersonic and hypersonic flow. For increased accuracy in resolution of details of strong aerodynamic shocks, code incorporates implicit, finite-volume, upwind numerical-integration scheme. Performs well in numerical simulations of flows around simple bodies.

Lawrence, Scott L.↗
NASA NTRS1990

Applications of an adaptive unstructured solution algorithm to the analysis of high speed flows

An upwind cell-centered scheme for the solution of steady laminar viscous high-speed flows is implemented on unstructured two-dimensional meshes. The first-order implementation employs Roe's (1981) approximate Riemann solver, and a higher-order extension is produced by using linear reconstruction with limiting. The procedure is applied to the solution of inviscid subsonic flow over an airfoil, inviscid supersonic flow past a cylinder, and viscous hypersonic flow past a double ellipse. A detailed study is then made of a hypersonic laminar viscous flow on a 24-deg compression corner. It is shown that good agreement is achieved with previous predictions using finite-difference and finite-volume schemes. However, these predictions do not agree with experimental observations. With refinement of the structured grid at the leading edge, good agreement with experimental observations for the distributions of wall pressure, heating rate and skin friction is obtained.

Thareja, R. R.↗
NASA NTRS1990

Computational study for passive control of supersonic asymmetric vortical flows around cones

The unsteady, compressible, double and single thin-layer Navier-Stokes equations are used to numerically study for the passive control of steady and unsteady, supersonic asymmetric flows around circular and noncircular cones. The main computational scheme of the present study is the implicit upwind, flux-difference splitting, finite-volume scheme. Passive control of flow asymmetry is studied by using a vertical fin in the leeward plane of geometric symmetry and side strakes with and without thickness at different orientations. The study is focused on circular-section cones since they are the most potential section-shapes for strong flow asymmetry. Side-strakes passive control has been shown to be more efficient and practical than the vertical-fin passive control.

Liu, C. H.↗