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

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma production and composition from hypervelocity impacts on solar cell cover glass

Satellites are subjected to detrimental effects from the space environment, including impacts by orbital debris and meteoroids that can cause mechanical or electrical damage. While mechanical damage is well studied, electrical effects resulting from hypervelocity impact (HVI) remain poorly understood. In HVIs, where the projectile speed exceeds the speed of sound in the target material, the impact energy ionizes the material near the surface, creating a plasma that can cause radio frequency emission. This emission is the main source of electrical anomalies in satellites, but the mechanism behind its creation is still unknown. Ground-based experiments have been used to characterize HVI plasma, including the empirical power law relations that describe the impact charge produced as a function of impactor mass and speed. Quantifying these power law relations gives us a better understanding of the species formed in the impact which can eventually lead to plasma evolution models and help us first, understand the mechanism behind radio frequency emission and second, ways to prevent its damage to the satellites. In this paper, we interpret the results of the HVI experiments on solar cell coverglass targets and identify the potential chemical or molecular reactions that form the impact ion species. This is done by looking simultaneously at multiple data sets for different species within the same impact configuration, and using a new technique for quantifying the power law relations for each emitted species. We use time of flight analysis (TOF) to detect the charged particle compositions and determine their initial speeds. The results suggest that the charge production power law is species specific, with different formation mechanisms involved in a single HVI event, some of which are associated with dissociation and some with ionization of a molecule. Furthermore, our results also show that electrons have different power laws that depend on the spacecraft surface charge. TOF measurements of ions from impact events on solar cell cover glass shows the presence of SiO + and SiO – 2 as the dominant species present in plasmas on positively and negatively charged targets, respectively, as well as the formation of FeO + .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics-Based Real-time Analysis and Control to Achieve Transient-Free Operations for the ITER Era (Final Report)

This is the final report for the early career award for Prof. Egemen Kolemen that aimed to develop the physics and technical basis of real-time transient – i.e. disruptions and Edge Localized Modes (ELMs) – avoidance and feedback control at the DIII-D nuclear fusion facility for the ITER era. During this project, we reconstructed high quality real-time DIII-D equilibria using kinetic information from the motional Stark effect (MSE), Thomson Scattering (TS), and Charge Exchange Recombination (CER) spectroscopy. These equilibria were analyzed in real-time to determine the low-n ideal MHD stability and the pedestal stability thresholds, and to obtain plasma response models. These calculations were used to understand which physical process evolves to lead to the onset of disruption in ITER-relevant scenarios at DIII-D, and how disruption can be predicted using stability analysis as opposed to the current methods of correlation of various signals, which are descriptive but lack predictive power. In the disruption studies, we focused on tearing modes, since these constitute the vast majority of disruptions in this regime. With this knowledge, real-time disruption avoidance algorithms were developed and tested. The enhanced situation awareness and physics insight were used in advanced real-time feedback control of pedestal and core profiles, which allowed the pursuit and sustainment of transient-free high-confinement regimes that cannot be achieved via feed-forward design based on in-between shot analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Helicon edge heating in Proto-MPEX: Azimuthally resolved RF modeling

The Material Plasma Exposure eXperiment (MPEX) is a linear device under construction at ORNL. MPEX employs a variety of heating systems, including a helicon antenna, which is the subject of this work. Lacking a resonant absorption mechanism, as tokamaks have, modeling the RF power absorption in such antennas requires detailed consideration of cold collisional physics. Previous modeling of the heating in this antenna was done exclusively in 2D, assuming axisymmetry and a single azimuthal mode number m = 1⁠. In this study, we introduce a fully 3D model. By including the full 3D physics, with its self-consistent geometry-imposed azimuthal spectrum, we are able to reproduce the experimentally observed transition from edge heating to core heating.

Tierens, Wouter [Oak Ridge National Laboratory (OR↗
DOE OSTI2023

Self-consistent integrated modeling of combined hybrid discharge-laser produced plasma devices for extreme ultraviolet metrology

Discharge- and laser-produced plasma (DLPP) devices are being used as light sources for extreme ultraviolet (EUV) generation. A key challenge for both, DPP and LPP, is achieving sufficient brightness to support the throughput requirements of nanometrology tools. To simulate the environment of a hybrid DLPP device and optimize the EUV output, we have developed an integrated HEIGHTS-DLPP computer simulation package. The package integrates simulation of two evolving plasmas (DPP and LPP) and includes modeling of a set of integrated self-consistent processes: external power source and plasma energy balance, plasma resistive magnetohydrodynamics, plasma heat conduction, detailed radiation transport (RT), and laser absorption and refraction. We simulated and optimized DLPP devices using Xe gas as a target material. We synchronized the external circuit parameters, chamber gas parameters, and laser beam temporal and spatial profiles to achieve maximum EUV output. The full 3D Monte Carlo scheme was integrated for detailed RT and EUV output calculations in Xe using more than 3600 spectral groups. The modeling results are in good agreement with Julich Forschungszentrum experimental data. Theoretical models, developed and integrated into the HEIGHTS package, showed wide capabilities and flexibility. In conclusion, the models and package can be used for optimization of the experimental parameters and settings, investigation of DLPP devices with complex design, analyzing the impact of integrated spatial effects and working timeline arrangement on the final EUV output, and EUV source size, shape, and angular distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effects of sodium and calcium chloride ionic stresses on model yeast membranes revealed by molecular dynamics simulation

As efforts to move a renewable economy grow, it will be necessary to make use of microbial conversion strategies for the production of novel materials or the upgrading of waste to high-value products. One critical technical challenge currently limiting waste upgrading remains the difficulty in obtaining single-pot conversion techniques where physical, chemical, and biological conversion are performed in a single step. To overcome this challenge, a detailed understanding of how different stresses impact microbial membrane stability will be necessary. In this work, using all-atom molecular dynamics simulations, we examine the impacts of moderate concentrations of NaCl and CaCl 2 on a model yeast plasma membrane. Weak, though statistically significant, changes in membrane morphology and dynamics functions are observed that are consistent with swelling and stiffening. Additionally, an examination of the ion-lipid contacts and the behavior of water at the water-membrane interface suggests that the impacts of these common salts may, in part, be mediated through changes to water-membrane hydrogen-bonding and hydration water dynamics.

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

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Extending the single-fluid solvability conditions for more general plasma systems

We extend the single-fluid solvability conditions to include plasma systems at arbitrary $\beta$ with arbitrary flows and external forcing terms. This treatment includes both the isotropic and the anisotropic cases. The generalized conditions that result can be used to generate certain classes of single-fluid equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Metal-induced oxidative stress and human plasma protein oxidation after SARS-CoV-2 infection

Pathogenesis of COVID-19 by SARS-CoV-2 resulted in a global pandemic and public health emergency in 2020. Viral infection can induce oxidative stress through reactive oxygen species (ROS). Inflammation and environmental stress are major sources of oxidative stress after infection. Micronutrients such as iron, copper, zinc, and manganese play various roles in human tissues and their imbalance in blood can impact immune responses against pathogens including SARS CoV-2. We hypothesized that alteration of free metal ions during infection and metal-catalyzed oxidation plays a critical role towards pathogenesis after infection. We analyzed convalescent and hospitalized COVID-19 patient plasma using orthogonal analytical techniques to determine redox active metal concentrations, overall protein oxidation, oxidative modifications, and protein levels via proteomics to understand the consequences of metal-induced oxidative stress in COVID-19 plasma proteins. Metal analysis using ICP-MS showed significantly greater concentrations of copper in COVID-19 plasma compared to healthy controls. We demonstrate significantly greater total protein carbonylation, other oxidative modifications, and deamidation of plasma proteins in COVID-19 plasma compared to healthy controls. Proteomics analysis showed that levels of redox active proteins including hemoglobulin were elevated in COVID-19 plasma. Molecular modeling concurred with potential interactions between iron binding proteins and SARS CoV-2 surface proteins. Overall, increased levels of redox active metals and protein oxidation indicate that oxidative stress-induced protein oxidation in COVID-19 may be a consequence of the interactions of SARS-CoV-2 proteins with host cell metal binding proteins resulting in altered cellular homeostasis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Extended particle absorber for efficient modeling of intense laser-solid interactions

Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation and diagnostic application of Kr K-shell emission in magnetized liner inertial fusion experiments at Z

In a series of Magnetized Liner Inertial Fusion (MagLIF) experiments performed at the Z pulsed power accelerator of Sandia National Laboratories, beryllium liners filled with deuterium gas pressures in the 4–8 atm range and a tracer amount of krypton were imploded. At the collapse of the cylindrical implosion, temperatures in the 1–3 keV range and atom number densities of ∼1023 cm−3 were expected. The plasma was magnetized with a 10 T axial magnetic field. Krypton was added to the fuel for diagnosing implosion plasma conditions. Krypton K-shell line emission was recorded with the CRITR time-integrated transmission crystal x-ray spectrometer. The observation shows n = 2 to n = 1 line emissions in B-, Be-, Li-, and He-like Kr ions and is characteristic of the highest electron temperatures achieved in the thermonuclear plasma. Detailed modeling of the krypton atomic kinetics and radiation physics permits us to interpret the composite spectral feature, and it demonstrates that the spectrum is temperature sensitive. We discuss temperatures extracted from the krypton data analysis for experiments performed with several filling pressures.

Clapp, J. T. (ORCID:0000000219928146)↗
DOE OSTI2022

Monte Carlo N-Particle forward modeling for density reconstruction of double shell capsule radiographs

In the Double Shell Inertial Confinement Fusion concept, characterizing the shape asymmetry of imploding metal shells is vital for understanding energy-efficient compression and radiative losses of the thermonuclear fuel. The Monte Carlo N-Particle MCNP® code forward models radiography of Double Shell capsule implosions using the Advanced Radiographic Capability at the National Ignition Facility. A procedure is developed for using MCNP to reconstruct density profiles from the radiograph image intensity. For a given Double Shell imploding target geometry, MCNP radiographs predict image contrast, which can help guide experimental design. Finally, in future work, the calculated MCNP synthetic radiographs will be compared with experimental radiographs to determine the radial and azimuthal density profiles of the Double Shell capsules.

47 OTHER INSTRUMENTATION↗