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 1,027 records · Page 57

Kinetic Simulations of Current-Sheet Formation and Reconnection at a Magnetic X Line

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.

Black, C.↗

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗

Thermal Gradient Effects on Redox Evolution and Volatility-Driven Fractionation in Ternary U/Ce/Cs Condensates

Understanding how thermal history influences redox evolution and chemical fractionation is essential for characterizing high-temperature condensation in complex materials, including nuclear debris. Here, we tested the hypothesis that distinct thermal regimes in a plasma flow reactor influence redox pathways and elemental partitioning in ternary U/Ce/Cs systems. A configurable plasma flow reactor was modified with an external tube furnace to impose two distinct thermal gradients: continuous ambient cooling and a furnace-assisted thermal hold-up near 1400 K followed by rapid cooling. Transmission electron microscopy characterized phase identity, morphology, and nanoscale element distributions, while inductively coupled plasma-mass spectrometry quantified bulk elemental ratios. Across both thermal regimes, uranium and cerium condensed as UO 2 and CeO 2 as dominant refractory oxide products. Uranium partially oxidized to α-UO 3 during extended ambient cooling, while furnace-assisted hold-up preserved UO 2 and produced partial reduction of cerium to Ce 2 O 3 . Cesium remained volatile upstream and condensed later in the reactor, forming Cs 2 O and Cs-uranate phases with the highest incorporation after thermal hold-up. Bulk ICP-MS measurements supported these observations. U/Ce ratios remained comparatively stable and Cs displayed delayed and apparent transient enrichment that matched the nanoscale measurements. This integrated approach provides a quantitative method for linking thermal gradients to redox evolution and volatility-driven fractionation. These results show how the plasma flow reactor can identify where equilibrium descriptions remain adequate and where kinetic effects from residence time and temperature history must be considered when interpreting condensation behavior in multicomponent systems.

and nuclear chemistry↗

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reply to Comment on ‘Deuterium supersaturated surface layer in tungsten: ion energy dependence’

Abstract We reply to the comment by Li et al (submitted to Nucl. Fusion with this response) on our recent paper Nishijima et al (2023 Nucl. Fusion 63 126003). In this response, we address the existence of an incident ion energy, E i , threshold for the deuterium (D) supersaturated surface layer (DSSL) formation with a newly conducted D plasma exposure experiment at E i ∼ 20 eV. It is also further demonstrated, based on new experiments where the ion flux, fluence, and sample temperature are scanned, that non-kinetic (ballistic) processes play a role in the DSSL formation and growth. In addition, the effect of impurities in our plasma is discussed also with a new analysis of the surface composition made after a D plasma exposure.

Nishijima, D. (ORCID:0000000231194827)↗

Particle-in-cell simulations of burning inertial confinement fusion capsule implosions

Anomalies observed in the neutron spectral shift of high-yield shots at the National Ignition Facility (NIF) suggest the presence of suprathermal ions [E. P. Hartouni et al., Nat. Phys. 19, 72 (2023)], implying that kinetic effects play a significant role in burning inertial confinement fusion (ICF) plasmas. Furthermore, recent measurements of reaction-in-flight (RIF) neutrons offer a direct probe of the stopping power in the burning fuel region of high-energy alpha particles and up-scattered fuel ions. We have developed the particle-in-cell code PICNIC, an exactly energy-conserving particle-in-cell Monte-Carlo collision (PIC-MCC) code to simulate the burn stage in ICF. We present results from 1D spherical simulations of NIF shot N210808, which was the first to exceed the Lawson criterion for laser fusion. We find that the suprathermal ions generated by large-angle Rutherford and nuclear elastic scattering (NES) with fusion alphas produce an alpha knock-on neutron (AKN) signal consistent with the extent and relative yield of the AKN spectrum identified in ignition experiments at the NIF. Furthermore, we also find that the inclusion of large-angle scattering physics does not explain the anomalously large spectral shift observed in the experiment.

High-energy-density plasmas↗

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Capabilities Development at the University of Texas at El Paso for Hydrogen Generation Research and Education

Gasification-based systems have recently received much attention due to their capability of converting wastes into useful fuels. The gasification of biomass, municipal solid waste (MSW), plastics, and other sustainable resources has shown promise in hydrogen production. When coupled to a CCU or CCS unit, these systems have the potential to achieve carbon neutrality or carbon-negative emissions. In addition, these systems have an edge over conventional incineration or landfill systems by reducing greenhouse gas emissions and recovering useful energy. However, the gasification of MSW and other wastes is still in its early stages. In particular, co-gasification of wastes with biomass has significant unknowns in optimizing the reaction kinetics, operability, design and performance improvements. Currently, Supercritical Water Gasifiers (SCWG) and Plasma Gasifiers are the major systems that are capable of producing high hydrogen amounts from Biomass and MSW, respectively. However, both systems have high capital and operating costs, adversely affecting efficiency and the hydrogen production cost. Additionally, due to lower temperatures, SCWGs are prone to tar formation and fouling in the heat exchangers. Hence, there is a need to look for alternative solutions for MSW and biomass gasifications, particularly in co-gasification. The current effort aims to develop a strategic plan to establish a sustainable gasification facility for hydrogen research at the University of Texas at El (UTEP). A major part of this effort involves an extensive literature survey to identify technological gaps and potential areas of interest. Several concepts were developed in accordance with the current demands from the literature survey. Based on the concepts, a center-wide capability assessment was conducted to measure the current capacity and feasibility of establishing hydrogen research. Afterward, a strategic research plan was developed, and a list of required resources was made for the expansion of hydrogen research at UTEP. In addition, successful partnerships with the local county and the city were developed to pursue hydrogen research. The strategic goal setting and planning of UTEP Aerospace Center resulted in securing $2.5 million in external support to expand the hydrogen research during the project performance period. In addition, during the project period, a course in Hydrogen Energy Systems was developed to expand the energy curriculum at the UTEP Aerospace and Mechanical Engineering Department. The course focused on hydrogen production, storage, supply and delivery and application. The cross-listed course was offered at both undergraduate and graduate levels during the Spring 2024 semester and had 27 enrolled students. Moreover, the students supported under this award were trained in gasification process modeling, CAD, CFD and FEA during the project period. The training enabled the students to move into new projects to support gasification design, integrated gasification combined cycle process plant analysis, gasifier structural and operational analysis, and digitally threading the systems.

08 HYDROGEN↗

The diurnal heat budget of the thermosphere.

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photo-dissociation of O2 is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons; (2) energy exchange from thermal plasma; (3) chemical reaction (ion-electron dissociative recombination) energy gain; (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are (1) thermal conduction at the lower boundary (120 km); and (2) radiative cooling of atomic oxygen. It is shown that the combined energy from processes 2 to 4 constitutes only a small fraction of the total energy available from photoelectrons and is in phase with the latter. These secondary sources (processes 2 to 4), therefore, do not constitute a significant energy source and their contribution can be simply incorporated into photoelectron energy (process 1) by defining an effective photoionization heating efficiency.

Chandra, S.↗

Energy released by the interaction of coronal magnetic fields

Comparisons between coronal spectroheliograms and photospheric magnetograms are presented to support the idea that as coronal magnetic fields interact, a process of field-line reconnection usually takes place as a natural way of preventing magnetic stresses from building up in the lower corona. This suggests that the energy which would have been stored in stressed fields is continuously released as kinetic energy of material being driven aside to make way for the reconnecting fields. However, this kinetic energy is negligible compared with the thermal energy of the coronal plasma. Therefore, it appears that these slow adjustments of coronal magnetic fields cannot account for even the normal heating of the corona, much less the energetic events associated with solar flares.

Sheeley, N. R., Jr.↗

An example of the association of X-ray and UV emission with H-alpha surges

From H-alpha UV, and X-ray data, the nature of a well-observed surge on November 12, 1980 has been studied to try to understand the relationships between the mechanical motions and the high-temperature emissions. The cool (H-alpha) and the hot (O v) parts of the surge plasma both show velocities in the range of 100-120 km/s. The kinetic and potential energies of the surge are two orders of magnitude higher than the radiative energy. The observations suggest that the surge occurs in open structures parallel to one leg of a large scale coronal magnetic loop. The energy is released principally in the open structure (surge) and only a small amount heats the coronal loop (X-ray signature). This energy division should probably not be treated as a general characteristic of surge events. Different events are expected to exhibit a different energy balance, and this can account for the lack of a consensus in the previous literature regarding surge/X-ray associations.

Schmieder, B.↗

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗

Self-Consistent Simulation of the Magnetotail

Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.

Peroomian, Vahe↗

Nonlinear Waves and Particle Heating in the Auroral Ionosphere

In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Goldman, Martin↗

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗

Wave Telescope Technique for MMS Magnetometer

Multipoint measurements are a powerful method in studying wavefields in space plasmas.The wave telescope technique is tested against magnetic field fluctuations in the terrestrial magnetosheath measured by the four Magnetospheric Multiscale (MMS) spacecraft on a spatial scale of about 20 km.The dispersion relation diagram and the wave vector distribution are determined for the first time in the ion-kinetic range. Moreover, the dispersion relation diagram is determined in a proxy plasma restframe by regarding the low-frequency dispersion relation as a Doppler relation and compensating for the apparent phase velocity. Fluctuations are highly compressible, and the wave vectors have an angle of about 60 from the mean magnetic field. We interpret that the measured fluctuations represent akinetic-drift mirror mode in the magnetosheath which is dispersive and in a turbulent state accompanied by a sideband formation.

Magnetometer↗