Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics”

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,441 records · Page 80

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Infrared compatible rapid mixer to probe millisecond chemical kinetics

Fast microfluidic mixers are a valuable tool for studying solution-phase chemical reaction kinetics and molecular processes with spectroscopy. However, microfluidic mixers that are compatible with infrared vibrational spectroscopy have seen only limited development due to the poor infrared transparency of the current microfabrication material. We describe the design, fabrication, and characterization of CaF2-based continuous flow turbulent mixers, which are capable of measuring kinetics in the millisecond time window with infrared spectroscopy, when integrated into an infrared microscope. Kinetics measurements demonstrate the ability to resolve relaxation processes with 1 millisecond time resolution, and straightforward improvements are described that should result in sub-100 µs time-resolution.

Instruments & Instrumentation↗

The kinetics of the PuO 2 to Pu 2 O 3 conversion

Here in an oxidizing environment, the oxide formed on plutonium (Pu) metal is composed of a plutonium dioxide (PuO 2 ) top layer and a thin cubic plutonium sesquioxide (Pu 2 O 3 ) middle layer. In a reducing environment, the PuO 2 layer auto-reduces to cubic Pu 2 O 3 . The speed and extent of this conversion depend on the combination of temperature and time. While PuO 2 provides a strong diffusion barrier against unwanted Pu corrosion by gaseous species (like hydrogen), Pu 2 O 3 does not, since its crystal structure has chains of oxygen vacancies. The kinetics of the PuO 2 reduction are, therefore, of fundamental interest and enable researchers to better protect Pu from corrosion. In this report, the oxygen-diffusion-limited kinetics of the dioxide to sesquioxide conversion were obtained by dynamically heating a PuO 2 -covered Pu sample from 294 to 418 K in a high-vacuum vessel equipped with an in situ spectroscopic ellipsometer. The physical/chemical constraints in the conversion process were combined with the ellipsometry method of multi-sample analysis to track the percentage of PuO 2 and to compute the extent of Pu 2 O 3 formation. The resulting diffusion coefficients were compared against and then combined with complementary literature data to produce a comprehensive set of kinetic parameters for reliably modeling oxide conversion over a larger temperature range than spanned by prior studies. The extracted thermal activation energy barrier (43.7 kJ/mol) and pre-exponential factor (5.0 × 10 -10 cm 2 /s) for the oxygen-diffusion-limited process can be used to accurately model the PuO 2 to Pu 2 O 3 transformation in vacuum and/or inert gas applications.

36 MATERIALS SCIENCE↗

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerated kinetic model for global macro stability studies of high-beta fusion reactors

The field reversed configuration (FRC), such as studied in the C-2W experiment at TAE Technologies, is an attractive candidate for realizing a nuclear fusion reactor. In an FRC, kinetic ion effects play the majority role in macroscopic stability, which allows global stability studies to make use of fluid-kinetic hybrid (also referred to as Ohm's law) models wherein ions are treated kinetically while electrons are treated as a fluid. The development and validation of such a hybrid particle-in-cell algorithm in the Exascale Computing Project code WarpX are reported here. Implementation of this model in the WarpX framework benefits from the numerical efficiency of WarpX as well as its scalability on large HPC systems and portability to different architectures. Performance benchmarks of the new algorithm for large, 3-dimensional, full device simulations from the Perlmutter supercomputer are presented. Results of a series of FRC simulations are discussed in which the impact of two-fluid effects on the tilt-mode growth rate was studied. It was observed that, in agreement with previous Hall-MHD studies, two-fluid effects have a stabilizing impact on the tilt mode.

Physics↗

Scale filtering analysis of kinetic reconnection and its associated turbulence

Previously, using an incompressible von Kármán–Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges [Adhikari et al., Phys. Plasmas 30, 082904 (2023)]. However, this description did not account for the energy transfer in the dissipation range for kinetic plasmas. In this Letter, we adopt a scale-filtering approach to investigate this previously unaccounted-for energy transfer channel in reconnection. Using kinetic particle-in-cell simulations of antiparallel and component reconnection, we show that the pressure–strain interaction becomes important at scales smaller than the ion inertial length, where the nonlinear energy transfer term drops off. Also, the presence of a guide field makes a significant difference in the morphology of the scale-filtered energy transfer. These results are consistent with kinetic turbulence simulations, suggesting that the pressure strain interaction is the dominant energy transfer channel between electron scales and ion scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗

KINETIC INSTABILITY OF DRIFT-ALFVEN WAVES IN SOLAR CORONA AND STOCHASTIC HEATING

The solar atmosphere is structured and inhomogeneous, both horizontally and vertically. The omnipresence of coronal magnetic loops implies gradients of the equilibrium plasma quantities such as the density, magnetic field, and temperature. These gradients are responsible for the excitation of drift waves that grow both within the two-component fluid description (both in the presence of collisions and without it) and within the two-component kinetic descriptions (due to purely kinetic effects). In this work, the effects of the density gradient in the direction perpendicular to the magnetic field vector are investigated within the kinetic theory, in both electrostatic (ES) and electromagnetic (EM) regimes. The EM regime implies the coupling of the gradient-driven drift wave with the Alfven wave. The growth rates for the two cases are calculated and compared. It is found that, in general, the ES regime is characterized by stronger growth rates, as compared with the EM perturbations. Also discussed is the stochastic heating associated with the drift wave. The released amount of energy density due to this heating should be more dependent on the magnitude of the background magnetic field than on the coupling of the drift and Alfven waves. The stochastic heating is expected to be much higher in regions with a stronger magnetic field. On the whole, the energy release rate caused by the stochastic heating can be several orders of magnitude above the value presently accepted as necessary for a sustainable coronal heating. The vertical stratification and the very long wavelengths along the magnetic loops imply that a drift-Alfven wave, propagating as a twisted structure along the loop, in fact occupies regions with different plasma-{beta} and, therefore, may have different (EM-ES) properties, resulting in different heating rates within just one or two wavelengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

KINETIC VERSUS MULTI-FLUID APPROACH FOR INTERSTELLAR NEUTRALS IN THE HELIOSPHERE: EXPLORATION OF THE INTERSTELLAR MAGNETIC FIELD EFFECTS

We present a new three-dimensional (3D) self-consistent two-component (plasma and neutral hydrogen) model of the solar wind interaction with the local interstellar medium (LISM). This model (K-MHD) combines the magnetohydrodynamic treatment of the solar wind and the ionized LISM component with a kinetic model of neutral interstellar hydrogen (LISH). The local interstellar magnetic field (B{sub LISM}) intensity and orientation are chosen based on an early analysis of the heliosheath flows. The properties of the plasma and neutrals obtained using the K-MHD model are compared to previous multi-fluid and kinetic models. The new treatment of LISH revealed important changes in the heliospheric properties not captured by the multi-fluid model. These include a decrease in the heliocentric distance to the termination shock (TS), a thinner heliosheath, and a reduced deflection angle ({theta}) of the heliosheath flows. The asymmetry of the TS, however, seems to be unchanged by the kinetic aspect of the LISH.

79 ASTRONOMY AND ASTROPHYSICS↗

SUPPRESSION OF REFLECTED ELECTRONS BY KINETIC ALFVEN TURBULENCE IN A QUASI-PERPENDICULAR SHOCK: PARTICLE-IN-CELL SIMULATIONS

Shock drift acceleration is one of the important mechanisms for electron acceleration associated with magnetic mirror reflection along the magnetic field in a quasi-perpendicular collisionless shock. We study the influence of a rippled shock surface in the in-plane magnetic field on the magnetic mirror reflection of electrons. Simulation results show that the number of reflected electrons reduces after generation of the rippled shock surface. Electric and magnetic wavenumber spectra of the generated fluctuations in the shock transition region indicate the existence of kinetic Alfven turbulence. The kinetic Alfven turbulence decreases the electron pitch angle by parallel scattering, which reduces the magnetic mirror force acting on the electrons. These results suggest that the shock-generated kinetic Alfven turbulence suppresses the magnetic mirror reflection of electrons during the shock drift acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

A CROSSED MOLECULAR BEAM, LOW-TEMPERATURE KINETICS, AND THEORETICAL INVESTIGATION OF THE REACTION OF THE CYANO RADICAL (CN) WITH 1,3-BUTADIENE (C{sub 4}H{sub 6}). A ROUTE TO COMPLEX NITROGEN-BEARING MOLECULES IN LOW-TEMPERATURE EXTRATERRESTRIAL ENVIRONMENTS

We present a joint crossed molecular beam and kinetics investigation combined with electronic structure and statistical calculations on the reaction of the ground-state cyano radical, CN(X {sup 2}{Sigma}{sup +}), with the 1,3-butadiene molecule, H{sub 2}CCHCHCH{sub 2}(X {sup 1} A{sub g}), and its partially deuterated counterparts, H{sub 2}CCDCDCH{sub 2}(X {sup 1} A{sub g}) and D{sub 2}CCHCHCD{sub 2}(X {sup 1} A{sub g}). The crossed beam studies indicate that the reaction proceeds via a long-lived C{sub 5}H{sub 6}N complex, yielding C{sub 5}H{sub 5}N isomer(s) plus atomic hydrogen under single collision conditions as the nascent product(s). Experiments with the partially deuterated 1,3-butadienes indicate that the atomic hydrogen loss originates from one of the terminal carbon atoms of 1,3-butadiene. A combination of the experimental data with electronic structure calculations suggests that the thermodynamically less favorable 1-cyano-1,3-butadiene isomer represents the dominant reaction product; possible minor contributions of less than a few percent from the aromatic pyridine molecule might be feasible. Low-temperature kinetics studies demonstrate that the overall reaction is very fast from room temperature down to 23 K with rate coefficients close to the gas kinetic limit. This finding, combined with theoretical calculations, indicates that the reaction proceeds on an entrance barrier-less potential energy surface (PES). This combined experimental and theoretical approach represents an important step toward a systematic understanding of the formation of complex, nitrogen-bearing molecules-here on the C{sub 5}H{sub 6}N PES-in low-temperature extraterrestrial environments. These results are compared to the reaction dynamics of D1-ethynyl radicals (C{sub 2}D; X {sup 2}{Sigma}{sup +}) with 1,3-butadiene accessing the isoelectronic C{sub 6}H{sub 7} surface as tackled earlier in our laboratories.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

KINETIC ALFVEN WAVE INSTABILITY DRIVEN BY FIELD-ALIGNED CURRENTS IN SOLAR CORONAL LOOPS

Magneto-plasma loops, which trace closed solar magnetic field lines, are the primary structural elements of the solar corona. Kinetic Alfven wave (KAW) can play an important role in inhomogeneous heating of these magneto-plasma structures in the corona. By the use of a low-frequency kinetic dispersion equation, which is presented in this paper and is valid in a finite-{beta} plasma with Q < {beta} < 1 plasma (where {beta} is the kinetic to magnetic pressure ratio and Q = m{sub e} /m{sub i} is the mass ratio of electrons to ions), KAW instability driven by a field-aligned current in the current-carrying loops in the solar corona is investigated. The results show that the KAW instability can occur in wave number regimes 0 < k{sub z} < k{sup c}{sub z} and 0 < k < k{sup c} , and that the critical wave numbers k{sup c}{sub z} and k{sup c} and the growth rate both considerably increase as the drift velocity V{sub D} of the current-carrying electrons increases in the loops. In particular, for typical parameters of the current-carrying loops in the solar corona this instability mechanism results in a high growth rate of KAWs, {omega}{sub i} {approx} 0.01-0.1{omega}{sub ci} {approx} 10{sup 3}-10{sup 4} s{sup -1}. The results are of importance in understanding the physics of the electric current dissipation and plasma heating of the current-carrying loops in the solar corona.

79 ASTRONOMY AND ASTROPHYSICS↗

A measurement of the kinetic Sunyaev-Zel'dovich signal toward MACS J0717.5+3745

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km s{sup –1} (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit, line-of-sight velocity of v{sub z} = +3450 ± 900 km s{sup –1}, with (1 – Prob[v{sub z} ≥ 0]) = 1.3 × 10{sup –5} (4.2σ away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v{sub z} = +2550 ± 1050 km s{sup –1}, with (1 – Prob[v{sub z} ≥ 0]) = 2.2 × 10{sup –3} (2.9σ).

79 ASTRONOMY AND ASTROPHYSICS↗

Active galactic nucleus feedback in an isolated elliptical galaxy: The effect of strong radiative feedback in the kinetic mode

Based on two-dimensional high-resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central active galactic nucleus (AGN) on the cosmological evolution of an isolated elliptical galaxy. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. It is well known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being ≅ 10{sup –3.5}), and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

79 ASTRONOMY AND ASTROPHYSICS↗