Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Space 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

Simulated signatures of ignition

Ignition on the National Ignition Facility (NIF) provides a novel opportunity to evaluate past data to identify signatures of capsule failure mechanisms. We have used new simulations of high-yield implosions as well as some from past studies in order to identify unique signatures of different ignition failure mechanisms: jetting due to the presence of voids or defects, jetting due to the capsule fill tube, interfacial mixing due to instabilities or due to plasma transport, radiative cooling due to the presence of contaminant in the hot spot, long-wavelength drive asymmetry, and preheat. Many of these failure mechanisms exhibit unique trajectories that can be distinguished through variations in experimental observables such as neutron yield, down-scattered ratio (DSR), and burn width. Our simulations include capsules using both plastic and high-density carbon ablators and span all high-yield designs considered since the beginning of the National Ignition Campaign in 2011. We observe that the variability in trajectories through the space of neutron yield, DSR, and burn width varies little across capsule design yet are unique to the failure mechanism. The experimental trajectories are most consistent with simulated preheat and jetting due to voids and defects, which are the only failure mechanisms that are indistinguishable in our analysis. This suggests that improvements to capsule compression due to improved capsule quality or reduced preheat have played a primary role in enabling high yields on NIF. Furthermore, our analysis suggests that further improvements have the potential to increase yields further.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of oxidation and impurities in lithium surfaces on the emitting wall plasma sheath

The use of lithium as a surface coating in fusion devices improves the plasma performance, but the change in the wall properties affects the secondary electron emission properties of the material. Lithium oxidizes easily, which drives the emission yield well above unity. We present here simulations demonstrating the change in the sheath structure from monotonic to the nonmonotonic space-charge limited sheath using an energy-dependent data-driven emission model, which self-consistently captures both secondary emission and backscattering populations. Increased secondary electron emission from the material has ramifications for the degradation and erosion of the wall. The results show that the oxidation leads to an increased electron energy flux into the wall and a reduced ion energy flux. The net transfer of energy to the surface is significantly greater for the oxidized case than for the pure lithium case. High backscattering rates of low-energy particles lead to a high re-emission rate at the wall.

Fusion reactors↗

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗

Extension of Virtual Test Bed Advanced Modeling and Simulation Capabilities for Fusion Energy

The National Reactor Innovation Center (NRIC) was established to accelerate the deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing reactor experiments. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test beds. It is a collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the mission to accelerate the deployment of advanced reactors by facilitating the adoption of advanced modeling and simulation (M&S) tools developed by the DOE. This mission has been carried out by the VTB since 2020 by hosting and featuring dozens of advanced fission nuclear reactor models developed by national laboratories and academia. The charter of the NRIC’s definition of advanced reactors also includes fusion nuclear reactors. As the tools developed by the NEAMS program are increasingly used for modeling fusion energy devices, there is an increasing need to host fusion reactor models on the VTB repository. The VTB will be extended in 2024 to support fusion energy modeling and simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probabilistic locked mode predictor in the presence of a resistive wall and finite island saturation in tokamaks

We present a framework for estimating the probability of locking to an error field in a rotating tokamak plasma. This leverages machine learning methods trained on data from a mode-locking model, including an error field, resistive magnetohydrodynamics modeling of the plasma, a resistive wall, and an external vacuum region, leading to a fifth-order ordinary differential equation (ODE) system. It is an extension of the model without a resistive wall introduced by Akçay et al. [Phys. Plasmas 28, 082106 (2021)]. Tearing mode saturation by a finite island width is also modeled. We vary three pairs of control parameters in our studies: the momentum source plus either the error field, the tearing stability index, or the island saturation term. The order parameters are the time-asymptotic values of the five ODE variables. Normalization of them reduces the system to 2D and facilitates the classification into locked (L) or unlocked (U) states, as illustrated by Akçay et al., [Phys. Plasmas 28, 082106 (2021)]. This classification splits the control space into three regions: L̂, with only L states; Û, with only U states; and a hysteresis (hysteretic) region Ĥ, with both L and U states. In regions L̂ and Û, the cubic equation of torque balance yields one real root. Region Ĥ has three roots, allowing bifurcations between the L and U states. The classification of the ODE solutions into L/U is used to estimate the locking probability, conditional on the pair of the control parameters, using a neural network. We also explore estimating the locking probability for a sparse dataset, using a transfer learning method based on a dense model dataset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗

Scalings of laser-driven MeV electron sources with relativistic multi-picosecond laser pulses

We investigate the generation of MeV electrons in relativistic multi-picosecond laser–solid interactions, using two-dimensional particle-in-cell simulations. Scanning a parameter space of intensities from 10 18 to 10 20 W/cm 2 and pulse durations between 100 fs and 10 ps, we discuss scalings of the effective hot electron temperature and MeV electron dose relevant to applications, as well as the direct laser acceleration mechanism that drives electrons to energies much higher than the ponderomotive potential of the laser pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conceptual design of HTS magnets for fusion nuclear science facility

Second-generation high temperature superconductors (HTS) are available for producing >20 T at the magnet bore compared to 13–16 T for lower temperature superconducting (LTS) toroidal field magnets proposed in recent fusion energy systems studies (FESS) of Fusion Nuclear Science Facility (FNSF). HTS may enable higher fusion power density and smaller device size. High current density cables of multi-layered REBCO tapes have achieved >10 kA at 4–20 K operation in short sample tests for fusion. High current density cables are required for engineering design of FNSF to allow space for interior plasma components. High current density HTS magnets are particularly attractive in reducing the size of a fusion device, beneficial for compact tokamaks, due to their space constraints. Successful HTS magnet development may enable the design of smaller and cheaper fusion pilot plants with a mission of demonstrating net electricity. It may also offer significant cost and performance advantages in non-fusion applicants such as nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI). Furthermore, we developed HTS magnet design concepts for a compact FNSF radial build in order to define the coil size, winding pack mechanical loading and engineering requirements. Partnering with vendors in the US, PPPL is also testing high current cable prototypes aiming at enabling low cost cable technology toward 100 A/mm 2 engineering current density over the winding pack desired in high field model coil development for compact fusion pilot plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High yield polar direct drive fusion neutron sources at the National Ignition Facility

Polar direct drive neutron source experiments were performed at the National Ignition Facility showing substantial improvement in total neutron yield and efficiency of conversion of laser energy to fusion output. Plastic capsules 3–4 mm in diameter were filled with 1.5 mg/cc of deuterium–tritium (DT) fuel and imploded with laser beam pointing and defocus designed to compensate for polar asymmetry introduced by the facility beam entrance angles. Radiation-hydrodynamics simulations were employed to optimize the multi-dimensional laser and target parameter space, within facility and target fabrication constraints. Ensembles of 1D simulations tuned to match the outputs of early shots in the series were used to design subsequent shots in the series. This allowed the later shots to be designed based on empirically motivated sensitivities to laser and target input parameters, while eliminating the need to explicitly model phenomena such as hydrodynamic instabilities and nonlinear laser–plasma interactions. Additionally, one experiment with a 3.0 mm diameter CH capsule produced 13.6 kJ (4.81 × 10 15 DT neutrons) from a laser input below the NIF optics damage threshold at 585 kJ, 328 TW. Two experiments with 4.0 mm capsules produced 31.3 and 33.6 kJ of fusion output (1.11 × 10 16 and 1.19 × 10 16 DT neutrons) with 1.10 MJ, 390 TW and 1.26 MJ, 425 TW of laser input, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Grassroots Network and Community Roadmap for Interconnected Autonomous Science Laboratories for Accelerated Discovery

Scientific discovery is being revolutionized by AI and autonomous systems, yet current autonomous laboratories remain isolated islands unable to collaborate across institutions. We present the Autonomous Interconnected Science Lab Ecosystem (AISLE), a grassroots network transforming fragmented capabilities into a unified system that shorten the path from ideation to innovation to impact and accelerates discovery from decades to months. AISLE addresses five critical dimensions: (1) cross-institutional equipment orchestration, (2) intelligent data management with FAIR compliance, (3) AI-agent driven orchestration grounded in scientific principles, (4) interoperable agent communication interfaces, and (5) AI/ML-integrated scientific education. By connecting autonomous agents across institutional boundaries, autonomous science can unlock research spaces inaccessible to traditional approaches while democratizing cutting-edge technologies. This paradigm shift toward collaborative autonomous science promises breakthroughs in sustainable energy, materials development, and public health.

Ferreira da Silva, Rafael [Oak Ridge National Labo↗

Evolution of an arched magnetized laboratory plasma in a sheared magnetic field

Arched magnetized structures are a common occurrence in space and laboratory plasmas. Results from a laboratory experiment on spatio-temporal evolution of an arched magnetized plasma ( $\beta \approx 10^{-3}$ , Lundquist number $\approx 10^{4}$ , plasma radius/ion gyroradius $\approx 20$ ) in a sheared magnetic configuration are presented. The experiment is designed to model conditions relevant to the formation and destabilization of similar structures in the solar atmosphere. The magnitude of a nearly horizontal overlying magnetic field was varied to study its effects on the writhe and twist of the arched plasma. In addition, the direction of the guiding magnetic field along the arch was varied to investigate its role in the formation of either forward- or reverse-S shaped plasma structures. The electrical current in the arched plasma was well below the current required to make it kink unstable. A significant increase in the writhe of the arched plasma was observed with larger magnitudes of overlying magnetic field. A forward-S shaped arched plasma was observed for a guiding magnetic field oriented nearly antiparallel to the initial arched plasma current, while the parallel orientation yielded the reverse-S shaped arched plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lasing efficiency of krypton ions in the (8 – 14)-nm band upon pulsed laser excitation

The emission spectra of krypton plasma in the range of 8 – 14 nm upon excitation of a pulsed gas jet by 1.06-mm Nd : YAG laser radiation with a pulse energy of 0.85 J, pulse duration of 5.2 ns, and repetition rate of 10 Hz are investigated. The krypton emission spectrum is a wide (8 – 14 nm) band, peaking at 10.3 nm, which is formed by a series of much narrower lines. The observed lines are identified, and the fraction of laser pulse energy converted into the (8 – 14)-nm emission band and emitted into half-space (2π sr) is determined. The maximum conversion efficiency is found to be 21 %. The expected throughputs of lithographic systems with sources based on Sn, Xe, and Kr ions for different wavelengths, corresponding to the emission peaks of ions of these materials, are compared. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energetic Electron Transport in Magnetized Plasma with Magnetic Islands (Final Technical Report)

The main scientific goal of this project was to investigate the interactions between energetic electrons and magnetic islands in magnetized plasmas from lab to space. The proposed plan of addressed the following specific objectives: (1) Determine if energetic electrons are trapped by magnetic islands. Specifically, test the hypothesis that energetic electrons are trapped near island O-points due to stable orbits and de-confined near island X-points due to X-point tangles. (2) Determine how electron transport changes with island width and location. Specifically, test the hypothesis that electrons are accelerated by contracting magnetic islands through a Fermi acceleration process. (3) Assess how electron transport is affected by island dynamics, including island rotation, bifurcation, and/or island overlap. Specifically, test the hypothesis that electrons are deconfined and possibly accelerated during island rotations, bifurcations, and overlap due to stochastization of the magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing the impact of oxygen negative ions on the self-organized pattern in 1 atm DC glow with liquid anode

In an atmospheric DC glow discharge with a liquid anode, the plasma anode glow attached to the grounded liquid surface under certain conditions self-organizes into coherent patterns. Optical emission spectroscopy revealed that the emission consists primarily of the second positive system of nitrogen, N 2 (C-B), whose excitation energy is low and sensitive to changes in the electron energy distribution. In addition to electrons, negative ions can accumulate in the anode sheath and affect the local space charge. It has been speculated that these negative ions play a role in pattern formation at the anode surface. In this work, the role of oxygen negative ions was explored. It was found that the formation of anode patterns requires at least a 7% volume fraction of oxygen in the ambient gas. Results showed that O 2 - is the dominant negative ion species in atmospheric DC glow discharge, with a density of ~10 12 cm -3 . While the presence of oxygen appears to be crucial for pattern formation, this study indicated that patterns still formed without geometric changes even when 62% of negative ions in the plasma were detached by a laser. This suggests that negative ions do not support the patterns, while oxygen's heating effect may induce instability at the anode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion-Exchange Kinetics of Alpha-Zirconium Phosphate Nanoplatelets for Application in Targeted Alpha Therapy

The loading and kinetics of uptake for the inorganic ion-exchange material, α-zirconium phosphate, Zr(HPO 4 ) 2 •H 2 O, were investigated with ions relevant to targeted alpha therapy including natural Cs + (as a surrogate for 221 Fr), natural Bi 3+ (as a surrogate for 213 Bi), and natural La 3+ (as a surrogate for 225 Ac). The narrow interlayer spacing of the α-phase (d = 7.6 Å) was found not to be advantageous with respect to rapid ion exchange, while converting the material to the hydrated Na-phase (d = 12.2 Å) allowed ion exchange to proceed more quickly, reaching equilibrium in under an hour. High loading was achieved for all three ions, at 72 ± 4% for Cs + and 95 ± 5% for La 3+ at pH 3 and 108 ± 10% for Bi 3+ at pH 0.5. Leaching profiles of the loaded materials in a simulated human plasma environment of 10 mM carbonate buffer at pH 7.4 and 37°C, showed negligible release, with less than 1% in all cases. Furthermore, the ion-exchange properties and overall benign nature of zirconium phosphate highlight these materials as a promising platform for a targeted alpha therapy drug-delivery vehicle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗