Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space plasma”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Commercialization of microcavity plasma devices and arrays: Systems for VUV photolithography and nanopatterning, disinfection of drinking water and air, and biofilm deactivation for medical therapeutics

Abstract A little more than two decades after the introduction of the first microcavity plasma devices, a growing body of commercial products based on the remarkable properties of these low‐temperature, atmospheric plasmas is now available. Following a brief review of early microplasma lamp development, this article describes microplasma‐based devices and systems currently being manufactured for applications in photolithography, photopatterning, and other nanofabrication processes, such as atomic layer deposition, spectroscopic identification and analysis of materials (including wide bandgap crystals), the disinfection of drinking water in off‐grid communities, deactivating COVID‐19 and other pathogens in air and on surfaces in public spaces, and vacuum ultraviolet lamps for driving atomic clocks. The microplasma‐based treatment of otitis media in the human ear will also be discussed.

42 ENGINEERING↗
DOE OSTI2022

Metal (hydr)oxide surface precipitates and their effects on potassium sorption

Surface precipitation has been shown to occur on rapid time scales in clay and metal oxide mineral systems. The formation of surface precipitates is hypothesized to present new potential sorption sites for potassium (K), where K can become incorporated into newly formed interlayer spaces (e.g., between tetrahedral–octahedral–tetrahedral stacked sheets). The objective of this study is to determine the effects of newly formed mineral surface precipitates on K sorption. Potassium adsorption experiments were conducted by utilizing Al 2 O 3 and SiO 2 sorbents in the presence of various cations (magnesium, zinc, and nickel) that helped to catalyze the formation of surface precipitates. Dissolved concentrations of elements were monitored via inductively coupled plasma optical emission spectrometry (ICP-OES). Solids were characterized via X-ray diffraction (XRD), and K surface complexation was analyzed via X-ray absorption near edge structure (XANES) spectroscopy. X-ray diffraction analysis indicated bayerite, layered double hydroxides (LDH), and silicated LDH were formed as reaction products, thus creating new surface sites for potential K adsorption. The presence of Si increased K adsorption perhaps due to its role in the formation of LDH surface precipitates. When the differences between observed and theoretical surface area normalized K sorption densities were averaged, a 31% increase in K adsorption was observed in the presence of Si. XANES analysis indicated that the binding mechanism of K to Si is different than that of K to Al, perhaps due to the presence of inner-sphere complexation of K to Al-oxide. Samples reacted for one month versus one week yielded more intense XANES post-edge peaks which indicated that the K sorption complex changes over time. Altogether, our findings provide novel insights into the mechanisms of K fixation in soil and has high implication in providing improved K fertilizer recommendation to growers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bifurcation-like transition of divertor conditions induced by X-point radiation in KSTAR L-mode plasmas *

Abstract Density ramps with ion grad B drift directed into lower single null KSTAR L-mode plasmas are associated with a simultaneous and abrupt reduction of the divertor particle flux on both low- and high-field-side targets when the mid-plane line averaged electron density reaches a given level. Target embedded Langmuir probe signals show a clear ‘cliff edge’ behavior similar to that observed in the divertor target electron temperature in DIII-D H-mode plasmas (Eldon et al 2017 Nucl. Fusion 57 066039; McLean et al 2015 J. Nucl. Mater. 463 533–6). The collapse of the particle flux is observed along the whole divertor target area (from private flux region to the far scrape-off layer (SOL)). The critical upstream density of this target flux cliff is invariant under fuel gas throughput modulation. The transition along the cliff occurs in tens of milliseconds. With the cliff, carbon impurities and deuterium neutrals transported through the X-point to the core produce a strong radiation spot near the X-point, seen on bolometric signals, and increase the upstream density. The experimental observations are consistent with time-dependent SOLPS-ITER simulations, which also demonstrate an abrupt transition of the target flux and upstream density with the increase in X-point radiation. The timescale of the cliff predicted by SOLPS-ITER is consistent with the experiment, although, it is influenced by gas throughput or time-dependent numerical methods. In the L-mode phase space of separatrix electron density and temperature, branches are divided based on target temperature, because the latter is strongly coupled to the radiation front and ionization front due to the monotonic characteristic of the parallel electron temperature distribution. Since the H-mode condition operates at a much higher upstream density and electron temperature in phase space, dissipation from sputtered carbon alone leads to the density limit before reaching the X-point radiation condition. This is therefore consistent with the fact that cliffs have never been observed in H-mode KSTAR experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Orbital and short-term variability of iron absorption

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe K band on second timescales. Using data from the Performance Verification phase in 2024 April, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser, and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve’s capability to probe plasma environments in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2020

Plasma-arc lamp high heat flux cycling exposure of neutron irradiated tungsten materials

Thick plate, unalloyed W was neutron irradiated in the High Flux Isotope Reactor (HFIR) at 550 °C to a fast fluence of 1.24 × 10 25 n m -2 E > 0.1 MeV (~0.24 dpa). Unirradiated and irradiated specimens of the material were high heat flux (HHF) tested in the Plasma Arc Lamp (PAL) facility. The PAL uses a high-power photon source to provide a broad and even heat distribution on the sample surface. To simulate on/off cycling of normal operating plasma, the samples were exposed to approximately 800 cycles at 4.73 MW m -2 absorbed heat flux (incident heat fluxes of 10.95 MW m -2 ). Additionally, after PAL exposure, slight changes were observed on the surfaces of the samples with SEM. The samples showed some annealing in the near surface polished region, but they were all below the damage threshold for cracking or other destructive features. The PAL has a large parameter space for future testing. The use of the HFIR and PAL to sequentially expose neutron irradiated samples to HHF will be a powerful tool for understanding materials behavior in a fusion-like environment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.

Zhang, Xin↗
DOE OSTI2022

Design of proton deflectometry with in situ x-ray fiducial for magnetized high-energy-density systems

We report a design and implementation of proton deflectometry with an in situ reference x-ray image of a mesh to precisely measure non-uniform magnetic fields in expanding plasmas at the OMEGA and OMEGA EP laser facilities. The technique has been developed with proton and x-ray sources generated from both directly driven capsule implosions and short pulse laser–solid interactions. The accuracy of the measurement depends on the contrast of both the proton and x-ray images. Here we present numerical and analytic studies to optimize the image contrast using a variety of mesh materials and grid spacings. Our results show clear enhancement of the image contrast by factors of four to six using a high Z mesh with large grid spacing. This leads to further improvement in the accuracy of the magnetic field measurement using this technique in comparison with its first demonstration at the OMEGA laser facility [Rev. Sci. Instrum. 93, 023502 (2022) [CrossRef]].

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Solitary Magnetic Structures Developed From Gyro-Resonance With Solar Wind Ions at Mars and Earth

We investigate solitary magnetic structures in the foreshock of Mars and Earth. The structures exhibit pulse-like magnetic field and density enhancements along with plasma heating and local solar wind (SW) ion reflection. The structures at Mars resemble the foreshock structures developed from ultra-low-frequency electromagnetic waves gyro-resonant with SW ions at Earth, and propagate toward the magnetosphere. We perform fully kinetic simulations to reproduce the solitary structures and the nonlinear evolution of ion distribution functions, illustrating their resonance with SW ions. The structures present self-induced foreshock turbulence that can have space weather effects. Our results advance the fundamental understanding of how SW interacts with planetary magnetospheres, and have potential impact on the current picture of planet-origin ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling Cosmic Reionization

The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the histories of star formation and black hole accretion. A number of massive data sets from new ground- and space-based instruments and facilities over the next decade are poised to revolutionize our understanding of primeval galaxies, the reionization photon budget, the physics of the intergalactic medium (IGM), and the fine-grained properties of hydrogen gas in the "cosmic web". In this review we survey the physics and key aspects of reionization-era modeling and describe the diverse range of computational techniques and tools currently available in this field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
DOE OSTI2022

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗
DOE OSTI2022

A process optimization framework for laser direct energy deposition: Densification, microstructure, and mechanical properties of an Fe-Cr alloy

Laser Direct Energy Deposition (DED) is a metal additive manufacturing technique with the ability to fabricate large and complex parts through deposition of metal powders. However, achieving high-density parts and targeted build heights using DED can be challenging due to the large number of highly sensitive process variables. This work proposes a robust fabrication parameter optimization framework to generate process maps for primary parameters in DED, including laser power, scan speed, mass flow rate, hatch spacing, and layer height. Simple single-track experiments were utilized to map out the parameter space, and a combination of geometric criteria for hatch spacing and layer height were proposed to determine parameter sets that achieve both targeted build heights and mitigate porosity formation. Using this framework, specimens with >99 % density and consistent mechanical properties were successfully fabricated over a wide range of process parameters for an Fe-9wt.%Cr (Fe9Cr) alloy, a surrogate for radiation damage-resistant reduced activation ferritic/martensitic (RAFM) steels. Processing these materials using DED is of particular interest in the development of plasma facing components for nuclear fusion applications. The microstructure and mechanical properties of as-printed Fe9Cr were characterized using optical and electron microscopy, X-ray diffraction, and uniaxial tensile tests. As-printed Fe9Cr displayed ~25 % elongation and ultimate tensile strengths of up to 475 MPa which is comparable to similar wrought alloys. Finally, the proposed framework will allow for accelerated DED parameter optimization for novel alloy systems, as well as open the possibility for local microstructure control while simultaneously mitigating defect formation.

42 ENGINEERING↗
DOE OSTI2020

A WKB based preconditioner for the 1D Helmholtz wave equation

Frequency-domain full-wave solutions to the cold-plasma problem have become ubiquitous in the study of radio frequency power in fusion plasmas. However, recent efforts at extreme levels of geometric fidelity have revealed fundamental limits in the problem size that can be solved by typical sparse direct solver based methods. These limits are of particular importance in the 3D study of RF launchers, where the number of degrees of freedom required can exceed 100 million. In such cases, it would be advantageous to solve the system via iterative means, but due to the large null space of the curl-curl operator, the convergence properties of algorithms like GMRES are poor. Here we present a physics-based preconditioner in the form of a WKB solution and demonstrate the iterative solution to the frequency-domain Helmholtz problem in 1D for several cases ranging from satisfying the WKB approximation to strongly violating it.

Green, David↗
DOE OSTI2023

Synthesis and Thermal Oxidation Resistance of Boron-Rich Boron–Carbide Material

A boron-rich boron–carbide material (B4+δC) was synthesized by spark plasma sintering of a ball-milled mixture of high-purity boron powder and graphitic carbon at a pressure of 7 MPa and a temperature of 1930 °C. This high-pressure, high-temperature synthesized material was recovered and characterized by X-ray diffraction, X-ray photoelectron spectroscopy, Raman spectroscopy, Vickers hardness measurements, and thermal oxidation studies. The X-ray diffraction studies revealed a single-phase rhombohedral structure (space group R-3m) with lattice parameters in hexagonal representation as a = 5.609 ± 0.007 Å and c = 12.082 ± 0.02 Å. The experimental lattice parameters result in a value of δ = 0.55, or the composition of the synthesized compound as B4.55C. The high-resolution scans of boron binding energy reveal the existence of a B-C bond at 188.5 eV. Raman spectroscopy reveals the existence of a 386 cm−1 vibrational mode representative of C-B-B linear chain formation due to excess boron in the lattice. The measured Vickers microhardness at a load of 200 gf shows a high hardness value of 33.8 ± 2.3 GPa. Thermal gravimetric studies on B4.55C were conducted at a temperature of 1300 °C in a compressed dry air environment, and its behavior is compared to other high-temperature ceramic materials such as high-entropy transition metal boride. The high neutron absorption cross section, high melting point, high mechanical strength, and thermal oxidation resistance make this material ideal for applications in extreme environments.

36 MATERIALS SCIENCE↗