SEARCH · Engineering Papers
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.
The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.
Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.
Solutions for the inviscid and viscous supersonic conical flow in a complete rectilinear inlet consisting of four planes intersecting at arbitrary wedge and sweep angles are obtained. To compute the flow on a specially constructed mesh, a three-dimensional flow solver ARC3D is used. It is shown that a single-pass mesh re-adaptation procedure can be used to obtain improved shock capture without the necessity of modifying the flow code, provided that the flow solver works in curvilinear coordinates and is restartable. Results computed by the method show good agreement with experimental measurements and previous calculations.
Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.
It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.
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.
Supersonic Implicit Marching Program (SIMP) applies numerical method, based on conservative form of full potential equation, to problem of threedimensional supersonic flows with embedded subsonic regions. Conservative formulation of problem provides ability to capture shocks and to assess accurately impact of sweep, thickness, and lift for conditions where linear theory unsatisfactory. Technique uses characteristic signal-propagation theory to control density biasing for treatment of shocks (including embedded shocks) and mixed elliptic/hyperbolic crossflow. SIMP written in FORTRAN 77.
A shock fitting method applied to the transonic small disturbance (TSD) potential equation is described. This method is then applied to a simple, two dimensional (2-D) rotating disturbance which is analogous to a shock radiating from the tip of a rotor blade in high speed hover. A comparison is made between the results of this method and the more standard shock capturing method. This comparison makes it clear that the effect of the results on the acoustic signature of the 2-D model is significant, and similar results can be expected when the method is extended to the three dimensional (3-D) case.
A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The lower-upper implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.
In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.
The finite volume explicit time marching method was refined and improved. Previously, extension had been made to the finite volume method to improve the accuracy of the calculation of total pressure in inviscid flow, extend the method to allow the calculation of laminar and turbulent boundary layers in internal flows, and improve the shock capturing properties of the method by introducing a Mach number dependent interpolation scheme for the pressure used in the calculating the density. The current work extends these developments by using the new pressure interpolation scheme in two dimensional viscous calculations, including a more complete description of the viscous stresses, introducing a criteria for the transverse upwind differencing which is a function of the ratio of transverse and streamwise mass fluxes, and allowing the calculation of internal flow where boundary layers are present on both walls of the duct. The manner in which the viscous stresses are evaluated in the nonorthogonal, nonuniform grid is detailed. The convergence is investigated and results for calculations of laminar flow in a converging duct are presented. Results for calculations of transonic flow in a converging-diverging nozzle are presented and the results are compared with Sajben's measurements and calculations by others.
Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.
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.
A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The LU implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.
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.