Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DYNAMICS”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗
DOE OSTI2018

Nonlinear Wave Structures and Plasma−Dust Effects in the Earth’s Atmosphere

The interaction of charged dust grains with nonlinear vortical structures in the Earth’s atmosphere is analyzed. Certain aspects of the atmosphere−ionosphere interaction, in particular, mechanisms for the appearance of dust grains at ionospheric altitudes, are discussed. It is shown that, at certain altitudes, there are regions in the wavenumber space in which conditions leading to the excitation of acoustic−gravity waves are satisfied. The interaction of nonlinear acoustic−gravity waves with dust grains of meteoric origin at ionospheric altitudes, which leads to the mixing and redistribution of dust grains over the region where vortices exist, is investigated. The possibility of formation of vertical and horizontal dust flows in dusty ionospheric plasma as a result of modulational instability is analyzed. The dynamics of dust grains in dust devils frequently arising in the atmosphere above well-heated surfaces is modeled. The vortical structure of such a dust devil is characterized by a reduced pressure in the center, which facilitates the lifting of small dust grains from the surface. The formulated model is used to calculate the trajectories of dust grains in dust devils with allowance for the influence of the electric field generated in the vortex by colliding dust grains. The calculations show that dust devils play an important role in the transport of dust grains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging (GPI) at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g. ELM filaments) in the scrape-off layer and the at the plasma edge with 2.5 us temporal and 10 mm optical resolution.A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics, and then extended to a database of 169 ELM events.Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average an ELM crash consists of a filament with a circular cross-section which is propelled outwards with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100us until the radial propulsion settles back to the pre-ELM level. The experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current-filament interaction model.

correlation↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of Coulomb collisions on lower hybrid drift waves inside a laboratory reconnection current sheet

We have developed a local, linear theoretical model for lower hybrid drift waves that can be used for plasmas in the weakly collisional regime. Two cases with typical plasma and field parameters for the current sheet of the magnetic reconnection experiment have been studied. For a case with a low electron beta (β e =0.25, high guide field case), the quasi-electrostatic lower hybrid drift wave is unstable, while the electromagnetic lower hybrid drift wave has a positive growth rate for a high-β e case (β e =8.9, low guide field case). For both cases, including the effects of Coulomb collisions reduces the growth rate but collisional impacts on the dispersion and growth rate are limited (≲20%).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗
DOE OSTI2024

Fundamental data for modeling electron-induced processes in plasma remediation of perfluoroalkyl substances

Plasma treatment of per- and polyfluoroalkyl substances (PFAS) contaminated water is a potentially energy efficient remediation method. In this treatment, an atmospheric pressure plasma interacts with surface-resident PFAS molecules. Developing a reaction mechanism and modeling of plasma–PFAS interactions requires fundamental data for electron–molecule reactions. In this paper, we present results of electron scattering calculations, potential energy landscapes and their implications for plasma modelling of a dielectric barrier discharge in PFAS contaminated gases, a first step towards modelling of plasma–water–PFAS intereactions. It is found that the plasma degradation of PFAS is dominated by dissociative electron attachment with the importance of other contributing processes varying depending on the molecule. All molecules posses a large number of shape resonances – transient negative ion states – from near-threshold up to ionization threshold. These states lie in the region of the most probable electron energies in the plasma (4–5 eV) and consequently are expected to further enhance the fragmentation dynamics in both dissociative attachment and dissociative excitation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Coincident 1.3-year Periodicities in the ap Geomagnetic Index and the Solar Wind

Recent observations show an approximately 1.3-year period in the speed of the solar wind detected by the IMP 8 and Voyager 2 spacecraft. A similar period is also seen in the north-south (GSE) component of the magnetic field observed by IMP 8. Since both parameters are commonly used as input to models of geomagnetic activity, the 'ap' index (a measure of geomagnetic disturbance) is examined to look for this periodicity. The Lomb-Scargle periodogram method is used on the ap, plasma, and magnetic field data during the 1973-1994 time range. A dynamic FFT periodogram method is also used to analyze the ap data during this time, as well as to look for periods present between 1932 and 1972. A clear 1.3-year periodicity is present in the post-1986 data when the same period is observed in the plasma and field data. The V(2)B(zsm) and V(2)B(s) proxies for geomagnetic activity also show this periodicity. However, the southward (GSM) component of the magnetic field does not have a 1.3-year period, and neither do solar wind or ap data from 1973-1985. This demonstrates that the ap geomagnetic index can act as a proxy for solar wind periodicities at this time scale. Historic ap data are examined, and show that a similar periodicity in ap exists around 1942. Since auroral data show a 1.4-year periodicity, all these similar periods may result from a common underlying solar mechanism.

Paularena, K. I.↗
DOE OSTI2017

Plasma-Generated OH Radical Production for Analyzing Three-Dimensional Structure in Protein Therapeutics

Protein three-dimensional structure dynamically changes in solution depending on the presence of ligands and interacting proteins. Methods for detecting these changes in protein conformation include ‘protein footprinting,’ using mass spectrometry. We describe herein a new technique, PLIMB (Plasma Induced Modification of Biomolecules), that generates µs bursts of hydroxyl radicals from water, to measure changes in protein structure via altered solvent accessibility of amino acid side chains. PLIMB was first benchmarked with model compounds, and then applied to a biological problem, i.e., ligand (EGF) induced changes in the conformation of the external (ecto) domain of Epidermal Growth Factor Receptor (EGFR). Regions in which oxidation decreased upon adding EGF fall along the dimerization interface, consistent with models derived from crystal structures. Furthermore these results demonstrate that plasma-generated hydroxyl radicals from water can be used to map protein conformational changes, and provide a readily accessible means of studying protein structure in solution.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2013

Prominence Mass Supply and the Cavity

A prevalent but untested paradigm is often used to describe the prominence-cavity system; the cavity is under-dense because it it evacuated by supplying mass to the condensed prominence. The thermal non-equilibrium (TNE) model of prominence formation offers a theoretical framework to predict the thermodynamic evolutin of the prominence and the surrounding corona. We examine the evidence for a prominence-cavity connection by comparing the TNE model and diagnostics of dynamic extreme ultraviolet (EUV) emission surrounding the prominence, specifically prominence horns. Horns are correlated extensions of prminence plasma and coronal plasma which appear to connect the prominence and cavity. The TNE model predicts that large-scale brightenings will occur in the Solar Dynamics Observatory Atmospheric Imaging Assembly 171 A badpass near he prominence that are associated with the cooling phase of condensation formation. In our simulations, variations in the magnitude of footpoint heating lead to variations in the duration, spatial scale, and temporal offset between emission enhancements in the other EUV bandpasses. While these predictions match well a subset of the horn observations, the range of variations in the observed structures is not captured by the model. We discuss the implications of one-dimensional loop simulations for the three-dimensional time-averaged equilibrium in the prominence and the cavity. Evidence suggests that horns are likely caused by condensing prominence plasma, but the larger question of whether this process produces a density-depleted cavity requires a more tightly constrained model of heating and better knowledge of the associated magnetic structure.

PROMINENCES↗
DOE OSTI2024

Physics-assisted Latent Space Dynamics Learning for Stiff Collisional-radiative Models

Collisional-radiative (CR) models describe the atomic processes in a plasma by tracking the population density in the ground and excited states for each charge state of the atom/ion. These models predict important plasma properties such as charge state distributions and radiative emissivity and opacity. Accurate descriptions of the CR balance of the plasma are essential in fusion whole device modeling, especially when significant impurities are introduced into the plasmas. In an integrated fusion plasma and CR simulation, a CR model, which is a high-dimensional stiff ODE, needs to be solved on each grid point in the configuration space, and can overwhelm the overall computational cost. In this work, we propose a machine-learning-based method that discovers a latent space and learns its corresponding latent dynamics, which can capture the essential physics to make accurate predictions at much lower online computational cost. The proposed approach is physics-assisted, due to its combination of a physical latent space and a data-driven latent space. It has been demonstrated that the proposed architecture can predict both the full-order CR dynamics and a physical quantity of interest accurately.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗