Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity equation”

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.

At least 253 records · Page 14

Linear and nonlinear interpretation of the direct strike lightning response of the NASA F106B thunderstorm research aircraft

The objective of the work reported here is to develop a methodology by which electromagnetic measurements of inflight lightning strike data can be understood and extended to other aircraft. A linear and time invariant approach based on a combination of Fourier transform and three dimensional finite difference techniques is demonstrated. This approach can obtain the lightning channel current in the absence of the aircraft for given channel characteristic impedance and resistive loading. The model is applied to several measurements from the NASA F106B lightning research program. A non-linear three dimensional finite difference code has also been developed to study the response of the F106B to a lightning leader attachment. This model includes three species air chemistry and fluid continuity equations and can incorporate an experimentally based streamer formulation. Calculated responses are presented for various attachment locations and leader parameters. The results are compared qualitatively with measured inflight data.

Rudolph, T. H.↗

The X-ray background and the evolution of quasars

The contribution of QSOs and active galaxies to the 2-10-keV X-ray background is calculated theoretically, assuming that the X-ray luminosity of each GSO has a power-law time evolution (dL/dt = AL exp alpha) and applying the continuity equation to derive the X-ray luminosity function at arbitrary redshift. The observed X-ray background is shown to require alpha greater than 1.2, ruling out pure luminosity evolution (alpha = 1).

Tucker, W. H.↗

Incompressible Navier-Stokes solutions for a sharp-edged double-delta wing

An implicit finite-difference scheme is used to compute the incompressible laminar vortical flow around a double-delta wing with an aspect ratio of 2.06. By adding a time derivative of the pressure to the continuity equation, the unsteady incompressible Navier-Stokes equations can be integrated like a conventional parabolic time-dependent system of equations. The flux-difference split scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The trajectory of the vortical core is well predicted in comparison with experimental data.

Hsu, Chung-Hao↗

Derivation of photochemical information near 1 mbar from ozone and temperature data

An analysis procedure has been developed for derivation of information about the photochemical behavior of ozone near 1 mbar by explicitly accounting for the dynamical transport terms in the continuity equation for perturbations from the zonal mean. The procedure is shown to be valid using data from a numerical transport model and is then applied to LIMS ozone and temperature data, using geostrophic winds to estimate the transport terms. The data study is restricted to March at 2, 1, and 0.7 mbar. Because the temperature deviations are dynamically produced, large temperature deviations are associated with significant ozone transport terms. The anticorrelation between the deviations of ozone and temperature disappears when the transport terms are small. The derived photochemical information is compared to photochemical theory. Although there is overall agreement in the magnitude as well as the latitude, altitude, and time dependencies, discrepancies are suggested which may be related to the long-standing failure of photochemical models to calculate ozone accurately near 1 mbar. The theory suggests that the addition of Cl(x) to the stratosphere will affect the relationship of ozone and temperature at 2 and 1 mbar. Comparison of the photochemical information derived from future ozone and temperature measurements with the results of the present analysis should provide a critical test of the photochemical scheme thought to describe the behavior of odd chlorine in the stratosphere.

Douglass, Anne R.↗

Feedback between neutral winds and auroral arc electrodynamics

The feedback between neutral atmospheric winds and the electrodynamics of a stable, discrete auroral arc is analyzed. The ionospheric current continuity equation and the equation for neutral gas acceleration by ion drag are solved simultaneously, as a function of time. The results show that, in general, the electric field in the ionosphere adjusts to neutral wind acceleration so as to keep auroral field-aligned currents and electron acceleration approximately independent of time. It is thus concluded that the neutral winds that develop as a result of the electrodynamical forcing associated with an arc do not significantly affect the intensity of the arc.

Lyons, L. R.↗

A computational study of the flowfield surrounding the Aeroassist Flight Experiment vehicle

A symmetric total variation diminishing (STVD) algorithm has been applied to the solution of the three-dimensional hypersonic flowfield surrounding the Aeroassist Flight Experiment (AFE) vehicle. Both perfect-gas and chemical nonequilibrium models have been used. The perfect-gas flows were computed at two different Reynolds numbers, including a flight trajectory point at maximum dynamic pressure, and on two different grids. Procedures for coupling the solution of the species continuity equations with the Navier-Stokes equations in the presence of chemical nonequilibrium are reviewed and tested on the forebody of the AFE and on the complete flowfield assuming noncatalytic wall and no species diffusion. Problems with the STVD algorithm unique to flows with variable thermodynamic properties (real gas) are identified and algorithm modifications are suggested. A potential heating problem caused by strong flow impingement on the nozzle lip in the near wake at 0-deg angle of attack has been identified.

Gnoffo, Peter A.↗

The normal modes of the thermosphere

The linearized momentum, energy, and continuity equations for the thermosphere can be reduced to a form that gives the vertical structure for each horizontal wave mode. The vertical structure equation can be described in terms of the normal modes, or eigenmodes, of the thermosphere. The latter are obtained by using a 27-layer model that includes a realistic temperature profile and the effects of the Lorentz force, viscosity, and heat conduction. The normal modes have one real eigenfrequency for every two complex conjugate eigenfrequency values. The real modes have a dominant rotational wind component and are nonpropagating. The complex modes have comparable divergent and rotational wind components. The complex eigenvalues give vertically propagating modes, primarily associated with the transient response to forcing, and are significantly affected by the dissipation in the upper E region and F region. Results show that the rotational wind component dominates in the steady state when the forcing is due to the two-cell convection pattern at high latitudes and that the normal modes explain the large shears and large winds speeds that are typically observed in the high-latitude E region. The vertical energy flux for the normal modes is also calculated. The results show that the flux is upward above 130 km but downward in the lower E region for the total solution. The downward energy flux is a contribution from the real eigenmode structure.

Larsen, M. F.↗

Application of a flux-split algorithm to chemically relaxing, hypervelocity blunt-body flows

Viscous, nonequilibrium, hypervelocity flow fields over two axisymmetric configurations are numerically simulated using a factored, implicit, flux-split algorithm. The governing gas-dynamic and species-continuity equations for laminar flow are presented. The gas-dynamics/nonequilibrium-chemistry coupling procedure is developed as part of the solution procedure and is described in detail. Numerical solutions are presented for hypervelocity flows over a hemisphere and over an axisymmetric aeroassisted orbital transfer vehicle using three different chemistry models. The gas models considered are those for an ideal gas, for a frozen gas, and for chemically relaxing air consisting of five species. The calculated results are compared with existing numerical solutions in the literature along the stagnation line of the hemisphere. The effects of free-stream Reynolds number on the nonequilibrium flow field are discussed.

Balakrishnan, A.↗

Efficient numerical techniques for complex fluid flows

The aim of the present research program is the development of more efficient and reliable calculation schemes for the coupled momentum and continuity equations. The resulting schemes would significantly reduce the expense of computing complex flows such as those in combustion chambers, gas turbines, and heat exchangers.

Patankar, Suhas V.↗

Mass loss

Observational evidence for mass loss from cool stars is reviewed. Spectra line profiles are used for the derivation of mass-loss rates with the aid of the equation of continuity. This equation implies steady mass loss with spherical symmetry. Data from binary stars, Mira variables, and red giants in globular clusters are examined. Silicate emission is discussed as a useful indicator of mass loss in the middle infrared spectra. The use of thermal millimeter-wave radiation, Very Large Array (VLA) measurement of radio emission, and OH/IR masers are discussed as a tool for mass loss measurement. Evidence for nonsteady mass loss is also reviewed.

Goldberg, Leo↗

Global atmospheric circulation statistics: Four year averages

Four year averages of the monthly mean global structure of the general circulation of the atmosphere are presented in the form of latitude-altitude, time-altitude, and time-latitude cross sections. The numerical values are given in tables. Basic parameters utilized include daily global maps of temperature and geopotential height for 18 pressure levels between 1000 and 0.4 mb for the period December 1, 1978 through November 30, 1982 supplied by NOAA/NMC. Geopotential heights and geostrophic winds are constructed using hydrostatic and geostrophic formulae. Meridional and vertical velocities are calculated using thermodynamic and continuity equations. Fields presented in this report are zonally averaged temperature, zonal, meridional, and vertical winds, and amplitude of the planetary waves in geopotential height with zonal wave numbers 1-3. The northward fluxes of sensible heat and eastward momentum by the standing and transient eddies along with their wavenumber decomposition and Eliassen-Palm flux propagation vectors and divergences by the standing and transient eddies along with their wavenumber decomposition are also given. Large interhemispheric differences and year-to-year variations are found to originate in the changes in the planetary wave activity.

Wu, M. F.↗

A simple, analytic 3-dimensional downburst model based on boundary layer stagnation flow

A simple downburst model is developed for use in batch and real-time piloted simulation studies of guidance strategies for terminal area transport aircraft operations in wind shear conditions. The model represents an axisymmetric stagnation point flow, based on velocity profiles from the Terminal Area Simulation System (TASS) model developed by Proctor and satisfies the mass continuity equation in cylindrical coordinates. Altitude dependence, including boundary layer effects near the ground, closely matches real-world measurements, as do the increase, peak, and decay of outflow and downflow with increasing distance from the downburst center. Equations for horizontal and vertical winds were derived, and found to be infinitely differentiable, with no singular points existent in the flow field. In addition, a simple relationship exists among the ratio of maximum horizontal to vertical velocities, the downdraft radius, depth of outflow, and altitude of maximum outflow. In use, a microburst can be modeled by specifying four characteristic parameters, velocity components in the x, y and z directions, and the corresponding nine partial derivatives are obtained easily from the velocity equations.

Oseguera, Rosa M.↗

Use of satellite radiometric imagery data for improvement in the analysis of divergent wind in the tropics

A scheme is proposed to incorporate satellite radiometric imagery data into the specification of initial conditions for the National Meteorological Center (NMC) operational global prediction model in order to improve the analysis of the divergent wind field in the tropics. The basic assumptions are that outgoing longwave radiation (OLR) data can provide 1) the division between convective (upward motion) and clear sky (downward motion) areas, and 2) the height of convection cells. The intensity of ascending motion in the convective areas is estimated based on OLR data. The intensity of descending motion is evaluated from the thermodynamic energy balance between radiative cooling and adiabatic warming, since the local time change of temperature is small in the tropics. Once the vertical motion field is determined, the horizontal divergence field can be calculated from the mass continuity equation. Then a divergent wind field is determined. The total wind is the sum of the new divergent wind and the rotational part, which is assumed to be unchanged. The proposed scheme is tested using the NMC analysis dataset of 21 January 1985 with satisfactory results.

Kasahara, Akira↗

An upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow around two-dimensional or axisymmetric bodies. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. The algorithm solves the gas dynamic and species continuity equations in a 'loosely' coupled manner. The new code has been validated by computing the laminar flow (at free stream Mach number 25) of chemically reacting air over a wedge and a cone. The results of these computations are compared with the results from a centrally-differenced, fully coupled, nonequilibrium PNS code. The agreement is excellent, except in the vicinity of the shock wave where the present code exhibits superior shock capturing capabilities.

Tannehill, John C.↗

Efficient numerical techniques for complex fluid flows

The central feature in any flow prediction method is the treatment of the coupling between the momentum and continuity equations. In natural-convection flows, the energy equation also becomes strongly coupled with the momentum equations. Because of the nonlinear nature of the coupling, these equations are solved iteratively. Iterative methods are often prone to slow convergence, divergence, and extreme sensitivity to underrelaxation factors. The aim of the present research is to develop more efficient and reliable solution schemes for the coupled flow equations. Such schemes will significantly reduce the expense of computing complex flows encountered in combustion chambers, gas turbines, heat exchangers, and other practical equipment. In the work completed so far, a technique employing norm reduction in conjunction with the successive-substitution and Newton-Raphson techniques was developed. Also, a block-correction procedure for the flow equations is currently being formulated and tested.

Patankar, Suhas V.↗

The effect of exhaust plume/afterbody interaction on installed Scramjet performance

Newly emerging aerospace technology points to the feasibility of sustained hypersonic flight. Designing a propulsion system capable of generating the necessary thrust is now the major obstacle. First-generation vehicles will be driven by air-breathing scramjet (supersonic combustion ramjet) engines. Because of engine size limitations, the exhaust gas leaving the nozzle will be highly underexpanded. Consequently, a significant amount of thrust and lift can be extracted by allowing the exhaust gases to expand along the underbody of the vehicle. Predicting how these forces influence overall vehicle thrust, lift, and moment is essential to a successful design. This work represents an important first step toward that objective. The UWIN code, an upwind, implicit Navier-Stokes computer program, has been applied to hypersonic exhaust plume/afterbody flow fields. The capability to solve entire vehicle geometries at hypersonic speeds, including an interacting exhaust plume, has been demonstrated for the first time. Comparison of the numerical results with available experimental data shows good agreement in all cases investigated. For moderately underexpanded jets, afterbody forces were found to vary linearly with the nozzle exit pressure, and increasing the exit pressure produced additional nose-down pitching moment. Coupling a species continuity equation to the UWIN code enabled calculations indicating that exhaust gases with low isentropic exponents (gamma) contribute larger afterbody forces than high-gamma exhaust gases. Moderately underexpanded jets, which remain attached to unswept afterbodies, underwent streamwise separation on upswept afterbodies. Highly underexpanded jets produced altogether different flow patterns, however. The highly underexpanded jet creates a strong plume shock, and the interaction of this shock with the afterbody was found to produce complicated patterns of crossflow separation. Finally, the effect of thrust vectoring on vehicle balance has been shown to alter dramatically the vehicle pitching moment.

Edwards, Thomas Alan↗

A multidimensional model for aerosols - Description of computational analogs

The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.

Toon, O. B.↗

Chemical nonequilibrium and viscous flow computation for conic aerobrake bodies

Three-dimensional analyses are presented for the viscous, reactive flow over a complete entry-body configuration with a wide-angle conic surface. The predictive method uses a split approach that solves iteratively the Navier-Stokes and the continuity equations of chemical species. The finite-difference formulation and the computational grid are adapted to the bow shock and the conformally mapped body such that the velocity components are in the computational spherical-polar space. Combinations of several conic forebody and afterbody configurations have been studied using wind-tunnel, Space Shuttle, and aerobraking orbital transfer vehicle (AOTV) entry conditions. The effects of the borebody bluntness and of finite-rate chemical reactions on the shock layer, the wall catalycity on the boundary layer, the shear-layer impingement on the afterbody, and the base-flow environment are discussed.

Li, C. P.↗