Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗
DOE OSTI2024

Temperature and density dependence of Kr L-shell spectrum in hot dense plasmas

Kr L-shell spectroscopy modeling results are discussed in this paper, focusing on the n = 4 to n = 2 line transitions of Be- and Li-like Kr ions. Collisional radiative atomic kinetic and Stark-broadened spectral line shape calculations show electron temperature Te and density ne sensitivity in the spectrum. The combination of the Te dependence due to the relative intensity of Be-like to Li-like line emissions in the range from 1.5 to 3 keV and the ne sensitivity from the Stark broadening effect on the line shapes in the range from 5×1023 to 2×1024/ cc results in a spectrum that can be employed to diagnose Te and ne. Two different collisional radiative atomic kinetic models i.e., Prismspect [J. J. MacFarlane, et al., Int. Fusion Sci. Appl. Conf. Proc. 457 (2003)] and ABAKO [Florido, et al., PRE, 80, 056402 (2009)] produce similar results in level populations and spectra. In x-ray spectroscopy of implosion cores, this Kr L-shell spectrum may prove useful in an intermediate Te range in which Ar is too ionized for its K-shell to be of diagnostic value and Kr is not ionized enough for its K-shell emission to be useful.

Physics↗
DOE OSTI2024

Inverse aspect-ratio expanded tokamak equilibria

Following Greene et al. [Phys. Fluids 14, 671 (1971)] and Connor et al. [Phys. Plasmas 31, 577 (1988); Plasma Phys. Control. Fusion 34, 161 (1992); and Nucl. Fusion 33, 1533 (1993)], the Grad-Shafranov equation for an axisymmetric tokamak plasma equilibrium is solved via an expansion in the, supposedly small, inverse aspect-ratio of the plasma, ϵ. The displacements of equilibrium magnetic flux-surfaces due to plasma shaping are assumed to be $\mathcal{O}$(ϵ) smaller than the minor radii of the surfaces, but no other restriction is placed on the nature of the shaping. The solution of the Grad-Shafranov equation is matched to a vacuum solution that extends to infinity, and consists of an expansion in toroidal functions. The external poloidal magnetic field generated by a finite set of discrete external poloidal magnetic field-coils is calculated, and incorporated into the toroidal function expansion. In this manner, the shape of a large aspect-ratio tokamak plasma is directly related to the currents flowing in the external poloidal field-coils. Finally, a pedestal in the plasma pressure, and the associated spike in the bootstrap current, are incorporated into the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

DOE Plasma Science Center - Predictive Control of Plasma Kinetics: Multi-Phase and Bounded Systems (Final Report DE-SC0001939)

Low temperature plasmas (LTPs) are the plasmas of electron-volt (eV) physics and eV technologies. LTPs have characteristic electron temperatures of a few eV and fractional ionizations that are typically small. Since LTPs have electron temperatures commensurate with the threshold energies of excited states in neutral atoms and molecules, power transfer from electrons to these atoms and molecules efficiently produces activated species (e.g., radicals, excited states, photons). Acceleration of ions in the sheaths of LTPs to energies of tens to hundreds of eV enable activation of surface modifying processes – sputtering, etching, deposition. With such properties, LTPs are often and beneficially used in technological devices, ranging from etching and deposition in microelectronics and solar cell fabrication, to hardening of surgical instruments. LTPs harbor fundamental science issues that are intellectually challenging and rewarding. At the same time, there are enormous societal benefits that are enabled by LTPs. The entire present-day and future information technology infrastructure owes its very existence to LTPs. Renewable energy sources, such as solar cell arrays, could not be economically produced in the absence of deposition and etching by LTPs. In acknowledgment of the importance of LTPs, the DOE Office of Fusion Energy Sciences supported the Plasma Science Center for Predictive Control of Plasma Kinetics: Multi-Phase and Bounded Systems from August 2009 to July 2021, consisting of an initial 5-year grant, 3 funded extensions and 2 no-cost extensions. This is the final report of the Center in which the productivity of the Center is discussed in terms of publications, impact and personnel. High-lights of the research performed in the Center are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗
DOE OSTI2022

Full fluid moment modeling of rotating spokes in Penning-type configuration

Rotating spokes are observed in a partially magnetized plasma using a two-dimensional full fluid moment (FFM) model. In the present setup, where the radial electric field and plasma density gradient exist in opposite directions, it is observed that the spokes propagate in the direction of the diamagnetic drift and not the E × B drift. This is contrary to the modified Simon–Hoh instability, and the results suggest that the spokes can be driven by a strong diamagnetic drift. Different parameters, including magnetic field amplitude and physical domain size, influence the growth of the rotational instability as well as the dominant wave modes that arise. Furthermore, the propagation speed of the rotating spokes obtained from the FFM simulation are in good agreement with the observations in experimental and other computational work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Resistive drift wave turbulence and anomalous transport of multi-species plasma

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.

Physics↗
DOE OSTI2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Application of physiologically based pharmacokinetic modeling for sertraline dosing recommendations in pregnancy

Pregnancy is a period of significant change that impacts physiological and metabolic status leading to alterations in the disposition of drugs. Uncertainty in drug dosing in pregnancy can lead to suboptimal therapy, which can contribute to disease exacerbation. A few studies show there are increased dosing requirements for antidepressants in late pregnancy; however, the quantitative data to guide dose adjustments are sparse. We aimed to develop a physiologically based pharmacokinetic (PBPK) model that allows gestational-age dependent prediction of sertraline dosing in pregnancy. A minimal physiological model with defined gut, liver, plasma, and lumped placental-fetal compartments was constructed using the ordinary differential equation solver package, ‘mrgsolve’, in R. We extracted data from the literature to parameterize the model, including sertraline physicochemical properties, in vitro metabolism studies, disposition in nonpregnant women, and physiological changes during pregnancy. The model predicted the pharmacokinetic parameters from a clinical study with eight subjects for the second trimester and six subjects for the third trimester. Based on the model, gestational-dependent changes in physiology and metabolism account for increased clearance of sertraline (up to 143% at 40 weeks gestational age), potentially leading to under-dosing of pregnant women when nonpregnancy doses are used. The PBPK model was converted to a prototype web-based interactive dosing tool to demonstrate how the output of a PBPK model may translate into optimal sertraline dosing in pregnancy. Quantitative prediction of drug exposure using PBPK modeling in pregnancy will support clinically appropriate dosing and increase the therapeutic benefit for pregnant women.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Off-harmonic optical probing of high intensity laser plasma expansion dynamics in solid density hydrogen jets

Due to the non-linear nature of relativistic laser induced plasma processes, the development of laser-plasma accelerators requires precise numerical modeling. Especially high intensity laser-solid interactions are sensitive to the temporal laser rising edge and the predictive capability of simulations suffers from incomplete information on the plasma state at the onset of the relativistic interaction. Experimental diagnostics utilizing ultra-fast optical backlighters can help to ease this challenge by providing temporally resolved inside into the plasma density evolution. We present the successful implementation of an off-harmonic optical probe laser setup to investigate the interaction of a high-intensity laser at 5.4 x 10 21 W/cm 2 peak intensity with a solid-density cylindrical cryogenic hydrogen jet target of 5μm diameter as a target test bed. The temporal synchronization of pump and probe laser, spectral filtering and spectrally resolved data of the parasitic plasma self-emission are discussed. The probing technique mitigates detector saturation by self-emission and allowed to record a temporal scan of shadowgraphy data revealing details of the target ionization and expansion dynamics that were so far not accessible for the given laser intensity. Plasma expansion speeds of up to (2.3 ± 0.4)x 10 7 m/s followed by full target transparency at 1.4ps after the high intensity laser peak are observed. A three dimensional particle-in-cell simulation initiated with the diagnosed target pre-expansion at -0.2ps and post processed by ray tracing simulations supports the experimental observations and demonstrates the capability of time resolved optical diagnostics to provide quantitative input and feedback to the numerical treatment within the time frame of the relativistic laser-plasma interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Demonstrating the impact of ro-vibrationally excited H 2 on divertor detachment via population modelling

EDGE2D-EIRENE profiles characteristic for low recycling, high recycling and detached plasmas are postprocessed applying a vibrationally resolved collisional radiative (CR) model based on the Yacora solver to predict excited state densities and thus the molecular Lyman-, Werner- and Fulcher band emission intensities. These emission intensities are compared to predictions from the AMJUEL database standardly used in EIRENE. Lower Fulcher band and higher Lyman- and Werner band emission is predicted by the Yacora model in comparison to AMJUEL suggesting that inconsistent results from model based evaluation of measured spectra are expected depending on which model and molecular emission band is considered. The implications on detachment rate predictions are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Influence of triangularity on the plasma response to resonant magnetic perturbations

The effect of the axisymmetric plasma shape on the non-axisymmetric plasma response to resonant magnetic perturbations is investigated in experiment and modeling for the DIII-D, EAST and ASDEX Upgrade tokamaks. Systematically modeling the effect of the triangularity, whilst keeping other equilibrium quantities largely unchanged, reveals that the plasma response is strongly suppressed at high triangularity compared to that at low triangularity. This is validated through targeted comparison with experiments at DIII-D, EAST, and ASDEX Upgrade. DIII-D and EAST magnetic measurements are used to validate simulations, while at ASDEX Upgrade the plasma edge displacement is measured. Both experiments and modeling find a reduced magnetic plasma response on the high-field side at high triangularity across devices. Multi-modal analysis of the simulation results extracts the mode structure and applied spectrum dependence of each mode. The amplitude of the dominant mode reveals similar trends with the edge resonance and radial displacement near the X-point, which suggests that the multi-mode plasma response provides another way to understand the edge localized mode (ELM) control physics. The plasma response is strongly reduced at high triangularity compared to that at low triangularity, which implies different ELM control effects as shaping is varied. Furthermore, these findings indicate that the plasma shape should be taken into consideration when designing an RMP-ELM control strategy in experiment, and that predictive plasma response calculations can be used to maximize access to RMP-ELM control in future devices by maximizing the coupling between coils and the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Validation of realistic Monte Carlo plasma gamma-ray source on JET discharges

A novel modelling methodology has been developed for the creation of a realistic plasma gamma-ray source for Monte Carlo transport simulations in the tokamak JET. The methodology couples the TRANSP code for plasma transport calculations with the MCNP Monte Carlo particle transport code, thus connecting plasma physics with gamma-ray transport. Here, we present the validation of the developed source methodology by comparing calculated gamma-ray spectra with measurements performed at JET. The validation focuses on gamma-ray spectra measured by the tangential gamma-ray spectrometer during two JET three ion RF scenario discharges, performed in the JET 2019 deuterium experimental campaign. For validation the calculated plasma gamma-ray spectrum was combined with the neutron induced prompt gamma-ray background, originating in the vacuum vessel, and scaled to absolute values calculating the total number of plasma gamma-ray and neutron emitting reactions. The comparison between calculated and measured gamma-ray spectra shows good agreement with the shape of the calculated gamma-ray spectra matching that of measurements for both studied discharges. Moreover, the calculated absolute values of the gamma-ray spectra were of the same order of magnitude at the position of the gamma-ray detector located at the end of a long line-of-sight in a biological shield. Finally, the comparison has validated the developed plasma gamma-ray source methodology for MCNP photon transport calculations at JET. The validation provides a basis for the developed plasma gamma-ray source to be used as a support for the development of future tokamaks such as DEMO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-Principles Theory of the Relativistic Magnetic Reconnection Rate in Astrophysical Pair Plasmas

In this paper, we develop a first-principles model for the relativistic magnetic reconnection rate in strongly magnetized pair plasmas. By considering the energy budget and required current density near the x-line, we analytically show that in the magnetically dominated relativistic regime, the x-line thermal pressure is significantly lower than the upstream magnetic pressure due to the extreme energy needed to sustain the current density, consistent with kinetic simulations. This causes the upstream magnetic field lines to collapse in, producing the open outflow geometry which enables fast reconnection. The result is important for understanding a wide range of extreme astrophysical environments, where fast reconnection has been evoked to explain observations such as transient flares and nonthermal particle signatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

K-shell radiation and bright spot characteristics of high-energy-density Fe-Cr-Ni plasmas influenced by X-pinch load geometry

K-shell radiation is diagnostically advantageous for the study of high-energy-density plasmas. Here, the current work investigates the influence of load geometry on K-shell x-rays in Z-pinch plasmas produced on the 1-MA Zebra generator. Stainless steel (Fe: 69%, Cr: 20%, Ni: 9%) X-pinch wire loads are fielded with an interwire angle of 31° (small-angle) or 62.5° (large-angle), respectively and studied using various x-ray diode, imaging, and spectroscopic diagnostics. The large-angle geometry produced large, individual x-ray-emitting bright spots (> 3 keV) at the wire cross-point averaging areas ≥ 1 mm 2 , with noteworthy soft x-ray (> 0.75 keV) bursts and radiation yields ≤ 14.6 kJ. The small-angle geometry produced multiple x-ray-emitting bright spots extending the pinch axis with smaller average sizes (≤ 0.5 mm 2 ), consistent with stronger current, increased radiation yield (≥ 15.5 kJ), and notable hard x-ray (> 9 keV) bursts, which peak with largest cumulative bright spot size. Spectroscopic analysis is performed with non-LTE collisional-radiative modeling, indicating a cold, nonthermal plasma region radiating from neutral to Ne-like ions, a hot, thermal region radiating from Fe and Cr He- and Li-like ions, and an intermediate region radiating Kα satellites from Li-like to O-like ions. Modeling of the nonthermal and satellite features imply a fractional hot electron abundance ~0.1% and ~0.5% for large- and small-angle geometries, respectively. Relative intensity analysis is performed on Heα and Kα lines, revealing intensity ratios that deviate from the original wire composition. Optical depth estimates allude to optically thick, thermal K-shell plasmas with densities consistent with non-LTE modeling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Plasma heating and current sheet structure in anti-parallel magnetic reconnection

A theoretical model and an analytic theory of current sheet structure are presented for understanding anti-parallel driven magnetic reconnection in 2-1/2 dimension in collisionless plasmas. The theoretical model provides formulation to compute the current sheet y-profiles by specifying the profiles of electron and ion flow velocities V ex (x,y) and V ix (x,y). The current sheet solutions depend on the plasma density n in , merging magnetic field B 0 , ion velocity v i , and electron velocity ve in the upstream and the S evz =V ez /V dz parameter where V ez is the electron velocity accelerated by the reconnection electric field E z in the electron orbit meandering region, V dz ≃cE y /B x is the $\overrightarrow{E}$ × $\overrightarrow{B}$ drift velocity as electrons enter the orbit meandering region, B x is the merging magnetic field, and E y is the electrostatic electric field. With simplifying assumptions on the y-profiles of V ex and V ix , we have also developed an analytic theory of the current sheet structure. Analytic expressions for the anomalous resistivity, the electrostatic potential drop, and the maximum E y amplitude E max are obtained. The analytic results agree reasonably well with both the particle-in-cell simulation results and the numerical solutions of the theoretical model. The ions energy gain due to the potential drop is $∝B^2_0/n_{in}$. The electron energy gain is $∝(B^2_0 /8πn_{in})S_{evz}$. The $B^2_0/n_{in}$ scaling of the average ion and electron energy gains are consistent with laboratory experiments and space plasma observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗