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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Analysis of EUV Dayglow Spectra of Triton, Titan and Earth

We have constructed a coupled ion and neutral chemical model of the ionosphere and thermosphere of Titan. Along with density profiles of ions and minor neutrals, we have computed the heating rates and heating efficiencies for the neutrals and chemical heating rates and efficiencies for the ions. We find that the neutral heating efficiency in our standard model varies from about 30% near 800 km to 22% near 2000 km. The most important heating processes are neutral-neutral reactions and photodissociation. The ion chemical heating rates maximize at about 10 eV cm(sup -3) s(sup -1), and the corresponding ion heating efficiencies peak near 1000 km at about 0.6%.

Fox, Jane L.↗

NASA NTRS1974

Neutral winds and electric fields from model studies using reduced ionograms

A relationship between the vertical component of the ion velocity and electron density profiles derived from reduced ionograms is developed. Methods for determining the horizontal components of the neutral winds and electric fields by using this relationship and making use of the variations of the inclinations and declinations of the earth's magnetic field are presented. The effects that electric fields have on the neutral wind calculations are estimated to be small but not second order. Seasonal and latitudinal variations of the calculated neutral winds are presented. From the calculated neutral winds a new set of neutral pressure gradients is determined. The new pressure gradients are compared with those generated from several static neutral atmospheric models. Sensitivity factors relating the pressure gradients and neutral winds are calculated and these indicate that mode coupling and harmonic generation are important to studies which assume linearized theories.

Baran, D. E.↗

DOE OSTI2025

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.

divertor↗

DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

PlasMod0 v.1.0

SAND2022-12609 O PlasMod0 is a 0-dimensional plasma modeling which approximates a plasma as a well-stirred reactor. The model uses two temperatures, one for electrons and one for ions and neutrals. This model was originally developed for laser produced plasma in single-particle aerosol mass spectrometers. As a result, there are some specialized capabilities that include spherical expansion into vacuum at a user-specified rate and photon reactions. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Lietz, Amanda↗

NASA NTRS1990

Determination of Comet Halley gas emission characteristics from mass loading of the solar wind

The velocity profile of the solar wind during Giotto's approach to Comet Halley is fitted with the mass loading produced by a simple model of the neutral particle distribution. The model is used to calculate the implanted ion flux at Giotto for any given time and position along the spacecraft trajectory. Comparing the added flux with Giotto's solar wind proton data from the inbound leg outside the bow shock, the ratio of the total mass-loaded ion flux to the solar wind flux at the spacecraft is computed. Hence, using equations given by Galeev et al. (1985) values of the solar wind velocity, u-infinity, far upstream from the comet are inferred. Imposing the condition that u-infinity should be as nearly as possible constant in time, values of gas production rate and the ratio between radial expansion velocity and ionization rate required to fit the model to data are derived. The values obtained are consistent with those derived by more direct methods.

Huddleston, D. E.↗

NASA NTRS1995

On Ba(+) production in the CRIT 2 experiment

Analysis of partical data from the CRIT 2 experiment, studying Alfven's critical ionization velocity (CIV) effect, shows that the density of newly created ions (presumably Ba(+) from the shaped-charge beam) is consistent with the increase in total plasma density measured by the independent RF plasma probe on board (Swenson et al., 1990) at the most active time period. We model this ion production using the measured electron flux data and the neutral barium model of Stenbaek-Nielsen et al. (1990a). To identify the main source mechanisms which may contribute most to the barium ionization, a simple model for the barium ion density at the payload location is developed based on Liouvilles theorem. We estimate that the electron impact ionization is responsible for 90% of the barium ion production observed by CRIT 2 in the first release and up to 45% in the second release. By employing a two-state approximation calculation (Rapp and Francis, 1962), the Ba-O(+) charge exchange cross section is found to range from about 2.0 X 10(exp -17) sq cm at a velocity of 4 km/s to 2.0 X 10(exp -15) sq cm at a velocity of 20 km/s. This result suggests that the Ba-O(+) charge exchange is probably dominant among all the non-CIV ionization processes. By considering the charge exchange process in our density model, the barrium ion densities are calculated for the two releases on CRIT II. The comparison between the model results and the observed data is found to be resonably consistent if the cross sections, as calculated above, are multiplied by 0.3 for the first release and 1.0 for the second release. Our result suggests that the charge exchange process could be the most important non-CIV ionization mechanism in the CRIT II experiment and it should be considered carefully case by case in CIV experiments.

Liou, K.↗

NASA NTRS2014

Titan's Plasma Environment : 3D Hybrid Kinetic Modeling of the TA Flyby and Comparison with CAPS-ELS and RPWS LP Observations

In this report we discuss the global plasma environment of the TA flyby from the perspective of 3D hybrid modeling. In our model the background, pickup, and ionospheric ions are considered as particles, whereas the electrons are described as a fluid. In homogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. Our modeling shows that mass loading of the background plasma (H(+), O(+)) by pick up ions H2(+), CH4(+) and N2(+) differs from theT9 encounter simulations when O(+) ions are not introduced into the background plasma. In our hybrid-modeling we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M=28 amu ions that were generated inside the ionosphere. Titan's interior is considered as a weakly conducting body. Special attention has been paid to comparing the simulated pickup ion density distribution with CAPS-ELS and with RPWS LP observations by the Cassini-Huygens spacecraft along theta trajectory. Our modeling shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures.

induced magnetospheres↗

NASA NTRS1991

On the position of the near-earth neutral sheet - A comparison of magnetic model predictions with empirical formulas

A detailed comparison of the near-earth neutral sheet position according to several approximations has been made. In particular, two empirical formulas for the neutral sheet location, recently derived from AMPTE/CCE data, have been compared with the corresponding positions derived from Tsyganenko's magnetic field models. Cylindrical coordinates referenced to the geodipole have been used, and the position of the neutral sheet has been identified by the reversal of the radial component of the magnetic field. The analysis reveals that, in the region X(GSM) = -9 to -5 earth radii, abs. value of Y(GSM) not greater than 5 earth radii, close agreement exists between the empirical formulas and the Tsyganenko models. Furthermore, the Tsyganenko models provide a representation for the neutral sheet position in the transition region between near-earth and distant tail models.

Peredo, Mauricio↗

DOE OSTI2022

Multi-qubit entanglement and algorithms on a neutral-atom quantum computer

Gate model quantum computers promise to solve currently intractable computational problems if they can be operated at scale with long coherence times and high fidelity logic. Neutral atom hyperfine qubits provide inherent scalability due to their identical characteristics, long coherence times, and ability to be trapped in dense multi-dimensional arrays. Combined with the strong entangling interactions provided by Rydberg states, all the necessary characteristics for quantum computation are available. Here we demonstrate several quantum algorithms on a programmable gate model neutral atom quantum computer in an architecture based on individual addressing of single atoms with tightly focused optical beams scanned across a two-dimensional array of qubits. Preparation of entangled Greenberger-Horne-Zeilinger (GHZ) states with up to 6 qubits, quantum phase estimation for a chemistry problem, and the Quantum Approximate Optimization Algorithm (QAOA) for the MaxCut graph problem are demonstrated. These results highlight the emergent capability of neutral atom qubit arrays for universal, programmable quantum computation, as well as preparation of non-classical states of use for quantum enhanced sensing.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2017

High-Latitude Thermosphere Neutral Density Response to Solar Wind Dynamic Pressure Enhancement

We examine the response of the thermosphere to the impact of solar wind dynamic pressure enhancements using observations and global magneto-hydrodynamics (MHD) simulations by the OpenGGCM model. Combining neutral density observations from the Challenging Minisatellite Payload (CHAMP) and the Gravity Recovery and Climate Experiment (GRACE) satellites with simultaneous Poynting flux measurements from the Defense Meteorological Satellite Program (DMSP) F16, we find that thermospheric density as well as downward Poynting flux intensified shortly after a sudden enhancement of the solar wind dynamic pressure. The intensification manifested mostly on the dayside high-latitude region with peak intensity in the vicinity of the noon and prenoon cusp. OpenGGCM modeling results show that the ionospheric Joule heating increased abruptly in response to the sudden enhancement of the dynamic pressure in the same region as the observed Poynting flux and neutral density enhancements. The modeling results show that the enhanced Joule heating coincides, both in time and location, with the appearance of a pair of high-latitude localized field-aligned currents (FACs) in the cusp region. The FACs intensified and extended azimuthally. Coincidental with the solar wind dynamic pressure enhancement, the y component of the interplanetary magnetic field (IMF) By became strongly positive and, in addition, had some large fluctuations. We explore the separate and combined effects of the dynamic pressure and IMF By perturbations, with specifically designed simulation experiments that isolate the effect of each solar wind parameter. We find that the dynamic pressure enhancement is the primary source for the Joule heating and neutral density enhancements, but the IMF By modulates the level of enhancement.

magneto-hydrodynamics↗

NASA NTRS1995

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗

NASA NTRS1989

The behavior of the electron density and temperatue at Millstone Hill during the equinox transition study September 1984

The ionospheric electron density and temperature variations is simulated during the equinox transition study in September 1984 and the results are compared with measurements made at Millstone Hill. The agreement between the modeled and measured electron density and temperature for the quiet day (18 September) is very good but there are large differences on the day of the storm (19 September). On the storm day, the measured electron density decreases by a factor of 1.7 over the previous day, while the model density actually increases slightly. The model failure is attributed to an inadequate increase in the ratio of atomic oxygen to molecular neutral densities in the MSIS neutral atmosphere model, for this particular storm. A factor of 3 to 5 increase in the molecular to atomic oxygen density ratio at 300 km is needed to explain the observed decrease in electron density. The effect of vibrationally excited N sub 2 on the electron density were studied and found to be small.

Richards, P. G.↗

NASA NTRS1989

The behavior of the electron density and temperature at Millstone Hill during the equinox transition study September 1984

Richards, P. G.↗

NASA NTRS1987

Effect of a localized minimum in equatorial field strength on resistive tearing instability in the geomagnetotail

A two-dimensional, resistive-MHD computer code is used to investigate the spontaneous reconnection of magnetotaillike configurations. The initial conditions adopted in the simulations are of two types: (1) in which the equatorial normal magnetic field component B(ze) declines monotonically down the tail, and (2) in which B(ze) exhibits a deep minimum in the near-earth plasma sheet. To represent the case where the earthward convection stops before the X line forms, zero-flow boundary conditions are imposed at the edges of the computational box. The initial configurations are in equilibrium and table within ideal MHD. The dynamic evolution of the system starts after the resistivity is turned on. The main results of these simulations basically support the neutral-line model of substorms and confirm Birn's (1980) computer studies. Spontaneous formation of an X-type neutral point and a single O-type plasmoid with strong tailward flow on the tailward side of the X point is found. The time interval from the turning on of the resistivity to the formation of a plasmoid is much shorter in the case where there is an initial deep minimum. A simple analytic calculation is also carried out to demonstrate why the configuration with a deep minimum is more susceptible to the development of the neutral point.

Hau, L.-N.↗

NASA NTRS2014

Anisotropic Solar Wind Sputtering of the Lunar Surface Induced by Crustal Magnetic Anomalies

The lunar exosphere is generated by several processes each of which generates neutral distributions with different spatial and temporal variability. Solar wind sputtering of the lunar surface is a major process for many regolith-derived species and typically generates neutral distributions with a cosine dependence on solar zenith angle. Complicating this picture are remanent crustal magnetic anomalies on the lunar surface, which decelerate and partially reflect the solar wind before it strikes the surface. We use Kaguya maps of solar wind reflection efficiencies, Lunar Prospector maps of crustal field strengths, and published neutral sputtering yields to calculate anisotropic solar wind sputtering maps. We feed these maps to a Monte Carlo neutral exospheric model to explore three-dimensional exospheric anisotropies and find that significant anisotropies should be present in the neutral exosphere depending on selenographic location and solar wind conditions. Better understanding of solar wind/crustal anomaly interactions could potentially improve our results.

exosphere↗

NASA NTRS1992

The simulated features of heliospheric cosmic-ray modulation with a time-dependent drift model. IV - The role of heliospheric neutral sheet deformation

Previous calculations with a time-dependent drift model revealed the model to be less successful in describing time-dependent modulation during periods of moderate to large solar activity. In this paper, it is argued that a major reason for this is that the previously used wavy heliospheric neutral sheet (HNS) description was based on an idealized HNS not subject to any spatial evolution while propagating radially outward. It is suggested that the deformation and compression of HNS wave peaks will lead to significant increases in the crossfield diffusion across these peaks (short-circuiting). The cosmic rays will effectively experience reduced tilt angles and therefore a reduction in the integrated HNS modulation effect between an observer and the heliospheric boundary. During periods of moderate to large solar activity these HNS deformation processes are progressively more frequent and should lead to a significant reduction in time-dependent modulation as predicted by drift models. Calculations done with radially propagating tilt angles that effectively decrease with radial distance give the expected reduction which improves the general description of modulation from 1987-1988.

Le Roux, J. A.↗

DOE OSTI2024

STREAM: Strategic Technology Roadmapping and Energy, Environmental, and Economic Analysis Model

This is the implementation of the framework for the planning of the technological makeup of the industrial sector. Motivated by the efforts to achieve carbon neutrality, the model of this framework modifies the portfolio of technologies over time for the sector of interest such that, - Net Present Value (NPV) is minimized - Constraint on Greenhouse Emissions (GHG), e.g. carbon dioxide, is satisfied - Demand of the underlying commodity is satisfied - The current implementation reflects a case study for the electric power sector. The for a given initial set of capacities of different vintages, the space of decisions include, - Retirement of the existing capacities. - Retrofitting the existing capacities to alternative characteristics. - Creation of new capacities from a technology portfolio. All the associated quantities with the deployments, e.g. CO2, heat requirement, etc.

Thiery, David↗