Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vlasov equation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Spacecraft-environment interaction: The environmental plasma aspect

The effects involved in the interaction between an obstacle and a space plasma can be divided into: (1) effects on the obstacle itselt (i.e., its charging); and (2) effects on the environmental plasma due to the motion of the obstacle (i.e., the creation of shocks ahead of the obstacle and complicated wakes behind the obstacle). In the wake (or antisolar direction), plasma oscillations are excited and instabilities, wave-particle interactions, turbulence, etc., are believed to take place. The effects on the obstacle and on the environmental space plasma are coupled. Hence, simultaneous solutions to the Vlasov (or Boltzmann) and Poisson equations are sought. To obtain realistic solutions of practical use, three-dimensional and time-dependent models of the interaction are needed. Achieving the latter is indeed not simple.

Samir, U.↗
DOE OSTI2022

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Required toroidal confinement for fusion and omnigeneity

Deuterium–tritium (DT) burning requires a long energy confinement times compared to collision times, so the particle distribution functions must approximate local Maxwellians. Non-equilibrium thermodynamics is applicable, which gives relations among transport, entropy production, the collision frequency, and the deviation from a Maxwellian. The distribution functions are given by the Fokker–Planck equation, which is an advection–diffusion equation. A large hyperbolic operator, the Vlasov operator with the particle trajectories as its characteristics, equals a small diffusive operator, the collision operator. The collisionless particle trajectories would be chaotic in stellarators without careful optimization. This would lead to rapid entropy production and transport—far beyond what is consistent with a self-sustaining DT burn. Omnigeneity is the weakest general condition that is consistent with a sufficiently small entropy production associated with the thermal particle trajectories. Omnigeneity requires that the contours of constant magnetic field strength be unbounded in at least one of the two angular coordinates in magnetic surfaces and that there be a symmetry in the field-strength wells along the field lines. Even in omnigenous plasmas, fluctuations due to microturbulence can produce chaotic particle trajectories and the gyro-Bohm transport is seen in many stellarator and tokamak experiments. The higher the plasma temperature above 10 keV, the smaller the transport must be compared to gyro-Bohm for a self-sustaining DT burn. The hot alphas of DT fusion heat the electrons. When the ion–electron equilibration time is long compared to the ion energy confinement time, a self-sustaining DT burn is not possible, which sets a limit on the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Conservation of quasiparticles in weakly turbulent plasmas.

An expansion making use of the eikonal is shown to yield a solution to the equations of a general weak turbulent plasma which is weakly dependent on space and time. The method is used to derive a quasi-particle conservation equation for quasi-static perturbations of Vlasov plasmas with general equilibrium field configurations.

Bloomberg, H. W.↗
NASA NTRS1977

Calculation of sheath and wake structure about a pillbox-shaped spacecraft in a flowing plasma

A computer program was used for studies of the disturbed zones around bodies in flowing plasmas, particularly spacecraft and their associated sheaths and wakes. The program solved a coupled Poisson-Vlasov system of nonlinear partial differential integral equations to obtain distributions of electric potential and ion and electron density about a finite length cylinder in a plasma flow at arbitrary ion Mach numbers. The approach was applicable to a larger range of parameters than other available approaches. In sample calculations, bodies up to 100 Debye lengths in radius were treated, that is, larger than any previously treated realistically. Applications were made to in-situ satellite experiments.

Parker, L. W.↗
NASA NTRS1985

Thermal deformation of concentrators in an axisymmetric temperature field

Axisymmetric thermal deformations of paraboloid mirrors, due to heating, are examined for a mirror with a optical axis oriented toward the Sun. A governing differential equation is derived using Mushtari-Donnel-Vlasov simplifications, and a solution is presented which makes it possible to determine the principal deformation characteristics.

Bairamov, R.↗
DOE OSTI2024

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Nondimensional Parameters and Equations for Nonlinear and Bifurcation Analyses of Thin Anisotropic Quasi-Shallow Shells

A comprehensive development of nondimensional parameters and equations for nonlinear and bifurcations analyses of quasi-shallow shells, based on the Donnell-Mushtari-Vlasov theory for thin anisotropic shells, is presented. A complete set of field equations for geometrically imperfect shells is presented in terms general of lines-of-curvature coordinates. A systematic nondimensionalization of these equations is developed, several new nondimensional parameters are defined, and a comprehensive stress-function formulation is presented that includes variational principles for equilibrium and compatibility. Bifurcation analysis is applied to the nondimensional nonlinear field equations and a comprehensive set of bifurcation equations are presented. An extensive collection of tables and figures are presented that show the effects of lamina material properties and stacking sequence on the nondimensional parameters.

Nemeth, Michael P.↗
DOE OSTI2022

LOKI

LOKI is a simulation code that models weakly-coupled plasma evolution using Eulerian kinetic models in four phase space dimensions. The code is useful to understand basic plasma nonlinear plasma physics phenomena - in particular, wave-particle and wave-wave resonances that occur in laser plasma interactions. LOKI can simulate fundamental physical process in low-density plasmas on a scale of hundreds of light wavelengths. It does so by solving the relativistic Vlasov and Valsov-Fokker-Planck differential transport equations describing the collision less and low-collisionality evolution of charged particles in electromagnetic fields in two configuration space and two velocity space dimensions. Two principle electromagnetic models are provided: electrostatic (Poisson's Equation) and electromagnetic (Maxwell's equations).

DIMITS, ANDRIS↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

Electromagnetic ion beam instabilities

The linear theory of electromagnetic instabilities driven by an energetic ion beam streaming parallel to a magnetic field in a homogeneous Vlasov plasma is considered. Numerical solutions of the full dispersion equation are presented. At propagation parallel to the magnetic field, there are four distinct instabilities. A sufficiently energetic beam gives rise to two unstable modes with right-hand polarization, one resonant with the beam, the other nonresonant. A beam with sufficiently large T (perpendicular to B)/T (parallel to B) gives rise to the left-hand ion cyclotron anisotropy instability at relatively small beam velocities, and a sufficiently hot beam drives unstable a left-hand beam resonant mode. The parametric dependences of the growth rates for the three high beam velocity instabilities are presented here. In addition, some properties at oblique propagation are examined. It is demonstrated that, as the beam drift velocity is increased, relative maxima in growth rates can arise at harmonics of the ion cyclotron resonance for both right and left elliptically polarized modes.

Gary, S. P.↗
DOE OSTI2021

Hamiltonian kinetic-Hall magnetohydrodynamics with fluid and kinetic ions in the current and pressure coupling schemes

We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Low-frequency waves in a high-beta collisionless plasma Polarization, compressibility and helicity

This paper considers the linear theory of waves near and below the ion cyclotron frequency in an isothermal electron-ion Vlasov plasma which is isotropic, homogeneous and magnetized. Numerical solutions of the full dispersion equation for the magnetosonic/whistler and Alfven/ion cyclotron modes at beta(i) = 1.0 are presented, and the polarizations, compressibilities, helicities, ion Alfven ratios and ion cross-helicities are exhibited and compared. At sufficiently large beta(i) and theta, the angle of propagation with respect to the magnetic field, the real part of the polarization of the Alfven/ion cyclotron wave changes sign, so that, for such parameters, this mode is no longer left-hand polarized. The Alfven/ion cyclotron mode becomes more compressive as the wavenumber increases, whereas the magnetosonic/whistler becomes more compressive with increasing theta.

Gary, S. P.↗
NASA NTRS1983

The nonlocal theory of periodic density drift instabilities

The nonlocal theory of electrostatic density drift instabilities is developed for an arbitrary, periodic variation in plasma density. Linear Vlasov theory is applied to Fourier series representations of the plasma density, and the resulting equations are solved numerically. The growth rates of the universal drift instability are compared by using local and nonlocal theories. It is demonstrated that the local theory is inadequate if the nonlocal eigenfunctions of the electrostatic modes have significant amplitude over spatial regions wider than the regions of large plasma density gradients. Nonlocal effects can enhance wave-particle diffusion of a triangular-shaped variation and reduce wave-particle dissipation of a square-shaped irregularity.

Bernhardt, P. A.↗