Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Influence of land-surface evapotranspiration on the earth's climate

Land-surface evapotranspiration is shown to strongly influence global fields of rainfall, temperature and motion by calculations using a numerical model of the atmosphere, confirming the general belief in the importance of evapotranspiration-producing surface vegetation for the earth's climate. The current version of the Goddard Laboratory atmospheric general circulation model is used in the present experiment, in which conservation equations for mass, momentum, moisture and energy are expressed in finite-difference form for a spherical grid to calculate (1) surface pressure field evolution, and (2) the wind, temperature, and water vapor fields at nine levels between the surface and a 20 km height.

Shukla, J.↗

Impact mechanics of the Cretaceous-Tertiary extinction bolide

An examination of the mechanics of asteroidal, cometary, and meteor swarm impact on the earth determined if the enrichment of projectile material in the K-T layer is consistent with melts and impact breccias on the earth and moon, the size of the impacters, the distribution of the kinetic energy, and the sequence of impacts that could give rise to observed extinction phenomena. Flows resulting from spherical projectile impacts onto layers of air, water, and silicates were modeled and Eulerian finite difference algorithms were employed to solve conservation equations and equations of state. A range of speeds and impacter densities were considered, along with sizes from 0.17 km, which would be consumed in the atmosphere, to a 10 km object, which would have had a diameter greater than a reference 7.1 km atmosphere depth. It is concluded that an impact of the K-T bolide could result in global biotic extinction and worldwide material deposition.

Okeefe, J. D.↗

Cloud/climate sensitivity experiments

A study of the relationships between large-scale cloud fields and large scale circulation patterns is presented. The basic tool is a multi-level numerical model comprising conservation equations for temperature, water vapor and cloud water and appropriate parameterizations for evaporation, condensation, precipitation and radiative feedbacks. Incorporating an equation for cloud water in a large-scale model is somewhat novel and allows the formation and advection of clouds to be treated explicitly. The model is run on a two-dimensional, vertical-horizontal grid with constant winds. It is shown that cloud cover increases with decreased eddy vertical velocity, decreased horizontal advection, decreased atmospheric temperature, increased surface temperature, and decreased precipitation efficiency. The cloud field is found to be well correlated with the relative humidity field except at the highest levels. When radiative feedbacks are incorporated and the temperature increased by increasing CO2 content, cloud amounts decrease at upper-levels or equivalently cloud top height falls. This reduces the temperature response, especially at upper levels, compared with an experiment in which cloud cover is fixed.

Roads, J. O.↗

Numerical viscosity and the entropy condition for conservative difference schemes

Consider a scalar, nonlinear conservative difference scheme satisfying the entropy condition. It is shown that difference schemes containing more numerical viscosity will necessarily converge to the unique, physically relevant weak solution of the approximated conservation equation. In particular, entropy satisfying convergence follows for E schemes - those containing more numerical viscosity than Godunov's scheme.

Tadmor, E.↗

Small particle transport across turbulent nonisothermal boundary layers

The interaction between turbulent diffusion, Brownian diffusion, and particle thermophoresis in the limit of vanishing particle inertial effects is quantitatively modeled for applications in gas turbines. The model is initiated with consideration of the particle phase mass conservation equation for a two-dimensional boundary layer, including the thermophoretic flux term directed toward the cold wall. A formalism of a turbulent flow near a flat plate in a heat transfer problem is adopted, and variable property effects are neglected. Attention is given to the limit of very large Schmidt numbers and the particle concentration depletion outside of the Brownian sublayer. It is concluded that, in the parameter range of interest, thermophoresis augments the high Schmidt number mass-transfer coefficient by a factor equal to the product of the outer sink and the thermophoretic suction.

Rosner, D. E.↗

Ignition of confined gaseous mixtures by hot surfaces and hot wires

Ignition times and spatial and temporal variations of temperature and concentration in gaseous mixtures confined between two infinite parallel walls or two infinite cylinders have been obtained by numerical integration of the appropriate conservation equations written in Lagrangian coordinates. Ignition times and ignition energies are presented for the case of an isothermal wall in terms of the initial mixture pressure and equivalence ratio for both one and two-step chemical reaction mechanisms. The numerical results indicate that there is a critical mixture pressure for which the ignition time is minimum. The values of this critical pressure are larger (smaller) than 1 atm for the one- (two-) step reaction mechanism. The critical pressure for the ignition time is not equal to the critical pressure for the ignition energy. The ignition time and energy decrease with the equivalence ratio within a certain range and then remain constant.

Ramos, J. I.↗

Transient flow analysis of the AEDC/HPDE MHD generator

A hybrid Lax-Wendroff/Method of Characteristics computer code has been developed for numerical simulation of flow transients associated with the operation of MHD generator facilities. The code employs the shock-fitting method, with an Eulerian formulation of the basic conservation equations and explicit tracking of shock waves. Pressure, temperature, and velocity are used as primary integration variables to simplify interfacing of the code with real-gas thermodynamic and transport property tables. Application of the code to the simulation of selected transients for the AEDC/HPDE MHD generator produced results that are in good agreement with experimental observations.

Wilson, D. R.↗

Numerical studies of the formation and destruction of vortices in a motored four-stroke piston-cylinder configuration

A finite-difference procedure which solves the conservation equations of mass, momentum, and energy is used to investigate the effects of the compression ratio, engine speed, bore-to-stroke ratio, and air intake flow angle on the turbulent flow field within an axisymmetric piston-cylinder configuration. It is shown that in a four-stroke piston-cylinder configuration, the intake stroke is characterized by the formation of a piston vortex. The piston vortex is stretched during the intake stroke, and the head vortex has an almost constant diameter. For a 0-deg air intake flow angle, both vortices disappear by the end of the compression stroke; for an air intake flow angle of 45 deg, the flow field within the cylinder shows three elongated vortices which persist into the compression stroke and then break up and merge. It is also shown that larger bore-to-stroke ratios give rise to lower turbulent levels than smaller bore-to-stroke ratios and that the turbulent intensity is almost independent of the rpm.

Schock, H. J.↗

A numerical method based on the Fourier-Fourier transform approach for modeling 1-D electron plasma evolution

A numerical method is presented for studying one-dimensional electron plasma evolution under typical interplanetary conditions. The method applies the Fourier-Fourier transform approach to a plasma model that is a generalization of the electrostatic Vlasov-Poisson system of equations. Conservation laws that are modified to include the plasma model generalization and also the boundary effects of nonperiodic solutions are given. A new conservation law for entropy in the transformed space is then introduced. These conservation laws are used to verify the numerical solutions. A discretization error analysis is presented. Two numerical instabilities and the methods used for their suppression are treated. It is shown that in interplanetary plasma conditions, the bump-on-tail instability produces significant excitation of plasma oscillations at the Bohm-Gross frequency and its second harmonic. An explanation of the second harmonic excitation is given in terms of wave-wave coupling during the growth phase of the instability.

Klimas, A. J.↗

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗

Steam chugging analysis in single-vent vapor injection

A complete cycle of the periodic steam chugging phenomenon is analyed. Steam velocity and pressure variations in the vent are described by one-dimensional conservation equations. This is coupled either to the water slug model when water is in the vent, or, the infinite pool spherical bubble model at the vent exit during bubble growth. An isolated spherical bubble model is used for computing the collapse pressures. Comparisons of the model predictions with the UCLA 1/12-scale and the Japan 1/6-scale data indicate that the vent-pipe model predicts the vent-clearing times and the bubble growth times well. In addition, the predicted maximum chugging heights compared well with those measured in the Japan data. On bubble collapse pressures, the comparison with the spherical bubble model predictions is only fair. The model generally overpredicts the magnitude of the spikes. On examining the effects of pool subcooling and steam mass flux, general agreement is found between the predicted trends and those measured.

Lee, C. K. B.↗

A computational method for viscous incompressible flows

An implicit, finite-difference procedure for numerically solving viscous incompressible flows is presented. The pressure-field solution is based on the pseudocompressibility method in which a time-derivative pressure term is introduced into the mass-conservation equation to form a set of hyperbolic equations. The pressure-wave propagation and the spreading of the viscous effect is investigated using simple test problems. Computed results for external and internal flows are presented to verify the present method which has proved to be very robust in simulating incompressible flows.

Kwak, D.↗

Dilution jets in accelerated cross flows

Results of flow visualization experiments and measurements of the temperature field produced by a single jet and a row of dilution jets issued into a reverse flow combustor are presented. The flow in such combustors is typified by transverse and longitudinal acceleration during the passage through its bending section. The flow visualization experiments are designed to examine the separate effects of longitudinal and transverse acceleration on the jet trajectory and spreading rate. A model describing a dense single jet in a lighter accelerating cross flow is developed. The model is based on integral conservation equations, including the pressure terms appropriate to accelerating flows. It uses a modified entrainment correlation obtained from previous experiments of a jet in a cross stream. The flow visualization results are compared with the model calculations in terms of trajectories and spreading rates. Each experiment is typified by a set of three parameters: momentum ratio, density ratio and the densimetric Froude number.

Lipshitz, A.↗

Embedded shear layer computations for increased drag reduction

One of the most promising methods of minimizing drag is the reduction of skin friction by injection of low momentum fluid into the near-wall region of turbulent boundary layer flows. This method could be made more effective by limiting the spread rate of the resulting mixing region. In order to achieve a better understanding of how this goal might be achieved, numerical investigations of the relevant fluid dynamic processes governing these regions have been conducted. A compact finite-difference algorithm has been applied to the complete form of the governing conservation equations for a two-dimensional laminar mixing layer. The ability of this computational approach to model successfully the formation and interaction of the large scale vortical structures which dominate such flow fields is verified in the present study. Parameters which affect the spread rate of the mixing region are also identified. In addition, the relative importance of viscous and momentum transport effects in the vortex interactions is determined.

Gatski, T. B.↗

On the theory of cosmic-ray-mediated shocks with variable compression ratio

Cosmic-ray-mediated shocks may accelerate enough cosmic rays to high enough energies that they escape the shock, carrying an appreciable amount of energy before being convected to downstream infinity. Under such conditions, it is noted, the overall compression ratio cannot be determined from the conservation equations as in conventional hydrodynamic treatments, and the standard equations for shock acceleration admit arbitrarily high compression ratios. A procedure is outlined for obtaining the structure of high Mach number, cosmic-ray-mediated shocks, including their overall compresion ratio, around a low Mach number viscous subshock. Analytic solutions are obtained by quardrature for an energy-dependent diffusion coefficient in the limit of extreme sensitivity to energy, which, unlike previous solutions, include the finite thermal pressure of the preshock gas.

Eichler, D.↗

The Scaling of Coronal Models from One Star to Another

The requirements that must be met in order that stationary numerical corona models can be scaled from one star to another are discussed. A corona model is a solution of the conservation equations for mass, momentum, and energy, subject to appropriate boundary conditions, and of the equation of state. In general, the mass M and radius R of the star enter these equations and boundary conditions as free parameters. A given solution can be scaled to other stars only if all equations can be rewritten in such a form that M and R do no longer appear explicitly as free parameters, but only implicitly as scaling factors of the variables. An adequate means to find these scaling factors is a homologous transformation: one multiplies all variables and parameters by separate constants (i.e., scaling factors) and requires that the equations and boundary conditions remain valid. This leads to a set of nonlinear relations between the transformation constants. Only if in this set the two constants associated with M and R can be chosen independently, can a given numerical corona model be scaled to arbitrary stars.

Hammer, R.↗

Numerical viscosity and the entropy condition for conservative difference schemes

Consider a scalar, nonlinear conservative difference scheme satisfying the entropy condition. It is shown that difference schemes containing more numerical viscosity will necessarily converge to the unique, physically relevant weak solution of the approximated conservation equation. In particular, entropy satisfying convergence follows for E schemes - those containing more numerical viscosity than Godunov's scheme.

Tadmor, E.↗

On convergence of computation of chemically reacting flows

The computational problems associated with high-temperature flows undergoing finite-rate ionization reactions is investigated. The conservation equations governing chemical species and vibrational and electron energies are solved simultaneously with those for overall mass, momentum, and energy for a one-dimensional subsonic flow, through a constant-area duct, originating behind a normal shock wave, using an implicit time-marching technique. Boundary conditions are imposed in the form of characteristic wave variables accounting for the effects of chemical reactions on the speed of sound. Converging solutions are obtained for cases in which chemical reactions are weak, but difficulty is encountered in other cases. The cause of the difficulty is investigated and shown to be the sharp pressure disturbances produced by such reactions.

Park, C.↗