SEARCH · Search NASA
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.
A divertor for system supporting plasma for fusion employs permanent magnets providing far weaker magnetic field strength than conventionally used in cryogenic magnetic systems through an ability to place the permanent magnets in close proximity to the plasma containment volume.
We have developed an experimental platform at the National Ignition Facility that employs colliding planar shocks to produce warm dense matter with uniform conditions and enable high-precision equation of state measurements. The platform uses simultaneous x-ray Thomson scattering and x-ray radiography to measure the density, electron temperature, and ionization state in warm dense matter. The experimental platform is designed to create a large volume of uniform plasma (approximately 700×700×150μm 3 ) at pressures approaching 100 Mbar and minimize the distribution of plasma conditions in the x-ray scattering volume, significantly improving the precision of the measurements. Here, in this study, we present the experimental design of the platform and compare hydrodynamic simulations to x-ray radiography data from initial experiments studying hydrocarbons, producing uniform densities within ±25% of the average probed condition. We show that the platform creates a homogeneous plasma that can be characterized using x-ray Thomson scattering. Thus, the new platform enables accurate measurements of plasma conditions necessary to test models for the equation of state and ionization potential depression in the warm dense matter regime.
This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.
Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.
Optical collective Thomson scattering is used to diagnose magnetized high energy density physics experiments at the Magpie pulsed-power generator at Imperial College London. The system uses an amplified pulse from the 2nd harmonic of a Nd:YAG laser (3 J, 8 ns, 532 nm) to probe a wide diversity of hightemperature plasma objects; with densities in the range of 10 17 -10 19 cm -3 and temperatures between 10 eV and a few keV. The scattered light is collected from 100 μm-scale volumes within the plasmas, which are imaged onto optical fiber arrays. Multiple collection systems observe these volumes from different directions, providing simultaneous probing with different scattering K-vectors (and different associated α-parameters, typically in the range 0.5 – 3) allowing independent measurements of separate velocity components of the bulk plasma flow. The fiber arrays are coupled to an imaging spectrometer with a gated ICCD. The spectrometer is configured to view the ion-acoustic waves (IAWs) of the collective Thomson scattered spectrum. Fits to the spectra with the theoretical spectral density function S(K,ω) yield measurements of the local plasma temperatures and velocities. Fitting is constrained by independent measurements of the electron density from laser interferometry, and the corresponding spectra for different scattering vectors. Furthermore, this TS diagnostic has been successfully implemented on a wide range of experiments, revealing temperature and flow velocity transitions across magnetized shocks, inside rotating plasma jets and imploding wire arrays, as well as providing direct measurements of drift velocities inside a magnetic reconnection current sheet.
We present two-dimensional (2D) optical Thomson scattering measurements of electron density and temperature in laser-produced plasmas. The novel instrument directly measures ne(x,y) and Te(x,y) in two dimensions over large spatial regions (cm2) with sub-mm spatial resolution, by automatically translating the scattering volume while the plasma is produced repeatedly by irradiating a solid target with a high-repetition-rate laser beam (10 J, ∼1012 W/cm2, 1 Hz). In this paper, we describe the design and motorized auto-alignment of the instrument and the computerized algorithm that autonomously fits the spectral distribution function to the tens-of-thousands of measured scattering spectra, and captures the transition from the collective to the non-collective regime with distance from the target. As an example, we present the first 2D scattering measurements in laser-driven shock waves in ambient nitrogen gas at a pressure of 0.13 mbar.
Activation of liquids with atmospheric pressure plasmas is being investigated for environmental and biomedical applications. When activating the liquid using gas plasma produced species (as opposed to plasmas sustained in the liquid), a rate limiting step is transport of these species into the liquid. To first order, the efficiency of activating the liquid is improved by increasing the ratio of the surface area of the water in contact with the plasma compared to its volume—often called the surface-to-volume ratio (SVR). Maximizing the SVR then motivates the plasma treatment of thin films of liquids. In this paper, results are discussed from a computational investigation using a global model of atmospheric pressure plasma treatment of thin water films by a dielectric barrier discharge (DBD) sustained in different gases (Ar, He, air, N 2 , O 2 ). The densities of reactive species in the plasma activated water (PAW) are evaluated. The residence time of the water in contact with the plasma is increased by recirculating the PAW in plasma reactor. Longer lived species such as H 2 O 2aq and NO 3 - aq accumulate over time (aq denotes an aqueous species). DBDs sustained in Ar and He are the most efficient at producing H 2 O 2aq , DBDs sustained in argon produces the largest density of NO 3 - aq with the lowest pH, and discharges sustained in O 2 and air produce the highest densities of O 3aq . Finally, comparisons to experiments by others show agreement in the trends in densities in PAW including O 3aq , OH aq , H 2 O 2aq and NO 3 - aq , and highlight the importance of controlling desolvation of species from the activated water.
The Gubser solution to inviscid relativistic fluid dynamics is used to examine the role of transverse expansion on the energy spectrum of photons radiated by quark-gluon plasma. Transverse flow is shown to be a modest effect on the energy spectrum of photons as a whole, despite its large effect on rare high-energy photons produced at low temperatures. An exact expression is derived for the volume of the plasma as a function of its temperature. Furthermore, a simple formula is obtained for the energy spectrum of high-energy thermal photons, which is used to relate the inverse slope T eff of the photon spectrum at energy E to the maximum temperature of the plasma T 0 , finding T eff ≈ T 0 / (1 + $\frac{5}{2}$ T 0 /E).
Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.
Exhausting power from the hot fusion core to the plasma-facing components is one fusion energy’s biggest challenges. The MAST Upgrade tokamak uniquely integrates strong containment of neutrals within the exhaust area (divertor) with extreme divertor shaping capability. By systematically altering the divertor shape, this study shows the strongest evidence to date to our knowledge that long-legged divertors with a high magnetic field gradient (total flux expansion) deliver key power exhaust benefits without adversely impacting the hot fusion core. These benefits are already achieved with relatively modest geometry adjustments that are more feasible to integrate in reactor designs. Benefits include reduced target heat loads and improved access to, and stability of, a neutral gas buffer that ‘shields’ the target and enhances power exhaust (detachment). Analysis and model comparisons shows these benefits are obtained by combining multiple shaping aspects: long-legged divertors have expanded plasma-neutral interaction volume that drive reductions in particle and power loads, while total flux expansion enhances detachment access and stability. Containing the neutrals in the exhaust area with physical structures further augments these shaping benefits. These results demonstrate strategic variation in the divertor geometry and magnetic topology is a potential solution to one of fusion’s power exhaust challenge.
Abstract The edge island domain in Wendelstein 7-X consists of divertor islands whose position coincide with the location of rational values of the rotational transform - ι ≈ ( 5 / 6 , 5 / 5 , 5 / 4 ) and surround the main confinement volume (the ‘main plasma’), providing the magnetic topology for an island divertor. The 5/5 edge consists of five individual islands that are unconnected. In contrast, a single island connects onto itself after 6 or 4 toroidal transits in the 5/6 or 5/4 edge, respectively. Many interesting phenomena are related to these islands and diagnostic analyses require a mapping from ‘laboratory’ or real space coordinates to the island coordinate system (ICS). A procedure is described here to calculate several scalar and vector quantities for closed island structures which can be utilized in fast interpolation schemes for inverse maps. For the 5/5 edge, a fixed-boundary vacuum (zero beta) magneto-hydrodynamic solution of the 5/5 island is found with the Variational Moments Equilibrium Code, VMEC. The solution is compatible with already existing routines which determine the radial and binormal vectors of VMEC solutions at arbitrary laboratory coordinates. VMEC does not support solutions for the 5/4 and 5/6 islands, but the radial and binormal vectors are available in a local 2D Fourier ICS.
In-plasma x-ray photoelectron spectroscopy (plasma-XPS) emerges as a powerful platform for real-time, in situ chemical analysis under conditions relevant to semiconductor processing and other plasma-enabled technologies. This study investigates the origins of binding energy (BE) shifts and the formation of “satellite” peaks observed during plasma-XPS measurements across conductive, dielectric, and gas-phase systems. Using a standard laboratory-based ambient pressure XPS apparatus coupled with an alternating current (AC)-driven capacitively coupled plasma source, we demonstrate that metastable surface species, such as transient Au oxides, can be detected during plasma exposure, revealing chemical states that are hardly accessible using conventional ultrahigh vacuum (UHV) XPS. In dielectric samples (e.g., undoped diamond, sapphire), we observe pressure- and plasma-type-dependent BE shifts of more than 50 eV, attributed to x-ray-induced and plasma-mediated surface charging. These shifts are mitigated at higher pressures/plasmas or in electronegative plasmas (e.g., O 2 ), the latter due to enhanced charge compensation mechanisms involving slow negative ions. For gas-phase species, AC-plasma excitation leads to spectral broadening and the emergence of “satellite” peaks with energy separations of a few electron volts, linked to oscillating local plasma potentials in the probing volume. Furthermore, these findings highlight the complex and important interplay between plasma parameters, surface charging, and local electric fields in shaping XPS spectra. Overall, plasma-XPS emerges as a critical metrological tool for probing transient surface chemistry, with implications for semiconductor processing, material synthesis, and plasma diagnostics.
A product having at least one plasma lamp that includes plates that are approximately parallel, with at least one array of microcavities formed in a surface of at least one plate. When desirable, the plates are separated a fixed distance by spacers with at least one spacer being placed near the plate's edge to form a hermetic seal therewith. A gas makes contact with the microcavity array. Electrodes capable of delivering a time-varying voltage are located such that the application of the time-varying voltage interacts with the gas to form a glow discharge plasma in the microcavities and the fixed volume between the plates. The glow discharge plasma efficiently and uniformly emits radiation that is predominantly in the UV/VUV spectral range with at least a portion of the radiation being emitted from the plasma lamp.
A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.
Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.
Abstract Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark–gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we determine the speed of sound in an extended volume of quark–gluon plasma using lead–lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb −1 . The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 ± 0.002 ( stat ) ± 0.016 ( syst ) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 ± 8 ( syst ) MeV . The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.
Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.
To date, 14 C tracer studies using accelerator mass spectrometry (AMS) have not yet resolved lipid-soluble analytes into individual lipoprotein density subclasses. The objective of this work was to develop a reliable method for lipoprotein separation and quantitative recovery for biokinetic modeling purposes. The novel method developed provides the means for use of small volumes (10–200 µL) of frozen plasma as a starting material for continuous isopycnic lipoprotein separation within a carbon- and pH-stable analyte matrix, which, following post-separation fraction clean up, created samples suitable for highly accurate 14 C/ 12 C isotope ratio determinations by AMS. Manual aspiration achieved 99.2 ± 0.41% recovery of [5- 14 CH 3 ]-(2R, 4'R, 8'R)-α-tocopherol contained within 25 µL plasma recovered in triacylglycerol rich lipoproteins (TRL = Chylomicrons + VLDL), LDL, HDL, and infranatant (INF) from each of 10 different sampling times for one male and one female subject, n = 20 total samples. Small sample volumes of previously frozen plasma and high analyte recoveries make this an attractive method for AMS studies using newer, smaller footprint AMS equipment to develop genuine tracer analyses of lipophilic nutrients or compounds in all human age ranges.