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 505 records · Page 28

Flare loop radiative hydrodynamics. I - Basic methods

The study presented here has two goals. The first is to examine the role of energy transport mechanisms in flare dynamics, for instance, thermal cone conduction in the corona and radiative transfer in the chromosphere. The second is to provide diagnostics of flare energization processes by predicting the evolutionary behavior of concurrent emissions from all regions of the flaring atmosphere. To attain these ends, a numerical method is developed for the simultaneous solution of thy continuity, momentum, and energy equations; the time-dependent atomic rate equations describing ionization and excitation; and the radiative transfer equations. With the physical model and computational methods used here, all components of the plasma move as a single fluid and in the initial studies are assumed to have a common kinetic temperature. The plasma is constrained to move in one dimension along the axis of a loop of nonuniform cross section, the geometry of which is defined by a sufficiently strong magnetic field.

Mcclymont, A. N.↗

Large-scale counterstreaming of H(+) and He(+) along plasmaspheric flux tubes

An interhemispheric plasma transport model is used to study the flow characteristics of H(+), He(+), and O(+) along closed geomagnetic field lines for solstice conditions. The model corresponds to a time-dependent solution of the coupled continuity, momentum, and energy equations for the ions and electrons. The equations are solved along an entire flux tube from 120 km in one hemisphere to 120 km in the other hemisphere. The calculations are carried out for noon conditions for the flux tube passsing through Millstone Hill, at L = 3.2. The main conclusion is that H(+)-He(+) counterstreaming can be expected along a large segment of a plasmaspheric flux tube at solstice. For both symmetric and asymmetric wind patterns, the He(+) flow is from the winter to the summer ionosphere not only in the steady state, but during flux tube refilling due to the winter helium bulge and the depletion of N2.

Richards, P. G.↗

Elasticity Solution of an Adhesively Bonded Cover Plate of Various Geometries

The plane strain of adhesively bonded structures consisting of two different isotropic adherends is considered. By expressing the x-y components of the displacements in terms of Fourier integrals and using the corresponding boundary and continuity conditions, the integral equations for the general problem are obtained and solved numerically by applying Gauss-Chebyshev integration scheme. The shear and the normal stresses in the adhesive are calculated for various geometries and material properties for a stiffened plate under uniaxial tension. Numerical results involving the stress intensity factors and the strain energy release rate are presented. The closed-form expressions for the Fredholm kernels are provided to obtain the solution for an arbitrary geometry and material properties. For the general geometry, the contribution of the normal stress is quite significant, while for symmetric geometries, the shear stress is dominant, the normal stress vanishes if the adherends are of the same material and the same thickness.

Aksel, G. N.↗

A time-dependent theoretical model of the polar wind Preliminary results

The coupled, time-dependent continuity, momentum, and energy equations of a two-ion O(+) and H(+) quasi-neutral plasma were solved in order to extend understanding of polar wind behavior. This numerical code allows studies of the time dependent behavior of polar wind-type flows into and out of the ionosphere. Initial studies indicate that the typical time constants for electron and ion temperature changes are of the order of minutes and tens of minutes, respectively. The response time of the minor high altitude ion O(+) is less than an hour, whereas that of the major ion, H(+), is many hours. The initial test runs also demonstrate the fact that temporary supersonic flows of both O(+) and H(+) are possible, especially in the presence of significant ion heating.

Gombosi, T. I.↗

Improved algorithms for circulation-control airfoils in transonic flow

A zonal model for aerodynamic analysis of two-dimensional transonic circulation control airfoils has been developed. The present approach combines a transonic full potential method for the global flow field and an integral boundary layer method for regions of the airfoil excluding the wall jet with a parabolized Navier-Stokes code for resolving the wall jet region beyond the slot. Existing methods suffer from two deficiencies: the insensitivity of the calculation to small changes in the Coanda surface geometry; and the inability to predict the shock structure of the underexpanded supersonic wall jets. The present wall jet procedure involves a pressure-split approach in the streamwise sense to enable noniterative solution of the coupled continuity and normal momentum equations for increased surface sensitivity and allows for expansion of applications to sonic slot exit conditions. Encouraging results are obtained in comparison with experimental data for two circulation airfoils with subsonic wall jets.

Dvorak, Frank A.↗

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗

Time-dependent numerical simulation of hot ion outflow from the polar ionosphere

The time-dependent hydrodynamic model of Gombosi et al. (1985) was used to investigate the effects of bulk plasma heating on ion flows in the polar ionosphere. In the analysis, coupled time-dependent continuity, momentum, and energy equations of a two-ion (H/+/ and O/+/) quasi-neutral plasma were solved between 200 and 8000 km for polar wind conditions. In spite of its several limitations, the model can serve to indicate the gross behavior of polar region plasma flows. Two specific cases, ion heating and electron heating, were considered. The calculated temperature and flow profiles suggest that ion heating, rather than electron heating, generates the upwelling ion events, for which large transient O(+) outflows accompanied by elevated ion temperatures were observed.

Gombosi, T. I.↗

A new model of cometary ionospheres

The coupled continuity, momentum, and energy equations were solved for ionospheric conditions appropriate for Comet Halley at 1 AU. The numerical scheme used is such that any shock transition appears naturally in the solution and no a priori assumptions are necessary. Solutions were obtained for a number of different assumptions concerning electron heating rates, but all showed that the electron temperatures increase rapidly and significantly at a distance from the nucleus where collisional electron-neutral cooling becomes unimportant. This temperature increase is accompanied by a sharp increase in both the plasma pressure and its associated polarization electric field, causing the supersonic plasma flow to go subsonic. It is not clear at this time whether or not this sonic transition is accompanied by a shock.

Korosmezey, A.↗

Formulation of a General Technique for Predicting Pneumatic Attenuation Errors in Airborne Pressure Sensing Devices

Presented is a mathematical model derived from the Navier-Stokes equations of momentum and continuity, which may be accurately used to predict the behavior of conventionally mounted pneumatic sensing systems subject to arbitrary pressure inputs. Numerical techniques for solving the general model are developed. Both step and frequency response lab tests were performed. These data are compared with solutions of the mathematical model and show excellent agreement. The procedures used to obtain the lab data are described. In-flight step and frequency response data were obtained. Comparisons with numerical solutions of the math model show good agreement. Procedures used to obtain the flight data are described. Difficulties encountered with obtaining the flight data are discussed.

Whitmore, Stephen A.↗

Formulaton of a general technique for predicting pneumatic attenuation errors in airborne pressure sensing devices

Presented is a mathematical model, derived from the Navier-Stokes equations of momentum and continuity, which may be accurately used to predict the behavior of conventionally mounted pneumatic sensing systems subject to arbitrary pressure inputs. Numerical techniques for solving the general model are developed. Both step and frequency response lab tests were performed. These data are compared against solutions of the mathematical model. The comparisons show excellent agreement. The procedures used to obtain the lab data are described. In-flight step and frequency response data were obtained. Comparisons with numerical solutions of the mathematical model show good agreement. Procedures used to obtain the flight data are described. Difficulties encountered with obtaining the flight data are discussed.

Whitmore, Stephen A.↗

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

A mathematical model of the middle and high latitude ionosphere

A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.

Schunk, R. W.↗

Compensating for pneumatic distortion in pressure sensing devices

A general numerical technique for obtaining unsteady pressure measurements using conventional pressure sensing technology has been developed. A pneumatic distortion model, based on the Navier-Stokes equations of momentum and continuity, was reduced to a low-order, state-variable model retaining most of the dynamic characteristics of the full model. The reduced-order model is coupled with results from minimum variance estimation theory to develop an algorithm to compensate the effects of pneumatic distortion. Both postflight and real-time algorithms were developed and evaluated using simulated and flight data.

Whitmore, Stephen A.↗

Vapor condensation at the free surface of an axisymmetric liquid mixed by a laminar jet

This paper presents numerical solutions of jet-induced mixing in a partially full cryogenic tank. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to slove the nondimensional form of steady-state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of previously stated parameters on the condensation Nusselt and Stanton numbers that characterize the steady-state interface condensation process are investigated. Detailed analysis is performed to gain a better understanding of the fundamentals of fluid mixing and interface condensation.

Lin, Chin-Shun↗

Supersonic propulsion simulation by incorporating component models in the large perturbation inlet (LAPIN) computer code

An approach to simulating the internal flows of supersonic propulsion systems is presented. The approach is based on a fairly simple modification of the Large Perturbation Inlet (LAPIN) computer code. LAPIN uses a quasi-one dimensional, inviscid, unsteady formulation of the continuity, momentum, and energy equations. The equations are solved using a shock capturing, finite difference algorithm. The original code, developed for simulating supersonic inlets, includes engineering models of unstart/restart, bleed, bypass, and variable duct geometry, by means of source terms in the equations. The source terms also provide a mechanism for incorporating, with the inlet, propulsion system components such as compressor stages, combustors, and turbine stages. This requires each component to be distributed axially over a number of grid points. Because of the distributed nature of such components, this representation should be more accurate than a lumped parameter model. Components can be modeled by performance map(s), which in turn are used to compute the source terms. The general approach is described. Then, simulation of a compressor/fan stage is discussed to show the approach in detail.

Cole, Gary L.↗

Numerical simulation of turbomachinery flows with advanced turbulence models

A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.

Lakshminarayana, B.↗

Helium in the polar wind

The coupled time-dependent continuity, momentum, and energy equations for four charged species (H(+), He(+), O(+) and electrons) moving through a stationary upper neutral atmosphere (composed of O2, N2, O, He, and H) were solved along open geomagnetic field lines using a modified version of the Gombosi et al. (1985) model. The thermospheric densities and temperatures, obtained from the MSIS-86 model (Hedin, 1987), were used to calculate neutral-ion sources, sinks, and collision frequencies. The various cross sections and conductivities employed by the Gombosi et al. model were also updated. The paper presents the first time-dependent solutions for the He(+) component of the polar wind. Special attention is paid to the solar cycle variation of the ion upflows.

Gombosi, T. I.↗