SEARCH · Engineering Papers

Results for “Conservation equations”

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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Diurnal variations in the thermosphere. I - Theoretical formulation

A nonlinear perturbation theory is formulated for the solution of the multicomponent equations of energy, mass, and momentum conservation in the atmosphere. The theory is three-dimensional and includes the effects of heat conduction and advection, viscosity, ion drag, and diffusion. The theory is described as a superposition of mathematical modes obtained by expanding the physical quantities into vector and spherical harmonics. The coupling between the various modes, both linear and nonlinear, is included. The theory provides a basis for the treatment of the thermosphere and its interaction with the lower atmosphere, where 'mode coupling' is most important. As an example, a comparison is presented between one-dimensional and three-dimensional calculations of the fundamental mode of the diurnal component in the thermosphere. Coupling between the lowest modes is considered to describe the physical conditions of the lower thermosphere where inertia and Coriolis forces become dominant over the ion-drag and viscous forces. In this region, the latitude structures of the temperature, wind field, and diffusively controlled oxygen are shown to change significantly.

Harris, I.↗

NASA NTRS1976

The multigrid method: Fast relaxation

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g. 4 x 2 mesh cells) to fine (e.g. 64 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with optimal iterations on the finest grid. To date the method was applied quite successfully to the solution of the transonic small-disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

NASA NTRS1979

Numerical modeling of two-dimensional confined flows

A numerical model of two-dimensional confined flows is presented. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying conservation principles directly to the individual streams. A listing of a computer code based on this approach in FORTRAN 4 language is presented. The code computes two dimensional compressible turbulent flows in ducts when the duct area along the flow is specified and the pressure gradient is unknown.

Greywall, M. S.↗

NASA NTRS1976

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse to fine. The coarser grids are used to diminish the magnitude of the smooth part of the residuals. The method was applied to the solution of the transonic small disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is studied with meshes of both constant and variable step size.

South, J. C., Jr.↗

NASA NTRS1981

Full potential solution of transonic quasi-3-D flow through a cascade using artificial compressability

The three-dimensional flow in a turbomachinery blade row was approximated by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on body fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗

NASA NTRS1981

Full potential solution of a transonic quasi-3-D flow through a cascade using artificial compressibility

A reliable method is presented for calculating the flowfield about a cascade of arbitrary 2-D airfoils. The method approximates the three-dimensional flow in a turbomachinery blade row by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to ensure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme for analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗

NASA NTRS1982

Computer program for calculating full potential transonic, quasi-three-dimensional flow through a rotating turbomachinery blade row

A fast, reliable computer code is described for calculating the flow field about a cascade of arbitrary two dimensional airfoils. The method approximates the three dimensional flow in a turbomachinery blade row by correcting for stream tube convergence and radius change in the throughflow direction. A fully conservative solution of the full potential equation is combined with the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. The instructions required to set up and use the code are included. The name of the code is QSONIC. A numerical example is also given to illustrate the output of the program.

Farrell, C. A.↗

NASA NTRS1983

A multi-grid method for the computation of viscid/inviscid interactions on airfoils

An improved version of the 'GRUMFOIL' code has been developed for the computation of airfoil flows. The method employs a conservative difference scheme for the potential equation, integral methods for the boundary layer, and viscous coupling conditions that fully account for the wake and strong interaction effects at trailing edges. The improved version uses Jameson's 'MAD' scheme to accelerate convergence of the inviscid solution, an improved 2nd order artificial viscosity and far field 'TARE' correction to reduce spatial truncation errors and Carter's semi-inverse method for the viscous solution. Results are presented which demonstrate a factor of ten reduction in computing cost.

Melnik, R. E.↗

NASA NTRS1984

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradient through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallel piped computational domain, and an explicit MacCormack (1972) explicit time-split predictor corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

NASA NTRS1984

An implicit form for the Osher upwind scheme

Conservative upwind schemes for the Euler equations, such as the Osher scheme, accurately resolve flow discontinuities and correctly model the physics of the problem. However, these schemes require many more arithmetic operations per integration step than simple central-difference schemes and hence result in large computing times. An implicit version of the first-order- and second-order-accurate Osher schemes in two spatial dimensions and generalized coordinates is developed in this study. Because implicit schemes permit the use of large integration steps, in many cases they require fewer integration steps to reach steady-state (especially in calculations on grids with widely varying mesh-cell sizes). The implicit scheme developed in this study accelerated convergence speeds by almost an order of magnitude in the problems considered. Test cases include quasi-one-dimensional nozzle flow and supersonic flow past a cylinder.

Rai, M. M.↗

NASA NTRS1984

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradients through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallelepiped computational domain, and an explicit MacCormack (1972) explicit time-split predictor-corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

NASA NTRS1984

A combined direct/inverse three-dimensional transonic wing design method for vector computers

A three-dimensional transonic-wing design algorithm for vector computers is developed, and the results of sample computations are presented graphically. The method incorporates the direct/inverse scheme of Carlson (1975), a Cartesian grid system with boundary conditions applied at a mean plane, and a potential-flow solver based on the conservative form of the full potential equation and using the ZEBRA II vectorizable solution algorithm of South et al. (1980). The accuracy and consistency of the method with regard to direct and inverse analysis and trailing-edge closure are verified in the test computations.

Weed, R. A.↗

NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗

NASA NTRS1985

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗

NASA NTRS1986

A generalization of the method of harmonic balance

A procedure is provided for generalizing the method of harmonic balance to obtain higher-order approximations to the periodic solutions of differential equations for two systems, one a conservative system and the other nonconservative. The procedure is currently being applied to investigate the possible solution behaviors of singular, nonlinear oscillators, where the usual perturbation methods do not work.

Mickens, R. E.↗

NASA NTRS1988

A stiffly-stable implicit Runge-Kutta algorithm for CFD applications

A stiffly-stable implicit Runge-Kutta integration algorithm is derived for CFD applications spanning the range of semidiscrete theories. The algorithm family contains the one-step 'theta' algorithms, including backwards Euler and the trapezoidal rule, and provides a versatile framework to identify expressions governing algorithm stability characteristics. Parameters of a Runge-Kutta optimal implicit algorithm, second-order accurate in time and stiffly-stable, are established. This algorithm is implemented within a weak statement finite element semidiscrete formulation for one- and two-dimensional conservation law systems. Numerical results are compared to theta-algorithm solutions, for unsteady quasi-one-dimensional Euler predictions with shocks, and for a specially derived two-dimensional conservation law system modeling the Euler equations.

Baker, A. J.↗

NASA NTRS1987

Computation of transonic potential flow on helicopter rotor blades

Two computer codes, the full-potential three-dimensional (FP3D) code and the full-potential rotor (FPR) code have recently been developed. Both of these codes solve the three-dimensional conservative formulation of the full potential equation. The FPR code was developed at the U.S. Army Aeroflightdynamics Directorate (AFDD) while the FP3D code was a joint development by ONERA and AFDD. Both of these codes were used to predict the nonlifting, unsteady flow over a rotor operating at high advance ratio and tip speed. Three different rotor tip planform shapes were studied: a rectangular tip, a 30 deg aft swept tip and a 30 deg forward swept tip. Results of these computations are compared to results obtained using an earlier small-disturbances code. Also, the lifting flow over a rectangular tip operating at a slightly different condition was computed. These results are also compared with the small disturbances computations and with experimental results.

Costes, M.↗

NASA NTRS1987

MHD aspects of magnetotail dynamics

The equations of ideal MHD represent the conservation of mass, momentum, and energy; thus, MHD must, in some broad sense, provide a correct description of large-scale structure of the earth's magnetosphere. However, steady-state, ideal MHD (the simplest description) allows only a static interaction between the solar wind and the magnetosphere. Hence an MHD theory of magnetospheric dynamics must involve unsteady flow turbulence, such as might be generated by Kelvin-Helmholtz turbulence at the magnetopause, and dissipative transport, such as reconnection, which violates ideal MHD. This brief review summarizes recent theoretical efforts to explore and extend the early MHD magnetospheric models of Axford and Petschek. Analytic calculations show that steady convection in the magnetotail is nearly impossible. Global magnetospheric simulations exhibit substorm reconnection phenomena for southward IMF and predict a novel convection pattern for northward IMF. Reconnection simulations show plasmoid formation and field-aligned current flows in the tail. The Kelvin-Helmholtz instability can provide a strong anomalous viscosity which can drive the viscous convection system.

Coroniti, F. V.↗