Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DRIFT”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A statistical study of gyro-averaging effects in a reduced model of drift-wave transport

Here, a statistical study of finite Larmor radius (FLR) effects on transport driven by electrostatic driftwaves is presented. The study is based on a reduced discrete Hamiltonian dynamical system known as the gyro-averaged standard map (GSM). In this system, FLR effects are incorporated through the gyro-averaging of a simplified weak-turbulence model of electrostatic fluctuations. Formally, the GSM is a modified version of the standard map in which the perturbation amplitude, K 0 , becomes K 0 J 0 ($\hat{p}$), where J 0 is the zeroth-order Bessel function and $\hat{p}$ s the Larmor radius. Assuming a Maxwellian probability density function (pdf) for $\hat{p}$ , we compute analytically and numerically the pdf and the cumulative distribution function of the effective drift-wave perturba- tion amplitude K 0 J 0 ($\hat{p}$). Using these results, we compute the probability of loss of confinement (i.e., global chaos), P c provides an upper bound for the escape rate, and that P t rovides a good estimate of the particle trapping rate. Lastly. the analytical results are compared with direct numerical Monte-Carlo simulations of particle transport.

42 ENGINEERING↗
DOE OSTI2016

Stable discrete representation of relativistically drifting plasmas

Representing the electrodynamics of relativistically drifting particle ensembles in discrete, co-propagating Galilean coordinates enables the derivation of a Particle-In-Cell algorithm that is intrinsically free of the numerical Cherenkov instability for plasmas flowing at a uniform velocity. Application of the method is shown by modeling plasma accelerators in a Lorentz-transformed optimal frame of reference.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Modelling enhanced confinement in drift-wave turbulence

Here, the results of modeling studies of an enhanced confinement in the drift wave turbulent plasma of the CSDX linear device are presented. The mechanism of enhanced confinement is investigated here using a reduced 1D, time-dependent model, which illustrates the exchange of enstrophy between two disparate scale structures: the mesoscale flow and profile, and the turbulence intensity fields. Mean density, mean vorticity, and turbulent potential enstrophy are the variables for this model. Total potential enstrophy is conserved in this model. Vorticity mixing occurs on a scale length related to an effective Rhines' scale of turbulence, and shrinks as both density and vorticity gradients steepen. Numerical results obtained from solution of the model agree well with the experimental data from CSDX showing: (i) a steepening of the mean density profile, indicating a radial transport barrier formation, (ii) the development of a radially sheared azimuthal flow velocity that coincides with the density steepening and initiates a turbulence quench, and (iii) negative Reynolds work values, indicating that fluctuations drive the shear flow. These observations as the magnitude of the magnetic field B increases are recovered using purely diffusive expressions for the vorticity and density fluxes. A new dimensionless turbulence parameter R DT -defined as the ratio of the integrated potential enstrophy transfer from turbulence to the flow, to the integrated potential enstrophy production due to relaxation of the density gradient is introduced as a turbulence collapse indicator that detects when the enhanced confinement state is triggered

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Transport dynamics of self-consistent, near-marginal drift-wave turbulence. I. Investigation of the ability of external flows to tune the non-diffusive dynamics

The reduction of turbulent transport across sheared flow regions has been known for a long time in magnetically confined toroidal plasmas. However, details of the dynamics are still unclear, in particular, in what refers to the changes caused by the flow on the nature of radial transport itself. In Paper II, we have shown in a simplified model of drift wave turbulence that, when the background profile is allowed to evolve self-consistently with fluctuations, a variety of transport regimes ranging from superdiffusive to subdiffusive open up depending on the properties of the underlying turbulence. In this study, we show that externally applied sheared flows can, under the proper conditions, cause the transport dynamics to be diffusive or subdiffusive

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas

A conservative scheme of drift kinetic electrons for gyrokinetic simulations of kinetic-magnetohydrodynamic processes in toroidal plasmas has been formulated and verified. Both vector potential and electron perturbed distribution function are decomposed into adiabatic part with analytic solution and non-adiabatic part solved numerically. The adiabatic parallel electric field is solved directly from the electron adiabatic response, resulting in a high degree of accuracy. Here, the consistency between electrostatic potential and parallel vector potential is enforced by using the electron continuity equation. Since particles are only used to calculate the non-adiabatic response, which is used to calculate the non-adiabatic vector potential through Ohm's law, the conservative scheme minimizes the electron particle noise and mitigates the cancellation problem. Linear dispersion relations of the kinetic Alfvén wave and the collisionless tearing mode in cylindrical geometry have been verified in gyrokinetic toroidal code simulations, which show that the perpendicular grid size can be larger than the electron collisionless skin depth when the mode wavelength is longer than the electron skin depth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Gyrokinetic theory of slab universal modes and the non-existence of the gradient drift coupling (GDC) instability

A gyrokinetic linear stability analysis of a collisionless slab geometry in the local approximation is introduced. We focus on k ∥ =0 universal (or entropy) modes driven by plasma gradients at small and large plasma β. These are small scale non-MHD instabilities with growth rates that typically peak near k ⊥ ρi~1 and vanish in the long wavelength k ⊥ →0 limit. This report also discusses a mode known as the Gradient Drift Coupling (GDC) instability previously reported in the gyrokinetic literature, which has a finite growth rate γ= $\sqrt{β/[2(1+β)]C_{s}/|L_{p}|}$ with C$^{2}_{s}$ =p 0 /ρ 0 for k ⊥ →0 and is universally unstable for 1/L p ≠0. Here, we show that the GDC instability is a spurious, unphysical artifact that erroneously arises due to the failure to respect the total equilibrium pressure balance p 0 +B$^{2}_{0}$ /(8π)=constant, which renders the assumption B$^{\prime}_{0}$ 0 =0 inconsistent if p$^{\prime}_{0}$ ≠0.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Design and implementation of a Thomson parabola for fluence dependent energy-loss measurements at the Neutralized Drift Compression eXperiment

The interaction of ion beams with matter includes the investigation of the basic principles of ion stopping in heated materials. An unsolved question is the effect of different, especially higher, ion beam fluences on ion stopping in solid targets. This is relevant in applications such as in fusion sciences. To address this question, a Thomson parabola was built for the Neutralized Drift Compression eXperiment (NDCX-II) for ion energy-loss measurements at different ion beam fluences. The linear induction accelerator NDCX-II delivers 2 ns short, intense ion pulses, up to several tens of nC/pulse, or 1010-1011 ions, with a peak kinetic energy of ~1.1 MeV and a minimal spot size of 2 mm FWHM. For this particular accelerator, the energy determination with conventional beam diagnostics, for example, time of flight measurements, is imprecise due to the non-trivial longitudinal phase space of the beam. In contrast, a Thomson parabola is well suited to reliably determine the beam energy distribution. The Thomson parabola differentiates charged particles by energy and charge-to-mass ratio, through deflection of charged particles by electric and magnetic fields. During first proof-of-principle experiments, we achieved to reproduce the average initial helium beam energy as predicted by computer simulations with a deviation of only 1.4%25. Successful energy-loss measurements with 1 μm thick silicon nitride foils show the suitability of the accelerator for such experiments. The initial ion energy was determined during a primary measurement without a target, while a second measurement, incorporating the target, was used to determine the transmitted energy. The energy-loss was then determined as the difference between the two energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Drift-ideal magnetohydrodynamic simulations of m = 0 modes in Z-pinch plasmas

In this paper, the effects of m = 0 modes on equilibrium Z-pinch plasmas are studied using a drift-ideal magnetohydrodynamic (MHD) model. The model equations are an extension of ideal MHD to include finite-ion-inertial-length/cyclotron-frequency (Ω i ) effects in Ohm's law and in the electron and ion heat transport equations. The linear modes contained in this model include the ideal interchange (sausage) mode and in the magnetized limit, Ω i τ i »1 with τ i the ion collision time, nonideal entropy modes. It is well known that these two modes are decoupled in the kρ s « 1 limit, where k is the axial mode number and ρ s = c s /Ω i is the gyro-Bohm scale with c s the sound speed [B. Kadomtsev, Sov. Phys. JETP-USSR 10, 780 (1960)]. For Bennett equilibrium profiles, it is shown that the regions of stability for both modes are completely governed by the adiabatic coefficient γ in these limits. Equilibria with Bennett profiles are stable to entropy modes for γ < 2 but unstable to ideal modes and vice versa for γ > 2. However, these modes are no longer decoupled when kρ s ≳1. The simulation results of the fully nonlinear set of equations in the magnetized limit show that seeded modes with kρ s ≳1 and γ = 5/3 display the characteristics of both ideal and entropy modes. The general heat flux for both ions and electrons as a function of the species magnetization is retained in the model. Both the linear and nonlinear behaviors of seeded modes for kρ s ≳1 display a strong dependence on the magnetization of the ions. Lastly, the growth rate increases linearly with k at large kρ s when the ions are magnetized but decreases with increasing k when Ω i τ i ≲1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Self-organized standing waves generated by AC-driven electron cyclotron drift instabilities

The formation of self-organized standing wave structures is observed due to the electron cyclotron drift instability (ECDI) driven by a time-varying external electric field and a crossed magnetic field. Using a particle-in-cell simulation, two standing wave mechanisms are identified: the linear mode and the beating mode. In the former, a standing wave emerges as a superposition of two counterpropagating ion acoustic waves predicted by the linear theory of ECDI. On the other hand, in the beating mode, the plasma wave is in resonance with the applied frequency. Nonlinear resonance of such standing waves results in a change in the dominant wavenumber and frequency. Such counterpropagating plasma waves are consistent with the experimental observations using coherent Thomson scattering in a crossed-field plasma discharge.

42 ENGINEERING↗
DOE OSTI2021

Drift waves enstrophy, zonal flow, and nonlinear evolution of the modulational instability

The interaction of the drift wave (DW) turbulence and zonal flow (ZF) is investigated with the modified Hasegawa–Mima equation taking into account the backreaction of ZF velocity on DW turbulence. It is shown that the y-averaged enstrophy of DW turbulence and the velocity of ZF are intrinsically related. By utilizing this feature, a nonlinear stage of DW modulational instability is considered within the framework of the wave kinetic equation. It is shown that in this approximation, the nonlinear stage of the modulational instability results in the collapsing solutions, accompanied by the “wave breaking” phenomenon. Numerical simulations based on the Hasegawa–Mima equation show that for a weak DW turbulence, Φ ̃ = ( e φ ̃ / T e ) ( L n / ρ s ) ⪝ 1 , the collapsing-like features on both ZF and y-averaged enstrophy of DW turbulence decay in time and then re-emerge again at different locations. For the case of a strong DW turbulence, Φ ̃ > 1 , where nonlinear interactions of DW harmonics dominate, stable spatial structures of ZF and y-averaged enstrophy of DW turbulence emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of magnetic angle on the E × B drift in a magnetic presheath

In this work, we provide experimental measurements of the E X B flow above a non-magnetic, absorbing boundary in a magnetized plasma. Measurements are taken as a function of the angle between the background magnetic field and the boundary normal. The measurements are compared to computational predictions for oblique magnetic fields based on Chodura's model. Ion flow measurements were obtained with laser induced fluorescence and the presheath potential structure was measured with an emissive probe. The ions were found to accelerate to nearly 30% of the sound speed parallel to the boundary at oblique angles of the magnetic field (~80°) and the ion speed at the electrostatic sheath edge was found to decrease with angle. The edge of the magnetic presheath has been experimentally determined based on the appearance of the E X B drift and is found to be independent of the angle of the magnetic field. We also demonstrate that laser induced fluorescence is suitable as a non-perturbative diagnostic to measure the electric field in a magnetized sheath.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Positive charging of grains in an afterglow plasma is enhanced by ions drifting in an electric field

In a plasma, the polarity of a dust grain's charge is typically negative, but it can reverse and become positive in an afterglow, when the power sustaining the plasma is switched off. This positive charging, which occurs in the afterglow's first few milliseconds, is studied for grains much larger than a few nm. It is hypothesized that the positive charging is enhanced by the presence of a dc electric field, which causes ions to drift through the neutral gas. A larger value of the reduced electric field E/N leads to a larger ion kinetic energy and thus a greater collection of positive charge on a grain. The maximum possible positive charge is attained if the grain's surface potential rises to match the ion kinetic energy, at a time before ions have departed and the grain's charge becomes frozen. Thereafter, when vacuum conditions prevail, the grain will retain its positive residual charge. In an experiment, dust grains were electrically levitated in a capacitively coupled plasma until the power was abruptly turned off. In the afterglow, grains fell faster than expected due to gravity alone, indicating a downward electric force, in the presence of a remaining dc electric field. Acceleration measurements yielded repeatable results for the residual charge's value, which was of the order +10 4 e and increased with E/N, supporting the hypothesis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

3D simulation of lower-hybrid drift waves in strong guide field asymmetric reconnection in laboratory experiments

The magnetic reconnection experiment has recently seen short wavelength ([Formula: see text]) lower-hybrid waves near the electron diffusion region in strong guide field reconnection. Based on plasma parameters from the experiment, we perform a three-dimensional fully kinetic simulation in order to investigate the generation of the lower-hybrid waves and their effects on the reconnection process. We find that the low-beta regions around the reconnection site are unstable to the lower-hybrid drift instability propagating in the outflow direction, driven by the difference between the electron and ion outflows. The waves modify the electron distributions, leading to periodic opening and closing of gaps in electron velocity space, and provide a small contribution to the anomalous resistivity. Finally, the simulation results are discussed in the context of space observations and laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2024

Drag exerted by drift-acoustic radiation on a magnetic island

The Island Equilibrium and Transport (IslET) model is extended to allow for the effects of ion compressibility albeit for a constant electron temperature. The model describes fully relaxed, saturated, steady-state magnetic islands arising from the evolution of a tearing mode caused either by an instability or an externally imposed resonant perturbation. It assumes that the collision frequency lies in the semi-collisional regime. A salient new feature of the extended model is an equation describing the transport of parallel momentum. The model is applied to evaluate the drag exerted on a thin island by the emission of drift-acoustic waves.

Physics↗
DOE OSTI2024

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↗