SEARCH · Engineering Papers

Results for “plasma boundary 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science↗

NASA NTRS1997

Computational Study of Electron-Molecule Collisions Related to Low-Temperature Plasmas

Computational study of electron-molecule collisions not only complements experimental measurements, but can also be used to investigate processes not readily accessible experimentally. A number of ab initio computational methods are available for these types of calculations. Here we describe a recently developed technique, the finite element Z-matrix method, Analogous to the R-matrix, method, it partitions the space into regions and employs real matrix elements. However, unlike the implementation of the R-matrix method commonly used in atomic and molecular physics, the Z-matrix method is fully variational. In the present implementation, a mixed basis of finite elements and Gaussians is used to represent the continuum electron, thus offering full flexibility without imposing fixed boundary conditions. Numerical examples include the electron-impact dissociation of N2 via the metastable A3Su+ state, a process which may be important in the lower thermosphere, and the dissociation of the CF radical, a process of interest to plasma etching. To understand the dissociation pathways, large scale quantum chemical calculations have been carried out for all target states which dissociate to the lowest five limits in the case of N2, and to the lowest two limits in the case of CF. For N2, the structural calculations clearly show the preference for predissociation if the initial state is the ground X1 Sg+ state, but direct dissociation appears to be preferable if the initial state is the A3Su+ state. Multi-configuration SCF target functions are used in the collisional calculation.

Huo, Winifred M.↗

NASA NTRS2017

The Mars Crustal Magnetic Field Control of Plasma Boundary Locations and Atmospheric Loss: MHD Prediction and Comparison with MAVEN

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

Global↗

DOE OSTI2026

Experimental validation of predictive plasma boundary simulations using BOUT++/Hermes-3 [Slides]

Predictive simulations of heat and particle fluxes in the plasma boundary are a critical missing capability for magnetic fusion energy power plant design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A magnetic diagnostic suite for the Pegasus-III experiment

Pegasus-III is an ultralow aspect ratio spherical tokamak providing a dedicated US experiment for comparative solenoid-free startup studies. A new magnetic diagnostic suite for equilibrium and low frequency (<200 kHz) magnetohydrodynamic mode analysis has been installed. These new diagnostics address the significant challenges of measuring magnetic field in a high noise environment with the majority constrained to fit in an 8 mm diagnostic gap on the high field side. Electrostatic switching noise generated by the 16 independent current feedback-controlled power supplies produces dV cm /dt ~ 1 kV/μs and volt level common mode noise on the magnetics. Immunity to this switching noise is accomplished through differential signal runs and signal processing, along with end-to-end electromagnetic interference shielding. The magnetic measurements are simultaneously digitized at 1 MHz and conditioned by precision 8 pole Butterworth filters with a corner frequency of 200 kHz to prevent aliasing down to the 16-bit level over the full passband. Ex-vessel calibrations of the B p coils were completed with a typical uncertainty of <0.5%. Stray toroidal field pickup from coil misalignment or positioning errors is corrected using a physics-based model. Comparisons of the corrected measurements to modeling agree to within 1.3% on average. This is within the 1.5% measurement uncertainty that a sensitivity analysis determined is needed for accurate fast boundary and equilibrium reconstruction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

Physics↗

DOE OSTI2022

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2000

Accuracy Improvement in Magnetic Field Modeling for an Axisymmetric Electromagnet

This paper examines the accuracy and calculation speed for the magnetic field computation in an axisymmetric electromagnet. Different numerical techniques, based on an adaptive nonuniform grid, high order finite difference approximations, and semi-analitical calculation of boundary conditions are considered. These techniques are being applied to the modeling of the Variable Specific Impulse Magnetoplasma Rocket. For high-accuracy calculations, a fourth-order scheme offers dramatic advantages over a second order scheme. For complex physical configurations of interest in plasma propulsion, a second-order scheme with nonuniform mesh gives the best results. Also, the relative advantages of various methods are described when the speed of computation is an important consideration.

Ilin, Andrew V.↗

NASA NTRSDate not supplied

What do we know about the reconnection electric field?

The reconnection electric field is at the core of the reconnection process. Its magnitude is directly related to the effectiveness of the magnetic flux, energy, and mass transport in general and across magnetic boundaries. The reconnection electric field has hence been the focus of intense research, fueled by the Magnetospheric Multiscale mission (MMS) and concurrent theory and modeling. In this presentation, we review the present state of knowledge pertaining to the reconnection electric fields and its plasma physical underpinnings. We will take a close look at reconnection in symmetric and asymmetric systems, and how they are related. The presentation will further relate relevant theory and modeling results to the ground truth provided by MMS observations, and how observations have driven theory and vice-versa. After summing up the current state of knowledge, we will identify a set of open questions, which call for future scientific investigations.

Michael Hesse↗

NASA NTRSDate not supplied

What Do We Know About the Reconnection Electric Field?

Michael Hesse↗

DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Free-boundary MRxMHD equilibrium calculations using the stepped-pressure equilibrium code

Here the stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 , 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', , is computed, and the virtual-casing principle is used iteratively to compute the normal field on so that the equilibrium is consistent with an externally produced magnetic field. Recent modifications to SPEC are described, such as the use of Chebyshev polynomials to describe the radial dependence of the magnetic vector potential, and a variety of free-boundary verification calculations are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Noninvasive THz-TDS measurements of plasma bounded and optically shielded by Hall thruster wall material

We experimentally demonstrate the capability of terahertz time-domain spectroscopy (THz-TDS) to noninvasively measure the electron density and collision frequency of plasma bounded and optically shielded by Hall thruster wall material. This paper augments the standard THz-TDS plasma diagnostic theory to account for plasma boundaries, presents THz optical property measurements of three different wall materials (grades M, M26, and HP boron nitride composite), and provides electron density and collision frequency measurements of an inductively coupled plasma bounded and optically shielded by each wall material. We find that the electron density measurement capability is weakly impacted by the boundaries, whereas the electron collision frequency measurement capability is strongly reduced by the boundaries. Here, the bounded plasma electron density trends deviate substantially from those of the unbounded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Combined plasma–coil optimization algorithms

Combined plasma–coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria for multiregion relaxed magnetohydrodynmics (MRxMHD) is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi-free-boundary optimization, and the free-boundary (coil) optimization. These are described using the MRxMHD energy functional, the Biot–Savart integral, the coil-penalty functional and the virtual casing integral and their derivatives. The proposed free-boundary equilibrium calculation differs from existing methods in how the boundary-value problem is posed, and for the new approach it seems that there is not an associated energy minimization principle because a non-symmetric functional arises. We propose to solve the weak formulation of this problem using a spectral-Galerkin method, and this will reduce the free-boundary equilibrium calculation to something comparable to a fixed-boundary calculation. In our discussion of combined plasma–coil optimization algorithms, we emphasize the importance of the stability matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1978

A theory of local and global processes which affect solar wind electrons. 1: The origin of typical 1 AU velocity distribution functions: Steady state theory

A detailed first principle kinetic theory for electrons which is neither a classical fluid treatment nor an exospheric calculation is presented. This theory illustrates the global and local properties of the solar wind expansion that shape the observed features of the electron distribution function, such as its bifurcation, its skewness and the differential temperatures of the thermal and suprathermal subpopulations. Coulomb collisions are substantial mediators of the interplanetary electron velocity distribution function and they place a zone for a bifurcation of the electron distribution function deep in the corona. The local cause and effect precept which permeates the physics of denser media is modified for electrons in the solar wind. The local form of transport laws and equations of state which apply to collision dominated plasmas are replaced with global relations that explicitly depend on the relative position of the observer to the boundaries of the system.

Scudder, J. D.↗

NASA NTRS1994

Comparison of the plasma tails of four comets: P/Halley, Okazaki-Levy-Rudenko, Austin, and Levy

Photographic and charge coupled device (CCD) plasma tail observations are compared for four comets: P/Halley (22 nights in 1985/1986), Okazaki-Levy-Rudenko 1989 XIX (1989 December 2), Austin 1990 V (nine nights in 1990), and Levy 1990 XX (two nights in 1991). We present a discussion of several image-processing techniques used to enhance the visibility of the plasma tail features in order to measure velocities, accelerations, and position angles. The data are used to assess the validity of various physical mechanisms proposed to explain plasma tail phenomena. Seven disconnection events were observed in the comet P/Halley data, two in the Austin data, and none for the other comets. Analysis of these data suggests that while the crossing of the solar neutral sheet (the sector boundary) is a prominent factor in the production of a disconnection event, it is likely that several mechanisms are at work. A sector boundary crossing has been ruled out as the cause of either the 1986 April 26 P/Halley disconnection or the 1990 May 5/6 Austin disconnection. The motions of the disconnection events, knots, and condensations in the tails were seen to increase from 30-60 km/s near the nucleus (within 10(exp 6) km) to 80-100 km/s at 10(exp 7) km, consistent with either bulk motion or Alfven waves. Distinguishing between the two cases is not possible with these data. It was found that although the tail ray rotation rate slows as the ray approaches the tail axis, it is not a good indicator of the solar wind speed. Historical plasma tail data are also used to look for clues as to why some comets form well-developed plasma tails and others do not.

Farnham, Tony L.↗

DOE OSTI2023

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗