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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Revealing the frequency-dependent oscillations in the nonlinear terahertz response induced by the Josephson current

ABSTRACT Nonlinear responses of superconductors to intense terahertz radiation has been an active research frontier. Using terahertz pump-terahertz probe spectroscopy, we investigate the c-axis nonlinear optical response of a high-temperature superconducting cuprate. After excitation by a single-cycle terahertz pump pulse, the reflectivity of the probe pulse oscillates as the pump-probe delay is varied. Interestingly, the oscillatory central frequency scales linearly with the probe frequency, a fact widely overlooked in pump-probe experiments. By theoretically solving the nonlinear optical reflection problem on the interface, we show that our observation is well explained by the Josephson-type third-order nonlinear electrodynamics, together with the emission coefficient from inside the material into free space. The latter results in a strong enhancement of the emitted signal whose physical frequency is around the Josephson plasma edge. Our result offers a benchmark for and new insights into strong-field terahertz spectroscopy of related quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTX-U liquid metal core-edge facility (LMCE)

NSTX-U/LMCE will provide a unique and world-leading research facility to address the primary challenge to delivering economic and timely magnetic fusion energy, namely the need to develop a power and particle exhaust and first-wall system that can withstand very high edge heat fluxes, maximize energy confinement, and avoid the production of large masses of solid eroded first-wall material. The NSTX-U/LMCE facility will assess the ability of liquid metals (LMs) – especially liquid lithium – to provide a new boundary condition for magnetic fusion systems, to extend the lifetime of the plasma facing components (PFCs) and improve core plasma confinement. Such capability is needed to establish the basis for next-step fusion facilities including fusion pilot plants, and to maintain U.S. world leadership in core-edge integration research. NSTX-U/LMCE will leverage the ability to generate very high divertor perpendicular heat flux q⊥ ~ 100MW/m 2 , extensive diagnostics, and liquid-metal-applicable infrastructure of NSTX-U. NSTX-U/LMCE will provide access to a high-confinement plasma core with majority self-driven plasma current, the flexibility to test a range of liquid metal divertor concepts, access to a range of separatrix collisionalities (from high to very low), and the ability to controllably vary the first-wall temperature to vary the plasma- wall interaction physics on liquid lithium components. Further, NSTX-U/LMCE will utilize more reactor-relevant high-Z refractory-metal PFC substrates. With these capabilities the NSTX-U/LMCE facility will explore the full continuum of core-edge solutions ranging from high core radiated power, to conditions with radiative losses concentrated in the scrape-off layer (SOL), and ultimately low recycling conditions. The low collisionality SOL that may be accessible in the low recycling regime is relatively unexplored and will require a kinetic treatment of the edge, which can be addressed theoretically, and with experiments in LTX-β. Additional smaller-scale preparatory R&D facilities will be required to reduce the risk of premature technical/engineering failure of liquid metal systems implemented in NSTX-U. The NSTX-U/LMCE facility aligns very well with recommendations in the FESAC Long-Range Plan and NASEM Pilot Plant reports and the Bold Decadal Vision, will be unique in the world program throughout the next decade, and is garnering private company interest in utilizing NSTX-U/LMCE for development of LM PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Overview of density pedestal structure: Role fueling versus transport

How much of the density pedestal structure profile is determined by fueling versus transport is still a valid question. In most current tokamaks, the Scrape-Off Layer is not opaque enough to screen the neutrals and limit their penetration past the separatrix. As such recycled neutrals as well as additional gas puffing can directly affect the pedestal structure through ionization. Yet, in ITER and other future burning plasma magnetic confinement devices, this ionization source will be pushed further out into the Scrape-Off Layer. This poses the question, how much of the pedestal density structure is governed by this edge ionization source versus plasma transport effects. Theoretically, several turbulent modes have been identified which could provide a 'pinch' like up-gradient transport mechanism. Up-gradient transport is necessary in a source-free region to obtain a peaked density profile, which is often observed in the core of tokamaks. In an opaque Scrape-Off Layer, without the existence of such a pinch, the pedestal density structure would eventually disappear. Not only would this limit our ability to fuel a fusion reactor, it would also affect the pedestal stability. In this paper we will give an overview of the state of the art on what sets the pedestal density structure based on experimental observations as well as theoretical models. We will complement these studies with new results from Alcator C-Mod and DIII-D where the opaqueness is increased to reach values similar to those in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

To investigate the reciprocal interactions between dielectric barrier discharges and different biological targets in cell culture medium

The objective of this project was to investigate the dynamic interactions between non-thermal plasma (NTP) and biological targets suspended in a physiological liquid medium. The experiments under this collaborative study were designed to combine theoretical and experimental approaches to test the hypothesis that when living cells (and the growth medium surrounding them) encounter plasma, they reciprocally and dynamically modify each other.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The theoretical study on intermittency and propagation of geodesic acoustic mode in L-mode discharge near tokamak edge

Abstract Through a systematically developed theory, we demonstrate that the motion of Instanton identified in Zhang et al (2017 Phys. Plasmas 24 122304) is highly correlated to the intermittent excitation and propagation of geodesic acoustic mode (GAM) that is observed in tokamaks. While many numerical simulations have observed the phenomena, it is the first theory that reveals the physical mechanism behind GAM intermittent excitation and propagation. The preceding work is based on the micro-turbulence associated with toroidal ion temperature gradient mode, and slab-based phenomenological model of zonal flow. When full toroidal effect is introduced into the system, two branches of zonal flow emerge: the torus-modified low frequency zonal flow (TLFZF), and GAM, necessitating a unified exploration of GAM and TLFZF. Indeed, we observe that the transition from the Caviton to Instanton is triggered by a rapid zero-crossing of radial group velocity of drift wave and is found to be strongly correlated with the GAM onset. Many features peculiar to intermittent GAMs, observed in real machines, are thus identified in the numerical experiment. The results will be displayed in figures and in a movie; first for single central rational surface, and then with coupled multiple central rational surfaces. The periodic bursting first shown disappears as being replaced by irregular one, more similar to the intermittent characteristics observed in GAM experiments.

Physics↗
DOE OSTI2020

Stabilization of Liner Implosions via a Dynamic Screw Pinch

Magnetically driven implosions are susceptible to magnetohydrodynamic instabilities, including the magneto-Rayleigh-Taylor instability (MRTI). To reduce MRTI growth in solid-metal liner implo-sions, the use of a dynamic screw pinch (DSP) has been proposed [P. F. Schmit et al., Phys. Rev. Lett. 117, 205001 (2016)]. In a DSP con guration, a helical return-current structure surrounds the liner, resulting in a helical magnetic eld that drives the implosion. Here, we present the rst experimental tests of a solid-metal liner implosion driven by a DSP. Using the 1-MA, 100{200-ns COBRA pulsed-power driver, we tested three DSP cases (with peak axial magnetic elds of 2 T, 14 T, and 20 T) and a standard z-pinch (SZP) case (with a straight return-current structure and thus zero axial eld). The liners had an initial radius of 3.2 mm and were made from 650-nm-thick aluminum foil. Images collected during the experiments reveal that helical MRTI modes developed in the DSP cases, while non-helical (azimuthally symmetric) MRTI modes developed in the SZP case. Additionally, the MRTI amplitudes for the 14-T and 20-T DSP cases were smaller than in the SZP case. Speci cally, when the liner had imploded to half of its initial radius, the MRTI amplitudes for the SZP case and for the 14-T and 20-T DSP cases were, respectively, 1.1-0.3 mm, 0.7-0.2 mm, and 0.3-0.1 mm. Relative to the SZP, the stabilization obtained using the DSP agrees reasonably well with theoretical estimates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probing Early-Time Dynamics and Quark-Gluon Plasma Transport Properties with Photons and Hadrons

In this work, we use a hybrid model which relies on QCD effective kinetic theory - KøMPøST - to dynamically bridge the gap between IP-Glasma initial states and viscous hydrodynamics, for the theoretical interpretation of heavy-ion collision results. The hydrodynamic phase is then followed by dynamical freeze-out handled by UrQMD. It is observed that the new pre-equilibrium/pre-hydro phase does influence the extraction of transport coefficients. We also show how the early dynamics is reflected in photon observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Nuclear modification factors for jet fragmentation

Using a recently-developed perturbative-QCD approach for jet evolution in a dense quark-gluon plasma, we study the nuclear modification factor for the jet fragmentation function. The qualitative behaviour that we find is in agreement with the respective experimental observations in Pb+Pb collisions at the LHC: a pronounced nuclear enhancement at both ends of the spectrum. Our Monte Carlo simulations are supplemented with analytic estimates which clarify the physical interpretation of the results. The main source of theoretical uncertainty is the sensitivity of our calculations to a low-momentum cutoff which mimics confinement. To reduce this sensitivity, we propose a new observable, which describes the jet fragmentation into subjets and is infrared-and-collinear safe by construction. We present Monte Carlo predictions for the associated nuclear modification factor together with their physical interpretation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electron heat flux and propagating fronts in plasma thermal quench via ambipolar transport

The thermal collapse of a nearly collisionless plasma interacting with a cooling spot, in which the electron parallel heat flux plays an essential role, is both theoretically and numerically investigated. Here, we show that such thermal collapse, which is known as thermal quench in tokamaks, comes about in the form of propagating fronts, originating from the cooling spot, along magnetic field lines. The slow fronts, propagating with local ion sound speed, limit the aggressive cooling of plasma, which is accompanied by a plasma cooling flow toward the cooling spot. The extraordinary physics underlying such a cooling flow is that the fundamental constraint of ambipolar transport along the field line limits the spatial gradient of electron thermal conduction flux to the much weaker convective scaling, as opposed to the free-streaming scaling, so that a large electron temperature and, hence, pressure gradient can be sustained. The last ion front for a radiative cooling spot is a shock front where cold but flowing ions meet hot ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inflationary stimulated Raman scattering in shock-ignition plasmas

In the shock-ignition inertial confinement fusion scheme, high-intensity lasers propagate through an inhomogeneous coronal plasma, driving a shock designed to cause fuel ignition. During the high-intensity ignitor laser pulse, in the long scale length coronal plasma, back-scattered stimulated Raman scattering (SRS) is likely to be in the kinetic regime. In this work, we use one-dimensional particle-in-cell simulations to show that there is a non-linear frequency shift caused by kinetic effects, resulting in the growth of SRS in an inhomogeneous plasma far exceeding the predictions of fluid theory, the so-called inflationary SRS or iSRS. We find that iSRS occurs over a wide range of density scale lengths relevant to shock-ignition and other directly-driven inertial confinement fusion schemes. The presence of iSRS in shock-ignition plasmas has implications for the theoretical gains from shock-ignition inertial confinement fusion. Here, we quantify the intensity threshold for the onset of iSRS for shock-ignition relevant parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lithium, a path to make fusion energy affordable

In this tutorial article, we review the technological, physics, and economic basis for a magnetic fusion device utilizing a flowing liquid lithium divertor (molten metal velocity in the range of cm/s) and operating in a low-recycling plasma regime. When extrapolated to magnetic fusion reactor scale, the observed effects of a liquid lithium boundary on recycling reduction, confinement increase, and anomalous heat transport mitigation may offer a fundamentally distinct and promising alternative route to fusion energy production. In addition, this lithium-driven low recycling regime could accelerate fusion’s commercial viability since such a device would be smaller, dramatically decreasing plant and electricity costs if all technological complexities are solved. First, the theoretical basis of the energy confinement and fusion performance as well as the related possibilities of low recycling regimes driven by flowing lithium plasma-facing components are reviewed. Then the paper emphasizes the technological obstacles that need to be overcome for developing the necessary systems for such a flowing liquid lithium solution at reactor scale and details how many of these have been overcome at laboratory and/or proof-of-concept scale. Finally, the current and planned scientific and engineering endeavors being performed at the University of Illinois at Urbana-Champaign regarding this alternative reactor option are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Close-Out Project Description for Koepke's Dept of Energy grant DE-SC0021405

Objectives: To establish, for lab & space conditions, EM-IEDDI’s (electromagnetic shear-driven instability's) dispersion relation, unstable range, instability threshold, and mode characteristics, we need LAPD’s Alfven-wave-favorable electromagnetic-style conditions, including higher "beta" (0.001 < beta ≤ 0.3) and low-collisionality. Also, we attempted to intentionally launch or spontaneously destabilize compressional and shear Alfven waves in the strong, localized, perpendicular-velocity-shear region at the interface between coaxial plasmas (one plasma cylinder inside an outer, otherwise hollow, tube, each having a different, controllable, value of plasma electrostatic potential, i.e., “space” potential (not to be confused with the temperature-dependent “floating” potential of an object immersed in the plasma). Nonlinear wave-wave interactions between same-family (EM-IEDD or Alfven) and cross-family (EM-IEDD-with-Alfven) fluctuations were targeted for documentation over a range of spectral overlap. Although laboratory experiments were conducted, the following theoretical work was left unfinished: Analytical non-modal prediction Computational non-modal prediction Check to see if Mikhailenko’s theory formulation leads to his published graphs

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Novel SOLPS-ITER simulations of X-point target and snowflake divertors

Abstract The design and understanding of alternative divertor configurations may be crucial for achieving acceptable steady-state heat and particle material loads for magnetic confinement fusion reactors. Multiple X-point alternative divertor geometries such as snowflakes and X-point targets have great potential in reducing power loads, but have not yet been simulated widely in codes with kinetic neutrals. This paper discusses recent changes made to the SOLPS-ITER code to allow for the simulation of X-point target and low-field side snowflake divertor geometries. Snowflake simulations using this method are presented, in addition to the first SOLPS-ITER simulation of the X-point target. Analysis of these results show reasonable consistency with the simple modelling and theoretical predictions, supporting the validity of the methodology implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗