Engineering Papers⌕ Search

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.

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Direct measurement of the electron turbulence-broadening edge transport barrier to facilitate core–edge integration in tokamak fusion plasmas

Abstract The integration of a high-performance core and a dissipative divertor, or the so-called ‘core–edge integration,’ has been widely identified as a critical gap in the design of future fusion reactors. In this letter, we report, for the first time, direct experimental evidence of electron turbulence at the DIII-D H-mode pedestal that correlates with the broadening of the pedestal and thus facilitates core–edge integration. In agreement with gyrokinetic simulations, this electron turbulence is enhanced by high η e ( η e = L n / L T e , where L n is the density scale length and L T e is the electron temperature scale length), which is due to a strong shift between the density and temperature pedestal profiles associated with a closed divertor. The modeled turbulence drives significant heat transport with a lower pressure gradient that may broaden the pedestal to a greater degree than the empirical and theoretically predicted pedestal width scalings. Such a wide pedestal, coupled with a closed divertor, enables us to achieve a good core–edge scenario that integrates a high-temperature low-collisionality pedestal (pedestal top temperature T e,ped > 0.8 keV and a pedestal top collisionality ν * ped < 1) under detached divertor conditions. This paves a new path toward solving the core–edge integration issue in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermally induced structural competitiveness and metastability of body-centered-cubic iron under nonequilibrium conditions

The structure and stability of iron near melting at multi-megabar pressures are of significant interest in high pressure physics and earth and planetary sciences. While the body-centered cubic (BCC) phase is generally recognized as unstable at lower temperatures, its stability relative to the hexagonal close-packed (HCP) phase at high temperatures (approximately 0.5 eV) in the Earth’s inner core (IC) remains a topic of ongoing theoretical and experimental debate. Motivated by a series of studies attempting to understand the state of the IC and resulting in contradictory findings, we aim to explore the physical factors that may contribute to this confusion in a broader context. Our ab initio calculations show a significant drop in energy, the emergence of a plateau and a local minimum in the potential energy surface, and stabilization of all phonon modes of BCC iron at elevated electron temperatures (>1–1.5 eV). These effects increase the competition among the BCC, HCP, and the face-centered cubic (FCC) phases and lead to the metastability of the BCC structure. Furthermore, at pressure-temperature conditions relevant to the IC, we can expect substantial coexistence of FCC and HCP iron due to their similar free energies and low transition barrier (∼0.2 eV/atom), and the thermodynamic stability of BCC iron is enhanced by its substantial lattice vibration entropy. In conclusion, these findings provide a clear theoretical framework for understanding iron phase relations and solidification processes, both experimentally and in the IC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Shock formation in flowing plasmas by temporally and spatially smoothed laser beams

The cumulative impact of multiple laser speckles on a supersonic plasma flow across optically smoothed laser beams is investigated. The bending of laser beams caused by ponderomotive laser–plasma coupling, together with flow, leads to plasma a momentum-conserving response that results in a deceleration of the flow. Once the flow velocity decreases to a subsonic level, the action of the laser beams can generate a shock within the plasma. This scenario has been predicted theoretically and confirmed by hydrodynamic simulations. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size, and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging only in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser–plasma coupling in the inertial confinement fusion experiments are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Uncertainties in Atomic Data for Modeling Astrophysical Charge Exchange Plasmas

Relevant uncertainties of theoretical atomic data are vital to determining the accuracy of plasma diagnostics in a number of areas, including, in particular, the astrophysical study. We present a new calculation of the uncertainties on the present theoretical ion-impact charge exchange atomic data and X-ray spectra, based on a set of comparisons with the existing laboratory data obtained in historical merged-beam, cold-target recoil-ion momentum spectroscopy, and electron beam ion traps experiments. The average systematic uncertainties are found to be 35–88% on the total cross sections, and 57–75% on the characteristic line ratios. The model deviation increases as the collision energy decreases. The errors on total cross sections further induce a significant uncertainty to the calculation of ionization balance for low-temperature collisional plasmas. Substantial improvements of the atomic database and dedicated laboratory measurements are needed to obtain the current models, ready for the X-ray spectra from the next X-ray spectroscopic mission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Preliminary design overview and performance assessment of the ITER low-field side reflectometer

The design of the ITER low-field side reflectometer (LFSR) has matured to a complete end-to-end preliminary design. LFSR will supply three key plasma measurements: (1) electron density profile, (2) electron density fluctuations, and (3) poloidal rotation. Simultaneous measurements of the three quantities are enabled by an array of six monostatic antennas which inject from an equatorial port on the outboard side of the ITER vessel. Low-loss transmission lines, consisting of corrugated, overmoded waveguide and miter bends, transmit the 30–165 GHz, O- and X-mode signals to and from the ITER plasma. Integrated transmission-line components serve a range of purposes, such as protection from high-power stray radiofrequency radiation, accommodation of transmission-line displacement, and simultaneous measurement of reference and plasma phases during the discharge. Broadband transmission signals are realized by full-band microwave transceivers combined with quasi-optical multiplexing. A field-programmable gate array (FPGA) processor demodulates the profile reflectometer signals, enabling real-time density profile measurements for plasma control system feedback. A full-scale transmission line test facility provides an integrated environment to assess the performance of critical LFSR components. Theoretical modeling together with insertion loss measurements provide the basis for a comprehensive power budget, which accounts for transmitted output power, transmission-line losses, antenna coupling, and plasma effects. Results indicate that high signal-to-noise ratios are achievable with the current design. A synthetic reflectometer model, using real design parameters and baseline ITER profiles, has been developed to estimate the return signal. With evolving microwave and data acquisition technologies, full-band, ultrafast sweeps (<1 μs) will be realizable for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of a gallium alloy liquid metal mirror for application in inertial fusion energy systems

The final optic in future inertial fusion energy plants must maintain high performance and survivability under gigashot laser irradiation while exposed to fusion byproducts. Although no practical solution has been demonstrated, a self-healing final turning mirror based on a grazing-incidence liquid metal mirror has been theoretically explored. This work investigates the suitability of a gallium alloy. The experimental system simulates key conditions, including vacuum exposure and 10-Hz, 355-nm, and 6-ns laser pulses. Experiments address reflectivity, damage, and response as a function of pressure. Results highlight the potential of liquid metal gallium alloys for robust final optics in laser-based fusion power plants and inform non-traditional optics for intense laser systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constraints on ion velocity distributions from fusion product spectroscopy

Abstract Recent inertial confinement fusion experiments have shown primary fusion spectral moments which are incompatible with a Maxwellian velocity distribution description. These results show that an ion kinetic description of the reacting ions is necessary. We develop a theoretical classification of non-Maxwellian ion velocity distributions using the spectral moments. At the mesoscopic level, a monoenergetic decomposition of the velocity distribution reveals there are constraints on the space of spectral moments accessible by isotropic distributions. General expressions for the directionally dependent spectral moments of anisotropic distributions are derived. At the macroscopic level, a distribution of fluid element velocities modifies the spectral moments in a constrained manner. Experimental observations can be compared to these constraints to identify the character and isotropy of the underlying reactant ion velocity distribution and determine if the plasma is hydrodynamic or kinetic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ballooning instability preventing the H-mode access in plasmas with negative triangularity shape on the DIII-D tokamak

Infinite toroidal mode number (n=∞ ) ballooning mode analysis of negative triangularity discharges on DIII-D shows that the access to 2nd stability becomes strongly restricted when the top triangularity decreases even modestly from -0.18 to -0.36. This is observed in experiment to coincide with the suppression of the L-H-transition. Further theoretical analysis with ballooning mode limited pedestals shows that the threshold for opening the 2nd stability access rises from a pedestal temperature of 0.3 keV to 1 keV when the top triangularity is decreased from -0.18 to -0.36 indicating that to access the 2nd stability with the more negative triangular shape would require unrealistically high pedestal temperature. Furthermore, the pedestal predicted by the EPED code agrees with the experimental H-mode profile for the negative triangularity case but in contrast to positive triangular shapes the prediction shows no increase in pedestal height with increasing core pressure when triangularity is negative. This work provides a first model to predict when negative triangularity plasmas can be expected to access the H-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗
DOE OSTI2022

Parametric studies of stream instability-induced higher harmonics in plasma ionization breakdown near an emissive dielectric surface

Here, in this work, we comprehensively investigate the generation of higher harmonic (HH) electric fields normal to the applied rf electric field in multipactor-coexisting plasma breakdown by fully kinetic particle-in-cell (PIC) simulations and a theoretical model. Firstly, a base case at driving frequency 1 GHz, transverse rf electric field amplitude 3 MV m –1 , and background gas pressure 0.2 Torr, is studied in detail. The enhanced harmonic frequency observed is around ten times the fundamental rf frequency, significantly lower than the Langmuir frequency. A theoretical model reveals that the fundamental mechanism of HHs generation is stream–plasma instability, which originates from stream-like secondary electron emission interacting with plasma. The resulting HH frequency and the growth rate of its oscillating amplitude from the theoretical model, agree well with the PIC simulations. With increasing pressure, the HH oscillations are found to be significantly reduced. This is because at higher pressure the gas ionization rate is higher, which causes a more rapidly increasing plasma density, leaving less time for the growth of instability. Furthermore, the parameter space in terms of background gas pressure and rf field amplitude within which the HHs appear is revealed. Finally, the effect of the driving rf frequency on HHs is also investigated, and it is found that the instability-induced oscillating HHs field remains when the driving frequency is increased to 2 GHz, however, it is highly reduced at higher driving frequency of 5 GHz, as oscillations at the fundamental frequency start playing a more important role.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗
DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Progress and challenges in small systems

We present a comprehensive review of the theoretical and experimental progress in the investigation of novel high-temperature quantum chromodynamics phenomena in small systems at both the Relativistic Heavy Ion Collider and the Large Hadron Collider. We highlight the challenges and opportunities associated with studying small systems, by which we generally mean collision systems that involve at least one light ion or even a photon projectile. Here, we discuss perspectives on possible future research directions to better understand the underlying physics at work in the collisions of small systems.

collective behavior↗
DOE OSTI2021

As a Matter of Tension: Kinetic Energy Spectra in MHD Turbulence

While magnetized turbulence is ubiquitous in many astrophysical and terrestrial systems, our understanding of even the simplest physical description of this phenomena, ideal magnetohydrodynamic (MHD) turbulence, remains substantially incomplete. As such, in this work we highlight the shortcomings of existing theoretical and phenomenological descriptions of MHD turbulence that focus on the joint (kinetic and magnetic) energy fluxes and spectra by demonstrating that treating these quantities separately enables fundamental insights into the dynamics of MHD turbulence. This is accomplished through the analysis of the scale-wise energy transfer over time within an implicit large eddy simulation of subsonic, super-Alfvénic MHD turbulence. Our key finding is that the kinetic energy spectrum develops a scaling of approximately k –4/3 in the stationary regime as magnetic tension mediates large-scale kinetic to magnetic energy conversion and significantly suppresses the kinetic energy cascade. This motivates a reevaluation of existing MHD turbulence theories with respect to a more differentiated modeling of the energy fluxes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Response of a charged particle in contact with a chaotic thermostat to an oscillating electric field

A numerical investigation is made of the response of a thermalized charge in contact with a chaotic thermostat to an externally applied, oscillating electric field, with and without a confining magnetic field. This study is based on the recent theoretical development of a chaotic thermostat [G. J. Morales, Phys. Rev. E 99, 062218 (2019)] for which the static properties of a thermalized charge (spatial diffusion and mobilities) were previously established. It is found here that in the unmagnetized case, to extract significant power from the oscillating electric field requires that the oscillatory velocity be larger than the thermal velocity. The spatial diffusion caused by the self-consistent, chaotic fluctuations is shown to be reduced, and even suppressed, by the coherent oscillations. The frequency-dependent mobilities, including cyclotron resonance, are determined, and the nonlinear modifications produced by large electric fields are explored. The value of the mobilities and the modifications to the diffusion coefficient exhibit non-monotonic dependencies on the control parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗