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 559 records · Page 31

Carbon neutrality in Malaysia and Kuala Lumpur: Insights from stakeholder-driven integrated assessment modeling

Several cities in Malaysia have established plans to reduce their CO2 emissions, in addition to Malaysia submitting a Nationally Determined Contribution to reduce its carbon intensity (against GDP) by 45% in 2030 compared to 2005. Meeting these emissions reduction goals will require a joint effort between governments, industries, and corporations at different scales and across sectors. In collaboration with national and sub-national stakeholders, we developed and used a global integrated assessment model to explore emissions mitigation pathways in Malaysia and Kuala Lumpur. Guided by current climate action plans, we created a suite of scenarios to reflect uncertainties in policy ambition, level of adoption, and implementation for reaching carbon neutrality. Through iterative engagement with all parties, we refined the scenarios and focus of the analysis to best meet the stakeholders’ needs. We found that Malaysia can reduce its carbon intensity and reach carbon neutrality by 2050, and that action in Kuala Lumpur can play a significant role. Decarbonization of the power sector paired with extensive electrification, energy efficiency improvements in buildings, transportation, and industry, and the use of advanced technologies such as hydrogen and carbon capture and storage will be major drivers to mitigate emissions, with carbon dioxide removal strategies being key to eliminate residual emissions. This study highlights the participatory process in which stakeholders contributed to the development of the model and guided the analysis, as well as insights into Malaysia’s decarbonization potential and the role of multilevel governance.

cities↗

Comparing event generator predictions and ab initio calculations of ν- 12 C neutral-current quasielastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs are used to predict kinematic observables for a range of neutrino energies. These simulations make use of simple models of nucleon dynamics. As such, they may fail to capture features of more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the adequacy of the implementation of two nuclear physics models in EG simulations. Neutral-current quasielastic (NCQE) scattering events for neutrinos and antineutrinos on a 12 C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the relativistic Fermi gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. The peak of the cross section dσ/(dΩ dω) is consistent in energy transfer, ω, for RFG and CBF simulations. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two. These results demonstrate the importance of benchmarking EGs so their nuclear physics implementations can be improved for the analysis of future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A model of Triton's role in Neptune's magnetosphere

Escape of neutral hydrogen (H) and nitrogen (N) from Triton's maintains a large neutral cloud, called the Triton torus, in Neptune's magnetosphere. We have developed the first detailed Monte Carlo simulation model of the Triton torus that includes the collisionality, the complex geometry, the injection of two neutral species from Triton (H and N), and the combined effects of photoionization, electron impact ionization, and charge exchange. Ionization in Neptune's plasma sheet was modeled using Voyager plasma observations. Collisions cause both the H and N neutral clouds to become more radially extended, both toward Neptune and out beyond the magnetopause, as well as more extended in latitude, when compared with collisionless models. Moreover, collisions of H with the much more massive N greatly enhance the collisional ejection of H from the system and into Neptune's atmosphere. This effect decreases the probability of H ionization within the magnetosphere relative to that for N, and furthermore causes model results for two-species injection from Triton to differ significantly from those for H injection alone. For a hydrogen escape rate from Triton of 5 x 10(exp 25)/s, as given by photo-chemical models of Triton's upper atmosphere, a nitrogen escape rate of 5 x 10(exp 24)/s gives proton and N(+) sources of 5.6 x 10(exp 24)/s and 3.3 x 10(exp 24)/s, respectively, whose ratio is close to the observed ratio of protons to heavies. A nitrogen escape rate of 2 x 10(exp 25)/s, yields an N(+) source more than twice that of protons, inconsistent with the Voyager data.

Decker, R. B.↗

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.

Lipatov, Alexander S.↗

The Global Mapping of Electron Precipitation and Ionospheric Conductance From Whistler-Mode Chorus Waves

Auroral precipitation is the second major energy source after solar irradiation that ionizes the Earth’s upper atmosphere. Diffuse electron aurora caused by wave-particle interaction in the inner magnetosphere (L < 8) takes over 60% of total auroral energy flux, strongly contributing to the ionospheric conductance and thus to the ionosphere-thermosphere dynamics. This paper quantifies the impact of chorus waves on the diffuse aurora and the ionospheric conductance during quiet, medium, and strong geomagnetic activities, parameterized by AE <100, 100 < AE < 300, and AE > 300, respectively. Using chorus wave statistics and inner-magnetosphere plasma conditions from Timed History Events and Macroscale Interactions during Substorms (THEMIS) observations, we directly derive the energy spectrum of diffuse electron precipitation under quasi-linear theory. We then calculate the height-integrated conductance from the wave-driven aurora spectrum using the electron impact ionization model of Fang et al. (Geophys. Res. Lett., 2010, 37) and the MSIS atmosphere model. By utilizing Fang’s ionization model, the US Naval Research Laboratory Mass Spectrometer and Incoherent Scattar Radar (NRLMSISE-00) model from 2000s for the neutral atmosphere components, and the University of California, Los Angeles (UCLA) Full Diffusion Code, we improve upon the standard generalization of Maxwellian diffuse electron precipitation patterns and their resulting ionosphere conductance. Our study of global auroral precipitation and ionospheric conductance from chorus wave statistics is the first statistical model of its kind. We show that the total electron flux and conductance pattern from our results agree with those of Ovation Prime model over the pre-midnight to post-dawn sector as geomagnetic activity increases. Our study examines the relative contributions of upper band chorus (UBC) and lower band chorus wave (LBC) driven conductance in the ionosphere. We found LBC waves drove diffuse electron precipitation significantly more than UBC waves, however it is possible that THEMIS data may have underestimated the upper chorus band wave observations for magnetic latitudes below 65 degrees. ° .

diffuse auroral precipitation↗

Constraints on the neutral hydrogen corona at Uranus from its interaction with magnetospheric plasma

Models of collisional reactions between the hot (several keV) protons and neutral atomic hydrogen corona observed at Uranus by Voyager 2 are compared with observations in order to deduce constraints on the coronal density. An upper limit is found for n sub H (5R sub U) of several tens of H/cu cm on the sunward side of Uranus and several H/cu cm on the dark side. The possibility that the hot plasma population might arise from photoionization and electron impact ionization of coronal and interstellar medium hydrogen in the magnetotail, followed by pickup and energization while convecting toward Uranus, is examined and rejected.

Herbert, Floyd↗

Tungsten and Barium Transport in the Internal Plasma of Hollow Cathodes

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the flow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushedback to the emitter surface by the electric field and drag from the xenon ion flow. Thisbarium ion flux is sufficient to maintain a barium surface coverage at the downstream endgreater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length,so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollowcathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

Hollow cathodes↗

Barium Depletion in Hollow Cathode Emitters

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the ow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushed back to the emitter surface by the electric field and drag from the xenon ion flow. This barium ion flux is sufficient to maintain a barium surface coverage at the downstream end greater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length, so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollow cathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

plasma discharges↗

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗

Implications of observed charge states of low-energy solar cosmic rays

Recent measurements of the charge states of low-energy (about 100 keV/nucleon) solar cosmic rays at 1 AU are discussed. The measurements are consistent with models involving charge equilibrium with neutral matter at the sun only if the particles lose about 90% of their energy owing to adiabatic deceleration in the solar wind. Such an energy loss is shown to be possible only if the diffusion coefficient for 1-MeV/nucleon particles is smaller than 100 quintillion sq cm per sec. The implications of these results for models of solar-cosmic-ray acceleration are discussed.

Jokipii, J. R.↗

Data analysis and interpretation related to space system/environment interactions at LEO altitude

Several studies made on the interaction of active systems with the LEO space environment experienced from orbital or suborbital platforms are covered. The issue of high voltage space interaction is covered by theoretical modeling studies of the interaction of charged solar cell arrays with the ionospheric plasma. The theoretical studies were complemented by experimental measurements made in a vacuum chamber. The other active system studied was the emission of effluent from a space platform. In one study the emission of plasma into the LEO environment was studied by using initially a 2-D model, and then extending this model to 3-D to correctly take account of plasma motion parallel to the geomagnetic field. The other effluent studies related to the releases of neutral gas from an orbiting platform. One model which was extended and used determined the density, velocity, and energy of both an effluent gas and the ambient upper atmospheric gases over a large volume around the platform. This model was adapted to study both ambient and contaminant distributions around smaller objects in the orbital frame of reference with scale sizes of 1 m. The other effluent studies related to the interaction of the released neutral gas with the ambient ionospheric plasma. An electrostatic model was used to help understand anomalously high plasma densities measured at times in the vicinity of the space shuttle orbiter.

Raitt, W. John↗

Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories

The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event by event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Non-equilibrium low-temperature plasma-assisted combustion of iso-octane: Perturbing pyrolysis and oxidation kinetics

Here, in this study, a plasma-coupled flow reactor facility is used to examine the effects of non-equilibrium low-temperature plasmas on perturbing the pyrolysis and oxidation kinetics of iso-octane. Experiments were performed in highly dilute reactive mixtures of nitrogen, at near isothermal conditions for temperatures ranging from 523 K to 1203 K. Experiments cumulatively demonstrated enhanced chemical reactivity with the plasma for temperatures below 900 K, and a lowering of the hot-ignition temperature. Detailed kinetic insight was derived from a 0D plasma-coupled kinetic model, utilizing a constructed mechanism that combined both plasma-specific chemistry and the neutral combustion chemistry. For pyrolysis conditions, the model displayed relatively good agreement with fuel consumption and the formation of most intermediates compared to the experimental data, demonstrating the model is able to accurately predict primary radical formation from the plasma directly interacting with the fuel. Enhanced reactivity was attributed to collisional quenching of excited-states of N 2 with fuel, which led to efficient fuel fragmentation and enhancement of the H-radical flux. For oxidation conditions, the model displayed satisfactory agreement with the experiments. Model predictions were able to accurately predict fuel consumption and most intermediate speciation data for T > 800 K, but most discrepancies were towards T < 800 K in particular with oxygenated intermediates. In the presence of oxygen, plasma effects were predominantly spent on efficient enhancement of O- and H-radical fluxes, leading to further fuel fragmentation and initiation of both the OH- and HO 2 -radical pools. Subsequent reactivity of iso-octane was then dictated by the response of the temperature-dependent neutral chemistry. At low-temperatures (T = 643 K), enhanced fuel radicals and O 2 -additon chemistry lead to the formation of oxygenated species, while at intermediate temperatures (T = 843 K) net decrease in OH-radical reactivity led to an increase in hydrocarbon speciation. Near the self-ignition threshold (T = 1163 K), radicals generated by high-temperature branching reactions dominate the oxidation process and effectively ignition. This study ultimately demonstrated that the enhancement of radicals afforded by the plasma causes a deviation in known understanding of iso-octane kinetics in some regards and warrants future studies to reconcile these discrepancies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Phase and amplitude relationships of wave structure observed in the lower thermosphere

Data from the Atmosphere Explorer-C satellite clearly exhibit wavelike variations in neutral composition, ion density, and electron temperature, which appear to be a general feature of the atmosphere. The neutral constituents do not exhibit uniform wave characteristics since the density variations of argon are approximately twice those of molecular nitrogen and helium structure has about one-half the amplitude of the N2 variation. The waves evident in the ion density are nearly in phase with the heavy neutrals, while the electron-temperature variations are predominantly out of phase with those in the ion density. A simple model is suggested to explain the neutral composition results, wherein the enhancements in the major gas densities are in phase with the vertical component of the perturbation velocity of the gas. The vertical velocity modifies the composition by transporting parcels of air to higher or lower regions where the composition is different. The phase relationship between density and velocity implies phase velocities (assuming that these are gravitational waves) of the order of 500 meters per sec.

Reber, C. A.↗

The global distribution of thermospheric odd nitrogen for solstice conditions during solar cycle minimum

A two-dimensional model of odd nitrogen in the thermosphere and upper mesosphere is described. The global distributions of nitric oxide and atomic nitrogen are calculated for the solstice period for quiet and moderate magnetic activity during the solar minimum period. The effect of thermospheric transport by winds is investigated along with the importance of particle-induced ionization in the auroral zones. The results are compared with rocket and satellite measurements, and the sensitivity of the model to eddy diffusion and neutral winds is investigated. Downward fluxes of NO into the mesosphere are given, and their importance for stratospheric ozone is discussed. The results show that the summer-to-winter pole meridional circulation transports both NO and N(S-4) across the solar terminator into the polar night region where there is a downward vertical transport toward the mesosphere. The model shows that odd nitrogen densities at high winter latitudes are entirely controlled by particle precipitation and transport processes.

Gerard, J.-C.↗

Storm-Time Variability of Terrestrial Hydrogen Exosphere: Kinetic Simulation Results

Recent studies of TWINS Lyman-α observations have reported an increase in geocoronal column brightness during geomagnetic storms, indicating enhanced exospheric hydrogen atom density (NH). This suggests a complex role of exospheric neutrals in determining storm-time magnetosphere dynamics and their energy release through charge-exchange processes. We developed a Model for Analyzing Terrestrial Exosphere (MATE) to investigate storm-time exospheric behaviors and their physical drivers. MATE traces test hydrogen atoms backward in time from locations in the exosphere to a nominal exobase altitude of 500 km, employing Newtonian mechanics with gravitational force. The model then calculates the phase-space densities (PSDs) of test hydrogen atoms at the exobase using the Maxwellian distribution with physics-based exobase conditions from the TIMEGCM upper atmosphere model. MATE maps PSDs at the exobase to the exosphere using Liouville’s Theorem under collisionless assumptions and derives NH by integrating the PSDs across velocity space. We conducted MATE simulation before, during, and after a minor geomagnetic storm from 12 to 18 June 2008, and compared the model results with NH estimates from the TWINS geocorona data. MATE reproduces storm-time density enhancements soon after the minimum Dst is reached, matching well with a general trend of TWINS NH estimates. The results suggest that upper atmospheric heating during a geomagnetic storm increases the number of ballistic and escaping hydrogen atoms entering the exosphere from the exobase, thereby boosting NH. However, the magnitude of modeled NH mismatches the TWINS NH estimates. The potential mechanisms of this density discrepancy include the physics excluded in the MATE model — such as neutralneutral collisions, neutral-plasma charge exchange, solar radiation pressure, and photoionization — as well as the higher exobase hydrogen density of TIMEGCM compared to typical empirical values, which will be addressed in future.

Hyunju K Connor↗

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.

Lipatov, Alexander S.↗