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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Three-Dimensional Multiscale MHD Model of Cometary Plasma Environments

First results of a three-dimensional multiscale MHD model of the interaction of an expanding cometary atmosphere with the magnetized solar wind are presented. The model starts with a supersonic and super-Alfvenic solar wind far upstream of the comet (25 Gm upstream of the nucleus) with arbitrary interplanetary magnetic field orientation. The solar wind is continuously mass loaded with cometary ions originating from a 10-km size nucleus. The effects of photoionization, electron impact ionization, recombination, and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our new multiscale upwind scalar conservation laws-type numerical technique (MUSCL). We have named this the multiscale adaptive upwind scheme for MHD (MAUS-MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. The main findings are the following: (1) Mass loading decelerates the solar wind flow upstream of the weak cometary shock wave (M approximately equals 2, M(sub A) approximately equals 2), which forms at a subsolar standoff distance of about 0.35 Gm. (2) A cometary plasma cavity is formed at around 3 x 10(exp 3) km from the nucleus. Inside this cavity the plasma expands outward due to the frictional interaction between ions and neutrals. On the nightside this plasma cavity considerably narrows and a relatively fast and dense cometary plasma beam is ejected into the tail. (3) Inside the plasma cavity a teardrop-shaped inner shock is formed, which is terminated by a Mach disk on the nightside. Only the region inside the inner shock is the 'true' diamagnetic cavity. (4) The model predicts four distinct current systems in the inner coma: the density peak current, the cavity boundary current, the inner shock current, and finally the cross-tail current. (5) The calculated plasma parameters (magnetic field, plasma density, speed, and temperature) are in very good agreement with published Giotto observations.

Gombosi, Tamas I.↗
NASA NTRS1987

A boundary layer model for magnetospheric substorms

An alternative framework for understanding magnetospheric substorm activity is presented. It is argued that observations of magnetic field and plasma flow variations in the magnetotail can be explained in terms of the passage of the plasma sheet boundary layer over the satellite detecting the tail signatures. It is shown that field-aligned currents and particle acceleration processes on magnetic field lines threading the ionospheric Harang discontinuity lead to the distinctive particle and field signatures observed in the magnetotail during substorms. It is demonstrated that edge effects of field-aligned currents associated with the westward traveling surge can lead to the negative B(z) perturbations observed in the tail that are presently attributed to observations made on the anti-earthward side of a near-earth neutral line. Finally, it is shown that the model can provide a physical explanation of both the driven system and the loading-unloading system whose combined effects provide the observed substorm perturbation pattern in the magnetosphere and ionosphere.

Rostoker, Gordon↗
NASA NTRS1972

The responses of the thermosphere due to a geomagnetic storm: A MHD model

A magnetohydrodynamics theory was used to study the dynamic response of the neutral atmosphere to a geomagnetic storm. A full set of magnetohydrodynamic equations appropriate for the present problem is derived and their various orders of approximation are discussed in some detail. In order to demonstrate the usefulness of this theoretical model, the May 1967 geomagnetic storm data were used in the resulting set of nonlinear, time dependent, partial differential magnetohydrodynamic equations to calculate variations of the thermosphere due to the storm. The numerical results are presented for wind speeds, electric field strength, and amount of joule heating at a constant altitude for the data recorded. Data show that the strongest thermospheric responses are at the polar region becoming weaker in the equatorial region. This may lead to the speculation that a thermospheric wave is generated in the polar region due to the geomagnetic storm which propagates towards the equator.

Wu, S. T.↗
NASA NTRS1989

Silicon chemistry in interstellar clouds

Interstellar SiO was discovered shortly after CO but it has been detected mainly in high density and high temperature regions associated with outflow sources. A new model of interstellar silicon chemistry that explains the lack of SiO detections in cold clouds is presented which contains an exponential temperature dependence for the SiO abundance. A key aspect of the model is the sensitivity of SiO production by neutral silicon reactions to density and temperature, which arises from the dependence of the rate coefficients on the population of the excited fine structure levels of the silicon atom. This effect was originally pointed out in the context of neutral reactions of carbon and oxygen by Graff, who noted that the leading term in neutral atom-molecule interactions involves the quadrupole moment of the atom. Similar to the case of carbon, the requirement that Si has a quadrupole moment requires population of the J = 1 level, which lies 111K above the J = 0 ground state and has a critical density n(cr) equal to or greater than 10(6)/cu cm. The SiO abundance then has a temperature dependence proportional to exp(-111/T) and a quadratic density dependence for n less than n(cr). As part of the explanation of the lack of SiO detections at low temperatures and densities, this model also emphasizes the small efficiencies of the production routes and the correspondingly long times needed to reach equilibrium. Measurements of the abundance of SiO, in conjunction with theory, can provide information on the physical properties of interstellar clouds such as the abundances of oxygen bearing molecules and the depletion of interstellar silicon.

Langer, William D.↗
NASA NTRS1980

A thermalized ion explosion model for high energy sputtering and track registration

A velocity spectrum of neutral sputtered particles as well as a low resolution mass spectrum of sputtered molecular ions was measured for 4.74 MeV F-19(+2) incident of UF4. The velocity spectrum is dramatically different from spectra taken with low energy (keV) bombarding ions, and is shown to be consistent with a hot plasma of atoms in thermal equilibrium inside the target. A thermalized ion explosion model is proposed for high energy sputtering which is expected to describe track formation in dielectric materials. The model is shown to be consistent with the observed total sputtering yield and the dependence of the yield on the primary ionization rate of the incident ion.

Seiberling, L. E.↗
NASA NTRS2011

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRS1973

A study of the electric field in an open magnetospheric model

Recently, Svalgaard and Heppner reported two separate features of the polar electromagnetic field that correlate with the dawn-dusk component of the interplanetary magnetic field. This work attempts to explain these findings in terms of properties of the open magnetosphere. The topology and qualitative properties of the open magnetosphere are first studied by means of a simple model, consisting of a dipole in a constant field. Many such properties are found to depend on the separation line, a curve connecting neutral points and separating different field line regimes. In the simple model it turns out that the electric field in the central polar cap tends to point from dawn to dusk for a wide variety of external fields, but, near the boundary of the polar cap, electric equipotentials are deformed into crescents.

Stern, D. P.↗
DOE OSTI2023

Two-photon double ionization with finite pulses: Application of the virtual sequential model to helium

As a step toward the full ab initio description of two-photon double ionization processes, we present a finite-pulse version of the virtual-sequential model for polyelectronic atoms. The model relies on the ab initio description of the single ionization scattering states of both the neutral and ionized target system. As a proof of principle and a benchmark, the model is applied to the helium atom using the newstock atomic photoionization code. The results of angularly integrated observables, which are in excellent agreement with existing time-dependent Schrödinger equation (TDSE) simulations, show how the model is able to capture the role of electron correlation in the nonsequential regime, and the influence of autoionizing states in the sequential regime, at a comparatively modest computational cost. The model also reproduces the two-particle interference with ultrashort pulses, which is within reach of current experimental technologies. Furthermore, the model shows the modulation of the joint energy distribution in the vicinity of autoionizing states, which can be probed with extreme-ultraviolet pulses of duration much longer than the characteristic lifetime of the resonance. Furthermore, the formalism discussed here applies also to polyelectronic atoms and molecules, thus opening a window on nonsequential and sequential double ionization in these more complex systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

DarkNews: a Python-based event generator for heavy neutral lepton production in neutrino-nucleus scattering

We introduce DarkNews, a lightweight Python-based Monte-Carlo generator for beyond-the-Standard-Model neutrino-nucleus scattering. The generator handles the production and decay of heavy neutral leptons via additional vector or scalar mediators, as well as through transition magnetic moments. DarkNews samples pre-computed neutrino-nucleus upscattering cross sections and heavy neutrino decay rates to produce dilepton and single-photon events in accelerator neutrino experiments. We present two case studies with differential distributions for models that can explain the MiniBooNE excess. The aim of this code is to aid the neutrino theory and experimental communities in performing searches and sensitivity studies for new particles produced in neutrino upscattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1990

A multigenerational Monte Carlo model of a comet coma

A steady state multigenerational Monte Carlo model was developed to describe the distribution of a neutral coma species which was produced after several photodissociation steps. This distribution can be compared with either line profiles derived from long slit spectroscopic observations or narrow band imaging, and be used to determine the elemental abundances present in the comet coma. A Monte Carlo model was chosen due to limitations in standard forms of more analytic models, such as those of Haser and Festou. In contrast to the analytical models, the Monte Carlo method can be modified to include new physical parameters.

Ferro, Anthony J.↗
DOE OSTI2022

Neutral-current neutrino cross section and expected supernova signals for 40 Ar from a three-fold increase in the magnetic dipole strength

In view of the great interest in liquid argon neutrino detectors, the 40 Ar(γ,γ') 40 Ar reaction was revisited to guide a calculation of the neutral current neutrino cross section at supernova energies. Using the nuclear resonance fluorescence technique with a monoenergetic, 99% linearly polarized photon beam, in this work we report a three-fold increase in magnetic dipole strength at around 10 MeV in 40 Ar. Based on shell-model calculations, and using the experimentally identified transitions, the neutral current neutrino cross sections for low-energy reactions on 40 Ar are calculated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Magnetic measurements in model experiments

Plasma parameters for modeling of the interaction between the solar wind and the earth's magnetic field are selected on the limited modeling principle (plasma velocity 3 x 10 to the 7th power cm/sec, density 10 to the 13th power cm/3, electron temperature approximately 15 eV). Measurements were made with magnetic probes capable of developing all 3 components inside the magnetic field. The resulting magnetic field charts indicate that a neutral layer forms on the night side in the model experiments. A magnetosphere free of plasma streams is established on the daytime side. A shock wave forms on the boundary of the magnetosphere, its width c/omega sub 0 in accordance with theory, where omega sub 0 is the total Langmuir width.

Dubinin, E. M.↗
DOE OSTI2023

Systems and methods for controlling an industrial asset in the presence of a cyber-attack

Systems and methods are provided for the control of an industrial asset, such as a power generating asset. Accordingly, a cyber-attack model predicts a plurality of operational impacts on the industrial asset resulting from a plurality of potential cyber-attacks. The cyber-attack model also predicts a corresponding plurality of potential mitigation responses. In operation, a cyber-attack impacting at least one component of the industrial asset is detected via the cyber-attack neutralization module and a protected operational impact of the cyber-attack is identified based on the cyber-attack model. The cyber-attack neutralization module selects at least one mitigation response of the plurality of mitigation responses based on the predicted operational impact and an operating state of the industrial asset is altered based on the selected mitigation response.

D'Amato, Fernando Javier↗
DOE OSTI2023

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Quantitative simulation of a magnetospheric substorm. I - Model logic and overview

Results of a comprehensive computer simulation of the behavior of the earth's inner magnetosphere during a substorm-type event are reported. It is pointed out that the computer model self-consistently computes electric fields, currents, and plasma distributions and velocities in the inner-magnetosphere/ionosphere system; parallel electric fields and ionospheric neutral winds, however, are not included. The basic equations of the model are derived, and the inputs are described. An overview of the results is also given. The first appendix contains derivations of general, useful laws of bounce-averaged gradient, curvature, and E x B drifts in a plasma with isotropic pitch angle distribution. The second appendix gives a description of the numerical method used in the simulation.

Harel, M.↗
NASA NTRS2000

Identification of Gas Phase PAHs in Absorption Towards Protostellar Sources

The infrared emission bands (also known as the UIR bands.) have recently been observed in absorption at 3.25 micrometers in the ices surrounding a few proto-stellar objects at 11.2 micrometers in MonR2, and at 6.2 micrometers towards two sources near the galactic center. The UIR bands have been observed in emission for many years, but identifying these bands has proven to be both difficult and contentious as no one has yet found a single material that provides a good match to the features. However, most investigators agree that some form of carbon-based material with aromatic bonds is the most likely candidate, and many arguments favor free molecules (polycyclic aromatic hydrocarbons, PAHs) as the carriers of at least the narrow emission bands. Since the emission arises not from a single molecule but from a family of molecules, identifying which PAHs are contributing to the infrared emission bands is difficult. The identification is further complicated by the fact that the emission at short wavelengths is dominated by small molecules while at long wavelengths it is dominated by large molecules. Thus, for example, the emission at 3.3 micrometers is from a different mix of molecules than those which produce the 11.2 micrometer band. To complicate matters further, the molecular mix includes both neutral and ionic species. In absorption, the same mixture of molecules contributes at all wavelengths and the molecules should be neutral, potentially simplifying comparisons with lab data. Also, absorption strengths measured in the lab are directly applicable to interstellar absorption bands without the need to model an emission spectrum of an unknown mixture of ionized and neutral PAHs. In this paper we show that a mixture of argon matrix isolated PAH molecules can reproduce the 3.25 micrometers absorption band seen in the ISO SWS spectra of four embedded Infrared sources, S140 IRS1, AFGL 2591, Elias 29, and AFGL 989. In section 2 we describe the ISO SWS data analysis and in section 3 discuss the results.

Bregman, Jesse D.↗
NASA NTRS2016

Investigating Global Ion and Neutral Atom Populations with IBEX and Voyager

The main objective of this project was to investigate pickup ion (PUI) production in the solar wind and heliosheath (the region between the termination shock and the heliopause) and compute the distributed energetic neutral atom fluxes throughout the helioshpere. The simulations were constrained by comparing the model output against observations from Ulysses, New Horizons, Voyager 1 and 2, and IBEX space probes. As evidenced by the number of peer reviewed journal publications resulting from the project (13 plus three submitted) and their citation rate (156 citations over three years), the project has made a lasting contribution to the field. The outcome is a significant improvement of our understanding of the pickup ion production and distribution in the distant heliosphere. The team has accomplished the entire set of tasks A-H set forth in the proposal. Namely, the transport modeling framework has been augmented with two populations of pickup ions (PUIs), the boundary conditions for the plasma and interstellar neutral hydrogen were verified against Ulysses and New Horizons PUI and an optimal set of velocity diffusion parameters established. The multi-component fluxes of PUIs were computed and isotropic velocity distributions generated for each cell in the computer simulation that covered the heliosphere from 1.5 AU to the heliopause. The distributions were carefully compared with in situ measurements at 3 AU (Ulysses), 12 AU (New Horizons), and 80-90 AU (Voyager 1 and 2) as well as those inferred from ENA fluxes measured by Cassini and IBEX (Wu et al., 2016). Some examples of modeldata comparison are shown in Figure 1. We have used coupled MHD-plasma and kinetic-neutral code to investigate the likely range of plasma and magnetic field parameters in the local interstellar medium (LISM), based on the assumption that the shape of the IBEX ribbon could be used to determine the orientation of the interstellar magnetic field. While the magnetic field is believed to be oriented toward the center of the ribbon, constraining its strength requires comparing the model-predicted angular diameter and circularity of the ribbon with the observations. The study, published in Heerikhuisen et al. (2014), found that the most likely range for the LISM magnetic field strength is between 0.2 and 0.3 nT, which is less than previously thought. Figure 2 shows the IBEX data (left) and compares it to the simulation with a 0.2 nT interstellar magnetic field (center) and a 0.4 nT (right).

Florinski, Vladimir↗
NASA NTRS1984

Ion thruster performance model

A model of ion thruster performance is developed for high flux density, cusped magnetic field thruster designs. This model is formulated in terms of the average energy required to produce an ion in the discharge chamber plasma and the fraction of these ions that are extracted to form the beam. The direct loss of high energy (primary) electrons from the plasma to the anode is shown to have a major effect on thruster performance. The model provides simple algebraic equations enabling one to calculate the beam ion energy cost, the average discharge chamber plasma ion energy cost, the primary electron density, the primary-to-Maxwellian electron density ratio and the Maxwellian electron temperature. Experiments indicate that the model correctly predicts the variation in plasma ion energy cost for changes in propellant gas (Ar, Kr and Xe), grid transparency to neutral atoms, beam extraction area, discharge voltage, and discharge chamber wall temperature. The model and experiments indicate that thruster performance may be described in terms of only four thruster configuration dependent parameters and two operating parameters. The model also suggests that improved performance should be exhibited by thruster designs which extract a large fraction of the ions produced in the discharge chamber, which have good primary electron and neutral atom containment and which operate at high propellant flow rates.

Brophy, J. R.↗