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,099 records · Page 61

Method for Impeding Degradation of Porous Silicon Structures

This invention relates to surface modification of porosified silicon (pSi) structures with poly(alkylene) glycols for the purpose of controlled degradation of the silicon matrix and tailored release of encapsulated substances for biomedical applications. The pSi structures are currently used in diverse biomedical applications including bio-molecular screening, optical bio-sensoring, and drug delivery by means of injectable/orally administered carriers and implantable devices. The size of the pores and the surface chemistry of the pSi structure can be controlled during the microfabrication process and thereafter. A fine regulation of the degradation kinetics of mesoporous silicon structures is of fundamental importance. Polyethylene glycols (PEGs) represent the major category of surface modifying agents used in classical drug delivery systems and in pharmaceutical dosage forms. PEGylation enables avoidance of RES uptake, thus prolonging circulation time of intravenously injectable nanovectors. PEG molecules demonstrate little toxicity and immunogenicity, and are cleared from the body through the urine (molecular weight, MW less than 30 kDa) or in the feces (MW greater than 30kDa). The invention focuses on the possibility of finely tuning the degradation kinetics of the pSi nanovectors and other structures through surface conjugation of PEGs with various backbone lengths/MWs. To prove the concept, pSi nanovectors were covalently conjugated to seven PEGs with MW from 245 to 5,000 Da and their degradation kinetics in physiologically relevant media (phosphate buffer saline, PBS pH7.4, and fetal bovine serum) was assessed by the elemental analysis of the Si using inductive coupled plasma atomic emission spectroscopy (ICP-AES). The conjugation of the PEG with lowest MW to the nanovectors surface did not induce any change in the degradation kinetics in serum, but inhibited degradation and consequently the release of orthosilicic acid into buffer. When PEGs with the longer chains were evaluated, Si mass loss from the nanovectors was slowed down, and the PEGylated structures were almost fully degraded within 18 24 hours in serum and within 48 hours in PBS. The most dramatic effect was observed for high MW PEGs 3,400 and 5,000 Da, which prominently inhibited the degradation of the systems, with complete degradation achieved only after four days. For these PEGs, during the early stages of the degradation, there was a lag period of little or no Si mass loss from the nanovector. The obtained profiles were in agreement with the erosion of the nanovector surface as observed by scanning electron microscopy.

Vilentchouk, Biana Godin↗

Bolt: A Fast Solver for Kinetic Theories Using a High-Resolution Constrained Transport Scheme [Slides]

Our understanding of collisionless and semi-collisional plasmas in the nonlinear regime is limited by the expense of computing solutions numerically. Bolt is a fast, GPU-accelerated code for rapidly computing such solutions with accurate transport and an approximate collision operator. Such calculations are relevant to both problems in astrophysics, such as heat conduction and magnetic reconnection in accretion disks around black holes, and also to programmatic interests at LANL. The bolt code paper, demonstrating accuracy via a suite of test problems calculated on kodiak, is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport coefficients of quasiparticle models within a new relaxation time approximation of the Boltzmann equation

We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak

The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψ N ) or large (>4% ψ N ) island width depending on the magnitude of the edge bootstrap current, where ψ N is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χ | /χ ⊥ ) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗

Kinetic theory for electrostatic waves due to transverse velocity shears

A kinetic theory in the form of an integral equation is provided to study the electrostatic oscillations in a collisionless plasma immersed in a uniform magnetic field and a nonuniform transverse electric field. In the low temperature limit the dispersion differential equation is recovered for the transverse Kelvin-Helmholtz modes for arbitrary values of K parallel, where K parallel is the component of the wave vector in the direction of the external magnetic field assumed in the z direction. For higher temperatures the ion-cyclotron-like modes described earlier in the literature by Ganguli, Lee and Plamadesso are recovered. In this article, the integral equation is reduced to a second-order differential equation and a study is made of the kinetic Kelvin-Helmholtz and ion-cyclotron-like modes that constitute the two branches of oscillation in a magnetized plasma including a transverse inhomogeneous dc electric field.

Ganguli, G.↗

Linear gyrokinetic stability of a high β non-inductive spherical tokamak

Spherical tokamaks (STs) have been shown to possess properties desirable for a fusion power plant such as achieving high plasma β and having increased vertical stability. To understand the confinement properties that might be expected in the conceptual design for a high β ST fusion reactor, a 1 GW ST plasma equilibrium was analyzed using local linear gyrokinetics to determine the type of micro-instabilities that arise. Kinetic ballooning modes and micro-tearing modes are found to be the dominant instabilities. The parametric dependence of these linear modes was determined and, from the insights gained, the equilibrium was tuned to find a regime marginally stable to all micro-instabilities at θ 0 = 0.0. This work identifies the most important micro-instabilities expected to generate turbulent transport in high β STs. The impact of such modes must be faithfully captured in first-principles-based reduced models of anomalous transport that are needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetogasdynamic Power Extraction and Flow Conditioning for a Gas Turbine

An extension of the Russian AJAX concept to a turbojet is being explored. This magnetohydrodynamic (MHD) energy bypass engine cycle incorporating conventional gas turbine technology has MHD flow conditioning at the inlet to electromagnetically extract part of the inlet air kinetic energy. The electrical power generated can be used for various on-board vehicle requirements including plasma flow control around the vehicle or it may be used for augmenting the expanding flow in the high speed nozzle by MHD forces to generate more thrust. In order to achieve this interaction, the air needs to be ionized by an external means even up to fairly high flight speeds, and the leading candidates may be classified as electrical discharge devices. The present kinetic modeling calculations suggest that the use of electron beams with characteristics close to the commercially available e-beam systems (electron energy approx. 60 keV, beam current approx. 0.2 mA/sq cm) to sustain ionization in intermediate pressure, low-temperature (P = 0.1 atm, T = 300 K) supersonic air flows allows considerable reduction of the flow kinetic energy (up to 10 to 20 percent in M = 3 flows). The calculations also suggest that this can be achieved at a reasonable electron beam efficiency (eta approx. 5), even if the e-beam window losses are taken into account. At these conditions, the exit NO and O atom concentrations due to e-beam initiated chemical reactions do not exceed 30 ppm. Increasing the beam current up to approx. 2 mA/sq cm, which corresponds to a maximum electrical conductivity of sigma(sub max) approx. 0.8 mho/m at the loading parameter of K = 0.5, would result in a much greater reduction of the flow kinetic energy (up to 30 to 40 percent). The MHD channel efficiency at these conditions would be greatly reduced (to eta approx. 1) due to increased electron recombination losses in the channel. At these conditions, partial energy conversion from kinetic energy to heat would result in a significant total pressure loss (P(sub 0)/P(sub 0i) approx. 0.3). The total pressure loss can be reduced operating at the loading parameter closer to unity, at the expense of the reduced electrical power output. Raising the beam current would also result in the increase of the exit O atom concentrations (up to 600 ppm) and NO (up to 150 ppm).

Adamovich, Igor V.↗

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗

Solar Wind Strahl Broadening by Self-Generated Plasma Waves

This Letter reports on the results of numerical simulations which may provide a possible explanation for the strahl broadening during quiet solar conditions. The relevant processes involved in the broadening are due to kinetic quasi-linear wave-particle interaction. Making use of static analytical electron distribution in an inhomogeneous field, it is found that self-generated electrostatic waves at the plasma frequency, i.e., Langmuir waves, are capable of scattering the strahl component, resulting in energy and pitch-angle diffusion that broadens its velocity distribution significantly. The present theoretical results provide an alternative or complementary explanation to the usual whistler diffusion scenario, suggesting that self-induced electrostatic waves at the plasma frequency might play a key role in broadening the solar wind strahl during quiet solar conditions.

strahl↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Low temperature hydrogen plasma permeation in palladium and its alloys for fuel recycling in fusion systems

Superpermeation of hydrogen isotopes through metal foils is a critical component for efficient fuel recycling in fusion power systems. In that context hydrogen permeation through foils of palladium and its alloys with silver and copper was studied at low temperature (60 - 200°C) under plasma exposure. These alloys differ significantly in both bulk and surface properties, and comparisons can provide mechanistic insights. Here, permeation was observed only during plasma operation, confirming the negligible contribution of molecular hydrogen to the observed flux. As-received foils required surface treatment to achieve top performance. For Pd and Pd75Ag25 an oxidation treatment increased permeation an order of magnitude, but proved unstable as this desirable surface was reduced under hydrogen plasma exposure. In contrast, an Ar plasma cleaning step provided both high and stable flux. As-received FCC phase Pd60Cu40 foils required annealing to transform it into the high permeability BCC phase that delivered top performance. All foils displayed similar temperature dependence with flux declining with temperature, suggesting that the primary rate-limiting step is absorption of superthermal hydrogen. Among these foils the hydrogen flux through PdCu was 3 - 5X greater than that of Pd or PdAg, which were similar. The superiority of BCC PdCu is attributed to its superior hydrogen desorption kinetics. Using PdCu 100% permeation of supplied hydrogen was achieved, and the flux saturated with increasing plasma power at values >10 –2 mol H•m –2 •s– 1 , and under these conditions permeation rates are equivalent with or without the membrane present. The fluxes achieved are the highest reported to date at these conditions, and the results highlight the important roles of both surface and bulk properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-Beam Energy Transfer Saturation by Ion Heating

Cross-beam energy transfer (CBET) saturation by ion heating was measured in a gas-jet plasma characterized using Thomson scattering. Here, a wavelength-tunable ultraviolet (UV) probe laser beam was interacted with four intense UV pump beams to drive large-amplitude ion-acoustic waves. For the highest-intensity interactions, the power transferred to the probe laser dropped, demonstrating ion-acoustic wave saturation. Over this time, the ion temperature was measured to increase by a factor of 7 during the 500-ps interaction. Particle-in-cell simulations show ion trapping and a subsequent ion heating consistent with measurements. Linear kinetic CBET models were found to agree well with the observed energy transfer when the measured plasma conditions were used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1MA compared to ~8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.4 and with kinetic energies ~0.2–15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 m 2 s -1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field's effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗

Time-resolved CO 2 , CO, and N 2 vibrational population measurements in Ns pulse discharge plasmas

Abstract Time-resolved CO 2 and N 2 vibrational populations and translational-rotational temperature are measured in a CO 2 –N 2 plasma sustained by a ns pulse discharge burst in plane-to-plane geometry. Time-resolved, absolute number density of CO generated in the plasma is also inferred from the experimental data. CO 2 and CO vibrational populations are measured by mid-IR, tunable quantum cascade laser absorption spectroscopy, and N 2 vibrational populations are measured by the ns broadband vibrational CARS. Transient excitation of N 2 and CO 2 asymmetric stretch vibrational energy modes is detected during the discharge burst. The time-resolved rate of CO generation does not correlate with N 2 or CO 2 ( ν 3 ) vibrational temperatures, indicating that CO 2 dissociation via the vibrational excitation is insignificant at the present conditions. The rate of CO generation decreases gradually during the discharge burst. The estimated specific energy cost of the CO product is close to that of N atoms in pure nitrogen, measured previously at similar operating conditions. Comparison of the experimental data with the kinetic modeling analysis indicates that CO 2 dissociation in collisions with electronically excited N 2 molecules is the dominant channel of CO generation at the present conditions, although the inferred CO yield in these processes is significantly lower than 1. The effect of vibrational energy transfer between N 2 and CO 2 on the plasma chemical processes is insignificant. The kinetic model underpredicts a rapid reduction of the N 2 and CO 2 ( ν 3 ) vibrational temperatures during the later half of the discharge burst and in the afterglow. V–T relaxation of N 2 by N and O atoms generated in the ns pulse discharge plasma does not affect the vibrational relaxation rate in a significant way. However, rapid V–T relaxation of CO 2 by O atoms has a significant effect on the relaxation rate. The difference between the experimental data and the modeling predictions may be due to the unknown scaling of the CO 2 –O V–T rates with the vibrational quantum number.

Physics↗

A physics-informed deep learning description of Knudsen layer reactivity reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗