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 163 records · Page 9

A unified kinetic model of the tangential magnetopause structure

In the self-consistent model of the tangential magnetopause, formulated in the present paper on the basis of the Vlasov-Maxwell equations, the plasmas on both sides are magnetized and the magnetic field is everywhere parallel to the magnetopause (i.e., the normal field component is zero) and rotates through an arbitrary angle across the magnetopause. It is shown that the thickness of the magnetopause is greater than the gyroradius of the plasma ions. The presence of a trapped particle population within the magnetopause is shown to be required to allow the magnetic field to rotate more than a certain critical angle (-90 degrees). The model proposed can reproduce the observed features of the tangential magnetopause structure by specifying boundary conditions on both sides of the magnetopause.

Lee, L. C.↗

Double layers on auroral field lines

Time-stationary solutions to the Vlasov-Poisson equation for ion holes and double layers were examined along with particle simulations which pertain to recent observations of small amplitude (e phi)/t sub e approx. 1 electric field structures on auroral field lines. Both the time-stationary analysis and the simulations suggest that double layers evolve from holes in ion phase space when their amplitude reaches (e phi)/t sub e approx. 1. Multiple small amplitude double layers which are seen in long simulation systems and are seen to propagate past spacecraft may account for the acceleration of plasma sheet electrons to produce the discrete aurora.

Hudson, M. K.↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

Solitary waves and double layers on auroral field lines

Time stationary solutions to the Vlasov-Poisson equations for ion holes and double layers are examined along with particle simulations that pertain to recent observations of small amplitude electric field structures on auroral field lines. Both the time stationary analysis and the simulations suggest that the observed double layers evolve from holes in ion phase space. Multiple small amplitude double layers, as seen in long simulation systems, are observed to propagate past the spacecraft and may account for the acceleration of plasma sheet electrons to produce inverted-V precipitation.

Hudson, M. K.↗

Nonlinear upper hybrid drift waves for a longitudinal electric field perpendicular to a uniform magnetic field in the Vlasov-Maxwell approximation

Upper hybrid drift waves are found as a special solution to a Vlasov-Maxwell plasma which has a longitudinal electric field and a perpendicular uniform magnetic field. A single-species plasma with a constant-density mobile neutralizing background supports spatially varying disturbances that oscillate at the upper hybrid frequency. The general functional dependences of the electric field, the plasma number density, and the one-particle distribution function for the special case are found from more general Vlasov-Maxwell equations invariant under a Lie group point transformation. The one-particle distribution function for the plasma is a function of the Liouville invariant, which is the energy in the generalized Bernstein-Greene-Kruskal (BGK) reference frame, and the momentum in the drift direction.

Abraham-Shrauner, B.↗

Theory of beat-resonant coupling of electrostatic modes

A general expression is derived for the beat-resonant coupling electrostatic modes in a Vlasov plasma. The result for the coupling of two modes has a simple structure: the appropriate momentum gradient of the equilibrium particle distribution is weighted by a positive coupling coefficient and averaged over the resonance surface in momentum space. The contributions of all the resonance surfaces are then summed. This basic structure had been previously exhibited only for specific homogeneous plasma models. The present theory, which unifies and greatly simplifies these individual treatments, is based on a variational formulation of the Vlasov-Poisson equations. Using Lie transforms, the variational principle is reexpressed in oscillation-center variables, and then the nonlinear wave dynamics are obtained from the independent variations of the wave phase and the wave amplitude. The power of the method is then applied to a strongly magnetized, strongly inhomogeneous, non-neutral plasma model.

Crawford, John David↗

Linearly polarized magnetic fluctuations at comet Giacobini-Zinner

The fully electromagnetic linear Vlasov dispersion equation for instabilities in a homogeneous, magnetized plasma is solved. Parameters appropriate to conditions observed during the ICE encounter with comet Giacobini-Zinner are used to study the polarization of electromagnetic instabilities driven by the relative streaming of two ion components. If the ion/ion relative drift is near the Alfven speed, the associated long-wavelength modes may be approximately linearly polarized at relatively small angles of propagation with respect to the magnetic field. The relatively inhomogeneous plasma background of the comet may impose such oblique propagation, so that this result provides an alternate explanation to the observations of linearly polarized waves near the comet.

Gary, S. Peter↗

Electron/ion whistler instabilities and magnetic noise bursts

Two whistler instabilities are investigated by means of the linear Vlasov dispersion equation. They are called the electron/ion parallel and oblique whistler instabilities, and are driven by electron/ion relative drifts along the magnetic field. It is demonstrated that the enhanced fluctuations from these instabilities can explain several properties of magnetic noise bursts in and near the plasma sheet in the presence of ion beams and/or field-aligned currents. At sufficiently high plasma beta, these instabilities may affect the current system in the magnetotail.

Akimoto, K.↗

Ion anisotropy instabilities in the magnetosheath

Recent observations in Earth's magnetosheath have delineated several different kinds of magnetic fluctuation spectra below the proton cyclotron frequency. This paper provides a theoretical interpretation for some of these observations describing solutions of the linear Vlasov dispersion equation for fully electromagnetic instabilities for particle distributions which model those observed in the magnetosheath. This model yields three growing modes: the proton cyclotron anisotropy the helium cyclotron anisotropy, and the mirror instabilities. The results show very good agreement with the observations of mirror-like and proton-cyclotron-like events. This agreement with observations implies that the transition between cyclotron and mirror fluctuation dominance is consistent with linear theory.

Gary, S. P.↗

Superposition of nonlinear plasma waves

We report results showing that spatially periodic Bernstein-Greene-Kruskal (BGK) waves, which are exact nonlinear traveling wave solutions of the Vlasov-Maxwell equations for collisionless plasmas, satisfy a nonlinear principle of superposition in the small amplitude limit. The analysis explicates the notion of superimposed BGK waves which, as recent numerical calculations suggest, is crucial in the proper description of the time-asymptotic state of a plasma when a large amplitude electrostatic wave undergoes nonlinear Landau damping.

Buchanan, Mark↗

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗

Mirror Instability in the Solar Wind: The Theory Revisited

"Magnetic holes", localized depressions in the interplanetary magnetic field, have been identified in Ulysses data over a range of several AU and as far as 23 degrees south in latitude by Winterhalter et al., who concluded that these structures are most likely the remnants of structures caused by occasional mirror-mode instability in the solar wind. However, these authors, like a number of previous investigators, used the mirror stability criterion derived from the kinetic theory under very special assumptions. On the other hand, theoretical investigations using the fully self-consistent kinetic theory (Vlasov-Maxwell equations) have shown that the mirror stability criterion is more complicated when electrons and ions have different anisotropies, as is normally the case in the solar wind. Winterhalter et al used an instability criterion of the form R is greater than 1, where R is a function of the thermal anisotropy; the correct criterion (for bi-Maxwellian distributions) is R R is greater than 1 - x(exp 2), where x is a real quantity that depends on both the proton anisotropy and electron anisotropy. So nonzero x would modify the Winterhalter et al results in the direction of reinforcing their conclusions. We have revisited the instability criterion in its most general form, allowing for (a) non-Maxwellian velocity distributions, (b) multiple ion species, and (c) interparticle streaming. These results should give sound theoretical grounding for future observational studies related to the mirror instability, by Ulysses and other spacecraft.

Barnes, A.↗

Higher-order space-charge stability in anisotropic beams: Vlasov-Poisson derivation, refined dispersion relations, and stability charts

The Hofmann stability chart is used to screen working points in space-charge-dominated linacs. We identify two errors in its published higher-order dispersion relations: missing $(1\mp2\hatη^2/α)$ factors in the third-order $S^4$ coupling residues, and a sign error in the stated isotropic reduction of the fourth-order relation. Both corrections follow from Hofmann's Vlasov-Poisson equations without fitted parameters. They reproduce coherent tune-shift coefficients in the author's later monograph that the printed forms miss by 24% and 127%. Mode-resolved figures from a published application agree with the corrected relations and reject the printed forms, indicating an inconsistency between the 1998 equations and the calculations underlying those tested figures. We quantify the effect on the non-oscillatory stability chart. Inside the adopted $S^2\le10$ comparison domain, printed and corrected forms disagree on 0.73-2.11% of cells, with no preferred direction. Among excluded cells, disagreement reaches 22%, and the printed relation over-predicts instability at every sampled anisotropy. This concentration may help explain why the errors persisted, although it does not establish their historical cause. For PIP-II, the corrected chart flags four of thirty-two evaluable periods, including one on a third-order odd branch missed by a second-order screen. This count covers non-oscillatory modes only and remains conditional on an unresolved factor-five disagreement between two codes on transverse emittance growth.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system

We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is an accurate model for studying hot plasma dynamics at a kinetic scale where small-angle Coulomb collisions dominate. Our scheme guarantees conservation of mass, momentum and energy as well as preservation of the monotonicity of entropy production in both the time-continuous and discrete systems. We employ the conservative integrator for both the Hamiltonian Vlasov-Poisson equations and the dissipative Landau equation using the PETSc library (www.mcs.anl.gov/petsc) to showcase structure-preserving properties.

Discrete gradients↗

Wave-particle transport from density drift instabilities - A comparison of local and nonlocal theories

Second-order Vlasov theory is used to compute the dissipation rates of plasma irregularities with a variety of shapes. A derivation of the nonlocal dispersion equation using linearized Vlasov theory is presented. Expressions for the normalized amplitudes of the first-order plasma density and electrostatic potential fluctuations are derived. Expressions are given for the saturation amplitudes of the electrostatic eigenmodes. The wave-particle transport and irregularity dissipation rate are computed by using formulas whose derivation is presented. Computational results for specific density variations are shown, and conclusions on the validity of the local theory as opposed to the nonlocal theory are given.

Bernhardt, P. A.↗

Implications of Navier-Stokes turbulence theory for plasma turbulence

The methodology of Navier-Stokes fluid turbulence theory is reviewed, with emphasis placed on the relevance of the Navier-Stokes concepts for understanding plasma turbulence. After a brief consideration of the three-dimensional case, two-dimensional problems are discussed. In addition, MHD turbulence and turbulence in Vlasov plasmas are treated. In particular, the direct interaction approximation developed by Kraichnan (1959) is generalized from Navier-Stokes turbulence theory to produce a computable set of differentio-integral equations for a Vlasov plasma.

Montgomery, D.↗

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↗

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↗