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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Seasonal and magnetic storm related changes in the thermosphere induced by eddy mixing

The possibility of explaining the seasonal and magnetic storm related changes in the thermosphere and the ionosphere through variation in eddy diffusion coefficient in the lower thermosphere is investigated theoretically by obtaining simultaneous numerical solutions of the relevant continuity, momentum, and energy balance equations in a self-consistent manner. It is shown that various important features of thermospheric seasonal behavior, including the winter helium bulge phenomenon, can be explained by assuming effective reduction in the winter-time eddy diffusion coefficient by about one order of magnitude. Storm related changes in the thermal structure and composition are seen to arise by assuming a downward shift of turbulence. Physical implications and relative significance of variation in thermospheric turbulence and large scale inter-hemispheric circulation are briefly discussed.

Sinha, A. K.↗

NASA NTRS1977

The influence of convection electric fields on thermal proton outflow from the ionosphere

The continuity, momentum and energy hydrodynamic equations for an O(+)-H(+) ionosphere have been solved self-consistently for steady state conditions when a perpendicular (convection) electric field is present. Comparison of the H(+) temperature profiles obtained with and without the electric field show that the effect of the electric field is to enhance the H(+) temperature at high altitudes from about 3600 to 6400 K. Due to ion heating by the electric field, there is a net reduction of O(+) in the F2-region as compared with the case of a nonconvecting ionosphere. When the reduction of O(+) is neglected, the electric field acts to increase the H(+) outward flux. However, when the reduction of O(+) is included, there is a net reduction in the outward H(+) flux. Nevertheless, the convection electric field still results in an increase in the rate of depletion of the F-region ionization due to H(+) outflow, by a factor of 2.2 for a 100 mV/m electric field.

Raitt, W. J.↗

NASA NTRS1978

A comprehensive model of the Venus ionosphere

Coupled time-dependent continuity-momentum and energy balance equations for the Venus ionosphere were simultaneously solved for CO2(+), O2(+), O(+), He(+), and H(+) densities and electron and ion temperatures for an altitude range of 120-500 km. Values of the solar zenith angle varied from 0 deg (subsolar point) to 90 deg (terminator). The calculations include the horizontal bulk transport of ions by neutral winds but not the horizontal diffusion. The two-stream photoelectron transport method was used to find the heating rates for the ambient electrons. Different boundary conditions were considered, and a nightside ionosphere was calculated. The results of these model calculations are in good agreement with measurements in the region of maximum electron density. Characteristics of topside and nightside densities and temperatures are discussed.

Chen, R. H.↗

NASA NTRS1979

Numerical analysis of turbulent wakes of axial flow turbomachinery rotor blades

The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.

Hah, C.↗

NASA NTRS1979

Equatorial wave-mean flow interaction - A numerical study of the role of latitudinal shear

A time-dependent primitive equation model for an equatorial channel is used to assess the interaction of equatorial Kelvin and mixed Rossby-gravity waves with the mean flow. The proposed model involves a semiimplicit time-differencing scheme and a finite-difference grid in the meridional plane. It is shown that forced equatorial waves interact with mean flow to produce equatorial jets characterized by downward-moving westerly (Kelvin wave forcing) and easterly (mixed Rossby-gravity wave forcing) shear zones, respectively. For parameters characteristic of the observed waves in the equatorial stratosphere, the wave-mean flow interaction process always reduces the amplitude of any initial cross-equatorial mean wind shear. The mean flow profile tends to become symmetric about the equator as the interaction process continues.

Holton, J. R.↗

NASA NTRS1980

A simulation of the atmospheric cloud physics laboratory to aid in its design and the design of the experiments within the laboratory

Using the finite difference method with overrelaxation, numerical solutions of the steady-state vorticity transport equation were obtained for a continuous flow diffusion chamber of the Hudson-Squires type. The calculation neglected the effects due to temperature, gravity, and saturation. The size and shape of the manifold used to inject the aerosol laden flow were varied to obtain a design which would improve the performance of the chamber from strictly low Reynolds number (less than 20) fluid dynamical considerations.

Winchester, L. W., Jr.↗

NASA NTRS1982

Numerical analysis and FORTRAN program for the computation of the turbulent wakes of turbomachinery rotor blades, isolated airfoils and cascade of airfoils

Turbulent wakes of turbomachinery rotor blades, isolated airfoils, and a cascade of airfoils were investigated both numerically and experimentally. Low subsonic and incompressible wake flows were examined. A finite difference procedure was employed in the numerical analysis utilizing the continuity, momentum, and turbulence closure equations in the rotating, curvilinear, and nonorthogonal coordinate system. A nonorthogonal curvilinear coordinate system was developed to improve the accuracy and efficiency of the numerical calculation. Three turbulence models were employed to obtain closure of the governing equations. The first model was comprised to transport equations for the turbulent kinetic energy and the rate of energy dissipation, and the second and third models were comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, respectively. The second model handles the convection and diffusion terms in the Reynolds stress transport equation collectively, while the third model handles them individually. The numerical results demonstrate that the second and third models provide accurate predictions, but the computer time and memory storage can be considerably saved with the second model.

Hah, C.↗

NASA NTRS1983

MHD stability of incompressible coronal loops with radiative energy loss

Previous studies of the magnetohydrodynamic (MHD) stability of solar coronal loops have not taken into account the effects of radiative energy loss in the energy equation. However, since coronal loops continuously lose energy by radiation and heat conduction, it is important to understand how these energy loss mechanisms affect MHD stability. We investigate the problem assuming that a magnetic loop has cylindrical geometry. As a first step, stability is studied for a localized mode, and the result is applied to a specific equilibrium. We find that the radiative energy loss effect not only changes the growth rate of ideally unstable modes, but also alters the stability boundary predicted by ideal MHD theory.

An, C.-H.↗

NASA NTRS1985

Noniterative cross-flow integration for the pressure-split analysis of subsonic mixing-layer problems

In utilizing pressure-split methodology for the two-dimensional analysis of curved wall jets and the subsonic regions of underexpanded free jets, Dash et al (1984) employed the cross flow procedure of Bradshaw et al (1979, 1984). In assessing this procedure, it was noted that the iterative steps required between the continuity and cross flow momentum equation solutions could be eliminated by combining these equations into a unified equation for the cross flow velocity, yielding a significant reduction in computer time and eliminating possible convergence problems which occur in iterative approaches. Attention is presently given to the application of this noniterative procedure to a two-dimensional curved wall jet problem.

Dash, S. M.↗

NASA NTRS1987

A one-dimensional time-dependent model of the magnetized ionosphere of Venus

The behavior and time evolution of the large-scale magnetic fields and plasma of the dayside Venus ionosphere are studied using a one-dimensional model. The coupled continuity, momentum, and Maxwell's equations are solved simultaneously for O(+), O2(+), and H(+), and the magnetic field. The calculated magnetic field profiles are in good agreement with Pioneer Venus orbiter magnetometer observations. The magnetic field structure is quasi-steady for slow changes of the solar wind dynamic pressure. The peak at 165 km is maintained by downward convection from higher altitudes. The time scale for the decay of the field by the pure one-dimensional vertical diffusion/convection process is several hours unless the flux is resupplied from the top of the ionosphere.

Shinagawa, H.↗

NASA NTRS1989

Analysis of regularized Navier-Stokes equations, 2

A practically important regularization of the Navier-Stokes equations was analyzed. As a continuation of the previous work, the structure of the attractors characterizing the solutins was studied. Local as well as global invariant manifolds were found. Regularity properties of these manifolds are analyzed.

Ou, Yuh-Roung↗

NASA NTRS1951

Behavior of fast moving flow of compressible gas in cylindrical pipe in presence of cooling

For compressible flow with friction in a cylindrical pipe the momentum, continuity, and heat-transfer equations are examined to determine whether an increase in Mach number ("thermal" Laval nozzle) is obtainable through heat conduction from the gas through the pipe walls. The analysis is based on the assumption that the wall temperature is negligibly small in comparison with the stagnation temperature of the gas. The analysis leads to a negative result. When the gas cooling is increased by also considering radiation to the wall, a limited region at high temperatures is obtained where Mach number increases were theoretically possible. Obtaining this condition practically is considered impossible.

FLOW, SUPERSONIC↗

NASA NTRS1993

The nonlinear evolution of modes on unstable stratified shear layers

The nonlinear development of disturbances in stratified shear flows (having a local Richardson number of value less than one quarter) is considered. Such modes are initially fast growing but, like related studies, we assume that the viscous, non-parallel spreading of the shear layer results in them evolving in a linear fashion until they reach a position where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. Four different basic integro-differential amplitude equations are possible, including one due to a novel mechanism; the relevant choice of amplitude equation, at a particular instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate of the disturbance, its wavenumber, and the viscosity of the fluid. This richness of choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential amplitude equations with a cubic nonlinearity but the resulting significant increase in the size of the disturbance's amplitude leads on to the next stage of the evolution process where the evolution of the mode is governed by an integro-differential amplitude equations with a quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which time the flow has become fully nonlinear.

Blackaby, Nicholas↗

NASA NTRS1993

Control of nonlinear systems represented in quasilinear form

Methods to synthesize controllers for nonlinear systems are developed by exploiting the fact that under mild differentiability conditions, systems of the form: x-dot = f(x) + G(x)u can be represented in quasilinear form, viz: x-dot = A(x)x + B(x)u. Two classes of control methods are investigated. The first is zero-look-ahead control, where the control input depends only on the current values of A(x) and B(x). For this case the control input is computed by continuously solving a matrix Riccati equation as the system progresses along a trajectory. The second is controllers with look-ahead, where the control input depends on the future behavior of A(x) and B(x). These controllers use the similarity between quasilinear systems and linear time varying systems to find approximate solutions to optimal control type problems. The methods that are developed are not guaranteed to be globally stable. However in simulation studies they were found to be useful alternatives for synthesizing control laws for a general class of nonlinear systems.

Coetsee, Josef A.↗

NASA NTRS1995

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

NASA NTRS1995

Three-fluid solar wind model with Alfven waves

We present a study of a three-fluid solar wind model. with continuity, momentum and separate energy equations for protons. alpha particles and electrons. Allowing separate coronal heat sources for all three species, we study the flow properties of the solar wind as a function of heat input, Alfven wave energy input, and alpha particle abundance.

Esser, Ruth↗

NASA NTRS1995

Conservative properties of finite difference schemes for incompressible flow

The purpose of this research is to construct accurate finite difference schemes for incompressible unsteady flow simulations such as LES (large-eddy simulation) or DNS (direct numerical simulation). In this report, conservation properties of the continuity, momentum, and kinetic energy equations for incompressible flow are specified as analytical requirements for a proper set of discretized equations. Existing finite difference schemes in staggered grid systems are checked for satisfaction of the requirements. Proper higher order accurate finite difference schemes in a staggered grid system are then proposed. Plane channel flow is simulated using the proposed fourth order accurate finite difference scheme and the results compared with those of the second order accurate Harlow and Welch algorithm.

Morinishi, Youhei↗

NASA NTRS1996

Analysis of Two-Phase Flow in Damper Seals for Cryogenic Turbopumps

Cryogenic damper seals operating close to the liquid-vapor region (near the critical point or slightly su-cooled) are likely to present two-phase flow conditions. Under single phase flow conditions the mechanical energy conveyed to the fluid increases its temperature and causes a phase change when the fluid temperature reaches the saturation value. A bulk-flow analysis for the prediction of the dynamic force response of damper seals operating under two-phase conditions is presented as: all-liquid, liquid-vapor, and all-vapor, i.e. a 'continuous vaporization' model. The two phase region is considered as a homogeneous saturated mixture in thermodynamic equilibrium. Th flow in each region is described by continuity, momentum and energy transport equations. The interdependency of fluid temperatures and pressure in the two-phase region (saturated mixture) does not allow the use of an energy equation in terms of fluid temperature. Instead, the energy transport is expressed in terms of fluid enthalpy. Temperature in the single phase regions, or mixture composition in the two phase region are determined based on the fluid enthalpy. The flow is also regarded as adiabatic since the large axial velocities typical of the seal application determine small levels of heat conduction to the walls as compared to the heat carried by fluid advection. Static and dynamic force characteristics for the seal are obtained from a perturbation analysis of the governing equations. The solution expressed in terms of zeroth and first order fields provide the static (leakage, torque, velocity, pressure, temperature, and mixture composition fields) and dynamic (rotordynamic force coefficients) seal parameters. Theoretical predictions show good agreement with experimental leakage pressure profiles, available from a Nitrogen at cryogenic temperatures. Force coefficient predictions for two phase flow conditions show significant fluid compressibility effects, particularly for mixtures with low mass content of vapor. Under these conditions, an increase on direct stiffness and reduction of whirl frequency ratio are shown to occur. Prediction of such important effects will motivate experimental studies as well as a more judicious selection of the operating conditions for seals used in cryogenic turbomachinery.

Arauz, Grigory L.↗