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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Dynamical features of moving double layers

Numerical simulations of the dynamics of double layers that form in response to an applied potential drop across a plasma are carried out for plasmas of lengths much greater than previously considered. The evolution of finite-size plasmas of lengths 100, 200 and 400 Debye lengths at the low-potential end of the simulation region was followed by the solution of the Vlasov and Poisson equations. Results show that the double layer moves towards the high-potential side of the layer, at a speed directly dependent on the length of the plasma. The moving double layer is accompanied by a moving density front which is similar to an ion acoustic shock created by the expansion of a high-density plasma into a low-density plasma, as well as plasma heating and evacuation, electron and ion current interruptions and recovery. Double layer formation, motion, and accompanying phenomena are found to repeat periodically. Results are in good agreement with laboratory experiments, and may be used to explain certain phenomena observed in auroral plasmas.

Singh, N.↗

NASA NTRS1983

Linear density drift instabilities in very low beta plasmas A different approach

This paper reports a study of the linear Vlasov electromagnetic dispersion equation for density drift instabilities in very low beta plasmas. A uniform magnetic field and a weak uniform density gradient are assumed. This paper differs from most other studies of this topic in three important ways: (1) no low frequency or long wavelength approximations are made, (2) no gauge condition is imposed and (3) a modification of the local approximation is used which is argued to be less arbitrary than the usual local approximation. Numerical solution of the resulting dispersion equation yields stabilization of the universal density drift instability at values of beta somewhat higher than those obtained from the local approximation, and near maximum growth shows no evidence of coupling between the universal and Alfven modes.

Gary, S. P.↗

NASA NTRS1990

High voltage plasma sheath analysis related to TSS-1

On the first mission of the Tethered Satellite System (TSS-1), a 1.8 m diameter spherical satellite will be deployed a distance of 20 km above the Space Shuttle Orbiter on an insulated conducting tether. The satellite will be held at electric potentials up to 5000 volts positive with respect to the ambient plasma. Due to the passage of the conducting tether through the Earth's magnetic field, an electromagnetic field (EMF) will be created, driving electrons down the tether to the Orbiter, out through an electron gun into the ionosphere and back into the positive-biased satellite. The main problem addressed here is the current-voltage characteristics of the ionospheric interaction with the satellite. The first problem is that while the satellite will be capable of measuring charged particle flow to the surface at several locations, the detectors have a limited range of acceptance angle. The second problem is that the angle of incidence of the incoming electrons will have to be relative to the local normal. This will be important in order to predict the magnitude of the detectable current at each detector location so the detector gain can be pre-set to the correct range. The plasma sheath was analyzed mathematically, and subroutines were written to solve relevant finite element, Taylor-Vlasov, and Poisson equations.

Sheldon, John W.↗

NASA NTRS1996

Formation of the Electric Double Layer and its Effects on Moving Bodies in a Space Plasma Environment

In this paper we solve the self-consistent Vlasov and Poisson equations by a numerical method to determine the local distribution function of the ion and the electron, within a thin layer near the moving body, respectively. Using these ion and electron distributions, the number density for the ions and electrons are determined, such that, the electric potential is obtained within this thin layer (i.e., measured by Debye length). Numerical results are presented for temporal evolution of the electron and ion density and its corresponding electric potential within the layer which shows the formation of electric double layer and its structures. From these numerical results, we are able to determine the maximum conditions of the electric potential, it may create satellite anomaly.

Yang, Qianli↗

NASA NTRS2009

Analysis of Electromagnetic Wave Propagation in a Magnetized Re-Entry Plasma Sheath Via the Kinetic Equation

Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically induced plasma, it has been demonstrated that the classical radiofrequency communications blackout that is experienced during atmospheric reentry can be mitigated through the appropriate control of an external magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and involves an analytical solution for the electric and magnetic fields within the plasma allowing for a description of the attendant transmission, reflection and absorption coefficients. The ability to transmit through the magnetized plasma is due to the magnetic windows that are created within the plasma via the well-known whistler modes of propagation. The case of 2 GHz transmission through a re-entry plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves through the re-entry phase.

Manning, Robert M.↗

DOE OSTI2025

Algebra of invariants for the Vlasov–Maxwell system

The algebra of invariants for both the relativistic and nonrelativistic multispecies Vlasov–Maxwell system is examined, including the case with a fixed ion background. Invariants and their associated fluxes are obtained directly from the Vlasov–Maxwell system. The invariants are shown to Poisson commute with the Hamiltonian and the rest of the Poisson bracket algebra of invariants is identified. Special attention is given to the role played by the monopole condition, ∇ · B.

Fundamental invariants↗

DOE OSTI2024

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗

DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Cosmological perturbation theory for large scale structure in phase space

We develop a framework for Large Scale Structure (LSS) perturbation theory, that solves the Vlasov-Poisson system of equations for the distribution function in full phase space. This approach relaxes the usual apriori assumption of negligible velocity dispersion underlying the Standard Perturbation Theory (SPT). We apply the new method to rederive the usual SPT kernels up to third order in the perturbative expansion. We also show that a counterterm, identical to the one introduced by standard Effective Field Theory (EFT) methods, naturally arises within our framework. We finish by making a precise connection to EFT techniques, which reveals the necessity of the EFTofLSS to self-consistently model the long-wavelength fluid, and illustrates the importance of having theoretical control over short distance fluctuations.

Cosmological perturbation theory in GR and beyond↗

DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗

DOE OSTI2020

Theoretical formulation of phase space microbunching instability in the presence of intrabeam scattering for single-pass or recirculation accelerators

Microbunching instability (MBI) has been one of the most challenging issues in designs of high-brightness beam transport lines for single-pass or recirculating accelerators. Although the intrabeam scattering (IBS) has long been studied in lepton or hadron storage rings as a slow diffusion process or in high-intensity proton linear accelerators as one mechanism for the beam halo, the effects of IBS on single-pass or recirculating electron accelerators have drawn attention only in the recent two decades due to emergence of linac-based or energy-recovery-linac-based fourth-generation light sources, which require high-quality electron beams during the beam transport. In this paper we develop a theoretical formulation of microbunching instability in the presence of IBS for single-pass or recirculation accelerators. To quantify MBI with inclusion of IBS, we start from the Vlasov-Fokker-Planck (VFP) equation, combining both collective interactions and incoherent IBS effects. The linearized VFP equation and the corresponding friction and diffusion coefficients are derived. The evolutions of the resultant density and energy modulations are formulated as a set of coupled integral equations. The theoretical formulation is then applied to a recirculating beamline design. The results from the semianalytical calculation are compared and show good agreement with massive particle tracking simulations.

43 PARTICLE ACCELERATORS↗

NASA NTRS1982

A mechanism for plasma waves at the harmonics of the plasma frequency foreshock boundary

A bump-on-tail unstable reduced velocity distribution, constructed from data obtained at the upstream boundary of the electron foreshock by the GSFC electron spectrometer experiment on the ISEE-1 satellite, is used as the initial plasma state for a numerical integration of the 1D-Vlasov-Maxwell system of equations. The integration is carried through the growth of the instability, beyond its saturation, and well into the stabilized plasma regime. A power spectrum computed for the electric field of the stabilized plasma is dominated by a narrow peak at the Bohm-Gross frequency of the unstable field mode but also contains significant power at the harmonics of the Bohm-Gross frequency. The harmonic power is in sharp peaks which are split into closely spaced doublets. The fundamental peak at the Bohm-Gross frequency is split into a closely spaced triplet. The mechanism for excitation of the second harmonic is shown to be second order wave-wave coupling.

Klimas, A. J.↗

NASA NTRS1983

A numerical method based on the Fourier-Fourier transform approach for modeling 1-D electron plasma evolution

A numerical method is presented for studying one-dimensional electron plasma evolution under typical interplanetary conditions. The method applies the Fourier-Fourier transform approach to a plasma model that is a generalization of the electrostatic Vlasov-Poisson system of equations. Conservation laws that are modified to include the plasma model generalization and also the boundary effects of nonperiodic solutions are given. A new conservation law for entropy in the transformed space is then introduced. These conservation laws are used to verify the numerical solutions. A discretization error analysis is presented. Two numerical instabilities and the methods used for their suppression are treated. It is shown that in interplanetary plasma conditions, the bump-on-tail instability produces significant excitation of plasma oscillations at the Bohm-Gross frequency and its second harmonic. An explanation of the second harmonic excitation is given in terms of wave-wave coupling during the growth phase of the instability.

Klimas, A. J.↗

NASA NTRS1983

A mechanism for plasma waves at the harmonics of the plasma frequency in the electron foreshock boundary

Klimas, A. J.↗

NASA NTRS1994

A splitting algorithm for Vlasov simulation with filamentation filtration

A Fourier-Fourier transformed version of the splitting algorithm for simulating solutions of the Vlasov-Poisson system of equations is introduced. It is shown that with the inclusion of filamentation filtration in this transformed algorithm it is both faster and more stable than the standard splitting algorithm. It is further shown that in a scalar computer environment this new algorithm is approximately equal in speed and far less noisy than its particle-in-cell counterpart. It is conjectured that in a multiprocessor environment the filtered splitting algorithm would be faster while producing more precise results.

Klimas, A. J.↗

DOE OSTI2020

Hot electrons between cold walls

We consider electron cooling in a collisionless plasma slab between two cold and freely emitting walls. Numerical calculations suggest a counterintuitive behavior of this system: the cooling rate slows down and eventually stops, leaving the system with a significant fraction of its initial thermal energy. Analytical treatment within the Vlasov–Poisson model reveals a set of steady-states with a two-component distribution of electrons: the primary electrons trapped within the potential wells and the secondary electrons forming the counterstreaming beams. We show that such steady-states are linearly stable with respect to one-dimensional perturbations. Establishment of a particular steady-state depends on initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1971

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗