Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 361 records · Page 20

The perturbed neutral circulation in the vicinity of a symmetric stable auroral arc

In the present investigation, a sophisticated high-resolution time-dependent numerical model is used to study the neutral response to changes in the E region electrodynamical forcing and particle heating associated with the appearance of a stable auroral arc. The initial condition is a nighttime ionosphere in the presence of a background precipitation of electrons associated with diffuse auroras. A simulation study is conducted of the perturbation in the neutral circulation. The simulations were performed on the basis of a model which describes nonlinear, nonhydrostatic viscous flow on a rotating atmosphere.

Walterscheid, R. L.↗

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗

Nonlinear triggered lightning models for use in finite difference calculations

Two nonlinear triggered lightning models have been developed for use in finite difference calculations. Both are based on three species of air chemistry physics and couple nonlinearly calculated air conductivity to Maxwell's equations. The first model is suitable for use in three-dimensional modeling and has been applied to the analysis of triggered lightning on the NASA F106B Thunderstorm Research Aircraft. The model calculates number densities of positive ions, negative ions, and electrons as a function of time and space through continuity equations, including convective derivative terms. The set of equations is closed by using experimentally determined mobilities, and the mobilities are also used to determine the air conductivity. Results from the model's application to the F106B are shown. The second model is two-dimensional and incorporates an enhanced air chemistry formulation. Momentum conservation equations replace the mobility assumption of the first model. Energy conservation equations for neutrals, heavy ions, and electrons are also used. Energy transfer into molecular vibrational modes is accounted for. The purpose for the enhanced model is to include the effects of temperature into the air breakdown, a necessary step if the model is to simulate more than the very earliest stages of breakdown. Therefore, the model also incorporates a temperature-dependent electron avalanche rate. Results from the model's application to breakdown around a conducting ellipsoid placed in an electric field are shown.

Rudolph, Terence↗

A Neutral-Point Expansion of the Ideal Magnetosphere

The idealized model of the geomagnetic field-solar wind interaction yields a singular (neutral) point on the magnetopause at which the magnetic field vanishes. If we expand the fields in a power series around this point, including quadratic terms, we can derive an approximate equation of the magnetopause in a small neighborhood of the neutral point which is consistent with the idealized boundary conditions to fourth order. We then consider the additional pressure due to multiple reflections of particles in this neighborhood and show that less than a 4 per cent correction to the single-reflection pressure condition is necessary.

MATHEMATICAL MODEL↗

Kinetic modeling of O(+) upflows resulting from E x B convection heating in the high-latitude F region ionosphere

We report here the results of modeling work aimed at understanding the development of ionospheric O(+) field-aligned upflows that develop in response to high-latitude E x B drift induced frictional heating. The model used is a collisional semikinetic model which includes ion-neutral resonant charge exchange and polarization collisions as well as Coulomb self-collisions. It also includes the process of chemical removal of O(+) as well as all of the macroscopic forces: ambipolar electric, gravity, magnetic mirror, and centripetal. Model results show the development of several types of non-Maxwellian velocity distributions including toroids at low altitude, distributions with large heat flow in the perpendicular component at intermediate altitudes, and distributions with a separate upflowing population or upward superthermal tail at high altitudes. Whenever the convection electric field increases from a small value (less than 25 mV/m) to a large value (100-200 mV/m) in 6 min or less large upflows develop with parallel drift speeds which peak (below 1000 km) at values between 500 m/s and 2 km/s, parallel fluxes which peak between 6.0 x 10(exp 8) and 3.2 x 10(exp 9)/sq cm/s, and parallel per particle heat flows which peak between 8.0 x 10(exp -9) and 8.0 x 10(exp -8) ergs cm/s. The higher values in these ranges occur for a cooler neutral atmosphere, with a larger convection electric field that is turned on quickly. The model produces field-aligned O(+) flow speeds that are larger than those produced by a 20-moment generalized transport model but smaller then those produced by an isotropic hydrodynamic model for comparable values of the convection turn on times. The model results compare favorably with some topside satellite and radar data.

Wilson, G. R.↗

Low-temperature Galactic fountains

Analytic and one- and two-dimensional hydrodynamic calculations are performed to determine whether or not the Galactic fountain model can explain observations of neutral gas in the Galactic halo. Galactic halo gasdynamics are examined, including thermal conduction, Galactic rotation, and Galactic fountains. Results are presented from numerical modeling of low-temperature Galactic fountains. It is found that the best reproduction of the observations is provided by a model near the transonic regime with a temperature and density at the base of the fountain of 300,000 K and 0.001/cu cm, respectively. The cooling time for the hot gas is 33 million yr and the time for newly formed clouds to return to the disk is 47 million yr, neglecting drag forces.

Houck, John C.↗

Development of boundary layer rolls from dynamic instabilities

Using three three-dimensional spectral models of flow in a neutral atmosphere, the development of atmospheric boundary layer rolls from the inflection point and parallel instabilities is examined. Using the Ekman profile, the ability of the various models to reproduce the modes arising from the parallel and inflection point instability mechanisms was tested by comparing the results from the spectral models with the results from larger models used by other investigators. It was found that the truncated spectral models used in this investigation yielded preferred values of orientation angle and aspect ratio that agreed with those produced by the larger models, as well as with the observations data taken on September 1981 during the West-German KonTur experiment.

Stensrud, David J.↗

The structure of cometary ionospheres. II - CO-rich comets

The structure of the ionosphere of a CO-rich comet is computed using two different models. The first one, the photochemical model, assumes that the dissociation and ionization of cometary neutrals and ions are due to photoionization, and photodissociation by solar UV radiation together with dissociative recombinations and ion-neutral reactions. The second one, the internal source model, also incorporates the ionization and dissociation effects of an electric current discharging through the inner coma. The generation of this current has been discussed in earlier papers. It is concluded that the internal source model can explain qualitatively the basic morphology of the ionospheres of CO-rich comets such as Humason (1962, VIII) and Morehouse (1908, III), whereas the photochemical model cannot. The main aim of this paper is not so much to provide accurate numerical estimates as to draw attention to a process which may very well dominate the structures of cometary ionospheres.

Ip, W.-H.↗

Ionospheric composition: The seasonal anomaly explained

The main photochemical processes of the ionosphere are reanalyzed in the light of laboratory measurements of rate coefficients, using the Atmosphere Explorer data. Major changes to the chemistry include the transfer of nearly all metastable 0(+) ions to N2(+) via charge exchange with N2. The N2(+) ions become vibrationally excited by resonant fluorescence of solar near UV and near infrared radiation, leading to a return transfer of N2(+) ionization to 0(+) by charge exchange or vibrationally excited N2(+) with atomic oxygen. With this chemistry the seasonal variations in the peak electron densities are then shown to be caused primarily by anomalous seasonal variations in neutral composition. The required neutral composition variations are empirically produced by the MSIS model atmosphere. The circulation derived from recent 3D models of the global thermosphere qualitatively accounts for the seasonal variations in neutral composition predicted by the MSIS model. In addition to the composition effect, vibrationally excited N2 is found to contribute a 20% effect to the anomalous seasonal behavior at solar maximum.

Torr, D. G.↗

Modeling the unidentified infrared emission with combinations of polycyclic aromatic hydrocarbons

The infrared emission band spectrum associated with many different interstellar objects can be modeled successfully by using combined laboratory spectra of neutral and positively charged polycyclic aromatic hydrocarbons (PAHs). These model spectra, shown here for the first time, alleviate the principal spectroscopic criticisms previously leveled at the PAH hypothesis and demonstrate that mixtures of free molecular PAHs can indeed account for the overall appearance of the widespread interstellar infrared emission spectrum. Furthermore, these models give us insight into the structures, stabilities, abundances, and ionization balance of the interstellar PAH population. These, in turn, reflect conditions in the emission zones and shed light on the microscopic processes involved in the carbon nucleation, growth, and evolution in circumstellar shells and the interstellar medium.

Spectrophotometry, Infrared↗

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗

Retarding potential analyzer research

Studies were continued in the areas of ionospheric electrodynamics, plasma irregularities, and ion-neutral coupling. The spectral characteristics of structure in the ion concentration and the electric field at high latitudes were examined. It was found that in the region of velocity shears the relationship between the slopes of the two spectra are dependent on the existence of large field-aligned currents and apparently large ionospheric conductivity gradients. The relationship between the Joule dissipation and the locally measured value of the ion temperature were studied by examining the horizontal ion drift velocity. The importance of ion-neutral coupling was investigated in additional studies involving the nature of the coupling at high latitude boundaries and during times of magnetic disturbances. Much of this behavior can be modeled using globular circulation models. In order to study the storm time response of the thermosphere the global ion convection pattern was derived from DE-2 measurements and used as input in a circulation model. The calculated variables of neutral velocity and composition and temperature were compared with measurements to validate the technique.

Hanson, W. B.↗

The auroral 6300 A emission - Observations and modeling

A tomographic inversion is used to analyze measurements of the auroral atomic oxygen emission line at 6300 A made by the atmosphere explorer visible airglow experiment. A comparison is made between emission altitude profiles and the results from an electron transport and chemical reaction model. Measurements of the energetic electron flux, neutral composition, ion composition, and electron density are incorporated in the model.

Solomon, Stanley C.↗

PATRAN 2.5/EAL interface guide

The PATRAN/EA interface guide describes two programs, EALPAT and PATEAL. EALPAT reads an EAL LO1 file and translates the model and results into a PATRAN 2.5 neutral file, element results file, and nodal results file. An EAL model can be color coded in PATRAN, and geometry, loads, boundary conditions, section and material properties, rigid masses, springs, and beam orientations can be plotted and debugged. EAL results can be brought into PATRAN as element or nodal quantities and displayed as deformed plots, animated shapes, color coded elements, or color filled contour plots. PATEAL converts a PATRAN 2.5 data base into an EAL runstream. Geometry, including all element types, alternate coordinate systems, material properties, section properties, loads and boundary conditions are all converted. EALPAT and PATEAL can also be used together with the PATNAS translator from PDA Engineering to convert an EAL runstream to an MSC/NASTRAN bulk data file.

Feldhaus, Winifred S.↗

An investigation of the solar zenith angle variation of D-region ionization

Model calculations are carried out with a view to interpreting the solar zenith angle variation of D-region ionization. A model is developed for the neutral chemistry including the transport terms relating to molecular and eddy diffusion. The diurnal behavior is described of the minor neutral constituents formed in an oxygen-hydrogen-nitrogen atmosphere, in the height interval between 30 and 120 km. Computations carried out for two cases of the eddy diffusion coefficients models indicate that the constituents which are important for the D-region positive-ion chemistry do not show a significant variation with zenith angle for values up to 75 deg over the D-region heights. In the ion chemistry model, ion-pair production rates are calculated for solar X-rays between 1 A and 100 A, EUV radiations from 100 A up to the Lyman-alpha line, precipitating electrons, and galactic cosmic rays. The solar zenith angle variation of the positive-ion composition, negative-ion composition, and the electron densities are described up to 75 deg zenith angle, in the height interval between 60 and 100 km.

Ratnasiri, P. A. J.↗

Numerical Investigation of Occupant Injury Risks During A Realistic Transport Aircraft Crash Conditions

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a full-scale crash test of a Fokker F28 MK1000 aircraft to investigate the performance of transport aircraft under realistic crash conditions. This crash test was computationally recreated using finite element (FE) human body models (HBMs) to further explore potential injury risks to occupants and analyze the utilization of HBMs in the aerospace crash environment. The Global Human Body Model Consortium (GHBMC) male 50th percentile occupant detailed model (v6.0) and the Toyota Human Model for Safety (THUMS) male 50th percentile occupant model (v6.1) were selected to be used in the crash simulations. The HBMs were simulated in conditions matching those of anthropomorphic test device (ATD) experiments included within the aircraft cabin during the crash test. Seven occupant locations within the cabin were simulated utilizing each of the models. The models were positioned in a neutral upright posture with hands resting on the legs and the feet contacting the floor. Head, brain, neck, and lumbar vertebra injury metrics were calculated for all trials. Both HBM models required minor modifications to stabilize these simulations. The GHBMC model required added erosion for six parts while the THUMS model only required one in order to complete the full simulation. The THUMS model, however, required a much smaller timestep for stability and therefore took significantly more computational time. In addition the GHBMC model includes integrated instrumentation while the THUMS model requires development and implementation of instrumentation. Both models predicted 100% injury risk for lumbar vertebra fracture in all test conditions. This prediction was in family with high lumbar load values measured by the ATDs during the crash test. The THUMS model consistently predicted lower injury risks than the GHBMC model in all three other metrics varying depending on the crash pulse. Overall, the THUMS model required less modifications to allow for this study. However, the GHBMC models significantly faster run time and integrated instrumentation make it a more intuitive model for this research.

Crashworthiness↗

A quantitative model of plasma in Neptune's magnetosphere

A model encompassing plasma transport and energy processes is applied to Neptune's magnetosphere. Starting with profiles of the neutral densities and the electron temperature, the model calculates the plasma density and ion temperature profiles. Good agreement between model results and observations is obtained for a neutral source of 5 x 10 exp 25/s if the diffusion coefficient is 10 exp -8 L3R(N)/2s, plasma is lost at a rate 1/3 that of the strong diffusion rate, and plasma subcorotates in the region outside Triton.

Richardson, John D.↗