Engineering Papers⌕ Search

Engineering topics

Warming, R. F.

Publications and source records attributed to Warming, R. F..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗
NASA NTRS1985

Application of TVD schemes for the Euler equations of gas dynamics

First-order, second-order, and implicit total variation diminishing (TVD) schemes are reviewed using the modified flux approach. Some transient and steady-state calculations are then carried out to illustrate the applicability of these schemes to the Euler equations. It is shown that the second-order explicit TVD schemes generate good shock resolution for both transient and steady-state one-dimensional and two-dimensional problems. Numerical experiments for a quasi-one-dimensional nozzle problem show that the second-order implicit TVD scheme produces a fairly rapid convergence rate and remains stable even when running with a Courant number of 10 to the 6th.

Yee, H. C.↗
NASA NTRS1985

Application of TVD schemes for the Euler equations of gas dynamics

Highly accurate and yet stable shock-capturing finite difference schemes have been designed for the computation of the Euler equations of gas dynamics. Four different principles for the construction of high resolution total variation diminishing (TVD) schemes are available, including hybrid schemes, a second-order extension of Godunov's scheme by van Leer (1979), the modified flux approach of Harten (1983, 1984), and the numerical fluctuation approach of Roe (1985). The present paper has the objective to review the class of second-order TVD schemes via the modified flux approach. Attention is given to first-order TVD schemes, a second-order accurate explicit TVD scheme, the global order of accuracy of the second-order TVD scheme, extensions to systems and two-dimensional conservation laws, numerical experiments with a second-order explicit TVD scheme, implicit TVD schemes, and second-order implicit TVD schemes.

Yee, H. C.↗
NASA NTRS1985

On a class of TVD schemes for gas dynamic calculations

The purpose of this paper is to review a class of explicit and implicit second-order accurate Total Variation Diminishing (TVD) schemes and to show by numerical experiments, the performance of these schemes to the Euler equations of gas dynamics. The method of constructing these second-order accurate TVD schemes is sometimes known as the modified flux approach.

Yee, H. C.↗
NASA NTRS1983

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme.

Yee, H. C.↗
NASA NTRS1983

Implicit total variation diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. Previously announced in STAR as N83-23085

Yee, H. C.↗
NASA NTRS1982

On the application and extension of Harten's high resolution scheme

Extensions of a second order high resolution explicit method for the numerical computation of weak solutions of one dimensonal hyperbolic conservation laws are discussed. The main objectives were (1) to examine the shock resoluton of Harten's method for a two dimensional shock reflection problem, (2) to study the use of a high resolution scheme as a post-processor to an approximate steady state solution, and (3) to construct an implicit in the delta-form using Harten's scheme for the explicit operator and a simplified iteration matrix for the implicit operator.

Yee, H. C.↗
NASA NTRS1982

A high-resolution numerical technique for inviscid gas-dynamic problems with weak solutions

The shock resolution of Harten's (1982) second-order explicit method for one-dimensional hyperbolic conservation laws is investigated for a two-dimensional gas-dynamic problem. The possible extension to a high resolution implicit method for both one- and two-dimensional problems is also investigated. Applications of Harten's method to the quasi-one-dimensional nozzle problem with two nozzle shapes (divergent and convergent-divergent) and the two-dimensional shock-reflection problem resulted in high shock resolution steady-state numerical solutions.

Yee, H. C.↗
NASA NTRS1981

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗
NASA NTRS1981

Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods

The conservation-law form of the inviscid gasdynamic equations has the remarkable property that the nonlinear flux vectors are homogeneous functions of degree one. This property readily permits the splitting of flux vectors into subvectors by similarity transformations so that each subvector has associated with it a specified eigenvalue spectrum. As a consequence of flux vector splitting, new explicit and implicit dissipative finite-difference schemes are developed for first-order hyperbolic systems of equations. Appropriate one-sided spatial differences for each split flux vector are used throughout the computational field even if the flow is locally subsonic. The results of some preliminary numerical computations are included.

Steger, J. L.↗
NASA NTRS1981

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗
NASA NTRS1981

Recent advances in the development of implicit schemes for the equations of fluid dynamics

Innovations and extensions of implicit schemes for equations of fluid dynamics are presented. The notation and theory for linear multistep methods are reviewed, and extensions of work by Beam and Warming (1979) include the implementation of one-leg methods, ADI methods for equations with mixed derivatives, flux vector splitting, the P-dimensional wave equation, and boundary conditions. Numerical experiments indicate that implicit treatment of the boundary conditions is necessary for unconditional stability, and the improvement and implementation of the boundary condition theory should improve the implicit algorithms for gas dynamic equations.

Warming, R. F.↗
NASA NTRS1980

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A(0)-stable linear two-step methods in conjunction with the method of approximate factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A(0)-stable method. The parameter space for which the resulting ADI schemes are second-order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗
NASA NTRS1979

Flux vector splitting of the inviscid equations with application to finite difference methods

The conservation-law form of the inviscid gasdynamic equations has the remarkable property that the nonlinear flux vectors are homogeneous functions of degree one. This property readily permits the splitting of flux vectors into subvectors by similarity transformations so that each subvector has associated with it a specified eigenvalue spectrum. As a consequence of flux vector splitting, new explicit and implicit dissipative finite-difference schemes are developed for first-order hyperbolic systems of equations. Appropriate one-sided spatial differences for each split flux vector are used throughout the computational field even if the flow is locally subsonic. The results of some preliminary numerical computations are included.

Steger, J. L.↗
NASA NTRS1979

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A sub 0-stable linear two-step methods in conjunction with the method of approximation factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A sub 0-stable method. The parameter space for which the resulting ADI schemes are second order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗
NASA NTRS1979

An implicit factored scheme for the compressible Navier-Stokes equations. II - The numerical ODE connection

An attempt is made to establish a connection between linear multistep methods for applications to ordinary differential equations and their extension (by approximate factorization) to alternating direction implicit methods for partial differential equations. An earlier implicit factored scheme for the compressible Navier-Stokes equations is generalized by innovations that (1) increase the class of temporal difference schemes to include all linear multistep methods, (2) optimize the class of unconditionally stable factored schemes by a new choice of unknown variable, and (3) improve the computational efficiency by the introduction of quasi-one-leg methods.

Beam, R. M.↗
NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗