Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cosmology of Sub-MeV Dark Matter Freeze-In

Dark matter (DM) could be a relic of freeze-in through a light mediator, where the DM is produced by extremely feeble, IR-dominated processes in the thermal standard model plasma. In the simplest viable models with DM lighter than 1 MeV, the DM has a small effective electric charge and is born with a nonthermal phase-space distribution. This DM candidate would cause observable departures from standard cosmological evolution. In this work, we combine data from the cosmic microwave background (CMB), Lyman-α forest, quasar lensing, stellar streams, and Milky Way satellite abundances to set limits on freeze-in DM masses up to ~20 keV, with the exact constraint depending on whether the DM thermalizes in its own sector. Further, we perform forecasts for the CMB-S4 experiment, the Hydrogen Epoch of Reionization Array, and the Vera Rubin Observatory, finding that freeze-in DM masses up to ~80 keV can be explored.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Effect of an arcjet plume on satellite reflector performance

The effect of an arcjet plume on the performance of satellite reflector antennas is studied. The arcjet plume is modeled as a weakly ionized plasma. The spatial permittivity distribution of the plume is approximated using the measured electron density profile and a cold plasma model. Geometrical optics is applied to determine the ray paths as well as the transmitted fields through the inhomogeneous plume. The ray optics results are compared against several exact solutions for scattering from inhomogeneous dielectrics, and good agreement is observed for sufficiently large scatterer size. The far-field antenna patterns of the reflector in the presence of the plume are calculated from the transmitted ray fields using a ray-tube integration scheme. For arcjet prototypes in the 1-kW class, the plume effect on the antenna performance is small. As the electron density increases, the main beam and sidelobe level gradually degrade. The main beam also tends to squint away from the plume region.

Ling, Hao↗
NASA NTRS2018

The Importance of Nightside Magnetometer Observations for Electromagnetic Sounding of the Moon

Understanding the structure and composition of the lunar interior is a fundamental goal in furthering our knowledge of the formation and subsequent evolution of the Earth-Moon system. Among various methods, electromagnetic sounding is a valuable approach to constraining lunar interior structure. Recent analyses of plasma and field observations provide a wealth of understanding about the dynamics of the lunar plasma environment. To perform Time Domain EM (TDEM) Sounding at the Moon, the first step is to characterize the dynamic plasma environment, and to be able to isolate geophysically induced currents from concurrently present plasma currents. The TDEM Sounding transfer function method focuses on analysis of the nightside observations when the Moon is immersed in the solar wind. This method requires two simultaneous observations: an upstream reference measuring the pristine solar wind, and one downstream at or near the lunar surface. This method was last performed during Apollo and assumed the induced fields on the nightside of the Moon expand as in an undisturbed vacuum within the wake cavity. Our results indicate that EM sounding of airless bodies in the solar wind must be interpreted via self-consistent plasma models in order to untangle plasma and induced field contributions, with implications not only at the Moon but at all airless bodies exposed to the solar wind. Nightside TDEM sounding has the capability to advance the state of knowledge of the field of lunar science. This requires magnetometer operations to withstand the harsh conditions of the lunar night.

Haviland, H. Fuqua↗
DOE OSTI2021

Formulation of 8-moment plasma transport with application to the Nernst effect

8-moment plasma models using two different distribution functions are used to study the Nernst effect and heat transport in dense plasma. These models are presented in hyperbolic form in contrast to traditional parabolic systems derived from perturbing the distribution function, as in Braginskii [Rev. Plasma Phys. 1, 205 (1965)]. The hyperbolic moment formulation can be solved implicitly in time with straightforward and fast local solvers. The numerical implementation of 8-moment models with the relaxation method in the PERSEUS code is also presented. To test 8-moment PERSEUS compared to Braginskii’s transport equations, a verification test for the Nernst thermomagnetic wave by Velikovich et al. [Phys. Plasmas 26, 112702 (2019)] is performed that confirms the presence of the same physics, but with slight differences in the transport coefficients, which are tabulated in the limits of high and low magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗
DOE OSTI2022

Development of Ion Stopping Models for HED Plasmas Using Unified Self-Consistent Field Models and Self-Consistent Electron Distributions

We have implemented several corrections to the electronic stopping power model combining the RPA dielectric response formalism and local density approximation with electronic density distribution calculated in an average atom model. These modifications include strong collision correction, local field correction, electron binding energy correction, and the Barkas effect. The combined results bring the RPA-LDA stopping power in cold targets to closer agreements with experiments for a wide range of materials. The same method is then applied to the stopping of ions in warm dense plasmas. The computational framework developed during this project is publicly available on GitHub (https://github.com/dedx-erpa/dedx). Tabulated data for protons in cold target for common materials are located in the data/ subdirectory of the repository.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Robe Development for Electrical Conductivity Analysis in an Electron Gun Produced Helium Plasma

The use of magnetohydrodynamic (MHD) power conversion systems, potentially coupled with a fission power source, is currently being investigated as a driver for an advanced propulsion system, such as a plasma thruster. The efficiency of a MHD generator is strongly dependent on the electrical conductivity of the fluid that passes through the generator; power density increases as fluid conductivity increases. Although traditional MHD flows depend on thermal ionization to enhance the electrical conductivity, ionization due to nuclear interactions may achieve a comparable or improved conductivity enhancement while avoiding many of the limitations inherent to thermal ionization. Calculations suggest that nuclear-enhanced electrical conductivity increases as the neutron flux increases; conductivity of pure He-3 greater than 10 mho/m may be achievable if exposed to a flux greater than 10(exp 12) neutrons/cm2/s.) However, this remains to be demonstrated experimentally. An experimental facility has been constructed at the Propulsion Research Center at the NASA Marshall Space Flight Center, using helium as the test fluid. High energy electrons will be used to simulate the effects of neutron-induced ionization of helium gas to produce a plasma. These experiments will be focused on diagnosis of the plasma in a virtually static system; results will be applied to future tests with a MHD system. Initial experiments will utilize a 50 keV electron gun that can operate at up to a current of 200 micro A. Spreading the electron beam over a four inch diameter window results in an electron flux of 1.5x 10(exp 13) e/sq cm/s. The equivalent neutron flux that would produce the same ionization fraction in helium is 1x10(exp 12) n/sq cm/s. Experiments will simulate the neutron generated plasma modeled by Bitteker, which takes into account the products of thermal neutron absorption in He-3, and includes various ion species in estimating the conductivity of the resulting plasma. Several different probes will be designed and implemented to verify the plasma kinetics model. System parameters and estimated operating ranges are summarized. The predicted ionization fraction, electron density, and conductivity levels are provided in for an equivalent neutron flux of 1x10(exp 12) n/cm2/s. Understanding the complex plasma kinetics throughout a MHD channel is necessary to design an optimal power conversion system for space propulsion applications. The proposed experiments seek to fully characterize the helium plasma and to determine the reliability of each measurement technique, such that they may be applied to more advanced MHD studies. The expected value of each plasma parameter determined from theoretical models will be verified experimentally by several independent techniques to determine the most reliable method of obtaining each parameter. The results of these experiments will be presented in the final paper.

Bragg-Sitton, Shannon M.↗
NASA NTRS2006

Analyzing Electric Field Morphology Through Data-Model Comparisons of the GEM IM/S Assessment Challenge Events

The storm-time inner magnetospheric electric field morphology and dynamics are assessed by comparing numerical modeling results of the plasmasphere and ring current with many in situ and remote sensing data sets. Two magnetic storms are analyzed, April 22,2001 and October 21-23,2001, which are the events selected for the Geospace Environment Modeling (GEM) Inner Magnetosphere/Storms (IM/S) Assessment Challenge (IMSAC). The IMSAC seeks to quantify the accuracy of inner magnetospheric models as well as synthesize our understanding of this region. For each storm, the ring current-atmosphere interaction model (RAM) and the dynamic global core plasma model (DGCPM) were run together with various settings for the large-scale convection electric field and the nightside ionospheric conductance. DGCPM plasmaspheric parameters were compared with IMAGE-EUV plasmapause extractions and LANL-MPA plume locations and velocities. RAM parameters were compared with Dst*, LANL-MPA fluxes and moments, IMAGE-MENA images, and IMAGE-HENA images. Both qualitative and quantitative comparisons were made to determine the electric field morphology that allows the model results to best fit the plasma data at various times during these events. The simulations with self-consistent electric fields were, in general, better than those with prescribed field choices. This indicates that the time-dependent modulation of the inner magnetospheric electric fields by the nightside ionosphere is quite significant for accurate determination of these fields (and their effects). It was determined that a shielded Volland-Stern field description driven by the 3-hour Kp index yields accurate results much of the time, but can be quite inconsistent. The modified Mcllwain field description clearly lagged in overall accuracy compared to the other fields, but matched some data sets (like Dst*) quite well. The rankings between the simulations varied depending on the storm and the individual data sets, indicating that each field description did well for some place, time, and energy range during the events, as well as doing less well in other places, times, and energies. Several unresolved issues regarding the storm-time inner magnetospheric electric field are discussed.

Liemohn, Michael W.↗
DOE OSTI2026

Power and particle exhaust in the ST-E1 fusion power plant

Power exhaust challenges and potential solutions for a 5 m major radius, low-aspect ratio burning tokamak have been explored. 1D edge plasma models have been used to screen for access to detachment using short and long outer divertor legs in double and single null configurations, using Ar as the primary impurity and assuming tungsten plasma-facing components (PFCs). These show that detachment access can be accessed for all but the most conservative assumptions on scrape-off layer (SOL) width and power, but that trade-offs will be required between magnet engineering and the size of the acceptable window of as-yet uncertain plasma parameters. SOLPS-ITER was used to further model selected plasma scenarios, confirming that Ar seeding can be used to achieve dissipative divertor scenarios with peak deposited heat fluxes below 15 MWm -2 . Initial scoping of first wall loads and positioning of limiters has been carried out, showing the feasibility of protecting the breeding blanket wall during steady state without impeding tritium breeding. Initial PFC technology selection is also presented, identifying this as a critical area where further work is needed to find an attractive solution for helium-cooled PFCs that can handle high heat fluxes without excessive power requirements. Key questions and trade-offs for concept development have been identified, including: how to achieve high radiation for reduction of SOL power without core performance degradation; whether power exhaust can be well-controlled in a double null plasma; mechanical design and materials challenges of high-heat flux PFCs; and control of material erosion, redeposition and tritium retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Understanding the cold plasma synthesis of ammonia with model metal catalysts through plasma diagnostics (Final Report)

The industrial synthesis of ammonia, which amounts to over 200 million tons annually, is the most energy-intensive chemical process. Therefore, there is a critical need and increased interest in exploring less energetic routes to produce ammonia. Not only does ammonia have a direct impact on the food market, but also it has the potential as a fuel and hydrogen carrier. Recently, plasma catalysis has emerged as a promising alternative for synthesizing ammonia at mild (pressure, temperature and power) conditions. The key to this catalytic process is the synergy between the plasma and the catalyst, where the non-equilibrium plasma allows the generation of excited species, which recombine at the catalyst surface to form ammonia. However, our current understanding of this process is in its infancy. In this respect, model metals are ideal candidates for gaining a basic understanding of this reaction. Moreover, a major roadblock to rationally designing novel effective catalysts for plasma-assisted ammonia production is the need for fundamental aspects of this process. Through a comprehensive plan that integrates model metals as catalysts and world-class diagnostics, the proposed work aims to provide fundamental knowledge about the nature of reactive processes occurring during plasma-assisted catalysis and to demonstrate the selective production of ammonia, catalyzed by employing selected metals under non-thermal plasma conditions. Toward this goal, the central thrust of this proposal was to demonstrate that the synergy between plasma and rationally selected model metals will boost ammonia yields during plasma-assisted ammonia synthesis by delaying hydrogen recombination. Specifically, we aimed to (1) understand the formation and role of gas-phase active species such as NH, N 2 , N 2 + and Ha during the plasma-enhanced synthesis of ammonia through OES and FTIRAS using different reaction configurations: a) only plasma (non-packed DBD reactor) and b) packed DBD reactor with metal nanoparticles. This proposal was awarded/recommended with a run time on the PCRF facility FY20.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Studies of Westward Electrojets and Field-Aligned Currents in the Magnetotail During Substorms: Implications for Magnetic Field Models

This section outlines those tasks undertaken in the final year that contribute integrally to the overarching project goals. Fast, during the final year, it is important to note that the project benefited greatly with the addition of a Boston University graduate student, Ms. Karen Hirsch. Jointly, we made substantial progress on the development of and improvements to magnetotail magnetic field and plasma models. The ultimate aim of this specific task was to assess critically the utility of such models for mapping low-altitude phenomena into the magnetotail (and vice-versa). The bulk of this effort centered around the finite-width- magnetotail convection model developed by and described by Spence and Kivelson (J. Geophys. Res., 98, 15,487, 1993). This analytic, theoretical model specifies the bulk plasma characteristics of the magnetotail plasma sheet (number density, temperature, pressure) across the full width of the tail from the inner edge of the plasma sheet to lunar distances. Model outputs are specified by boundary conditions of the source particle populations as well as the magnetic and electric field configuration. During the reporting period, we modified this code such that it can be interfaced with the auroral particle precipitation model developed by Dr. Terry Onsager. Together, our models provide a simple analytic specification of the equatorial distribution of fields and plasma along with their low-altitude consequences. Specifically, we have built a simple, yet powerful tool which allows us to indirectly 'map' auroral precipitation signatures (VDIS, inverted-V's, etc.) measured by polar orbiting spacecraft in the ionosphere, to the magnetospheric equatorial plane. The combined models allow us to associate latitudinal gradients measured in the ion energy fluxes at low-altitudes with the large-scale pressure gradients in the equatorial plane. Given this global, quasi-static association, we can then make fairly strong statements regarding the location of discrete features in the context of the global picture. We reported on our initial study at national and international meetings and published the results of our predictions of the low-altitude signatures of the plasma sheet. In addition, the PI was invited to contribute a publication to the so-called 'Great Debate in Space Physics' series that is a feature of EOS. The topic was on the nature of magnetospheric substorms. Specific questions of the when and where a substorm occurs and the connection between the auroral and magnetospheric components were discussed in that paper. This paper therefore was derived exclusively from the research supported by this grant. Attachment: Empirical modeling of the quite time nightside magnetosphere.' 'CRRES observations of particle flux dropout event.' The what, where, when, and why of magnetospheric substorm triggers'. and 'Low altitude signature of the plasma sheet: model prediction of local time dependence'.

Spence, Harlan E.↗
NASA NTRS1982

A review of studies on ion thruster beam and charge-exchange plasmas

Various experimental and analytical studies of the primary beam and charge-exchange plasmas of ion thrusters are reviewed. The history of plasma beam research is recounted, emphasizing experiments on beam neutralization, expansion of the beam, and determination of beam parameters such as electron temperature, plasma density, and plasma potential. The development of modern electron bombardment ion thrusters is treated, detailing experimental results. Studies on charge-exchange plasma are discussed, showing results such as the relationship between neutralizer emission current and plasma beam potential, ion energies as a function of neutralizer bias, charge-exchange ion current collected by an axially moving Faraday cup-RPA for 8-cm and 30-cm ion thrusters, beam density and potential data from a 15-cm ion thruster, and charge-exchange ion flow around a 30-cm thruster. A 20-cm thruster electrical configuration is depicted and facility effects are discussed. Finally, plasma modeling is covered in detail for plasma beam and charge-exchange plasma.

Carruth, M. R., Jr.↗
NASA NTRS2006

Comparison of Two IRI plasmasphere Extensions with GPS-TEC Observations

Comparisons of two model results with Global Positioning System GPS-TEC measurements have been carried out for different latitudinal, solar activity, magnetic activity, diurnal and seasonal conditions. The models evaluated are the Global Core Plasma Model (GCPM-2000) and the IRI extension with Russian plasmasphere model (IRI*).Data of 23 observatories providing GPS-TEC and ionosonde data have been used. It is shown that IRI* plasmasphere electron density is greater than GCPM results by an order of magnitude at 6370 km altitude (one Earth's radius) with this excess growing to 2-3 orders of magnitude towards the GPS satellite altitude of 20000 km. Another source of model and GPS-TEC differences is a way of selection of the F2 layer peak parameters driving the models either with ITU-R (former CCIR) maps or ionosonde observations. Plasmasphere amendment to IRI improves accuracy of TEC model predictions because the plasmasphere contribution to the total TEC varies from 10% by daytime under quiet magnetic conditions to more than 50% by night under stormy conditions.

Gulyaeva, T. L.↗