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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

In this paper, we consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Oblique and rippled heliosphere structures from the Interstellar Boundary Explorer

Abstract Past analysis has shown that the heliosphere structure can be deduced from correlations between long-scale solar wind pressure evolution and energetic neutral atom emissions. However, this required spatial and temporal averaging that smoothed out small or dynamic features of the heliosphere. In late 2014, the solar wind dynamic pressure increased by roughly 50% over a period of 6 months, causing a time and directional-dependent rise in around 2–6 keV energetic neutral atom fluxes from the heliosphere observed by the Interstellar Boundary Explorer. Here, we use the 2014 pressure enhancement to provide a simultaneous derivation of the three-dimensional heliospheric termination shock (HTS) and heliopause (HP) distances at high resolution from Interstellar Boundary Explorer measurements. The analysis reveals rippled HTS and HP surfaces that are oblique with respect to the local interstellar medium upwind direction, with significant asymmetries in the heliosphere structure compared to steady-state heliosphere models. We estimate that the heliosphere boundaries contain roughly ten astronomical unit-sized spatial variations, with slightly larger variations on the HTS surface than the HP and a large-scale, southwards-directed obliquity of the surfaces in the meridional plane. Comparisons of the derived HTS and HP distances with Voyager observations indicate substantial differences in the heliosphere boundaries in the northern versus southern hemispheres and their motion over time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Control of Large-Amplitude Combustion Oscillations Using Nanosecond Repetitively Pulsed Plasmas

This paper details the use of nanosecond repetitively pulsed discharges to attenuate combustion instabilities in a 14 kW swirl-stabilized methane/air combustor. The combustor exhibits large-amplitude pressure oscillations ranging from 1 to 4% of the mean pressure during which the flame exhibits bulk motion in each instability cycle, upstream and downstream, as revealed by high-speed chemiluminescence. Control is accomplished with an electrode comprising a pin anode at the centerline of the combustor, allowing a nanosecond spark to be generated in a region spanning close to the flame base, through the shear layers of the swirling flow and ending at the metallic combustor wall. The discharges are generated using 20 kV, 9 kHz pulses; and they correspond to about 120 W of mean power. This results in a suppression of the peak amplitude of the pressure oscillations by a factor of two to four, and 5 dB in the rms value. Using phase-averaged visualizations of the flame with and without plasma, we detail the sequence of flame motion in the course of the instability. With the plasma active, this reveals significant interactions between the flame and the plasma during the suppression. Finally, we present a state-space model of the thermoacoustic system, and we demonstrate open-loop control of the instabilities.

Engineering↗
DOE OSTI2022

System-on-chip integrated circuit technology applications on the DIII-D tokamak for multi-field measurements

Several mm-wave diagnostics on the DIII-D tokamak provide multi-scale and multi-dimensional measurements of plasma profile evolution and turbulence fluctuations. Mm-wave fusion plasma diagnostics that adopt system-on-chip integrated circuit technology can provide better space utilization, flexible installation, and improved sensitivity. In order to further extend this technology for additional fusion facilities with a higher toroidal magnetic field, V-band (55–75 GHz) and F-band (90–140 GHz) chips for Microwave Imaging Reflectometer (MIR) and Electron Cyclotron Emission Imaging (ECEI) instruments are developed and tested in the Davis Millimeter Wave Research Center (DMRC). Current measurement data show that correlation between these SoC-based diagnostic instruments with other state-of-the-art diagnostics enables co-located multi-field turbulence fluctuation measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Laboratory Demonstration of Collisionless Blob Formation via Laser-Produced Plasma Self-Focusing

Strongly localized, propagating plasma density structures that are capable of crossing magnetic field lines are known as “blobs.” Here we demonstrate a novel mechanism for the formation and propagation of an ion gyroradius-scale blob-cavity structure at the interface between a super-Alfvénic laser-produced plasma (LPP) and an ambient magnetized plasma. The LPP self-focuses along the edge of the diamagnetic cavity which results in a dense, jetlike structure as compared to ballistic motion. This collimated flow couples momentum to the ambient plasma through a collisionless process known as Larmor coupling. The Larmor electric fields locally displace the ambient ions forming a blob above the LPP flow. In the region between a gyrating blob and collimated LPP flow, a secondary cavity of expelled magnetic field forms. In conclusion, these findings are supported by particle-in-cell simulations that replicate the blob formation mechanism and provide insight to similar processes in space, astrophysical, and laboratory settings characterized by ion kinetic scales.

Fluorescence spectroscopy↗
DOE OSTI2020

Reversed Field Pinch Research in MST (Final Technical Report)

This is the Final Technical Report for the cooperative agreement DE-FC02-05ER54814 titled “Reversed Field Pinch Research in MST” spanning the funding period 4/1/2005 to 12/14/2019. It summarizes key results from MST research, major publications and presentations, and the mentoring of graduate students and postdocs. The Madison Symmetric Torus (MST) is unique within the U.S. and world fusion plasma program. As one of just five operating reversed field pinch (RFP) experiments, and the only one in the U.S., MST has had an important role in advancing the fundamental understanding of the RFP plasma configuration. The MST research program targets topics that intersect with all three goals. While MST provides a unique opportunity to advance the science and fusion potential of the RFP, as a cousin to the tokamak and stellarator plasma configurations the RFP is a valuable partner in establishing validated predictive capability for fusion science more generally. Specific combinations of the major variables in toroidal confinement, like magnetic field strength, plasma current, and shaping, define the different configurations. By exploring adjacent regions in this major variable parameter space, the RFP exposes dependencies not otherwise accessible in the tokamak and stellarator. This diversity enlarges the arena for scientific discovery, both for fusion and plasma physics. The basic science emphasis for MST research has been the self-organizing behavior of RFP plasmas. This inspired MST’s participation in the NSF Physics Frontier Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas (CMSO), which ran from 2004-2016 and motivated a substantial fraction of MST’s research goals and capabilities. By all measures the MST Program has been remarkably successful in publications, community engagement, and mentoring of young scientists. While the DOE-FES’s proof-of-principle program for the RFP now languishes through the cessation of this cooperative agreement, operation of MST continues as one of two devices included in the new Wisconsin Plasma Physics Laboratory (WiPPL) basic science user facility. The scope of present research for MST in WiPPL is narrowed to the basic science associated with the RFP magnetic configuration and the advancement of specific fusion topics, which now includes studies of disruptions in tokamak plasmas. The world’s effort in RFP fusion research continues, despite a lack of U.S. federal support to investigate the RFP configuration for fusion application, most notably through a major upgrade to the RFX-mod2 facility in Italy [1] and the recent construction of the KTX facility at USTC in China. First plasmas for RFX-mod2 are anticipated in 2021.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗
DOE OSTI2023

Conceptual design of the LHCD system on CFETR

With the capability of saving flux consumption in the current ramp up phase, controlling the safety factor (q) profile, providing the required off-axis current drive (CD), a lower hybrid current drive (LHCD) system was eventually determined to be used on the China Fusion Engineering Test Reactor (CFETR). Here, in this paper, a 20 MW/4.6 GHz system composed of 40 units of 500 kW continuous wave (CW) klystron amplifier is preliminarily designed. The feasibility of the frequency choice of 4.6 GHz in the physical aspect is discussed. In order to minimize the transmission loss, a TE01 over-mode of circular waveguide will be used in the transmission system. Although previous calculations indicates that high field side (HFS) launcher location results in waves damping at inner region than low field side (LFS) launcher, which is favorable for plasma stability, HFS launcher will be extremely difficult in engineering due to the limited space. Calculations show that the tritium breeding ratio (TBR) will be decreased by ~ 1.5% with HFS launcher, while it will be decreased by ~ 0.22% only with LFS launcher. As a result, the LH power will be coupled to plasma from the top port at LFS. The designed passive active multi-junction (PAM) launcher is arranged in an array of 5 rows and 8 columns with the height of 1181 mm and the width of 794 mm. A shielding block with a thickness of 20 cm will be equipped around the feeding waveguides to protect the feeding waveguides and to prevent neutron leakage. A preliminary scheme of the remote maintenance for the LFS launcher is given. The modeling results shows that the density for optimum coupling is ~ 1.8 × 10 17 /m 3 , lower than the 4.6 GHz cut-off density n e_co = 2.6 × 10 17 /m 3 . Around this density, the power directivity (D p ) can be high as 71%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Constructing the space of quasisymmetric stellarators through near-axis expansion

Abstract A simplified view of the space of optimised stellarators has the potential to guide and aid the design efforts of magnetic confinement configurations suitable for future fusion reactors. We present one such view for the class of quasisymmetric stellarators based on their approximate description near their centre (magnetic axis). The result is a space that captures existing designs and presents new ones, providing a common framework to study them. Such a simplified construction offers a basic topological approach, guided by certain theoretical and physical choices, which this paper presents in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics of high energy plasmas. Final report

Various papers and articles are highlighted. The main conclusions of the work in this area, which are of broad interest to fusion research, space physics and astrophysics, are presented in our forthcoming Physics Letters A publication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Visualization of phase-space orbit topological boundary using imaging neutral particle analyzer

A newly-developed imaging neutral particle analyzer (INPA) in the DIII-D tokamak interrogates phase space occupied by fast ions on multiple different orbit topologies, including passing, stagnation, trapped and potato orbits. Depending on plasma parameters and beam injection geometries, this new INPA system is capable of visualizing distributions of fast ions on the selected orbit topology and its associated orbit topology boundaries. More importantly, the system is able to directly visualize the pitch angle scattering process in phase space by measuring fast ions that are scattered across the trapped-passing orbit topology boundaries. It also enables visualization of fast ion confined-loss boundaries, which are determined by particle energy, plasma current and magnetic field.The key goal of this new INPA system is to directly measure pitch angle scattering across phase space induced by drift waves and its interaction with Alfv'en eigenmodes, i.e., a key issue towards \redit{a} future fusion power plant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Synthesis of Hf 6 Ta 2 O 17 superstructure via spark plasma sintering for improved oxidation resistance of multi-component ultra-high temperature ceramics

Ultra-high temperature ceramics (UHTCs) have shown aspiration to overcome challenges in the thermal protection system (TPS) by designing new materials referred to as multi-component UHTCs (MC-UHTCs) in the compositional space. MC-UHTCs have shown remarkable improvement in oxidation resistance due to the formation of the Hf6Ta2O17 superstructure during plasma exposure. Herein, the Hf 6 Ta 2 O 17 superstructure is synthesized via a solid-state reaction between HfO 2 and Ta 2 O 5 powder mixtures during spark plasma sintering (SPS). The compositions chosen are 50 vol% of HfO 2 -50 vol% of Ta 2 O 5 (50HO-50TO) and 70 vol% of HfO 2 -30 vol% of Ta 2 O 5 (70HO-30TO). The phase quantification via Rietveld analysis showed Hf 6 Ta 2 O 17 as a principal phase with some residual Ta 2 O 5 phase in both the samples. The high-temperature thermal stability of the samples was evaluated using high-velocity plasma jet exposure for up to 3 min. 50HO-50TO was able to withstand the intense plasma condition, which is attributed to the higher content of the Hf 6 Ta 2 O 17 phase (~84%) and lower strain in the Ta 2 O 5 phase. The augmentation in the Hf 6 Ta 2 O 17 phase to 94.7% (in 50HO-50TO) post plasma exposure has been attributed to the invariant transformation from a liquid state to Hf 6 Ta 2 O 17 at temperatures >2500 °C during testing. The mechanical integrity is elucidated from the insignificant change in the hardness ~13.3 GPa before and 11.2 GPa after plasma exposure of the 50HO-50TO sample. As a result, the Hf 6 Ta 2 O 17 superstructure's thermo-mechanical stability suggests developing novel oxidation-resistant MC-UHTCs in compositional space for reusable space vehicle applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Model-free stabilization via Extremum Seeking using a cost neural estimator

In this paper, a fully model-free architecture for vertical stabilization of thermonuclear plasmas in tokamak experimental reactors is presented. For the first time, an Extremum Seeking control algorithm is combined with neural networks to estimate the Lyapunov function to be minimized, resulting in a fully data-driven control architecture. The performance of different neural networks are compared. Specifically, Multilayer Perceptrons and Extreme Learning Machines are considered. The proposed architecture is tested in simulation to show that it can counteract relevant plasma disturbances, resulting in a significant improvement in terms of the achievable operative space compared to the Extremum Seeking algorithm, which still relies on model-based cost estimator.

42 ENGINEERING↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plasma-Bubble Imaging and Spectroscopy (Final Report)

The aim of this research was to capture time- and space-resolved behavior of breakdown initiation and development in water with gas bubbles injected between the electrodes when excited with microsecond voltage pulses. The results were directly compared to our previous results from a PCRF collaboration on the development of an electrical discharge in a gas bubble immersed in liquids with different conductivities when excited by nanosecond voltage pulses. For nanosecond regimes, cavitation seems to play a major role in plasma breakdown for low conductivities and discharge initiates at an electrode. On microsecond timescales, where cavitation does not play a role, it has been hypothesized that breakdown in water occurs first in bubbles on or near the electrodes, pre-existing or formed thermally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗