SEARCH · Engineering Papers

Results for “total variation”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1989

On some numerical schemes for transonic flow problems

New second-order-accurate finite-difference approximations for a class of nonlinear PDEs of mixed type, which includes the two-dimensional low-frequency transonic small-disturbance equation (TSD) and full-potential equation (FP), are presented. For the TSD equation, the scheme is implemented via a time-splitting algorithm; the inclusion of flux limiters keeps the total variation nonincreasing and eliminates spurious oscillations near shocks. Global-linear-stability, TVD, and entropy-stability results are proved. Numerical results for the flow over a thin airfoil are presented. Current techniques used to solve the TSD equation may easily be extended to second-order accuracy by this method. For the FP equation, the new scheme requires no subsonic/supersonic switching and no numerical flux biasing. Global linear stability for all values of the Mach number is proved.

Mostrel, Marco Mosche↗

NASA NTRS1989

Analysis of a second-order-accurate finite-volume method for temporally-growing compressible shear layers

A finite-volume method for solving the compressible laminar Navier-Stokes equation in two dimensions is assessed for use in the simulation of transitional flows. The method is second-order-accurate, uses alternating-direction-implicit time integration, and a total-variation-diminishing smoothing operator in the inviscid flux terms. The test problem was a free shear layer with a forced periodicity in x at a specified wavelength and with the parallel mean flow specified as a hyperbolic-tangent profile. A grid refinement investigation verified this method to be second-order-accurate with less than 4 percent error in the growth rate of isolated modes on a 32 by 64 grid. The interaction between modes of similar amplitude was small, less than 1 percent. However, in cases where there was a dominant mode, numerical phase error pumped energy into all other modes. This error results in nonphysical growth rate for modes whose energy content is several orders of magnitude below the dominant mode. In the case when the dominant mode saturates, the pumping action stops and the growth rate of the next largest growing mode regains physical significance.

Atkins, H. L.↗

NASA NTRS1990

Triangle based TVD schemes for hyperbolic conservation laws

A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented.

Durlofsky, Louis J.↗

NASA NTRS1990

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

Gnoffo, Peter A.↗

NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

NASA NTRS1990

Some Aspects of Essentially Nonoscillatory (ENO) Formulations for the Euler Equations, Part 3

An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.

Chakravarthy, Sukumar R.↗

NASA NTRS1990

On central-difference and upwind schemes

A class of numerical dissipation models for central-difference schemes constructed with second- and fourth-difference terms is considered. The notion of matrix dissipation associated with upwind schemes is used to establish improved shock capturing capability for these models. In addition, conditions are given that guarantee that such dissipation models produce a Total Variation Diminishing (TVD) scheme. Appropriate switches for this type of model to ensure satisfaction of the TVD property are presented. Significant improvements in the accuracy of a central-difference scheme are demonstrated by computing both inviscid and viscous transonic airfoil flows.

Swanson, R. C.↗

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 NTRS1991

High-Resolution Computations Of Hypersonic Flows

Report discusses extension, to hypersonic flows, of class of implicit, total-variation-diminishing (TVD) algorithms suitable for numerical simulation of transonic and supersonic flows obeying Euler and Navier-Stokes equations. Stabilities and rates of convergence, relative to those of other algorithms discussed. Study complemented by variety of compuations of steady and unsteady viscous and inviscid hypersonic flows about blunt bodies.

Yee, H. C.↗

NASA NTRS1989

Conservation equations and physical models for hypersonic air flows over the aeroassist flight experiment vehicle

The code development and application program for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the Aeroassist Flight Experiment (AFE) in the near term and Aeroassisted Space Transfer Vehicle (ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3-D, viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE wind tunnel model. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium over the full scale AFE configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

NASA NTRS1990

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗

NASA NTRS1991

An upwind-biased space marching algorithm for supersonic viscous flow

The modifications are documented which were made to the Langley Aerothermodynamic Upwind Relaxation Algorithm which allow it to compute solutions in a space marching manner. The space marching flux is formulated to be either first- or second-order accurate, using Roe's upwind differencing or Symmetric Total Variation Diminishing differencing, respectively. The algorithm solves the thin layer Navier-Stokes equations, and is subject to the same flow restrictions as a parabolized Navier-Stokes solver. Each cross flow plane is locally iterated in time until converged, then marched in space to the next station. The algorithm is tested on a sphere-cone geometry and a geometry which models the windward side of the Space Shuttle Orbiter. Computational results for surface heating are compared to ground based experimental data. In addition, space marching predictions for surface pressure are compared against values from the original algorithm.

Greene, Francis A.↗

NASA NTRS1991

Computation of viscous blast wave flowfields

A method to determine unsteady solutions of the Navier-Stokes equations was developed and applied. The structural finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using MacCormack's modified Steger-Warming flux or Roe flux difference splittings with total variation diminishing limiters, while the viscous flux is computed using central differences. The use of implicit boundary conditions in conjunction with a telescoping in time and space method permitted solutions to this strongly unsteady class of problems. Comparisons of numerical, analytical, and experimental results were made in two and three dimensions. These comparisons revealed accurate wave speed resolution with nonoscillatory discontinuity capturing. The purpose of this effort was to address the three-dimensional, viscous blast-wave problem. Test cases were undertaken to reveal these methods' weaknesses in three regimes: (1) viscous-dominated flow; (2) complex unsteady flow; and (3) three-dimensional flow. Comparisons of these computations to analytic and experimental results provided initial validation of the resultant code. Addition details on the numerical method and on the validation can be found in the appendix. Presently, the code is capable of single zone computations with selection of any permutation of solid wall or flow-through boundaries.

Atwood, Christopher A.↗

NASA NTRS1990

Three-dimensional Euler time accurate simulations of fan rotor-stator interactions

A numerical method useful to describe unsteady 3-D flow fields within turbomachinery stages is presented. The method solves the compressible, time dependent, Euler conservation equations with a finite volume, flux splitting, total variation diminishing, approximately factored, implicit scheme. Multiblock composite gridding is used to partition the flow field into a specified arrangement of blocks with static and dynamic interfaces. The code is optimized to take full advantage of the processing power and speed of the Cray Y/MP supercomputer. The method is applied to the computation of the flow field within a single stage, axial flow fan, thus reproducing the unsteady 3-D rotor-stator interaction.

Boretti, A. A.↗

NASA NTRS1991

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

Yungster, Shaye↗

NASA NTRS1990

A numerical study of hypersonic stagnation heat transfer predictions at a coordinate singularity

The problem of grid induced errors associated with a coordinate singularity on heating predictions in the stagnation region of a three-dimensional body in hypersonic flow is examined. The test problem is for Mach 10 flow over an Aeroassist Flight Experiment configuration. This configuration is composed of an elliptic nose, a raked elliptic cone, and a circular shoulder. Irregularities in the heating predictions in the vicinity of the coordinate singularity, located at the axis of the elliptic nose near the stagnation point, are examined with respect to grid refinement and grid restructuring. The algorithm is derived using a finite-volume formulation. An upwind-biased total-variation diminishing scheme is employed for the inviscid flux contribution, and central differences are used for the viscous terms.

Grasso, Francesco↗

NASA NTRS1991

Efficient real gas upwinded Navier-Stokes computations of high speed flows

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three-dimensional Navier-Stokes calculations has been demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature, and transport properties are directly related to the flow variables. The method is efficient since no subiterations are required to deduce the pressure and temperature from the flux variables and is modular by allowing different equations of state to be easily supplied to the code. The flexibility of the EOS approach is shown by its implementation into a high-order total variation diminishing upwinding scheme as well as a standard central-differencing scheme. The EOS approach is then demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates using both the upwind and central-differencing schemes.

Coirier, William J.↗

NASA NTRS1990

A user guide for the EMTAC-MZ CFD code

The computer code (EMTAC-MZ) was applied to investigate the flow field over a variety of very complex three-dimensional (3-D) configurations across the Mach number range (subsonic, transonic, supersonic, and hypersonic flow). In the code, a finite volume, multizone implementation of high accuracy, total variation diminishing (TVD) formulation (based on Roe's scheme) is used to solve the unsteady Euler equations. In the supersonic regions of the flow, an infinitely large time step and a space-marching scheme is employed. A finite time step and a relaxation or 3-D approximate factorization method is used in subsonic flow regions. The multizone technique allows very complicated configurations to be modeled without geometry modifications, and can easily handle combined internal and external flow problems. An elliptic grid generation package is built into the EMTAC-MZ code. To generate the computational grid, only the surface geometry data are required. Results obtained for a variety of configurations, such as fighter-like configurations (F-14, AVSTOL), flow through inlet, multi-bodies (shuttle with external tank and SRBs), are reported and shown to be in good agreement with available experimental data.

Szema, Kuo-Yen↗