SEARCH · Search NASA
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.
A global MHD model of the earth's magnetosphere is defined. An introduction to numerical methods for solving the MHD equations is given with emphasis on the shock-capturing technique. Finally, results concerning the shape of the magnetosphere and the plasma flows inside the magnetosphere are presented.
The interactive phenomena that occur in supersonic jet mixing flowfields, and numerical modeling techniques developed to analyze such phenomena are discussed. A spatial marching procedure based on solving the parabolized Navier-Stokes jet mixing equations is presented. This procedure combines shock-capturing methodology for the analysis of supersonic mixing regions with pressure-split methodology for the analysis of subsonic mixing regions. The two regions are coupled at viscous sonic lines utilizing a viscous-characteristic coupling procedure. Specialized techniques for the treatment of jet boundary growth, strong discontinuties (Mach disks), and small embedded subsonic zones (behind Mach disks) are presented. Turbulent processes are represented by two-equation turbulence model formulations. In Part II of this article, numerical studies are presented for a variety of supersonic jet interactive phenomena.
A new computational model (SCIPVIS) is described which predicts the multiple-cell wave/shock structure in underexpanded or overexpanded turbulent jets. SCIPVIS solves the parabolized Navier-Stokes jet-mixing equations utilizing a shock-capturing approach in supersonic regions of the jet and a pressure-split approach in subsonic regions. Turbulence processes are represented by the solution of compressibility-corrected two-equation turbulence models. SCIPVIS presently analyzes jets exhausting into a quiescent or supersonic external stream for which a single-pass spatial-marching solution can be obtained. The features of SCIPVIS are reviewed, and calculations are described exhibiting the influence of turbulence modelling, jet temperature, and flight velocity on the jet shock structure.
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.
The present numerical method for simulating the formation and propagation of interplanetary shocks is based on the shock-capturing finite difference scheme of Lax (1950) and Lax and Wendroff (1960), as well as the recent method of NEAR characteristics of Nakagawa (1980, 1981). Attention is given to examples which strongly suggest that all the shocked solar wind plasma parameters due to given physical perturbations, such as flare-generated shocks, can be predicted through the use of this method; the method is, however, limited to the supersonic and super-Alfvenic flow.
A brief review of 2D PNS methodology is first presented which describes the specialized features of supersonic shock-capturing and subsonic pressure-split models required for the analysis of aircraft, rocket and scramjet jet mixing problems. These features include techniques for dealing with various types of embedded and interfacing subsonic regions, the inclusion of finite-rate chemistry and the direct-coupling with potential flow solutions. Preliminary 3D extensions of this PNS methodology geared to supersonic and subsonic rectangular free jet mixing problems are also reviewed. New 3D PNS work will be described which includes the development of a hybrid supersonic/subsonic free jet mixing model, and, a supersonic model geared to the analysis of turbulent mixing and combustion processes occurring in scramjet combustor/nozzle flowfields.
An overview of computational models developed for the complete, design-oriented analysis of a scramjet propulsion system is provided. The modular approach taken involves the use of different PNS models to analyze the individual propulsion system components. The external compression and internal inlet flowfields are analyzed by the SCRAMP and SCRINT components discussed in Part II of this paper. The combustor is analyzed by the SCORCH code which is based upon SPLITP PNS pressure-split methodology formulated by Dash and Sinha. The nozzle is analyzed by the SCHNOZ code which is based upon SCIPVIS PNS shock-capturing methodology formulated by Dash and Wolf. The current status of these models, previous developments leading to this status, and, progress towards future hybrid and 3D versions are discussed in this paper.
Flows calculated for realistic engine-inlet conditions. Computer code LAPIN, large-perturbation inlet, developed to analyze large-perturbation, transient-flow fields in supersonic inlets. Robust, quick-running code capable of solving unsteady quasi-one-dimensional, inviscid-flow problems in mixed subsonic and supersonic regimes for inlets. Approach based upon quasi-one-dimensional, inviscid, unsteady formulation including engineering models of unstart/restart, bleed, bypass, and geometrical effects. Numerical solution of governing time-dependent equations of motion accomplished through shock-capturing, finite-difference algorithm. Program written in FORTRAN IV.
A modified version of the multiscale turbulence model of Hanjalic has been applied to the problem of supersonic jets exhausting into still air. In particular, the problem of shock-cell decay through turbulent interaction with the mixing layer has been studied for both mildly interacting and strongly resonant jet conditions. The modified Hanjalic model takes into account the nonequilibrium energy transfer between two different turbulent spectral scales. The turbulence model was incorporated into an existing shock-capturing, parabolized Navier-Stokes computational model in order to perform numerical experiments. The results show that the two-scale turbulence model provides significant improvement over one-scale models in the prediction of plume shock structure for underexpanded supersonic (Mach 2) and sonic (Mach 1) jets. For the supersonic jet, excellent agreement with experiment was obtained for the centerline shock-cell pressure decay up to 40 jet radii. For the sonic jet, the agreement with experiment was not so good, but the two-scale model still showed significant improvement over the one-scale model. It is shown that by relating some of the coefficients in the turbulent-transport equations to the relative time scale for transfer of energy between scales the two-scale model can provide predictions that bound the measured shock-cell decay rate for the sonic jet.
A two-dimensional, high resolution shock-capturing algorithm was used on a supercomputer to solve Eulerian gasdynamic equations in order to simulate nonstationary strong shock diffraction by a circular arc model in a shock tube. The hypersonic Mach shock wave was assumed to arrive at a high angle of incidence, and attention was given to the effect of varying values of the ratio of specific heats on the shock diffraction process. Details of the conservation equations of the numerical algorithm, written in curvilinear coordinates, are provided, and model output is illustrated with the results generated for a Mach shock encountering a 15 deg circular arc. The sample graphics include isopycnics, a shock surface density profile, and pressure and Mach number contours.
The shock-capturing capability of total variation diminishing (TVD) schemes is demonstrated for a more realistic complex shock-diffraction problem for which the experimental data are available. Second-order explicit upwind and symmetric TVD schemes are used to solve the time-dependent Euler equations of gas dynamics for the interaction of a blast wave with an airfoil at high angle-of-attack. The test cases considered are a time-dependent moving curved-shock wave and a contant moving planar-shock wave impinging at an angle-of-attack 30 deg on a NACA 0018 airfoil. Good agreement is obtained between isopycnic contours computed by the TVD schemes and those from experimental interferograms. No drastic difference in flow-field structure is found between the curved- and planar-shock wave cases, except for a difference in density level near the lower surface of the airfoil. Computation for cases with higher shock Mach numbers is also possible. Numerical experiments show that the symmetric TVD scheme is less sensitive to the boundary conditions treatment than the upwind scheme.
A modified version of the multiscale turbulence model of Hanjalic has been applied to the problem of underexpanded supersonic jets. In particular, the shock-cell decay resulting from shock-mixing layer interactions has been studied for both mildly interacting and strongly resonant jet conditions. A version of the Hanjalic model which accounts for nonequilibrium energy transfer between two spectral scales of turbulence was incorporated into an existing shock-capturing, parabolized Navier-Stokes computational model in order to perform numerical experiments. Results are presented for nominal initial jet Mach numbers of 2.0, 1.4, and 1.0 and are compared to experiments and to predictions made using single-scale models. The results show significant effects of multi-spectral turbulent energy transfer on the predicted shock-cell decay particularly for the lower jet Mach numbers.
Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.
An implicit, finite-difference, shock-capturing algorithm that calculates inviscid, hypersonic flows in chemical equilibrium is presented. The flux vectors and flux Jacobians are differenced using a first-order, flux-split technique. The equilibrium composition of the gas is determined by minimizing the Gibbs free energy at every node point. The code is validated by comparing results over an axisymmetric hemisphere against previously published results. The algorithm is also applied to more practical configurations. The accuracy, stability, and versatility of the algorithm have been promising.
In the present continuation of development and analysis efforts towards essentially nonoscillatory shock-capturing approximations of hyperbolic conservation laws, the Godunov (1959) scheme and its MUSCL second-order-accurate extension are generalized to an arbitrary order of accuracy by means of a hierarchy of uniformly high-order-accurate schemes. The design method, which involves the solution's essentially nonoscillatory piecewise polynomial reconstruction on the basis of its cell averages, followed by time-evolution through an approximate solution of the resulting initial value problem and its averaging over each cell, employs an adaptive stencil of grid points to yield schemes that are highly nonlinear.
A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicated that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.
The Osher upwind scheme, a shock-capturing algorithm based on an approximate Riemann solver, which has been developed for the solution of hyperbolic conservation laws, is briefly described, and the need for modification of the scheme for use in three dimensions is demonstrated. Two forms of a scheme modification for use on a three-dimensional arbitrary grid are then presented.
An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.