Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data

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

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

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

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

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

Boundary Approximations for Implicit Schemes for .One-Dimensional Inviscid Equations of Gasdynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations is examined for initial boundary-value problems of the hyperbolic type. For the numerical experiments the one-dimensional inviscid gasdynamic equations in conservation law form are 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 vs implicit boundary schemes, and unconditionally stable schemes

Stability Analysis of Numerical Boundary Conditions and Implicit Difference Approximations for Hyperbolic Equations

Implicit noniterative finite-difference schemes have recently been 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 has been that unconditional stability requires implicit numerical boundary conditions. An apparent counterexample was the space-time extrapolation considered by Gustafsson, Kreiss, and Sundstrom. In this paper we examine space (implicit) and space-time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals. The results indicate that for unconditional stability with a finite number of spatial mesh intervals the numerical boundary conditions must be implicit

On the Application and Extension of Harten\u27s High-Resolution Scheme

Xost first-order upstream conservative differencing methods can capture shocks quite well for one-dimensional problems. A direct application of these first-order methods to two-dimensional problems does not necessarily produce the same type of accuracy unless the shocks are locally aligned with the mesh. Harten has recently developed a second-order high-resolution explicit method for the numerical computation of weak solutions of one-dimensional hyperbolic conservation laws. The main objectives of this paper are (a) to examine the shock resolution of Harten\u27s method for a two-dimensional shock reflection problem, (b) to study the use of a high-resolution scheme as a post-processor to an approximate steady-state solution, and (c) to construct an implicit method in the delta-form using Harten\u27s scheme for the explicit operator and a simplified iteration matrix for the implicit operator