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 523 records · Page 29

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.↗

Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements

Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dark matter in the Alternative Left Right model

The Alternative Left-Right Model is an attractive variation of the usual Left-Right Symmetric Model because it avoids flavour-changing neutral currents, thus allowing the additional Higgs bosons in the model to be light. We show here that the model predicts several dark matter candidates naturally, through introduction of an R-parity similar to the one in supersymmetry, under which some of the new particles are odd, while all the SM particles are even. Dark matter candidates can be fermionic or bosonic. We present a comprehensive investigation of all possibilities. We analyze and restrict the parameter space where relic density, direct and indirect detection bounds are satisfied, and investigate the possibility of observing fermionic and bosonic dark matter signals at the LHC. Both the bosonic and fermionic candidates provide promising signals, the first in LHC at 300 fb -1 , the second at higher luminosity, 3000 fb -1 . Signals from bosonic candidates are indicative of the presence of exotic d' quarks, while fermionic candidates imply the existence of charged Higgs bosons, all with masses in the TeV region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Two-color scattering for the measurement of neutrals at the edge of fusion devices

Laser two-color scattering (TCS) is proposed to detect the neutral species in the edge of fusion devices, namely, tokamaks. Here, TCS uses two wavelengths to probe both the laser Rayleigh scattering and Thomson scattering of the neutral-electron bath, with emphasis on neutral density measurements such as that of hydrogen and deuterium. Modeling of the Rayleigh scattering of tokamak neutral species under various plasma conditions (electron density and temperature) shows that, with an appropriate filtering of the Thomson signal and by going to ultraviolet-region wavelengths, identification of the Rayleigh signal can be achieved. Photon count and signal fractions are calculated in two test cases, one in the midplane region of the National Spherical Torus Experiment and one in the divertor region of DIII-D. An uncertainty analysis and discussion of the feasibility of the TCS diagnostic is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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.↗

Validating the simulation of beam-ion charge exchange in MAST Upgrade

Abstract Simulation of the impact of charge-exchange (CX) reactions on beam ions in the Mega Amp Spherical Tokamak (MAST) Upgrade was compared to measurements carried out with a fission chamber (neutron fluxes) and a fast ion deuterium-alpha (FIDA) diagnostic. A simple model was developed to reconstruct the outer-midplane neutral density based on measurements of deuterium-alpha emission from edge neutrals, and on Thomson scattering measurements of electron density and temperature. The main computational tools used were the ASCOT orbit-following code and the FIDASIM code for producing synthetic FIDA signals. The neutral density reconstruction agrees qualitatively with SOLPS-ITER modelling and yields a synthetic passive FIDA signal that is consistent with measurement. When CX losses of beam ions are accounted for, predicted neutron emission rates are quantitatively more consistent with measurement. It was necessary to account for CX losses of beam ions in simulations to reproduce the measured passive FIDA signal quantitatively and qualitatively. The results suggest that the neutral density reconstruction is a good approximation, that CX with edge neutrals causes significant beam-ion losses in MAST Upgrade, typically 20% of beam power, and that the ASCOT fast-ion CX model can be used to accurately predict the redistribution and loss of beam ions due to CX.

Physics↗

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.↗

Density pedestal prediction model for tokamak plasmas

Abstract A model for the pedestal density prediction based on neutral penetration combined with pedestal transport is presented. The model is tested against a pedestal database of JET-ILW Type I ELMy H-modes showing good agreement over a wide range of parameters both in standalone modelling (using the experimental temperature profile) and in full Europed modelling that predicts both density and temperature pedestals simultaneously. The model is further tested for ASDEX Upgrade and MAST-U Type I ELMy H-modes and both are found to agree with the same model parameters as for JET-ILW. The JET-ILW experiment where the isotope of the main ion is varied in a D/T scan at constant gas rate and constant β N is successfully modelled as long as the separatrix density ( n e,sep ) and pedestal transport coefficient ratio ( D / χ ) are varied in accordance with the experimentally observed variation of n e,sep and the isotope dependence of D / χ found in gyrokinetic simulations. The predictions are found to be sensitive to n e,sep which is why the model is combined with an n e,sep model to predict the pedestal for the STEP fusion reactor.

Physics↗

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.↗