Engineering Papers⌕ Search

Engineering topics

Pulliam, T. H.

Publications and source records attributed to Pulliam, T. H..

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Recent improvements in efficiency, accuracy, and convergence for implicit approximate factorization algorithms

In 1977 and 1978, general purpose centrally space differenced implicit finite difference codes in two and three dimensions have been introduced. These codes, now called ARC2D and ARC3D, can run either in inviscid or viscous mode for steady or unsteady flow. Since the introduction of the ARC2D and ARC3D codes, overall computational efficiency could be improved by making use of a number of algorithmic changes. These changes are related to the use of a spatially varying time step, the use of a sequence of mesh refinements to establish approximate solutions, implementation of various ways to reduce inversion work, improved numerical dissipation terms, and more implicit treatment of terms. The present investigation has the objective to describe the considered improvements and to quantify advantages and disadvantages. It is found that using established and simple procedures, a computer code can be maintained which is competitive with specialized codes.

Pulliam, T. H.↗
NASA NTRS1985

Artificial dissipation models for the Euler equations

Various artificial dissipation models which are used with central difference algorithms for the Euler equations are analyzed for their effect on accuracy, stability and convergence rates. In particular, linear and nonlinear models are investigated using an implicit approximate factorization code (ARC2D) for transonic airfoils. Fully implicit application of the dissipation models is shown to improve robustness and convergence rates. The treatment of dissipation models at boundaries will be examined. It will be shown that accurate, error free solutions with sharp shocks can be obtained using a central difference algorithm coupled with an appropriate nonlinear artificial dissipation model.

Pulliam, T. H.↗
NASA NTRS1985

Implicit finite-difference methods for the Euler equations

The present paper is concerned with two-dimensional Euler equations and with schemes which are in use of the time of this writing. Most of the development presented carries over directly to three dimensions. The characteristics of the two-dimensional Euler equations in Cartesian coordinates are considered along with generalized curvilinear coordinate transformations, metric relations, invariants of the transformation, flux Jacobian matrices and eigensystems, numerical algorithms, flux split algorithms, implicit and explicit nonlinear control (smoothing), upwind differencing in supersonic regions, unsteady and steady-state computation, the diagonal form of implicit algorithm, metric differencing and invariants, boundary conditions, geometry and mesh generation, and sample solutions.

Pulliam, T. H.↗
NASA NTRS1984

A Newton multigrid method for the Euler equations

A multigrid method is used to apply Newton's method to the Euler equations in a two dimensional curvilinear coordinate system. The objective is to obtain rapid convergence for steady state problems. Solutions computed with the method evolve in a non-time-like manner. Stable pressure distributions typically develop in eight to ten Newton-multigrid steps, which is equivalent to the computational work of about 70 iterations with a factored implicit algorithm.

Childs, R. E.↗
NASA NTRS1984

Airfoil computation at high angles of attack, inviscid and viscous phenomena

An implicit central difference code is used to calculate two dimensional inviscid and thin-layer Navier-Stokes solutions for flow about an NACA0012 airfoil at high angles of attack. Among the issues addressed are whether separation can occur in an inviscid calculation and what the causes would be of such separation. Examples are shown of inviscid shocked flow with and without separation and shock-free flow with separation. An Euler solution with self-induced oscillation and separation driven by a strong shock is contrasted with a shock-free solution whose separation is caused by numerical error. Computed solutions to the Euler equations are compared to those of the potential equations. Comparisons are also made between experimental data from wind tunnel tests and viscous calculations at similar conditions.

Barton, J. T.↗
NASA NTRS1983

An enhanced version of an implicit code for the Euler equations

A two-dimensional implicit finite-difference code is applied to the inviscid Euler equations to compute transonic flow past airfoils in order to provide well-documented standard test cases for the general user community. The code is an improved version of Steger's 1976 implicit code. Enhancements include the use of up-wind differencing in supersonic regions before shocks and a variable time step to accelerate convergence. An airfoil grid generation routine based on algebraic techniques is employed. The grids are clustered near shocks to improve resolution. Computed results are compared with other numerical results from the literature.

Pulliam, T. H.↗
NASA NTRS1983

Flux vector splitting and approximate Newton methods

In the present investigation, the basic approach is employed to view an iterative scheme as Newton's method or as a modified Newton's method. Attention is given to various modified Newton methods which can arise from differencing schemes for the Euler equations. Flux vector splitting is considered as the basic spatial differencing technique. This technique is based on the partition of a flux vector into groups which have certain properties. The Euler equations fluxes can be split into two groups, the first group having a flux Jacobian with all positive eigenvalues, and the second group having a flux Jacobian with all negative eigenvalues. Flux vector splitting based on a velocity-sound speed split is considered along with the use of numerical techniques to analyze nonlinear systems, and the steady Euler equations for quasi-one-dimensional flow in a nozzle. Results are given for steady flows with shocks.

Jespersen, D. C.↗
NASA NTRS1983

A general perturbation approach for the equations of fluid dynamics

An efficient numerical technique to produce accurate solutions to the equations of fluid dynamics is presented where the governing equations are perturbed about an approximate solution and solved by finite-difference methods on a coarsened grid. The result is a scheme which substantially reduces the number of grid points necessary to accurately resolve the flow. Applications are presented for the two-dimensional Euler equations perturbed about a solution of the transonic full potential equation. However, the concept is applicable to arbitrary equation sets, higher dimensions and for a wide variety of applications.

Chow, L. J.↗
NASA NTRS1982

A zonal approach to solution of the Euler equations

A technique for the solution of the one- and two-dimensional Euler equations in a partitioned flow field is presented. The field is divided into distinct 'zones', each of which is computed separately. An implicit boundary procedure based on the characteristic propagaion of information and flux splitting methods is applied at the zonal interfaces. Numerical results are presented for both quasi-one-dimensional nozzle flows with shock waves and the unsteady shock-tube problem. These calculations demonstrate the capability of shock propagation through arbitrarily located zonal boundaries in a stable, conservative, and accurate manner. Two-dimensional results include the zonal computation of flows over blunt bodies and airfoils.

Hessenius, K. A.↗
NASA NTRS1982

Computation of the steady viscous flow over a tri-element 'augmentor wing' airfoil

The augmentor wing consists of a main airfoil with a slotted trailing edge for blowing, and two smaller aft airfoils which shroud the jet. This configuration has been modeled for numerical simulation by a novel discretization procedure which generates four separate grids: three surface-oriented airfoil grids and one outer free-stream grid. Grid lines and slopes are continuous across boundaries, so grid overlap at common boundaries provides boundary information without interpolation. A two-dimensional unsteady thin-layer Navier-Stokes code is used to calculate the flow for the no-blowing case at freestream Mach number = 0.7, Re = 12,600.000, and angles-of-incidence = 1.05 deg. Qualitative agreement with experimental data indicates the utility of this procedure in the analysis of multi-element configurations.

Lasinski, T. A.↗
NASA NTRS1981

A diagonal form of an implicit approximate-factorization algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to partial differential equations is presented. This algorithm is applied to the two- and three-dimensional Euler equations in general curvilinear coordinates. The modification transforms the coupled system of equations into an uncoupled diagonal form that requires less computational work. For steady-state applications, the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. The diagonal algorithm reduces the storage requirement of the implicit solution process and therefore has an important effect on the application of implicit finite-difference schemes to vector processors. Results are presented for realistic two-dimensional transonic flow fields about airfoils. Computation costs are reduced to 24-34%.

Pulliam, T. H.↗
NASA NTRS1981

Two-Dimensional Inlet Simulation Using a Diagonal Implicit Algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two-dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic freestream flow about and through inlets. The modification transforms the coupled system of equations Into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two-dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D.S.↗
NASA NTRS1981

A Diagonal Form of an Implicit Approximate-Factorization Algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to partial differential equations is presented. This algorithm is applied to the two- and three-dimensional Euler equations in general curvilinear coordinates. The modification transforms the coupled system of equations into an uncoupled diagonal form that requires less computational work. For steady-state applications, the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. The diagonal algorithm reduces the storage requirement of the implicit solution process and therefore has an important effect on the application of implicit finite-difference schemes to vector processors. Results are presented for realistic two-dimensional transonic flow fields about airfoils. Computation costs are reduced 24-34%.

Pulliam, T. H.↗
NASA NTRS1981

A numerical simulation of hypersonic viscous flow over arbitrary geometries at angle of attack

An implicit conservative, noniterative, finite-difference algorithm that predicts the supersonic, laminar or turbulent viscous flow about arbitrary geometries at large angles of attack is presented. The three-dimensional parabolized form of the thin-layer Navier-Stokes equations are written in generalized coordinates. These equations are solved using the delta form of the Beam-Warming implicit algorithm. Flow field simulations have been obtained for a blunt biconic with windward and leeward cuts and an X-24C lifting body for both laminar and turbulent flow at various Mach numbers and angles of attack. When compared with experiment or with previous theories, these computational predictions show good agreement.

Chaussee, D. S.↗
NASA NTRS1980

An implicit finite-difference code for inviscid and viscous cascade flow

An implicit finite-difference code is developed to solve either inviscid or viscous flow about two-dimensional cascade blade elements. General coordinate transformations are used so that boundaries can coincide with coordinate lines, and an automatic grid generation routine based on elliptic partial differential equations is employed to mesh arbitrary cascade elements. Characteristic combinations of the differential equations are used at inflow and outflow boundaries. Computed results for both inviscid and viscous flow are compared with other existing cascade solutions and experimental data.

Steger, J. L.↗
NASA NTRS1980

Supersonic flow over three-dimensional ablated nosetips using an unsteady implicit numerical procedure

The three-dimensional supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical domain, containing multiple zones of embedded subsonic flow and separated flow regions into a rectangular computational volume at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing numerical solutions and experimental data for simple spheres which validates the program. Results are also presented for analytically defined indented bodies for both laminar and turbulent flow conditions that demonstrate the program's capability for computing such flows.

Kutler, P.↗
NASA NTRS1980

A diagonal form of an implicit approximate-factorization algorithm with application to a two dimensional inlet

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic free stream flow about and through inlets. The modification transforms the coupled system of equations into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D. S.↗