Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma sheaths”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-energy ion beams generated with high efficiency using laser-driven 3D microstructures

Laser-driven ion acceleration in plasma is being proposed as a source of ion beams with a high peak current that can be useful in many fields of science and medicine. Using this method, high proton energies have been achieved by increasing the laser power and by using ultrathin (≤ 200 nm) foils. However, this approach is limited by survivability of the nanotargets to laser prepulses and by difficulty in controlling the plasma acceleration properties. Here, we introduce a new target platform using two-photon polymerization, 3D laser-printed “clone” microstructures with average densities lower than solid that are relatively insensitive to the laser prepulse. Two types of microstructured targets consisting of either a multilayered log-pile or a stochastic arrangement of one micron diameter wires are used. Both demonstrate a higher energy and higher yield proton acceleration compared to thin solid-density foil targets by the robust target normal sheath acceleration (TNSA) mechanism. We find that when such 10–20 μm thick structures are irradiated with a petawatt laser, protons with energies up to 110 MeV and a laser-to-proton conversion efficiency of ~ 10% are obtained. Our work suggests that such microstructures optimized for 60–200 MeV compact proton accelerators are promising for future radiotherapy and other applications.

Physics - Plasma physics↗
DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor

Type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ~65% in the peak and ~35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ~20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ~0.2 km s –1 in attached to ~0.8 km s –1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filaments fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantification of Lanthanides in a Molten Salt Reactor Surrogate Off-Gas Stream Using Laser-Induced Breakdown Spectroscopy

To enable the deployment of molten salt reactor technology, the development of off-gas treatment systems and advanced monitoring tools capable of operating with high temperatures and radiation fields while delivering near real-time information is necessary. Here, we aim to fulfill this requirement and proposes laser-induced breakdown spectroscopy (LIBS) for monitoring molten salt aerosol streams. A sheath gas measuring method was developed to protect optical elements from aerosol particles and to ensure a relatively constant aerosol stream for measurement. An aqueous system was studied to demonstrate the utility of LIBS for monitoring possible fission products in an aerosol stream: Gd, Nd, and Sm up to 2000 parts per million (ppm). A calibration model was built using partial least squares (PLS) regression with five, six, and nine latent variables for Gd, Nd, and Sm, respectively. This calibration model successfully estimated the concentrations of three test samples, which were validated with inductively charged plasma optical emission spectroscopy measurements at a 99.9% confidence interval. To enhance these models, a genetic algorithm was used to filter the spectra before entering the PLS model, thereby limiting the spectral features being regressed to those with greater correlations to concentration. This allowed for the number of latent variables used in the PLS models to be reduced to four, three, and three for Gd, Nd, and Sm, respectively. Finally, the genetic algorithm-filtered PLS models were used to predict the concentrations of the aerosol stream on a real-time dataset and resulted in a 73%, 18%, and 25% improvement in root mean squared error of prediction compared to the original PLS models developed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

ICRF near-field effects sensitivity to magnetic tilt-angle

This paper complements the study reported in G. Urbanczyk et al 2025 Nucl. Fusion 65 046018 by investigating how near-field effects from Ion Cyclotron Range of Frequency (ICRF) antennas depend on the magnetic geometry, particularly focusing on the sensitivity of rectified sheath potentials and associated impurity production to variations in magnetic field tilt-angle. Simulations were conducted using Petra-M, employing a flat model of the ASDEX Upgrade (AUG) three-strap antenna, and COMSOL, utilizing a curved model of the WEST antenna. Results demonstrate indications of sensitivity of the rectified sheath potentials and impurity sputtering to the alignment of Faraday screen bars with the magnetic field lines, emphasizing the necessity for careful consideration of geometric alignment in antenna design for fusion devices. This trend is consistent with experimental data measured with various diagnostics.

Urbanczyk, Guillaume [Universite de Lorraine, Nanc↗
DOE OSTI2024

Comparison of detachment in Ohmic plasmas with positive and negative triangularity

Abstract In recent years, negative triangularity (NT) has emerged as a potential high-confinement L-mode reactor solution. In this work, detachment is investigated using core density ramps in lower single null Ohmic L-mode plasmas across a wide range of upper, lower, and average triangularity (the mean of upper and lower triangularity: δ ) in the TCV tokamak. It is universally found that detachment is more difficult to access for NT shaping. The outer divertor leg of discharges with δ ≈ − 0.3 could not be cooled to below 5 eV through core density ramps alone. The behavior of the upstream plasma and geometrical divertor effects (e.g. a reduced connection length with negative lower triangularity) do not fully explain the challenges in detaching NT plasmas. Langmuir probe measurements of the target heat flux widths ( λ q ) were constant to within 30% across an upper triangularity scan, while the spreading factor S was lower by up to 50% for NT, indicating a generally lower integral scrape-off layer width, λ int . The line-averaged core density was typically higher for NT discharges for a given fuelling rate, possibly linked to higher particle confinement in NT. Conversely, the divertor neutral pressure and integrated particle fluxes to the targets were typically lower for the same line-averaged density, indicating that NT configurations may be closer to the sheath-limited regime than their PT counterparts, which may explain why NT is more challenging to detach.

Physics↗
DOE OSTI2023

Investigation of boosted proton energies through proton radiography of target normal sheath acceleration fields in the multi-ps regime

Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A laser parameter study on enhancing proton generation from microtube foil targets

Abstract The interaction of an intense laser with a solid foil target can drive $$\sim$$ ∼ TV/m electric fields, accelerating ions to MeV energies. In this study, we experimentally observe that structured targets can dramatically enhance proton acceleration in the target normal sheath acceleration regime. At the Texas Petawatt Laser facility, we compared proton acceleration from a $$1\, {\upmu }\hbox {m}$$ 1 μ m flat Ag foil, to a fixed microtube structure 3D printed on the front side of the same foil type. A pulse length (140–450 fs) and intensity ((4–10) $$\times 10^{20}$$ × 10 20 W/cm $$^2$$ 2 ) study found an optimum laser configuration (140 fs, 4 $$\times 10^{20}$$ × 10 20 W/cm $$^2$$ 2 ), in which microtube targets increase the proton cutoff energy by 50% and the yield of highly energetic protons ( $$>10$$ > 10 MeV) by a factor of 8 $$\times$$ × . When the laser intensity reaches $$10^{21}$$ 10 21 W/cm $$^2$$ 2 , the prepulse shutters the microtubes with an overcritical plasma, damping their performance. 2D particle-in-cell simulations are performed, with and without the preplasma profile imported, to better understand the coupling of laser energy to the microtube targets. The simulations are in qualitative agreement with the experimental results, and show that the prepulse is necessary to account for when the laser intensity is sufficiently high.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Kinetic beaming in radiative relativistic magnetic reconnection: a mechanism for rapid gamma-ray flares in jets

ABSTRACT Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other hand, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high energy particles – and their emitted photons – into beams, which may create rapid blips in flux as they cross a telescope’s line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative (specifically inverse Compton) cooling in mediating the observable signatures of this ‘kinetic beaming’ effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) set-ups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

SuperCIC: Enhanced Winding Current Density for Hybrid Windings of Tokamaks

The design is presented for a new approach to cable-in-conduit (SuperCIC) for use in the high-field windings of tokamaks. Two layers of high-field superconductor wires are cabled onto a thin-wall perforated center tube. An overwrap is applied and the cable is inserted as a loose fit into a sheath tube. The sheath tube is drawn down onto the cable to compress the wires onto the center tube and immobilize them. The SuperCIC is then co-wound with a high-strength armor extrusion, which is kerf-cut so that the co-winding onto a coil mandrel can be made without deformation within the armor or the CIC. SuperCIC facilitates hybrid windings, in which subwindings of NbTi are used where B < 7 T, Nb 3 Sn where B < 14 T, and Bi-2212where B > 14 T to minimize superconductor cost. Demount-able splices are used to interconnect layers. A conceptual design has been prepared for a tokamak with aspect ratio A ~ 2.0, magnetic field of 17.4Tat the inner windings and 6.7 T at the plasma,and winding current density ~140 A/mm 2 -sufficient for optimizing the fusion power density in a compact spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulations of divertor heat flux width using transport code with cross-field drifts under the BOUT++ framework

The fluid transport code [trans-electric field (Er) module] under the BOUT++ framework has been used to simulatedivertor heat flux width and boundary Er with all drifts and the sheath potential in the scrape-off layer. The calculated steady state radial Er in the pedestal region has been compared with that of experimental measurements from the Alcator C-Mod tokamak. The magnitude and shape of Er are similar to those of the experimental data. In order to understand the relative role of cross-field drifts vs turbulent transport in setting the heat flux width, four C-Mod enhanced Dα H-mode discharges with a lower single null divertor configuration should be simulated. BOUT++ transport simulations with cross-field drifts included yield similar heat flux width λq to that of experimental measurements (within a factor of 2) from both the probe and the surface thermocouple diagnostics and show a similar trend with plasma current to that of the Eich experimental scaling. The simulations show that both drifts and turbulent transport compete to determine the heat flux width. The magnetic drifts play a dominant role in setting the divertor heat-flux width, while the E × B drift decreases the heat flux width by 10%–25%, leading to improved agreement with the experiment relative to Goldston’s model. A turbulence diffusivity scan (χ) identifies two distinct regimes: a drift dominant regime when χ is small and a turbulence dominant regime when χ is large. The Goldston heuristic drift model yields a lower limit of the width λq.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Single-shot measurements of pulse-front tilt in intense ps laser pulses and its effect on accelerated electron and ion beam characteristics (invited)

Here we report recent single-shot spatiotemporal measurements of laser pulses, including pulse-front tilt (PFT) and spatial chirp, taken at the Compact Multipulse Terawatt laser at the Jupiter Laser Facility in Livermore, CA. STRIPED FISH, a device that measures the complete 3D electric field of fs to ps laser pulses on a single shot, was adapted to near infrared for these measurements. We present the design of the instrument used for these experiments, the on-shot measurements of systematic high-order PFT, and shot-to-shot variations in the measurements of spatiotemporal couplings. Finally, we simulate the effect of PFT in target normal sheath acceleration experiments. These simulations showed that pulse front tilt can steer hot electrons, shape the distribution of the accelerating sheath field, and increase the variability of cutoff energy in the resulting proton spectra. While these effects may be detrimental to experimental accuracy if the pulse front tilt is left unmeasured, hot electron steering shows promise for precision manipulation of the particle source for a range of applications, including irradiation of secondary targets for opacity measurements, radiography, or neutron generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional particle-in-cell simulations of laser-driven multiradiation sources based on double-layer targets

Double-layer targets (DLTs), made of a low-density foam on top of a solid substrate, can efficiently convert the energy of a high-intensity laser to provide sources of photons and protons. We investigate a 30-fs pulse with a peak intensity of I ~ 8.7×10 20 W/cm 2 and a peak power of ~120 TW interacting with a DLT using three-dimensional (3D) particle-in-cell simulations. We focus on providing quantitative results in full 3D geometry on the foam thickness dependence; on the competition between two photon-generating processes in DLTs, i.e., nonlinear inverse Compton scattering (NICS) and bremsstrahlung (BS); and on the acceleration of protons via enhanced target-normal sheath acceleration. Here, we discuss conversion efficiency, average energy, and angular distributions of such multiradiation sources. We find that NICS can prevail over BS if the DLT's substrate is thin enough (~µ⁢m) and that the optimal foam thickness that maximizes the conversion efficiency in NICS and BS photons and the proton cutoff energy, among those considered, is the same (15µ⁢m). These results show that DLTs constitute an excellent tool for developing relatively compact and optimized laser-driven multicomponent radiation sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterization of laser-accelerated proton beams from a 0.5 kJ sub-picosecond laser for radiography applications

Laser-accelerated ion beams show promise for many applications, including high-resolution flash imaging of static or dynamic objects in next-generation radiography to probe materials and plasmas in extreme environments and inertial confinement fusion. To scale up ion beam production for radiography applications, we conducted experiments using sub-picosecond lasers up to 0.5 kJ at the OMEGA-EP facility to characterize proton beams from solid targets, primarily CH/CD sub-micron thin films from which ion beams were also used for static and dynamic radiography for the first time. For standalone sub-micron thin CH films, the highest detected proton energy is in the range of 72–97 MeV. Proton beams with highest energy near or above 60 MeV at full laser energy and similar beam profiles are also measured from low-density CD foams and flat CH foil target of micrometer-scale thickness. The ~ 700–800 nm CH/CD foils achieve the highest ion yield among the targets tested. For sub-micron thin films, the laser prepulse can expand the target and lead to complex interactions, which is simulated using coupled hydrodynamic and two-step kinetic models. Simulations suggest the presence of a micrometer-scale preplasma plateau with near-critical density and further indicate that target normal sheath acceleration, electron heating from Relativistic transparency in the preplasma plateau, and background proton reflection from carbon ion front at the rear side contribute to the resulting proton spectrum from these sub-micron thin targets at various stages. These proton beams show strong potential for radiography and for production of secondary sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of Multi-Filament Textured-Powder Bi-2212/Ag Wire With Enhanced Local Area Ratio

The performance and cost of Bi-2212/Ag wire is limited by the large fraction of Ag matrix (~3:1) that is required in the oxide-powder-in-tube fabrication process. An alternative fabrication process is being developed in which fine-powder Bi-2212 is uniaxially compressed to form bars with a thin Ag foil sheath. The fine powder naturally textures (aligns the a-b planes perpendicular to the direction of compaction) with texture >80% using 200 MPa compression. A billet is formed by stacking trapezoidal- cross-section bars in a symmetric 8-12-16 pattern around an Ag rod and enclosed in an Ag-wall extrusion can. Here, the billet is extruded and drawn to fine wire. Results are presented on present status of the development and testing.

17 WIND ENERGY↗