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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRS2023

Equilibrium Selection via Current Sheet Relaxation and Guide Field Amplification

Although there is a continuous spectrum of current sheet equilibria, how a particular equilibrium is selected by a given system remains a mystery. Yet, only a limited number of equilibrium solutions are used for analyses of magnetized plasma phenomena. Here we present the exact process of equilibrium selection, by analyzing the relaxation process of a disequilibrated current sheet under a finite guide field. It is shown via phase-space analyses and particle-in-cell simulations that the current sheet relaxes in such a way that the guide field is locally amplified, yielding a mixed equilibrium from the spectrum. Comparisons to spacecraft observations and solar wind current sheet statistics demonstrate that such mixed equilibria are ubiquitous and exist as underlying local structures in various physical environments.

Young-Dae Yoon↗
DOE OSTI2023

GPLaSDI: Gaussian Process Latent Space Dynamics Identification

GPLaSDI is a software that introduces GPLaSDI, a novel LaSDI-based framework that relies on Gaussian process (GP) for latent space ODE interpolations. Using GPs offers two significant advantages. First, it enables the quantification of uncertainty over the reduced order model (ROM) predictions. Second, leveraging this prediction uncertainty allows for efficient adaptive training through a greedy selection of additional training data points. This approach does not require prior knowledge of the underlying PDEs. Consequently, GPLaSDI is inherently non-intrusive and can be applied to problems without a known PDE or its residual. We demonstrate the effectiveness of our approach on the Burgers equation, Vlasov equation for plasma physics, and a rising thermal bubble problem. Our proposed method achieves between 200 and 100,000 times speed-up, with up to 7% relative error.

Choi, Youngsoo↗
NASA NTRS1996

Studies for the loss of atomic and molecular species for Io

A summary is presented for research undertaken, physical insight gained, and new directives identified in this project. Significant enhancements for the SO2 neutral cloud model and its dissociative products (SO, O2, O, S) as well as its application to investigating the amount of SO2+ measured by the Voyager PLS instrument in the plasma torus inside of Io's orbit are discussed. New excitation rates for UV and visible emission lines were incorporated in the O and S neutral cloud models, and the very preliminary analysis of UV observations near Io acquired in 1992 by the Hubble Space telescope is discussed. Observations for O('D) 6300 A emission near Io, preliminary reduction of these observations, and an initial assessment of these observations are presented. The analysis of 1985 sodium eclipse and emission data for Io has been completed and is contained in a paper in the Appendix. A larger data set for 1987 sodium emission observations, which will provide a new base for more detailed future studies, is described. A preliminary discussion is also presented for the likely nature of neutral gas clouds for the outer three Galilean satellites.

Smyth, William H.↗
NASA NTRS2001

The Electromagnetic Compatibility (EMC) Design Challenge for Scientific Spacecraft Powered by a Stirling Power Converter

A 55 We free-piston Stirling Technology Demonstration Convertor (TDC) has been tested as part of an evaluation to determine its feasibility as a means for significantly reducing the amount of radioactive material required compared to Radioisotope Thermoelectric Generators (RTGs) to support long-term space science missions. Measurements were made to quantify the low frequency magnetic and electric fields radiated from the Stirling's 80 Hertz (Hz) linear alternator and control electronics in order to determine the magnitude of reduction that will be required to protect sensitive field sensors aboard some science missions. One identified "Solar Probe" mission requires a 100 dB reduction in the low frequency magnetic field over typical military standard design limits, to protect its plasma wave sensor. This paper discusses the electromagnetic interference (EMI) control options relative to the physical design impacts for this power system, composed of 3 basic electrical elements. They are (1) the Stirling Power Convertor with its linear alternator, (2) the power switching and control electronics to convert the 90 V, 80 Hz alternator output to DC for the use of the spacecraft, and (3) the interconnecting wiring including any instrumentation to monitor and control items 1 and 2.

Sargent, Noel B.↗
DOE OSTI2024

Avoiding fusion plasma tearing instability with deep reinforcement learning

For stable and efficient fusion energy production using a tokamak reactor, it is essential to maintain a high-pressure hydrogenic plasma without plasma disruption. Therefore, it is necessary to actively control the tokamak based on the observed plasma state, to manoeuvre high-pressure plasma while avoiding tearing instability, the leading cause of disruptions. This presents an obstacle-avoidance problem for which artificial intelligence based on reinforcement learning has recently shown remarkable performance. However, the obstacle here, the tearing instability, is difficult to forecast and is highly prone to terminating plasma operations, especially in the ITER baseline scenario. Previously, we developed a multimodal dynamic model that estimates the likelihood of future tearing instability based on signals from multiple diagnostics and actuators. Here we harness this dynamic model as a training environment for reinforcement-learning artificial intelligence, facilitating automated instability prevention. We demonstrate artificial intelligence control to lower the possibility of disruptive tearing instabilities in DIII-D, the largest magnetic fusion facility in the United States. The controller maintained the tearing likelihood under a given threshold, even under relatively unfavourable conditions of low safety factor and low torque. In particular, it allowed the plasma to actively track the stable path within the time-varying operational space while maintaining H-mode performance, which was challenging with traditional preprogrammed control. This controller paves the path to developing stable high-performance operational scenarios for future use in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Data-driven profile prediction for DIII-D

A new, fully data-driven algorithm has been developed that uses neural networks to predict plasma pro les on a scale of τ E into the future given actuators and the present plasma state. The model was trained and tested on DIII-D data from the 2013-2018 experimental campaigns. The model is accurate on average, with q predictions the worst and pressure predictions the best. The model can run in milliseconds and is very simple to use. This makes it a potentially useful tool for operators and physicists when planning plasma scenarios. It also is a candidate for doing phase-space exploration without going through the DIIID database or complicated and computationally expensive simulation codes. Here, a reduced model using only realtime diagnostics has also been developed and formed the basis for a model-predictive control algorithm implemented and successfully tested on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Extended particle absorber for efficient modeling of intense laser-solid interactions

Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A massively parallel time-domain coupled electrodynamics–micromagnetics solver

We present a high-performance coupled electrodynamics–micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain approach for Maxwell’s equations to a magnetization model described by the Landau–Lifshitz–Gilbert equation. The algorithm is implemented in the Exascale Computing Project software framework, AMReX, which provides effective scalability on manycore and GPU-based supercomputing architectures. Furthermore, the code leverages ongoing developments of the Exascale Application Code, WarpX, which is primarily being developed for plasma wakefield accelerator modeling. Our temporal coupling scheme provides second-order accuracy in space and time by combining the integration steps for the magnetic field and magnetization into an iterative sub-step that includes a trapezoidal temporal discretization for the magnetization. The performance of the algorithm is demonstrated by the excellent scaling results on NERSC multicore and GPU systems, with a significant (59×) speedup on the GPU using a node-by-node comparison. We demonstrate the utility of our code by performing simulations of an electromagnetic waveguide and a magnetically tunable filter.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Plasma motions in planetary magnetospheres

Interplanetary space is pervaded by a supersonic 'solar wind' plasma; five planets, in addition to the earth, have magnetic fields of sufficient strength to form the cometlike cavities called 'magnetospheres'. Comparative studies of these structures have indicated the specific environmental factor that can result in dramatic differences in the behavior of any pair of magnetospheres. Although planetary magnetospheres are large enough to serve as laboratories for in situ study of cosmic plasma and magnetic field behavior effects on particle acceleration and EM emission, much work remains to be done toward relating magnetospheric physics results to the study of remote astrophysical plasmas.

Hill, T. W.↗
NASA NTRS1992

Solar cell arcing: The role of outgassing and contamination

The effect of outgassing, venting, and thruster firing events on spacecraft system performance has been a long standing issue. Recent laboratory measurements on negatively biased high voltage solar cells at Physical Sciences Inc. (PSI) suggest that some currently designed and certainly future space power systems must address/re-evaluate this issue. Our observations show that exposing these cells to moderate levels (10 exp -3 torr-min) of H2O vapor enhances the arcing frequency, while heating to 85 C to remove water vapor significantly reduces the arc frequency. The interaction of the adhesive used to attach the cover glass to the solar cell with ambient water vapor is the key factor in determining arcing rates. Elimination of adhesive exposed to the environmental plasma reduces the arc frequency more than two orders of magnitude, and eliminates any sensitivity to H2O exposure. The adhesive may also become a source of spacecraft contamination. Macroscopic amounts were observed to blow off some arc events, and (we assume) electrostatically precipitate at other negatively biased locations. Data, analysis, and potential impact for future space platforms are discussed.

Marinelli, W. J.↗
NASA NTRSDate not supplied

Wavelength Calibration of the Full-sun Ultraviolet Rocket SpecTrograph (FURST)

The Sun has a well-known periodicity in sunspot number and magnetic field variation. The underlying cause of this 11-year cycle is not fully understood and has yet to be connected with those processes in other stellar objects. The Full-sun Ultraviolet Rocket SpecTrograph (FURST) is a sounding rocket payload being developed by Montana State University (MSU) alongside the Marshall Space Flight Center(MSFC) solar physics group. Scheduled to launch from White Sands Missile Range (WSMR) in 2022, this instrument is unique in that it will provide the connection between stellar observatories with measurements of our Sun. It will achieve this through measuring high-resolution full-disk spectral irradiance. We aim to obtain a wavelength resolution R>10,000 in the 120 - 181 nm UltraViolet (UV) range, on par with that of the Hubble (HST) Space Telescope Imaging Spectrograph (STIS). This resolution goal will allow us to study the relatively low-temperature plasma in the chromosphere and lower corona with spectral accuracy down to 0.1 Å (a Doppler-shift of about±30 km/s). In addition, the Lyman Alpha(121 nm) line is known to saturate most CCD electronics. These factors illustrate the particular challenge of precise wavelength calibration for this spectral range. We are building a collimator in order to calibrate the FURST instrument under these strict spectral requirements. This paper will present the results of our simulation of the diagnostic lamp signal to be used for wavelength calibration. The simulation allows us to begin to account for photon noise, electronic readout noise, and statistical error. These in turn lead to the development of our pre- and post-launch calibration plans. Future work includes absolute radiometric and wavelength calibration with this new collimator. In addition, the ability of FURST to measure small Doppler-shifts will provide capabilities for planetary atmospheric scientists. This impact is coupled with the diverse international partnership created by the closely-knit Sounding Rocket teams around the globe. Sounding Rockets like FURST have an even broader impact, as they encourage future satellite missions under the prospect of long-term observations

FURST↗
DOE OSTI2021

Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Physics↗
DOE OSTI2020

Self-acceleration and energy channeling in the saturation of the ion-sound instability in a bounded plasma

A novel regime of the saturation of the Pierce-type ion-sound instability in bounded ion-beam-plasma system is revealed in 1D PIC simulations. It is found that the saturation of the instability is mediated by the oscillating virtual anode potential structure. The periodically oscillating potential barrier separates the incoming beam ions into two groups. One component forms a supersonic beam which is accelerated to an energy exceeding the energy of the initial cold ion beam. The other component is organized as a self-consistent phase space structure of trapped ions with a wide energy spread – the ion hole. The effective temperature (energy spread) of the ions trapped in the hole is lower than the initial beam energy. In the final stage the ion hole expands over the whole system length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Counter-propagating toroidal Alfvén eigenmodes in tokamaks

Mechanisms of destabilization of toroidal Alfvén eigenmodes (TAEs) in tokamaks are analyzed with the aim to reveal those leading to modes propagating in the direction opposite to plasma current, i.e. counter-propagating modes (ctr-TAE). Plasmas with fast-ions sources [such as neutral beam injection (NBI), ion cyclotron resonance heating, fusion reactions] and without them (Ohmic discharges) are considered. A particular NSTX-U experiment with NBI, where co- and counter-propagating TAEs were observed simultaneously (Podestà et al 2018 Nucl. Fusion 58 082023), is considered. It is concluded that both types of TAEs occurred because their destabilization was caused by the velocity anisotropy of beam ions, which overrode effects of spatial inhomogeneity of these ions.

Alfvénic instabilities↗