Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma kinetic simulation”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A predicted energetic ion signature of a neutral line in the geomagnetic tail

Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.

Martin, R. F., Jr.↗
NASA NTRS2015

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗
DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Ion Sources and Acceleration Mechanisms Inferred from Local Distribution Functions

This study investigates the sources of the ions up the complex and nonisotropic H(+) velocity distribution functions observed by the Geotail spacecraft on May 23, 1995, in the near-Earth magnetotail region and recently reported by Frank et al. [1996]. A distribution function observed by Geotail at -10 R(sub E) downtail is used as input for the large scale kinetic (LSK) technique to follow the trajectories of approximately 90,000 H(+) ions backward in time. Time-dependent magnetic and electric fields are taken from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interactions with appropriate solar wind and IMF conditions. The ion population described by the Geotail distribution function was found to consist of a mixture of particles originating from three distinct sources: the ionosphere, the low latitude boundary layer (LLBL), and the high latitude plasma mantle. Ionospheric particles had direct access along field lines to Geotail, and LLBL ions convected adiabatically to the Geotail location. Plasma mantle ions, on the other hand, exhibited two distinct types of behavior. Most near-Earth mantle ions reached Geotail on adiabatic orbits, while distant mantle ions interacted with the current sheet tailward of Geotail and had mostly nonadiabatic orbits. Ions from the ionosphere, the LLBL, and the near-Earth mantle were directly responsible for the well-separated, low energy structures easily discernible in the observed and modeled distribution functions. Distant mantle ions formed the higher energy portion of the Geotail distribution. Thus, we have been successful in extracting useful information about particle sources, their relative contribution to the measured distribution and the acceleration processes that affected particle transport during this time.

Ashour-Abdalla, M.↗
NASA NTRS2002

Human muscle sympathetic nerve activity and plasma noradrenaline kinetics in space

Astronauts returning from space have reduced red blood cell masses, hypovolaemia and orthostatic intolerance, marked by greater cardio-acceleration during standing than before spaceflight, and in some, orthostatic hypotension and presyncope. Adaptation of the sympathetic nervous system occurring during spaceflight may be responsible for these postflight alterations. We tested the hypotheses that exposure to microgravity reduces sympathetic neural outflow and impairs sympathetic neural responses to orthostatic stress. We measured heart rate, photoplethysmographic finger arterial pressure, peroneal nerve muscle sympathetic activity and plasma noradrenaline spillover and clearance, in male astronauts before, during (flight day 12 or 13) and after the 16 day Neurolab space shuttle mission. Measurements were made during supine rest and orthostatic stress, as simulated on Earth and in space by 7 min periods of 15 and 30 mmHg lower body suction. Mean (+/- S.E.M.) heart rates before lower body suction were similar pre-flight and in flight. Heart rate responses to -30 mmHg were greater in flight (from 56 +/- 4 to 72 +/- 4 beats min(-1)) than pre-flight (from 56 +/- 4 at rest to 62 +/- 4 beats min(-1), P < 0.05). Noradrenaline spillover and clearance were increased from pre-flight levels during baseline periods and during lower body suction, both in flight (n = 3) and on post-flight days 1 or 2 (n = 5, P < 0.05). In-flight baseline sympathetic nerve activity was increased above pre-flight levels (by 10-33 %) in the same three subjects in whom noradrenaline spillover and clearance were increased. The sympathetic response to 30 mmHg lower body suction was at pre-flight levels or higher in each subject (35 pre-flight vs. 40 bursts min(-1) in flight). No astronaut experienced presyncope during lower body suction in space (or during upright tilt following the Neurolab mission). We conclude that in space, baseline sympathetic neural outflow is increased moderately and sympathetic responses to lower body suction are exaggerated. Therefore, notwithstanding hypovolaemia, astronauts respond normally to simulated orthostatic stress and are able to maintain their arterial pressures at normal levels.

STS-90 Shuttle Project↗
DOE OSTI2025

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗
NASA NTRS1984

Theory and modeling of atmospheric turbulence, part 2

Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.

Chen, C. M.↗
NASA NTRS1991

A simulation study of impulsive penetration of solar wind irregularities into the magnetosphere at the dayside magnetopause

The impulsive penetration processes that occur when a plasma irregularity in the magnetosheath, modeled as a field-aligned filament, impinges on the dayside magnetopause, are investigated via a 2D magnetohydrodynamic code. It is found that, when the magnetic fields of the magnetosheath and magnetosphere are parallel or antiparallel, even plasma filaments with small initial velocities can eventually penetrate through the magnetopause and enter the magnetosphere. If the magnetic fields are not aligned, penetration will occur only if the initial kinetic energy density of the plasma filament is more than 50 times larger than the magnetic energy density due to the transverse component of the magnetic field.

Ma, Z. W.↗
NASA NTRS1993

Self-consistent current sheet structures in the quiet-time magnetotail

The structure of the quiet-time magnetotail is studied using a test particle simulation. Vlasov equilibria are obtained in the regime where v(D) = E(y) c/B(z) is much less than the ion thermal velocity and are self-consistent in that the current and magnetic field satisfy Ampere's law. Force balance between the plasma and magnetic field is satisfied everywhere. The global structure of the current sheet is found to be critically dependent on the source distribution function. The pressure tensor is nondiagonal in the current sheet with anisotropic temperature. A kinetic mechanism is proposed whereby changes in the source distribution results in a thinning of the current sheet.

Holland, Daniel L.↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Low-frequency mobility response functions for the central plasma sheet with application to tearing modes

Consideration is given to the effect of constant cross-tail magnetic field By on the collisionless conductivity produced by chaotic scattering and stochastic diffusion of particles in the current sheet for a parabolic geometry. It is shown that the correlation time scales as (By/Bz)-squared, and from this strong By scaling a strong tendency toward stabilization of the linear tearing modes with increasing values of By is inferred. This effect of increased dawn-dusk mobility is particularly dramatic when electrons are introduced in the calculation, and is in agreement with the results of kinetic particle simulations. The collisionless conductivity is expressed in terms of the ensemble-averaged power spectrum of the single particle trajectories, which makes it possible to calculate directly the linear conductivity instead of deriving it from the calculation of the irreversible heating rates.

Hernandez, J.↗
DOE OSTI2024

Kinetics and Mechanism of the Singlet Oxygen Atom Reaction with Dimethyl Ether

Here, we combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H 2 O, HO 2 , and CH 2 O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH 3 OOCH 3 ) generated via O-insertion into the C–O bond undergoes subsequent dissociation to CH 3 O + CH 3 O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Laboratory Observations of Electron Heating and Non-Maxwellian Distributions at the Kinetic Scale during Electron-Only Magnetic Reconnection

Non-Maxwellian electron velocity distribution functions composed of a warm bulk population and a cold beam are directly measured during electron-only reconnection with a strong out-of-plane (guide) magnetic field in a laboratory plasma. Electron heating is localized to the separatrix, and the electron temperature increases continuously along the separatrix. The measured gain in enthalpy flux is 70% of the incoming Poynting flux. The electron beams are oppositely directed on either side of the X point, and their velocities are comparable to, and scale with, the electron Alfven speed. Particle-in-cell simulations are consistent with the measurements. Furthermore, the experimental results are consistent with, and go beyond, recent observations in the magnetosheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Effect of wave-particle interaction on the outflow of ions at high latitudes

The objective of the research is to use the low-frequency wave spectrum measured by the Plasma Wave Instrument on the DE-1 spacecraft to include the wave-particle interaction (WPI) in the different polar wind models that are available at Utah State University. A Monte Carlo technique was used to simulate the ion diffusion in the velocity space due to scattering by the waves. This enabled us to study the effect of WPI on the magnitude as well as the composition of the outflow of the ionospheric ions. In particular, in the first year the 1-D semi-kinetic code developed by Barakat and Schunk (1983) was modified in order to include the effect of the WPI.

Barakat, Abdallah R.↗
NASA NTRS1996

Progress Toward an Efficient and General CFD Tool for Propulsion Design/Analysis

The simulation of propulsive flows inherently involves chemical activity. Recent years have seen substantial strides made in the development of numerical schemes for reacting flowfields, in particular those involving finite-rate chemistry. However, finite-rate calculations are computationally intensive and require knowledge of the actual kinetics, which are not always known with sufficient accuracy. Alternatively, flow simulations based on the assumption of local chemical equilibrium are capable of obtaining physically reasonable results at far less computational cost. The present study summarizes the development of efficient numerical techniques for the simulation of flows in local chemical equilibrium, whereby a 'Black Box' chemical equilibrium solver is coupled to the usual gasdynamic equations. The generalization of the methods enables the modelling of any arbitrary mixture of thermally perfect gases, including air, combustion mixtures and plasmas. As demonstration of the potential of the methodologies, several solutions, involving reacting and perfect gas flows, will be presented. Included is a preliminary simulation of the SSME startup transient. Future enhancements to the proposed techniques will be discussed, including more efficient finite-rate and hybrid (partial equilibrium) schemes. The algorithms that have been developed and are being optimized provide for an efficient and general tool for the design and analysis of propulsion systems.

Cox, C. F.↗
DOE OSTI2022

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗