SEARCH · Search NASA

Results for “Upwinding”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1985

On the implementation of a class of upwind schemes for system of hyperbolic conservation laws

The relative computational effort among the spatially five point numerical flux functions of Harten, van Leer, and Osher and Chakravarthy is explored. These three methods typify the design principles most often used in constructing higher than first order upwind total variation diminishing (TVD) schemes. For the scalar case the difference in operation count between any two algorithms may be very small and yet the operation count for their system counterparts might be vastly different. The situation occurs even though one starts with two different yet equivalent representations for the scalar case.

Yee, H. C.↗

NASA NTRS1986

Upwind Navier-Stokes solutions for separated periodic flows

The application of an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using both the thin-layer and complete forms of the Navier-Stokes equations, the low Reynolds number laminar flow around a circular cylinder with periodic shedding is solved. The effect of grid density, grid extent, and time step on the Strouhal number is shown. Results from both sets of equations agree within the experimental data band. Unsteady, laminar flow computations around inclined plates and separated airfoils are also described. Strouhal numbers agree to within 5 percent of experiments for inclined plates. Differences between the complete equations and the thin-layer approximation for separated periodic flows are discussed. Computations of an impulsively started circular cylinder and airfoil yield time-accurate flowfield shapes in good agreement with experimental flow visualizations. The turbulent computation of an airfoil at a high angle-of-attack is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗

NASA NTRS1986

Computation of turbulent pipe and duct flow using third order upwind scheme

The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.

Kawamura, T.↗

NASA NTRS1986

An upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method does not require the addition of user specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming scheme in terms of accuracy, stability, computer time and storage, and programming effort. The new algorithm has been validated by applying it to three laminar test cases including flat plate boundary-layer flow, hypersonic flow past a 15 deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with the results obtained using the conventional Beam-Warming algorithm.

Lawrence, S. L.↗

NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗

NASA NTRS1987

Application of an upwind Navier-Stokes code to two-dimensional transonic airfoil flow

An upwind-biased implicit approximate factorization Navier-Stokes algorithm is applied to a variety of steady transonic airfoil cases, using the NACA 0012, RAE 2822, and Jones supercritical airfoils. The thin-layer form of the compressible Navier-Stokes equations is used. Both the CYBER 205 and CRAY 2 supercomputers are utilized, with average computational speeds of about 18 and 16 microsec/gridpoint/iteration, respectively. Lift curves, drag polars, and variations in drag coefficient with Mach number are determined for the NACA 0012 and Jones supercritical airfoils. Also, several cases are computed for comparison with experiment. The effect of grid density and grid extent on a typical turbulent airfoil solution is shown. An algebraic eddy-viscosity turbulence model is used for all of the computations.

Rumsey, Christopher L.↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

NASA NTRS1987

Application of an upwind algorithm to the three-dimensional parabolized Navier-Stokes equations

A new computer code for the solution of the three-dimensional parabolized Navier-Stokes equations has been developed. The code employs a state-of-the-art upwind algorithm to capture strong shock waves. The algorithm is implicit, uses finite volumes, and is second-order accurate in the crossflow directions. The new code is validated through application to laminar hypersonic flows past two simple body shapes: a circular cone of 10 deg half-angle, and a generic all-body hypersonic vehicle. Cone flow solutions were computed at angles of attack of 12, 20, and 24 deg and results are in agreement with experimental data. Results are also presented for the flow past the all-body vehicle at angles of incidence of 0 and 10 deg.

Lawrence, Scott L.↗

NASA NTRS1987

Three dimensional hypersonic flow simulations with the CSCM implicit upwind Navier-Stokes method

The Bardina and Lombard (1985) three-dimensional CSCM Navier-Stokes method is presently extended to the simulation of complex hypersonic reentry vehicle external flows at angle of attack. The robust stability of the method derives from the combination of conservative implicit upwind flux difference splitting with a three-dimensional diagonally-dominant approximate factorization and relaxation scheme and characteristic-based implicit boundary approximations. The method's efficiency derives from an implicit symmetric Gauss-Seidel 'method of planes' relaxation scheme with alternating directional space marching sweeps along the flow coordinate direction.

Bardina, Jorge↗

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 NTRS1987

Supercomputing of supersonic flows using upwind relaxation and MacCormack schemes

The performance of two numerical solution schemes, (1) an implicit upwind relaxation with a finite-volume discretization (Thomas and Walters, 1985) and (2) an explicit-implicit MacCormack (1981) scheme with a finite-difference discretization, is compared in two-dimensional simulations of supersonic flow past a flat plate with leading edge, a rearward-facing step, a 10-deg compression corner, a NACA 0012 airfoil at high angle of attack, and a cavity. The algebraic turbulence model, the solution methods, and the boundary conditions and SIMD coding are explained, and the results are presented in tables and graphs and characterized with reference to published experimental data. Scheme (1) is found to converge more rapidly and to give more accurate results than (2) in a wide range of problem types.

Baysal, O.↗

NASA NTRS1987

Modification of the Osher upwind scheme for use in three dimensions

The Osher upwind scheme, a shock-capturing algorithm based on an approximate Riemann solver, which has been developed for the solution of hyperbolic conservation laws, is briefly described, and the need for modification of the scheme for use in three dimensions is demonstrated. Two forms of a scheme modification for use on a three-dimensional arbitrary grid are then presented.

Hessenius, Kristin A.↗

NASA NTRS1988

Viscous flow computations using a second-order upwind differencing scheme

In the present computations of a wide range of fluid flow problems by means of the primitive variables-incorporating Navier-Stokes equations, a mixed second-order upwinding scheme approximates the convective terms of the transport equations and the scheme's accuracy is verified for convection-dominated high Re number flow problems. An adaptive dissipation scheme is used as a monotonic supersonic shock flow capture mechanism. Many benchmark fluid flow problems, including the compressible and incompressible, laminar and turbulent, over a wide range of M and Re numbers, are presently studied to verify the accuracy and robustness of this numerical method.

Chen, Y. S.↗

NASA NTRS1988

Explicit upwind algorithm for the parabolized Navier-Stokes equations

A new explicit upwind algorithm based on Roe's flux-difference splitting (FDS) method has been developed for the three-dimensional Parabolized Navier-Stokes (PNS) equations. For three-dimensional flows, FDS's are determined separately for the two nonmarching directions and modified to account for the calculated shock angle in the crossflow plane. Second-order FDS is applied to the pressure and convection terms with the streamwise pressure gradient limited in the subsonic region to maintain a hyperbolic inviscid equation set. Second-order central differencing is obtained in the two-step algorithm for the shear and heat flux terms. The new algorithm is demonstrated for three laminar flow test cases: supersonic flow over a flat plate, hypersonic flow over a 15 deg ramp, and hypersonic flow past a 10 deg cone at a 24 deg angle of attack. The computed results agree well with experimental measurements.

Korte, John J.↗

NASA NTRS1987

Streamline-upwind/Petrov-Galerkin procedures for the vorticity-stream function form of the Navier-Stokes equations

The paper presents procedures for the solution of the Navier-Stokes equations in the vorticity-stream function form. The difficulties involved are related to the convection term in the vorticity transport equation and to the lack of boundary conditions for voritcity at no-slip surfaces. The use of a streamline-upwind/Petrov-Galerkin finite element formulation for the solution of the vorticity transport equation. In the present scheme, the weighting functions are dependent on both spatial and temporal discretizations. A proper numerical treatment of the boundary conditions leading to an implicit treatment of the vorticity at no-slip surfaces is presented. These procedures have successfully been employed to simulate various flows of engineering interest.

Tezduyar, T. E.↗

NASA NTRS1988

Numerical simulation of underexpanded plumes using upwind algorithms

The solution of the flow characteristics of highly underexpanded plumes is presently sought by means of a high-resolution total variation-diminishing upwind scheme featuring finite-volume metrics, Roe averaging with entropy correction, and monotonic upstream schemes for conservation laws. The existing approximate factorization schemes are used both as preconditioners and in order to apply the conjugate residual algorithm. The plumes treated are axisymmetric and laminar, and for an ideal gas that is allowed to be thermally incomplete. Simulations of the SSME and SRB plumes have been successfully conducted.

Obayashi, Shigeru↗

NASA NTRS1988

A comparative study and validation of upwind and central-difference Navier-Stokes codes for high-speed flows

A comparative study was made using 4 different computer codes for solving the compressible Navier-Stokes equations. Three different test problems were used, each of which has features typical of high speed internal flow problems of practical importance in the design and analysis of propulsion systems for advanced hypersonic vehicles. These problems are the supersonic flow between two walls, one of which contains a 10 deg compression ramp, the flow through a hypersonic inlet, and the flow in a 3-D corner formed by the intersection of two symmetric wedges. Three of the computer codes use similar recently developed implicit upwind differencing technology, while the fourth uses a well established explicit method. The computed results were compared with experimental data where available.

Rudy, David H.↗

NASA NTRS1989

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗