Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Space 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Tomographic imaging of atmospheric pressure plasma on complex surfaces

Many plasma types and behaviors such as streamer, arcs, cathode spots, anode spots, ionization waves, and magnetic field interactions create non-symmetric, fully 3D plasma structures. The plasma distribution in 3D space is heavily influenced by complex surfaces and the coupling interactions between plasma properties and the interfacing material properties. For example, ionization waves propagate in directions where ionization rates are highest, leading to complex configurations that are not fully understood or well characterized. Recent advances in laser diagnostics and models have been able to investigate well-controlled idealized plasmas in 2D fashion, but the complex structure in actual plasmas requires a technique than can provide a more complete 3D picture. However, 3D plasma diagnostics do not currently exist. To address this limitation, this activity will leverage available equipment to build a new tomographic optical imaging capability and advance the state-of-the-art in plasma diagnostics to investigate 3D phenomena on complex surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cost function for low-dimensional manifold topology assessment

Abstract In reduced-order modeling, complex systems that exhibit high state-space dimensionality are described and evolved using a small number of parameters. These parameters can be obtained in a data-driven way, where a high-dimensional dataset is projected onto a lower-dimensional basis. A complex system is then restricted to states on a low-dimensional manifold where it can be efficiently modeled. While this approach brings computational benefits, obtaining a good quality of the manifold topology becomes a crucial aspect when models, such as nonlinear regression, are built on top of the manifold. Here, we present a quantitative metric for characterizing manifold topologies. Our metric pays attention to non-uniqueness and spatial gradients in physical quantities of interest, and can be applied to manifolds of arbitrary dimensionality. Using the metric as a cost function in optimization algorithms, we show that optimized low-dimensional projections can be found. We delineate a few applications of the cost function to datasets representing argon plasma, reacting flows and atmospheric pollutant dispersion. We demonstrate how the cost function can assess various dimensionality reduction and manifold learning techniques as well as data preprocessing strategies in their capacity to yield quality low-dimensional projections. We show that improved manifold topologies can facilitate building nonlinear regression models.

42 ENGINEERING↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optical shock-enhanced self-photon acceleration

Photon accelerators can spectrally broaden laser pulses with high efficiency in moving electron density gradients driven in a rapidly ionizing plasma. When driven by a conventional laser pulse, the group velocity walk-off experienced by the accelerated photons and deterioration of the gradient from diffraction and plasma-refraction limit the extent of spectral broadening. In this work, we show that a laser pulse with a shaped space-time and transverse intensity profile overcomes these limitations by creating a guiding density pro le at a tunable velocity. Self-photon acceleration in this profile leads to dramatic spectral broadening and intensity steepening, forming an optical shock that further enhances the rate of spectral broadening. In this new regime, multi-octave spectra extending from 400nm to 60nm wavelengths, which support near-transform limited < 400 as pulses, are generated over < 100 μm of interaction length.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Counting topological interface modes using simplicial characteristic classes

A computational approach for predicting the number of topological interface modes (TIMs) in hermitian systems using the spectral flow—monopole correspondence is presented. The number of TIMs is determined by calculating the Chern number of a complex line bundle of local polarisation vectors over a phase space sphere surrounding a Weyl point. The Chern number is computed by constructing the simplicial first Chern class of a discrete vector bundle on a simplicial mesh. This approach is gauge invariant, derivative free, structure preserving, and robust to noise. The algorithm is shown to reproduce the expected number of TIMs for the case of equatorial fluid waves and the topological Langmuir cyclotron wave. The possibility of using this algorithm to analyse experimental measurements of bulk wave polarisations and predict the associated number of TIMs is explored in a synthetic example.

discrete vector bundles↗
DOE OSTI2023

Porous Salts as Platforms for Heterogeneous Catalysis

In this study, the preparation of a new class of reactive porous solids, prepared via straightforward salt metathesis reactions, is described here. Reaction of the dimethylammonium salt of a magnesium-based porous coordination cage with the chloride salt of [Cr II Cl(Me 4 cyclam)] + affords a porous solid with concomitant removal of dimethylammonium chloride. The salt consists of the ions combined in the expected ratio based on their charge as confirmed by UV–vis and X-ray photoelectron spectroscopies, ion chromatography (IC), and inductively coupled plasma mass spectrometry (ICP-MS). The porous salt boasts a Brunauer-Emmett-Teller (BET) surface area of 213 m 2 g -1 . Single crystal X-ray diffraction reveals the chromium(II) cations in the structure reside in the interstitial space between porous cages. Importantly, the chromium(II) centers, previously shown to react with O 2 to afford reactive chromium(III)-superoxide adducts, are still accessible in the solid state as confirmed by UV–vis spectroscopy. The site-isolated reactive centers have competence toward hydrogen atom abstraction chemistry and display significantly increased stability and reactivity as compared to dissolved ions.

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

Improved stellarator permanent magnet designs through combined discrete and continuous optimizations

A common optimization problem in the areas of magnetized plasmas and fusion energy is the design of magnets to produce a given three-dimensional magnetic field distribution to high precision. When designing arrays of permanent magnets for stellarator plasma confinement, such problems have tens of thousands of degrees of freedom whose solutions, for practical reasons, should be constrained to discrete spaces. We perform a direct comparison between two algorithms that have been developed previously for this purpose, and demonstrate that composite procedures that apply both algorithms in sequence can produce substantially improved results. One approach uses a continuous, quasi-Newton procedure to optimize the dipole moments of a set of magnets and then projects the solution onto a discrete space. The second uses an inherently discrete greedy optimization procedure that has been enhanced and generalized for this work. Further, the approaches are both applied to design arrays cubic rare-Earth permanent magnets to confine a quasi-axisymmetric plasma with a magnetic field on axis of 0.5 T. The first approach tends to find solutions with higher field accuracy, whereas the second can find solutions with substantially (up to 30%) fewer magnets. When the approaches are combined, they can obtain solutions with magnet quantities comparable to the second approach while matching the field accuracy of the first.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2022

The Motional Stark Effect Diagnostic for NSTX-U

This report describes the development of the motional Stark Effect diagnostic for NSTX. The MSE concept relies upon observation of the Balmer-alpha (n=3 to n=2) emission from a neutral hydrogen beam traversing a plasma. As a beam at high velocity, $\vec{ν}$, passes through the magnetic field in the plasma, $\vec{B}$, it experiences in its reference frame a Lorentz electric field, $\vec{E}$ = $\vec{ν}$ x $\vec{B}$. This electric field causes the spectral emission to be split and polarized as described by the Stark effect. The transition is separated into a set of lines in which the line spacing is linearly proportional to the electric field magnitude, and the line polarization is dependent on the electric field direction. Traditionally, MSE has used polarimetry to determine the polarization angle of the light, and related that to the magnetic field pitch angle in the plasma. The MSE diagnostic was developed under high magnetic field ($>$ 1~T) conditions, and the lower fields of recent high-beta experiments such as NSTX posed challenges to the technique. At low magnetic fields Doppler broadening effects mask the Stark effect. With the use of a narrow bandwidth filter MSE measurements at very low magnetic fields have been obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Orbitron: A crossed-field device for co-confinement of high energy ions and electrons

To explore the confinement of high-energy ions above the space charge limit, we have developed a hybrid magnetic and electrostatic confinement device called an Orbitron. The Orbitron is a crossed-field device combining aspects of magnetic mirrors, magnetrons, and orbital ion traps. Ions are confined in orbits around a high-voltage cathode with co-rotating electrons confined by a relatively weak magnetic field. Experimental and computational investigations focus on reaching ion densities above the space charge limit through the co-confinement of electrons. The experimental apparatus and suite of diagnostics are being developed to measure the critical parameters, such as plasma density, particle energy, and fusion rate for high-energy, non-thermal plasma conditions in the Orbitron. Initial results from experimental and computational efforts have revealed the need for cathode voltages on the order of 100–300 kV, leading to the development of a custom high voltage, ultra-high vacuum bushing rated for 300 kV.

Materials Science↗
DOE OSTI2022

Meso-Scale Electrodynamic Coupling of the Earth Magnetosphere-Ionosphere System

Within the fully integrated magnetosphere-ionosphere system, many electrodynamic processes interact with each other. We review recent advances in understanding three major meso-scale coupling processes within the system: the transient field-aligned currents (FACs), mid-latitude plasma convection, and auroral particle precipitation. (1) Transient FACs arise due to disturbances from either dayside or nightside magnetosphere. As the interplanetary shocks suddenly compress the dayside magnetosphere, short-lived FACs are induced at high latitudes with their polarity successively changing. Magnetotail dynamics, such as substorm injections, can also disturb the current structures, leading to the formation of substorm current wedges and ring current disruption. (2) The mid-latitude plasma convection is closely associated with electric fields in the system. Recent studies have unraveled some important features and mechanisms of subauroral fast flows. (3) Charged particles, while drifting around the Earth, often experience precipitating loss down to the upper atmosphere, enhancing the auroral conductivity. Recent studies have been devoted to developing more self-consistent geospace circulation models by including a better representation of the auroral conductance. It is expected that including these new advances in geospace circulation models could promisingly strengthen their forecasting capability in space weather applications. The remaining challenges especially in the global modeling of the circulation system are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Self-consistent calculations of the electric charge, ion drag force, and the drift velocity of spherical grains using Langevin dynamics and comparisons against canonical experiments

We present trajectory simulation-based modeling to capture the interactions between ions and charged grains in dusty or complex plasmas. Our study is motivated by the need for a self-consistent and experimentally validated approach for accurately calculating the ion drag force and grain charge that determine grain collective behavior in plasmas. We implement Langevin dynamics in a computationally efficient predictor–corrector approach to capture multiscale ion and grain dynamics. Predictions of grain velocity, grain charge, and ion drag force are compared with prior measurements to assess our approach. The comparisons reveal excellent agreement to within ±20% between predicted and measured grain velocities [Yaroshenko et al., Phys. Plasmas 12, 093503 (2005) and Khrapak et al., Europhys. Lett. 97, 35001 (2012)] for 0.64, 1.25 μm grains at ∼20−500 Pa. Comparisons with the measured grain charge [Khrapak et al., Phys. Rev. E 72, 016406 (2005)] under similar conditions reveal agreement to within ∼20% as well. Measurements of the ion drag force [Hirt et al., Phys. Plasmas 11, 5690 (2004); IEEE Trans. Plasma Sci. 32, 582 (2004)] are used to assess the viability of the presented approach to calculate the ion drag force experienced by grains exposed to ion beams of well-defined energy. Excellent agreement between calculations and measurements is obtained for beam energies >10 eV, and the overprediction below 10 eV is attributed to the neglect of charge exchange collisions in our modeling. Along with critical assessments of our approach, suggestions for future experimental design to probe charging of and momentum transfer onto grains that capture the effect of space charge concentration and external fields are outlined.

Physics↗
DOE OSTI2021

A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments

The interaction of radio-frequency (RF) waves with edge turbulence modifies the incident wave spectrum, and can significantly affect RF heating and current drive in tokamaks. Previous lower hybrid (LH) scattering models have either used the weak-turbulence approximation, or treated more realistic, filamentary turbulence in the ray tracing limit. In this study, a new model is introduced which retains full-wave effects of RF scattering in filamentary turbulence. First, a Mie-scattering technique models the interaction of an incident wave with a single Gaussian filament. Next, an effective differential scattering width is derived for a statistical ensemble of filaments. Lastly, a Markov chain solves for the transmitted wave spectrum in slab geometry. This model is applied to LH launching for current drive. The resulting wave spectrum is asymmetrically broadened in angular wavenumber space. This asymmetry is not accounted for in previous LH scattering models. The modified wave spectrum is coupled to a ray tracing/Fokker–Planck solver (GENRAY/CQL3D) to study its impact on current drive. The resulting current profile is greatly altered, and there is significant increase in the on-axis current and decrease in the off-axis peaks. This is attributed to a portion of the modified wave spectrum that is strongly dampened on-axis during the first pass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A dimensionality reduction algorithm for mapping tokamak operational regimes using a variational autoencoder (VAE) neural network

Abstract A variational autoencoder (VAE) is a type of unsupervised neural network which is able to learn meaningful data representations in a reduced dimensional space. We present an application of VAE in identifying the operational stability boundary of tokamak plasma discharges. This model was implemented using a dataset of over 3000 discharges from the high beta tokamak-extended pulse (HBT-EP) device. We found the VAE model to be capable of forming a continuous low-dimensional operational space map and identifying the operational boundaries using a specified warning time window. By projecting the operational parameters onto the same reduced space, this provides an intuitive way for the machine operator or an automated control system to perform disruption avoidance using a relevant control actuator as a discharge approaches a boundary. Pre-programmed GPU control experiments were conducted to demonstrate this control technique using HBT-EP’s saddle control coils as a horizontal position actuator, showing the ability to avoid the oncoming disruptive event and extend the duration of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Self-Generated Plasma Rotation in a Z-Pinch Implosion with Preembedded Axial Magnetic Field

Using detailed spectroscopic measurements, highly resolved in both time and space, a self-generated plasma rotation is demonstrated using a cylindrical implosion with a preembedded axial magnetic field (B z0 ). The rotation direction is found to depend on the direction of B z0 and its velocity is found comparable to the peak implosion velocity, considerably affecting the force and energy balance throughout the implosion. Moreover, the evolution of the rotation is consistent with magnetic flux surface isorotation, a novel observation in a Z pinch, which is a prototypical time dependent system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗