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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Development and evaluation of improved numerical schemes for recirculating flows

The paper examines the performance of the flux-spline scheme for convection-diffusion. Computations are presented for a number of test cases, both linear and nonlinear. It is shown that in all cases the flux-spline scheme yields results which are superior to those obtained with the lower-order formulations such as hybrid differencing. In order to improve the computational efficiency, the flux-spline scheme has been combined with a direct solution algorithm for the continuity and momentum equations. Such an approach eliminates the need for an equation for pressure or pressure correction and is found to be rapidly convergent.

Patankar, S. V.↗
NASA NTRS1986

Temporal features of the refilling of a plasmaspheric flux tube

The refilling of plasmaspheric flux tubes was studied by assuming that the protonosphere provides an ionospheric boundary where the H(+) density can be assumed; the supersonic flow in the flux tube is driven by the depletion of the plasma from the flux tube, while the base density and pressure in the protonosphere remain constant. The time-dependent continuity and momentum equations for the H(+) ions were solved. The electron gas was assumed to obey the Boltzmann law, and the proton gas was assumed to be isothermal. In agreement with the postulate of Banks et al. (1971), it was found that an important feature of the refilling is the formation of a shock pair at the equator; as the shocks propagate toward the ionosphere, the refilling occurs. Depending on the density at the ionospheric boundaries, a fair agreement was found between the refilling rates obtained for L = 6.6 and those from the GEOS 2 observations.

Singh, Nagendra↗
NASA NTRS1988

A unidimensional model of comet ionosphere structure

The one-dimensional continuity and momentum equations for the plasma formed in the expanding coma of a comet near the sun are studied. An analytic expression for the magnetic field configuration in the presence of outflow, photoionization, dissociative recombination, plasma fluid pressure, and friction between the ions and neutrals is obtained. It is suggested that for a Halley-type comet there will be a region sunward of the nucleus from which the magnetic field is excluded, consistent with Giotto observations. Calculations have been performed for Halley and Giacobini-Zinner type comets, and in the field-free regions it is shown that the dominant terms in the momentum equation balancing the magnetic pressure gradient are the ion neutral friction and the net mass loading momentum gain.

Eviatar, Aharon↗
NASA NTRS1987

Deterministic representation of chaos with application to turbulence

Chaotic motions of nonlinear dynamical systems are decomposed into mean components and fluctuations. The approach is based upon the concept that the fluctuations driven by the instability of the original (unperturbed) motion grow until a new stable state is approached. The Reynolds-type equations written for continuous as well as for finite-degrees-of-freedom dynamical systems are closed by using this stabilization principle. The theory is applied to conservative systems, to strange attractors and to turbulent motions.

Zak, M.↗
NASA NTRS1989

Effects of transverse shear on large deflection random response of symmetric composite laminates with mixed boundary conditions

Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived considering von Karman strains and transverse shear effects. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The effects of transverse shear on large deflection vibration of laminated plates with mixed boundary support conditions under random excitation are studied. Mean-square deflections and mean-square inplane strains are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for plates with two opposing edges simply supported and the remaining edges clamped.

Prasad, C. B.↗
NASA NTRS1989

Precise solution for a finite set of spherical coefficients from equiangular gridded data

An important goal of geodesy is to determine the anomalous potential and its derivatives outside of the earth. Representing the surface anomalies by a series of spherical harmonics is useful since it is then possible to do a term by term solution of Laplace's equation and upward continuation. The problem of finding such a spherical harmonic series for anomaly values given on an equiangular surface grid is addressed. (This is a first step toward the more complicated problem of finding a function such that locally averaged values fit a grid of mean anomalies.) Three approaches to this fitting problem are discussed and compared: the discrete Fourier technique, the discrete integral technique, and a new approach. The peculiar nature of the equiangular grid, with its increasing density of (noisy) data toward the poles, causes each method to exhibit a different type of difficulty. The new method is shown to be practical as well as precise since the numerical conditioning problems which appear can be successfully handled by such well-known techniques as a (simple) Kalman filter.

Zucker, Paul A.↗
NASA NTRS1989

Near contact phenomena and transient effects in far infrared photoconductors

A combination of experimental and modeling work is summarized in two areas: first, the calculation of excess free carrier and space charge distributions near contacts and their effects on device resistivity, and second, the characterization of a slow transient response (tau approx. 1 sec) in Ge:Be detectors which is due to trapping associated with Be(+) formation. In both cases, analytical models, based on continuity and rate equations, have been developed to enable the application of these findings to a wide variety of photoconductor materials.

Haegel, Nancy M.↗
NASA NTRS1990

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary tasks during January 1990 to June 1990 have been the development and evaluation of various electron and electron-electronic energy equation models, the continued development of improved nonequilibrium radiation models for molecules and atoms, and the continued development and investigation of precursor models and their effects. In addition, work was initiated to develop a vibrational model for the viscous shock layer (VSL) nonequilibrium chemistry blunt body engineering code. Also, an effort was started associated with the effects of including carbon species, say from an ablator, in the flowfield.

Carlson, Leland A.↗
NASA NTRS1991

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗
NASA NTRS1990

Numerical modeling of the atmosphere with an isentropic vertical coordinate

A theta-coordinate model simulating the nonlinear evolution of a baroclinic wave is presented. In the model, vertical discretization maintains important integral constraints such as conservation of the angular momentum and total energy. A massless-layer approach is used in the treatment of the intersections of coordinate surfaces with the lower boundary. This formally eliminates the intersection problem, but raises other computational problems. Horizontal discretization of the continuity and momentum equations in the model are designed to overcome these problems. Selected results from a 10-day integration with the 25-layer, beta-plane version of the model are presented. It is concluded that the model can simulate the nonlinear evolution of a baroclinic wave and associated dynamical processes without major computational difficulties.

Hsu, Yueh-Jiuan G.↗
NASA NTRS1991

Plasma expansion characteristics of ionized clouds in the ionosphere - Macroscopic formulation

A macroscopic plasma expansion model, based on a numerical solution of the type-dependent nonlinear coupled continuity and momentum equations for background O(+) ions and several released ion species, was used to examine plasma expansion characteristics of ionized clouds in the ionosphere. Information is obtained on Ba(+), Li(+), and Ba(+) - Li(+) clouds; different cloud sizes (Gaussian half widths of 0.1, 0.5, and 1 km); cloud/background ion density ratios covering two orders of magnitude; electron/ion temperature ratios of 1, 5, and 10; and several cloud-background relative velocities along B (0, 1, 2, 4, 6, and 8 km/sec). The macroscopic expansion features obtained were found to be in general agreement with those obtained from the small-scale numerical simulations.

Schunk, R. W.↗
NASA NTRS1992

Observations of the structure and vertical transport of the polar upper ionosphere with the EISCAT VHF radar. I - Is EISCAT able to determine O(+) and H(+) polar wind characteristic? A simulation study

A method presented by Wu et al. (1992) for computing the H(+) vertical velocity from the main ionospheric parameters measured by the EISCAT VHF radar is tested in a fully controlled sequence which consists of generating an ideal ionospheric model by solving the coupled continuity and momentum equations for a two-ion plasma (O(+) and H(+)). Synthetic autocorrelation functions are generated from this model with the radar characteristics and used as actual measurements to compute the H(+) vertical velocities. Results of these simulations are shown and discussed for three cases of typical and low SNR and for low and increased mixing ratios. In most cases general agreement is found between computed H(+) velocities and generic ones with the altitude range considered, i.e., 200-1000 km. The method is shown to be reliable.

Blelly, Pierre-Louis↗
NASA NTRS1992

Numerical simulations of drainage flows on Mars

Data collected by Viking Landers have shown that the meteorology of the near surface Martian environment is analogous to desertlike terrestrial conditions. Geological evidence such as dunes and frost streaks indicate that the surface wind is a potentially important factor in scouring of the martian landscape. In particular, the north polar basin shows erosional features that suggest katabatic wind convergence into broad valleys near the margin of the polar cap. The pattern of katabatic wind drainage off the north polar cap is similar to that observed on Earth over Antarctica or Greenland. The sensitivity is explored of Martian drainage flows to variations in terrain slope and diurnal heating using a numerical modeling approach. The model used is a 2-D sigma coordinate primitive equation system that has been used for simulations of Antarctic drainage flows. Prognostic equations include the flux forms of the horizontal scalar momentum equations, temperature, and continuity. Parameterization of both longwave (terrestrial) and shortwave (solar) radiation is included. Turbulent transfer of heat and momentum in the Martian atmosphere remains uncertain since relevant measurements are essentially nonexistent.

Parish, Thomas R.↗
NASA NTRS1992

Dynamical polar wind and its response to kinetic ion heating

The effect of the kinetic ion heating of the polar wind on the dynamic expansion of the polar wind outflow is investigated by incorporating the kinetic ion heating effect into a time-dependent hydrodynamic polar wind model in which the polar wind O(+) and H(+) continuity and momentum equations are solved in a flux tube from ionospheric to magnetospheric altitudes. Two limiting cases were considered: (1) preferential O(+) escape and (2) preferential H(+) heating. It was found that ion heating does not have much effect on the escape of H(+) ions because of the limiting nature of the H(+) escape fluxes. However, O(+) heating can substantially increase the upward O(+) pressure gradient in the heating region, indicating that O(+) heating is a viable explanation for the unexpectedly large O(+) fluxes frequently observed in the polar magnetosphere.

Chen, Margaret W.↗
NASA NTRS1994

The energetics of Titan's ionosphere

We have developed a comprehensive model to study the dynamics and energetics of the ionosphere of Titan. We solved the one-dimensional, time-dependent, coupled continuity and momentum equations for several ion species, together with single ion and electron energy equations, in order to calculate density, velocity, and temperature profiles. Calculations were carried out for several cases corresponding to different local times and configurations of the Titan-Saturn system. In our model the effects of horizontal magnetic fields were assumed to be negligible, except for their effect on reducing the electron and ion thermal conductivities and inhibiting vertical transport in the subram region. The ionospheric density peak was found to be at an altitude of about 1100 km, in accordance with earlier model calculations. The ionosphere is chemically controlled below an altitude of about 1500 km. Above this level, ion densities differ significantly from their chemical equilibrium values due to strong upward ion velocities. Heat is deposited in a narrow region around the ionospheric peak, resulting in temperature profiles increasing sharply and reaching nearly constant values of 800-1000 deg K for electrons and 300 deg K for ions in the topside, assuming conditions appropriate for the wake region. In the subram region magnetic correction factors make the electron heat conductivities negligible, resulting in electron temperatures increasing strongly with altitude and reaching values in the order of 5000 deg K at our upper boundary located at 2200 km. Ion chemical heating is found to play an important role in shaping the ion energy balance in Titan's ionosphere.

Roboz, A.↗
NASA NTRS1995

Large eddy simulations of compressible turbulent flows

An evaluation of existing models for Large Eddy Simulations (LES) of incompressible turbulent flows has been completed. LES is a computation in which the large, energy-carrying structures to momentum and energy transfer is computed exactly, and only the effect of the smallest scales of turbulence is modeled. That is, the large eddies are computed and the smaller eddies are modeled. The dynamics of the largest eddies are believed to account for most of sound generation and transport properties in a turbulent flow. LES analysis is based on an observation that pressure, velocity, temperature, and other variables are the sum of their large-scale and small-scale parts. For instance, u(i) (velocity) can be written as the sum of bar-u(i) and u(i)-prime, where bar-u(i) is the large-scale and u(i)-prime is the subgrid-scale (SGS). The governing equations for large eddies in compressible flows are obtained after filtering the continuity, momentum, and energy equations, and recasting in terms of Favre averages. The filtering operation maintains only large scales. The effects of the small-scales are present in the governing equations through the SGS stress tensor tau(ij) and SGS heat flux q(i). The mathematical formulation of the Favre-averaged equations of motion for LES is complete.

Porter-Locklear, Freda↗
NASA NTRS1999

Evolution of Plasmaspheric Refilling: A Comparison of Measurements with an Interhemispheric Plasmasphere Model

The evolution of plasmasphere is closely coupled to processes governing flux tube refilling and transport. In this study we follow this evolution through a series of sequential ion observations made by RIMS on the DE 1 satellite and with the intervening time intervals filled by simulations made with the FLIP model. The FLIP model solves the continuity and momentum equations for the major ion species as well as the energy equations for ions and electrons along entire flux tubes rom 100 k altitude in both hemispheres. Convection has recently been included, driven a dawn-dusk electric field model. Since the observations are necessarily restricted to two local time passes through the plasmasphere per orbit, the FLIP code will follow the observed flux tubes as they corotate and convect through the other local times. The period of DE 1 is such that nearly the same flux tubes are observed every second day or every sixth orbit. The composite of observations and simulations allows assessment of the degree to which the degree to which the FLIP simulations accurately describe the processes which govern flux tube refilling over the two-day period between observations of the same flux tubes.

Comfort, Richard H.↗
NASA NTRS1998

Effects of Convection Electric Fields on Modeled Plasmaspheric Densities and ccc Temperatures

This paper examines the effects of convection electric fields on plasmaspheric H+, O+, He+, and N+ densities and electron and ion temperatures. These effects are studied with the aid of the Field Line Interhemispheric Plasma (FLIP) model, which has recently been extended to include the effects of ExB drifts. The FLIP model solves the continuity and momentum equations for the major ion species as well as the energy equations for ions and electrons along entire drifting flux tubes from 100 km altitude in the northern hemisphere to 100 km altitude in the southern hemisphere. Electron heating in the ionosphere and plasmasphere is provided by the solution of two-stream equations for photoelectrons. The dawn-dusk electric field imposed by the solar wind causes changes in plasmaspheric density and temperature as the plasma drifts onto flux tubes having different volumes. In an idealized convection model, outward drifts in the afternoon cause decreases in the plasmasphere density and temperature while inward drifts in the evening cause increases in plasmasphere density and temperature. In this paper we examine the effects of convection electric fields on the rate of refilling of flux tubes and investigate the hypothesis that convection electric fields are responsible for the unusually high evening electron temperatures and the post-midnight density maxima often observed in the winter ionosphere above Millstone Hill.

Comfort, Richard H.↗