SEARCH · Engineering Papers

Results for “Theoretical plasma physics”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Universality of scaled particle spectra in ultrarelativistic heavy-ion collisions

We study the transverse momentum spectra of identified particles in ultrarelativistic collisions of large and small collision systems. In order to isolate information contained in the momentum dependence, we propose to scale the spectra by the total particle number and mean transverse momentum, global quantities which are already well studied. We observe an interesting, nearly universal, centrality-independent shape in the scaled spectra, similar to scalings that have been studied previously. This scaling behavior breaks down at large transverse momentum and for very small systems, such as those produced in p − p collisions. We perform hybrid hydrodynamic simulations and show that, in these simulations, a centrality-independent shape is a consequence of an event-by-event independence. Here, our results motivate further theoretical and experimental investigations of the regime of validity of this scaling phenomenon and their physical interpretation at different collision energies and systems.

Collective flow↗

DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Internal measurements of electromagnetic geodesic acoustic mode (GAM) in EAST plasmas

Velocity, density, and magnetic fluctuations of the geodesic acoustic mode (GAM) have been measured using the Doppler backscattering system, Faraday-effect polarimeter-interferometer, and external pick-up coils in the Experimental Advanced Superconducting Tokamak. Simultaneous measurements of density and velocity fluctuations at the midplane and top of plasmas demonstrate that m = 1 density fluctuations are quantitatively balanced by the compression of perpendicular flow fluctuations. Furthermore, internal magnetic fluctuations associated with GAM have now been directly measured by laser-based Faraday-effect polarimetry for the first time. Line-averaged magnetic fluctuations (up to 16 Gauss, B̃¯R,GAMBT∼0.066%) are significantly larger than those extrapolated from edge coils (a few Gauss) and that magnetic fluctuations increase with β. The observed discrepancy between finite β theory and experimental data indicates the need for further theoretical investigations.

Physics↗

DOE OSTI2024

Effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability

Abstract The lack of understanding of the azimuthal instability and the resulting electron anomalous transport limits further improvement of Hall thrusters. Compared to theoretical and experimental approaches, the numerical particle-in-cell (PIC) simulation is a suitable and powerful tool, which has been widely applied to investigate the azimuthal instability, and great progress has been made in the past decades. However, PIC simulations are intrinsically computationally expensive, and it is realized that the Hall thruster azimuthal instability has a three dimensional nature. Therefore, massive 3D PIC simulation must be carried out to completely reveal the mechanism of the instability. In this paper, the effect of plasma initialization on 3D PIC simulation of Hall thruster azimuthal instability is studied as a starting point. It is found that by initializing with ion density and velocity fitting functions to the steady-state simulation results, a faster convergence can be obtained and the computational time can be reduced by about 1.5 times. Typical fitting functions of ion density, drifting velocity, and temperature are given, and the influence of different initialization profiles is presented.

Physics↗

DOE OSTI2021

On the violation of the zeroth law of turbulence in space plasmas

The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. This behaviour, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation – no matter how small – can thermalise the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetised plasma turbulence such as that recently observed by the Parker Solar Probe spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is ‘imbalanced’ – when the large-scale energy input is dominated by Alfvénic perturbations propagating in one direction (the most common situation in space plasmas) – nonlinear conservation laws imply the existence of a ‘barrier’ at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong resemblance to those observed in situ , exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales. The effect thus offers a possible solution to the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The existence of the ‘barrier’ also suggests that, how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.

Physics↗

DOE OSTI2024

Mining experimental magnetized liner inertial fusion data: Trends in stagnation morphology

In magnetized liner inertial fusion (MagLIF), a cylindrical liner filled with fusion fuel is imploded with the goal of producing a one-dimensional plasma column at thermonuclear conditions. However, structures attributed to three-dimensional effects are observed in self-emission x-ray images. Despite this, the impact of many experimental inputs on the column morphology has not been characterized. We demonstrate the use of a linear regression analysis to explore correlations between morphology and a wide variety of experimental inputs across 57 MagLIF experiments. Results indicate the possibility of several unexplored effects. For example, we demonstrate that increasing the initial magnetic field correlates with improved stability. Although intuitively expected, this has never been quantitatively assessed in integrated MagLIF experiments. We also demonstrate that azimuthal drive asymmetries resulting from the geometry of the “current return can” appear to measurably impact the morphology. In conjunction with several counterintuitive null results, we expect the observed correlations will encourage further experimental, theoretical, and simulation-based studies. Finally, we note that the method used in this work is general and may be applied to explore not only correlations between input conditions and morphology but also with other experimentally measured quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Highly selective Si 3 N 4 /SiO 2 etching using an NF 3 /N 2 /O 2 /H 2 remote plasma. II. Surface reaction mechanism

Developing processes for highly selective etching of silicon nitride (Si 3 N 4 ) with respect to silicon dioxide (SiO 2 ) is a major priority for semiconductor fabrication processing. In this paper and in Paper I [Volynets et al., J. Vac. Sci. Technol. A 38, 023007 (2020)], mechanisms are discussed for highly selective Si3N4 etching in a remote plasma based on experimental and theoretical investigations. The Si 3 N 4 /SiO 2 etch selectivity of up to 380 was experimentally produced using a remote plasma sustained in NF 3 /N 2 /O 2 /H 2 mixtures. A selectivity strongly depends on the flow rate of H 2 , an effect attributed to the formation of HF molecules in vibrationally excited states that accelerate etching reactions. Based on experimental measurements and zero-dimensional plasma simulations, an analytical etching model was developed for etch rates as a function of process parameters. Reaction rates and sticking coefficients were provided by quantum chemistry models and also fitted to the experimental results. Etch rates from the analytical model show good agreement with the experimental results and demonstrate why certain etchants accelerate or inhibit the etch process. In particular, the modeling shows the important role of HF molecules in the first vibrationally excited state [HF( v = 1)] in achieving high Si 3 N 4 /SiO 2 selectivity.

42 ENGINEERING↗

DOE OSTI2026

Mixed stochastic-deterministic density functional theoretic decomposition of Kubo–Greenwood conductivities in the projector augmented wave formalism

Pairing the accuracy of the Kohn–Sham density-functional framework with the efficiency of a stochastic algorithmic approach, mixed stochastic-deterministic density functional theory (mDFT) achieves a favorable computational scaling with system sizes and electronic temperatures. We employ the recently developed mDFT formalism to investigate the dynamic charge-transport properties of systems in the warm dense matter regime. The optical conductivity spectra are computed for single- and multi- component mixtures of carbon, hydrogen, and beryllium using two complementary approaches: Kubo–Greenwood in the mDFT picture and real-time time-dependent mDFT. We further devise a decomposition of the Onsager coefficients leading up to the Kubo–Greenwood spectra to exhibit contributions from the deterministic, stochastic, and mixed electronic state transitions at different incident photon energies.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Pressure probe study of the water relations of Phycomyces blakesleeanus sporangiophores

The physical characteristics which govern the water relations of the giant-celled sporangiophore of Phycomyces blakesleeanus were measured with the pressure probe technique and with nanoliter osmometry. These properties are important because they govern water uptake associated with cell growth and because they may influence expansion of the sporangiophore wall. Turgor pressure ranged from 1.1 to 6.6 bars (mean = 4.1 bars), and was the same for stage I and stage IV sporangiophores. Sporangiophore osmotic pressure averaged 11.5 bars. From the difference between cell osmotic pressure and turgor pressure, the average water potential of the sporangiophore was calculated to be about -7.4 bars. When sporangiophores were submerged under water, turgor remained nearly constant. We propose that the low cell turgor pressure is due to solutes in the cell wall solution, i.e., between the cuticle and the plasma membrane. Membrane hydraulic conductivity averaged 4.6 x 10(-6) cm s-1 bar-1, and was significantly greater in stage I sporangiophores than in stage IV sporangiophores. Contrary to previous reports, the sporangiophore is separated from the supporting mycelium by septa which prevent bulk volume flow between the two regions. The presence of a wall compartment between the cuticle and the plasma membrane results in anomalous osmosis during pressure clamp measurements. This behavior arises because of changes in solute concentration as water moves into or out of the wall compartment surrounding the sporangiophore. Theoretical analysis shows how the equations governing transient water flow are altered by the characteristics of the cell wall compartment.

NASA Discipline Number 40-30↗

DOE OSTI2024

Multimessenger measurements of the static structure of shock-compressed liquid silicon at 100 GPa

The ionic structure of high-pressure, high-temperature fluids is a challenging theoretical problem with applications to planetary interiors and fusion capsules. Here we report a multimessenger platform using velocimetry and angularly and spectrally resolved x-ray scattering to measure the thermodynamic conditions and ion structure factor of materials at extreme pressures. We document the pressure, density, and temperature of shocked silicon near 100 GPa with uncertainties of 6%, 2%, and 20%, respectively. The measurements are sufficient to distinguish between and rule out some ion screening models. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Thinking Bayesian for plasma physicists

Bayesian statistics offers a powerful technique for plasma physicists to infer knowledge from the heterogeneous data types encountered. To explain this power, a simple example, Gaussian Process Regression, and the application of Bayesian statistics to inverse problems are explained. The likelihood is the key distribution because it contains the data model, or theoretic predictions, of the desired quantities. By using prior knowledge, the distribution of the inferred quantities of interest based on the data given can be inferred. Because it is a distribution of inferred quantities given the data and not a single prediction, uncertainty quantification is a natural consequence of Bayesian statistics. The benefits of machine learning in developing surrogate models for solving inverse problems are discussed, as well as progress in quantitatively understanding the errors that such a model introduces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Diethylenetriamine-functionalized graphene oxide: Insights into ion adsorption and applications in rare earth element separation

The growing demand for critical minerals and materials requires atom- and energy-efficient, selective separations to overcome the challenges posed by the similar chemical and physical properties of the rare earth elements (REEs) and their low concentrations in unconventional domestic feedstocks. Here, in this study, we developed diethylenetriamine-functionalized graphene oxide (DETA-GO) as a membrane material for REE adsorption and separation. Synthesis conditions were optimized to maximize nitrogen incorporation while also preserving GO dispersibility for facile membrane fabrication. We investigated the mechanism of amine functionalization, the nitrogen-bonding configurations, the organization of the interlayer transport channels, and the resulting effects on ion and water transport for REE separations. Neat-GO and DETA-GO multilayer laminate membranes were fabricated by vacuum filtration onto polymer supports. To investigate the effects of amine functionalization, the membranes were characterized using scanning electron microscopy, Raman, Fourier transform infrared, and X-ray photoelectron spectroscopy, as well as grazing-incidence X-ray diffraction measurements. Ion permeation experiments with representative lanthanum (La 3+ ) and ytterbium (Yb 3+ ) solutions revealed enhanced ion adsorption and water transport through DETA-GO membranes compared to neat-GO. The strong affinity of the membranes for multivalent REEs was also validated with conductivity and inductively coupled plasma mass spectrometry measurements. Atomistic insight into the role of amine functionalization in modulating nanochannel architecture and long-term stability, optimizing adsorption sites, and regulating REE and water transport was obtained using classical molecular dynamics simulations. Collectively, our joint experimental and theoretical study demonstrates the potential of DETA-GO membranes for selective REE separations, offering insights into ion-binding mechanisms, water-transport properties, and nanochannel optimization for the recovery of critical materials from aqueous feedstocks.

Adsorbent↗

DOE OSTI2020

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Observation of impurity accumulation and its compatibility with high plasma performance in W7-X

At the W7-X stellarator, the bolometer system has measured an intensive radiation zone in the inner plasma region (at a normalized radius $ρ$ ~ 0.3–0.4) in the hydrogen plasma generated by electron cyclotron resonance heating; it differs from the normal plasma radiation distribution with an edge-localized emission zone. Spectroscopic diagnostics have recorded high-Z elements such as iron. This phenomenon happens in the plasma phases after gas supply turn-off, which results in all impurity relevant diagnostic signals increasing for several seconds. Despite the enhancement of the core radiation, the plasma energy confinement is improved. A transport analysis shows that this impurity radiation behavior is associated with a low diffusion coefficient ($D$~ 0.02 m 2 s -1 ) and a reversal of the convection around the radial position of the emission peak, which, under normal conditions, separates the zones of outward convection in the central (|$V$| ~ 0.1 m s -1 ) and inward convection in the outer region (|$V$| ~ 0.3 m s -1 ). An impurity accumulation around this radial position has been identified. The transport coefficients obtained are comparable with the theoretical predictions of collisional impurity transport. In the plasma phases studied, both impurity and energy confinement are enhanced. The mechanism responsible for the improvement is believed to be a reduction of micro-instabilities associated with the observed steepening of the density profile, initiated by a low edge plasma density (<1.0 × 10 19 m -3 ) after switching off the gas fueling. The normalized temperature and density gradients fulfil the condition for the suppression of ITG turbulence.

suppression of ITG turbulence↗

DOE OSTI2021

Azimuthal dependence of two-particle transverse momentum current correlations

Two-particle transverse momentum correlation functions are a powerful technique for understanding the dynamics of relativistic heavy-ion collisions. Among these, the transverse momentum correlator G 2 (Δη,Δφ) is of particular interest for its potential sensitivity to the shear viscosity per unit of entropy density η/s of the quark-gluon plasma formed in heavy-ion collisions. We use the UrQMD, AMPT, and EPOS models for Au–Au at √ s NN = 200 GeV and Pb–Pb at √ s NN = 2760 GeV to investigate the long range azimuthal dependence of G 2 (Δη,Δφ), and explore its utility to constrain η/s based on charged particle correlations. We find that the three models yield quantitatively distinct transverse momentum Fourier harmonics coefficients $a$ $^{pT}_{n}$. We also observe these coefficients exhibit a significant dependence on η/s in the context of the AMPT model. These observations suggest that exhaustive measurements of the dependence of G 2 (Δφ) with collision energy, system size, collision centrality, in particular, offer the potential to distinguish between different theoretical models and their underlying assumptions. Exhaustive analyses of G 2 (Δφ) obtained in large and small systems should also be instrumental in establishing new constraints for precise extraction of η/s.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Accommodating a hexagonal ζ- phase Mn 2 N film on a cubic MgO (001) substrate

The thin films of c -plane dominated hexagonal ζ- phase Mn 2 N were successfully grown on a cubic MgO (001) substrate directly using plasma-assisted molecular beam epitaxy. The surface was comprehensively studied through experimental and theoretical approaches. Reflection high energy electron diffraction revealed two pseudo-cubic domains along [100] MgO and [110] MgO , and one hexagonal domain, which is 30° apart from the [100] MgO direction. Scanning tunneling microscopy was used to image and resolve the high-quality surface, revealing a distorted hexagonal surface structure. Furthermore, atomic resolution of the hexagonal domain with a 2 × 2 reconstructed surface is presented. For c-plane surface, theoretical investigations were carried out using first principles studies to determine the surface formation energy for various reconstructed surfaces. The theoretical study shows that the nitrogen terminated 2 × 2 structure with a manganese adatom on the surface as the most stable reconstruction that reproduces the experimental results. A corresponding simulated scanning tunneling microscopy model is also presented, providing strong support for the experimentally observed in-plane lattice structures. Furthermore, the manganese: nitrogen stoichiometry within the bulk and surface is in good agreement with the expected ratio of 2:1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Observation of the Electromagnetic Field Effect via Charge-Dependent Directed Flow in Heavy-Ion Collisions at the Relativistic Heavy Ion Collider

The deconfined quark-gluon plasma (QGP) created in relativistic heavy-ion collisions enables the exploration of the fundamental properties of matter under extreme conditions. Noncentral collisions can produce strong magnetic fields on the order of 10 18 G , which offers a probe into the electrical conductivity of the QGP. In particular, quarks and antiquarks carry opposite charges and receive contrary electromagnetic forces that alter their momenta. This phenomenon can be manifested in the collective motion of final-state particles, specifically in the rapidity-odd directed flow, denoted as v 1 ( y ) . Here, we present the charge-dependent measurements of d v 1 / d y near midrapidities for π ± , K ± , and p ( p ¯ ) in Au + Au and isobar ( Ru 44 96 + Ru 44 96 and Zr 40 96 + Zr 40 96 ) collisions at s NN = 200 GeV , and in Au + Au collisions at 27 GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. The combined dependence of the v 1 signal on collision system, particle species, and collision centrality can be qualitatively and semiquantitatively understood as several effects on constituent quarks. While the results in central events can be explained by the u and d quarks transported from initial-state nuclei, those in peripheral events reveal the impacts of the electromagnetic field on the QGP. Our data put valuable constraints on the electrical conductivity of the QGP in theoretical calculations. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗