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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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 NTRS1993

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗
NASA NTRS1975

Upwind second-order difference schemes and applications in unsteady aerodynamic flows

Explicit second-order upwind difference schemes in combination with spatially symmetric schemes can produce larger stability bounds and better numerical resolution than symmetric schemes alone. However, if conservation form is essential, a special operator is required for transition between schemes. An operational approach has been devised for deriving transition operators so that strict conservation and local consistency are maintained. Various aspects of hybrid schemes are studied numerically for model linear and nonlinear equations. To demonstrate the utility of combining two different algorithms, MacCormack's explicit, noncentered, second-order method is combined with a completely upwind version, and numerical solutions of the Euler equations are obtained for two-dimensional, transonic flows with embedded supersonic regions and shock waves. The general utility of the operational approach for combining schemes is emphasized by deriving a second-order conservative scheme for the steady transonic small-disturbance potential equation.

Warming, R. F.↗
NASA NTRS1994

Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure

The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.

Jameson, Antony↗
NASA NTRS1986

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗
NASA NTRS1985

Finite difference methods for reducing numerical diffusion in TEACH-type calculations

A methodological evaluation for two-finite differencing schemes for computer-aided gas turbine design is presented. The two computational schemes include; a Bounded Skewed Finite Differencing Scheme (BSUDS); and a Quadratic Upwind Differencing Scheme (QSDS). In the evaluation, the derivations of the schemes were incorporated into two-dimensional and three-dimensional versions of the Teaching Axisymmetric Characteristics Heuristically (TEACH) computer code. Assessments were made according to performance criteria for the solution of problems of turbulent, laminar, and coannular turbulent flow. The specific performance criteria used in the evaluation were simplicity, accuracy, and computational economy. It is found that the BSUDS scheme performed better with respect to the criteria than the QUDS. Some of the reasons for the more successful performance BSUDS are discussed.

Syed, S. A.↗
NASA NTRS1985

A relaxation approach to patched-grid calculations with the Euler equations

A conservative zonal-boundary condition that was used with explicit integration schemes is extended to implicit, upwind, relaxation schemes; in particular to the Osher scheme. The rate of convergence was found to increase considerably with the use of the implicit, relaxation-zonal scheme when compared to the explicit scheme. The relaxation-zonal scheme has also been used in a time-accurate mode. Results demonstrating the time accuracy of the scheme and the feasibility of performing calculations in cases where some parts of the given system move relative to others (for example, rotor-stator configurations) are presented.

Rai, M. M.↗
NASA NTRS1985

A Relaxation Approach to Patched-Grid Calculations with the Euler Equations

A conservative zonal-boundary condition that was used with explicit integration schemes is extended to implicit, upwind, relaxation schemes; in particular to the Osher scheme, The rate of convergence was found to increase considerably with the use of the implicit, relaxation-zonal-scheme when compared to the explicit scheme. The relaxation-zonal scheme has also been used in a time-accurate mode. Results demonstrating the time accuracy of the scheme and the feasibility of performing calculations in cases where some parts of the given system move relative to others (e.g., rotor-stator configurations) are presented.

Rai, Man Mohan↗
NASA NTRS1989

Two-dimensional Euler computations on a triangular mesh using an upwind, finite-volume scheme

A numerical procedure was developed for the finite-volume solution of the Euler equations on unstructured triangular meshes based on a flux-difference split upwind method. Techniques for implementing Roe's (1985) approximate Reimann solver together with the preprocessing MUSCL differencing on unstructured grids are presented. Applications and comparisons with structured grid problems are carried out for a supersonic shock reflection problem, the supersonic flow over a blunt body, the transonic flow over NACA 0012 and RAE 2822 airfoils, and the flow about a double element Karman-Trefftz airfoil.

Whitaker, D. L.↗
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.↗
NASA NTRS1989

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

A scheme of solving the two-dimensional Euler equations is developed. 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.↗
NASA NTRS1989

Solution techniques for incompressible flow problems

A three-step Petrov-Galerkin (PG)/operator spliting scheme for the time-dependent incompressible Navier-Stokes equations is proposed. Each time step is split into two Stokes problems and one nonlinear convection-diffusion problem. Using a PG technique on the two outer Stokes problems ensures a stable scheme despite equal-order interpolation, while using a streamline upwind PG scheme on the inner convection-diffusion problem ensures a numerically stable solution at high Reynolds numbers. Numerical tests of this method have been carried out.

Tezduyar, T. E.↗
NASA NTRS2007

A Time-Accurate Upwind Unstructured Finite Volume Method for Compressible Flow with Cure of Pathological Behaviors

A time-accurate, upwind, finite volume method for computing compressible flows on unstructured grids is presented. The method is second order accurate in space and time and yields high resolution in the presence of discontinuities. For efficiency, the Roe approximate Riemann solver with an entropy correction is employed. In the basic Euler/Navier-Stokes scheme, many concepts of high order upwind schemes are adopted: the surface flux integrals are carefully treated, a Cauchy-Kowalewski time-stepping scheme is used in the time-marching stage, and a multidimensional limiter is applied in the reconstruction stage. However even with these up-to-date improvements, the basic upwind scheme is still plagued by the so-called "pathological behaviors," e.g., the carbuncle phenomenon, the expansion shock, etc. A solution to these limitations is presented which uses a very simple dissipation model while still preserving second order accuracy. This scheme is referred to as the enhanced time-accurate upwind (ETAU) scheme in this paper. The unstructured grid capability renders flexibility for use in complex geometry; and the present ETAU Euler/Navier-Stokes scheme is capable of handling a broad spectrum of flow regimes from high supersonic to subsonic at very low Mach number, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics). Numerous examples are included to demonstrate the robustness of the methods.

Loh, Ching Y.↗
NASA NTRS1995

Unification of some advection schemes in two dimensions

The relationship between two approaches towards construction of genuinely two-dimensional upwind advection schemes is established. One of these approaches is of the control volume type applicable on structured cartesian meshes. It resulted in the compact high resolution schemes capable of maintaining second order accuracy in both homogeneous and inhomogeneous cases. Another one is the fluctuation splitting approach, which is well suited for triangular (and possibly) unstructured meshes. Understanding the relationship between these two approaches allows us to formulate here a new fluctuation splitting high resolution (i.e. possible use of artificial compression, while maintaining positivity property) scheme. This scheme is shown to be linearity preserving in inhomogeneous as well as homogeneous cases.

Sidilkover, D.↗
NASA NTRS1989

Comparisons of TVD schemes applied to the Navier-Stokes equations

In this study, the following total variation diminishing (TVD) schemes for solving the Navier-Stokes equations have been tested: the Chakravarthy and Szema (1985) upwind biased TVD scheme, the Harten's upwind TVD scheme described by Yee et al. (1983), and the Yee's (1985) symmetric TVD scheme. The schemes have been compared using three test cases. The first case was the one-dimensional shock tube problem which tested the shock-capturing abilities of the schemes. Chakravarthy's and Harten's schemes gave similar results which were found to be more accurate than the results from Yee's scheme. The second case was a compressible boundary layer which tested the schemes's abilities to solve fiscous flows. In this case, the three schemes yielded almost identical results. Finally, the shock/boundary-layer interaction case studied experimentally by Hakkinen et al. (1959) was computed. Here, Chakravarthy's and Yee's schemes compared most favorably with the published data, with Yee's scheme giving slightly better results.

Buelow, Philip E.↗
NASA NTRS1989

Three-dimensional simulation of underexpanded plumes using upwind algorithms

Different numerical techniques have been examined for simulating the underexpanded plumes in supersonic flows. The results revealed that the monotone differentiable limiter is superior to the TVD (total variation diminishing) minmod limiter when the MUSCL (monotonic upstream schemes for conservation laws) approach is used for Roe's upwind scheme. Roe's upwind scheme is found to be superior to the central difference scheme not only for capturing shock waves, but also for capturing shear flows.

Obayashi, Shigeru↗
NASA NTRS1996

Assessment of an Unstructured-Grid Method for Predicting 3-D Turbulent Viscous Flows

A method Is presented for solving turbulent flow problems on three-dimensional unstructured grids. Spatial discretization Is accomplished by a cell-centered finite-volume formulation using an accurate lin- ear reconstruction scheme and upwind flux differencing. Time is advanced by an implicit backward- Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the sublayer region of the boundary layer. A systematic assessment of the method is presented to devise guidelines for more strategic application of the technology to complex problems. The assessment includes the accuracy In predictions of skin-friction coefficient, law-of-the-wall behavior, and surface pressure for a flat-plate turbulent boundary layer, and for the ONERA M6 wing under a high Reynolds number, transonic, separated flow condition.

Frink, Neal T.↗
NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗