Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 649 records · Page 36

Self-consistent electrostatic potential due to trapped plasma in the magnetosphere

A steady state solution for the self-consistent electrostatic potential due to a plasma confined in a magnetic flux tube is considered. A steady state distribution function is constructed for the trapped particles from the constants of the motion, in the absence of waves and collisions. Using Liouville's theorem, the particle density along the geomagnetic field is determined and found to depend on the local magnetic field, self-consistent electric potential, and the equatorial plasma distribution function. A hot anisotropic magnetospheric plasma in steady state is modeled by a bi-Maxwellian at the equator. The self-consistent electric potential along the magnetic field is calculated assuming quasineutrality, and the potential drop is found to be approximately equal to the average kinetic energy of the equatorially trapped plasma. The potential is compared with that obtained by Alfven and Faelthammar (1963).

Miller, Ronald H.↗

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sintering and Interface Strain Tolerance of Plasma-Sprayed Thermal and Environmental Barrier Coatings

Ceramic thermal and environmental barrier coatings will be more aggressively designed to protect gas turbine engine hot section SiC/SiC Ceramic Matrix Composite (CMC) components in order to meet future engine higher fuel efficiency and lower emission goals. A coating system consisting of a zirconia-based oxide topcoat (thermal barrier) and a mullite/BSAS silicate inner coat (environmental barrier) is often considered a model system for the CMC applications. However, the coating sintering, and thermal expansion mismatch between the zirconia oxide layer and the silicate environmental barrier/CMC substrate will be of major concern at high temperature and under thermal cycling conditions. In this study, the sintering behavior of plasma-sprayed freestanding zirconia-yttria-based thermal barrier coatings and mullite (and/or barium-strontium-aluminosilicate, i.e., BSAS) environmental barrier coatings was determined using a dilatometer in the temperature range of 1200-1500 C. The effects of test temperature on the coating sintering kinetics were systematically investigated. The plasma-sprayed zirconia-8wt.%yttria and mullite (BSAS) two-layer composite coating systems were also prepared to quantitatively evaluate the interface strain tolerance of the coating system under thermal cycling conditions based on the dilatomentry. The cyclic response of the coating strain tolerance behavior and interface degradation as a function of cycle number will also be discussed.

Zhu, Dongming↗

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗

Composition/Structure/Dynamics of comet and planetary satellite atmospheres

This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.

Combi, Michael R.↗

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

Furthering our understanding of many of today’s interesting problems in plasma physics – including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects – requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or ‘‘tiles,’’ along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. Here, we tested this implementation in the code Osiris and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.

97 MATHEMATICS AND COMPUTING↗

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗

Hybrid simulations of intermediate shocks - Coplanar and noncoplanar solutions

A hybrid code is used to investigate the kinetic structure and stability of subfast intermediate shocks (IS) formed dynamically by the interaction between a flowing plasma and a stationary piston. Results of the kinetic simulation of noncoplanar ISs are compared with predictions of the MHD theory. The relevance of the results of the study to observations of the magnetopause is discussed.

Karimabadi, H.↗

Plasma-Catalytic Ammonia Synthesis beyond the Equilibrium Limit

In this work, we explore the consequences of nonthermal plasma-activation on product yields in catalytic ammonia synthesis, a reaction that is equilibrium-limited at elevated temperatures. We employ a minimal microkinetic model that incorporates the influence of plasma-activation on N 2 dissociation rates to predict NH 3 yields into and across the equilibrium-limited regime. NH 3 yields are predicted to exceed bulk thermodynamic equilibrium limits on materials that are thermal-rate-limited by N 2 dissociation. In all cases, yields revert to bulk equilibrium at temperatures at which thermal reaction rates exceed plasma-activated ones. Beyond-equilibrium NH 3 yields are observed in a packed bed dielectric barrier discharge reactor and exhibit sensitivity to catalytic material choice in a way consistent with model predictions. The approach and results highlight the opportunity to exploit synergies between nonthermal plasmas and catalysts to affect transformations at conditions inaccessible through thermal routes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ringing After a High-Energy Collision: Ambipolar Oscillations During Impact Plasma Expansion

High-velocity impacts on the Moon and other airless bodies deliver energy and material to the lunar surface and exosphere. The target and i mpactor material may become vaporized and ionized to form a collision al plasma that expands outward and eventually becomes collisionless. In the present work, kinetic simulations of the later collision less stage of impact plasma expansion are performed. Attention is paid to characterizing "ambipolar oscillations" in which thermodynamic distur bances propagate outward to generate "ringing" within the expanding e lectron cloud, which could radiate an electromagnetic signature of lo cal plasma conditions. The process is not unlike a beam-plasma intera ction, with the perturbing electron population in the present case ac ting as a highly thermal "beam" that resonates along the expanding de nsity gradient. Understanding the electromagnetic aspects of impact p lasma expansion could provide insight into the lasting effects of nat ural, impact-generated currents on airless surfaces and charging haza rds to human exploration infrastructure and instrumentation.

Zimmerman, M. I.↗

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗

Continuum kinetic investigation of the impact of bias potentials in the current saturation regime on sheath formation

In this work, we examine sheath formation in the presence of bias potentials in the current saturation regime for pulsed power fusion experiments. It is important to understand how the particle and heat fluxes at the wall may impact the wall material and affect electrode degradation. Here, simulations are performed using the 1X-1V Boltzmann–Poisson system for a proton–electron plasma in the presence of bias potentials ranging from 0 to 10 kV. The results indicate that the sheath near the anode remains generally the same as that of a classical sheath without the presence of a bias potential. However, the sheath near the cathode becomes more prominent with a larger potential drop, a significant decrease in the electron density, and larger sheath lengths. The spatially constant current density increases to a saturation value with the increasing bias potential. For high bias potentials, the material choice needs to consider that the anode faces significantly larger particle and heat fluxes compared to the cathode. In general, the results trend with theory with differences attributed to the simplified assumptions in the theory and the kinetic effects considered in the simulations. Due to the significant computational cost of a well resolved 1X-2V simulation, only one such simulation is performed for the 5 kV case showing higher current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent resistivity, diffusion and heating

Experimental and theoretical studies are reported on ion acoustic and ion cyclotron turbulence and their roles in anomalous resistivity, viscosity, diffusion and heating and in the structure of collisionless electrostatic shocks. Resistance due to ion acoustic turbulence has been observed in experiments with a streaming cesium plasma in which electron current, potential rise due to turbulent resistivity, spectrum of unstable ion acoustic waves, and associated electron heating were all measured directly. Kinetic theory calculations for an expanding, unstable plasma, give results in agreement with the experiment. In a strong magnetic field, with T sub e/T sub i approximately 1 and current densities typical for present Tokomaks, the plasma is stable to ion acoustic but unstable to current driven electrostatic ion cyclotron waves. Relevant characteristics of these waves are calculated and it is shown that for ion, beta greater than m sub e/m sub i, the electromagnetic ion cyclotron wave has a lower instability threshold than the electrostatic one. However, when ion acoustic turbulence is present experiments with double plasma devices show rapid anomalous heating of an ion beam streaming through a plasma.

Fried, B. D.↗

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic instability in high-beta plasmas with generalized distribution functions

The electrostatic instabilities in an anisotropic plasma have been studied quite extensively, but the electromagnetic instabilities in high-beta plasmas, where they play an important role, have not been thoroughly investigated. Recently, Davidson and Wu (1970) looked into the ordinary mode electromagnetic instability which can arise in high-beta bi-Maxwellian plasmas. Here, the magnetic instability is discussed in (A, beta) space (A is the anisotropy in the temperature of the plasma species and beta is the ratio of the kinetic pressure to the magnetic pressure), which can occur in generalized non-Maxwellian plasmas with an inverted population of different species.

Buti, B.↗

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

Flame Chemistry Workshop: a perspective on challenges and strategic actions in combustion experiments and chemical kinetics modeling

Continued progress in the development of predictive models for combustion chemistry—including ignition and flame behavior, species evolution, and combustion system performance—relies on overcoming enduring and emerging challenges in experimental measurements, theoretical formulations, and chemical kinetics mechanism construction. As combustion science continues to coincide with advances in sustainable fuels development, plasma technologies, and automated modeling capabilities, the need for coordinated, community-driven strategies is essential. The Flame Chemistry Workshop (FCWS), held biennially before the International Symposium on Combustion, serves as a dedicated platform to identify, consolidate, and address these challenges in a structured and collaborative manner. This perspective arises from discussions at the 7th FCWS in Milan, Italy (2024), and presents a collective view of the critical barriers currently limiting progress. Across the five technical domains discussed during the 7th FCWS – sustainable fuels combustion, advanced diagnostics for combustion measurements, experiments and modeling in plasma combustion, artificial intelligence and automated methods for theory and mechanisms generation, and chemical kinetic models—a series of persistent and emerging scientific challenges were identified, highlighting the need for deeper integration between three areas: theory, experiments, and modeling. In conclusion, the present article concisely describes present challenges that were identified in each of the technical domains in an effort to streamline and coordinate solutions to accelerate progress in combustion science.

Chemical kinetics↗