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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1997

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRS2006

A Computational and Experimental Investigation of Shear Coaxial Jet Atomization

The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.

Ibrahim, Essam A.↗
NASA NTRS2008

High-Order Energy Stable WENO Schemes

A new third-order Energy Stable Weighted Essentially NonOscillatory (ESWENO) finite difference scheme for scalar and vector linear hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. A rigorous truncation error analysis is presented showing that the accuracy of the 3rd-order ESWENO scheme is drastically improved if the tuning parameters of the weight functions satisfy certain criteria. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO finite difference scheme of Jiang and Shu in terms of accuracy, while providing essentially nonoscillatory solutions near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS2010

The Beginner's Guide to Wind Tunnels with TunnelSim and TunnelSys

The Beginner's Guide to Wind Tunnels is a Web-based, on-line textbook that explains and demonstrates the history, physics, and mathematics involved with wind tunnels and wind tunnel testing. The Web site contains several interactive computer programs to demonstrate scientific principles. TunnelSim is an interactive, educational computer program that demonstrates basic wind tunnel design and operation. TunnelSim is a Java (Sun Microsystems Inc.) applet that solves the continuity and Bernoulli equations to determine the velocity and pressure throughout a tunnel design. TunnelSys is a group of Java applications that mimic wind tunnel testing techniques. Using TunnelSys, a team of students designs, tests, and post-processes the data for a virtual, low speed, and aircraft wing.

Benson, Thomas J.↗
NASA NTRS2013

The Aqua-planet Experiment (APE): Response to Changed Meridional SST Profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea- ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

AGCM↗
NASA NTRS2003

Numerical Simulations of STOVL Hot Gas Ingestion in Ground Proximity Using a Multigrid Solution Procedure

A multi grid solution procedure for the numerical simulation of turbulent flows in complex geometries has been developed. A Full Multigrid-Full Approximation Scheme (FMG-FAS) is incorporated into the continuity and momentum equations, while the scalars are decoupled from the multi grid V-cycle. A standard kappa-Epsilon turbulence model with wall functions has been used to close the governing equations. The numerical solution is accomplished by solving for the Cartesian velocity components either with a traditional grid staggering arrangement or with a multiple velocity grid staggering arrangement. The two solution methodologies are evaluated for relative computational efficiency. The solution procedure with traditional staggering arrangement is subsequently applied to calculate the flow and temperature fields around a model Short Take-off and Vertical Landing (STOVL) aircraft hovering in ground proximity.

Wang, Gang↗