SEARCH · Engineering Papers
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.
Long-time nonlinear behavior of single-mode longitudinal plasma waves is studied on the basis of the Vlasov equation with Fokker-Planck collision terms. The resonant layer, trapped island, collisional sublayers, and X-point neighborhoods are analyzed. A nonlinear evolution equation is obtained and the wave-generated current is also evaluated.
The analysis of plasma beam and shear problems in magnetic fields is usually based on a hydromagnetic fluid model. In a low-density collisionless plasma, however, the kinetic effects of the plasma, such as finite Larmor radius effects, are not yet clearly understood. In this paper, the kinetic equilibria of plasma shears in a uniform and fixed magnetic field, with full ion motion, are discussed by solving the Vlasov equation with a given electric field and drift velocity. In this model, the ion density profile through the plasma shear layer is quite different from the one predicted by a hydromagnetic model. As a result of a complicated ion gyromotion through the shear layer, single- and double-humped ion density profiles are obtained. The dependence on the temperature and the strength of the shear will be discussed. The results show a significant difference between positive and negative shears.
Observations are presented which show that plasma parameters vary on a scale length smaller than the ion gyroradius at the interface between the plasma sheet boundary layer and the lobe. The Vlasov equation is used to investigate the properties of such a boundary layer. The existence, at the interface, of a density gradient whose scale length is smaller than the ion gyroradius implies that an electrostatic potential is established in order to maintain quasi-neutrality. Strongly sheared (scale lengths smaller than the ion gyroradius) perpendicular and parallel (to the ambient magnetic field) electron flows develop whose peak velocities are on the order of the electron thermal speed and which carry a net current. The free energy of the sheared flows can give rise to a broadband spectrum of electrostatic instabilities starting near the electron plasma frequency and extending below the lower hybrid frequency.
The time dependent behavior of a plasma which surrounds a highly biased conducting sphere is considered. This plasma is treated as either a cold two component fluid or as a warm plasma whose time development can be found by solving the Vlasov equation. Both models demonstrate oscillatory behavior, but the electric fields predicted by the models are quantitatively quite different in regions close to the surface of the sphere and very similar otherwise. A broadening of the electron distribution function with time is observed indicating local heating of the plasma near the surface of the sphere.
The study of relativistic plasma double layers is described through the solution of the one-dimensional, unmagnetized, steady-state Poisson-Vlasov equations and by means of one-dimensional, unmagnetized, particle-in-cell simulations. The thickness vs potential-drop scaling law is extended to relativistic potential drops and relativistic plasma temperatures. The transition in the scaling law for 'strong' double layers suggested by analytical two-beam models by Carlqvist (1982) is confirmed, and causality problems of standard double-layer simulation techniques applied to relativistic plasma systems are discussed.
The paper derives an expression for the parallel current density for a plasma characterized by negligible bulk flow (magnetostatic velocity and a two-component (anisotropic) pressure tensor by expanding the equilibrium Vlasov equation for each species in the adiabatic parameter until such point as a nonvanishing moment j-parallel = integral d-cubed vv-parallel f is identified. The result is a nonlocal one: it relates j-parallel at one point s along a field line to j-parallel at another (reference) point s0 plus an integral function of the pressure and magnetic field between them. The equation derived by Vasyliunas (1970) is generalized and extended to a plasma in which the pressure tensor is comprised of two separate elements, P-perpendicular and P-parallel. This equation follows when P-perpendicular is set equal to P-parallel and s and s0 are taken to be at the ionosphere and the equator. These results are compared to others in the literature.
We present a general solution to the collisionless Boltzmann (Vlasov) equation for a free-flowing plasma along a magnetic field line using Liouville's theorem, allowing for an arbitrary potential structure including non-monotonicities. The constraints of the existing collisionless kinetic transport models are explored, and the need for a more general approach to the problem of self- consistent potential energy calculations is described. Then a technique that handles an arbitrary potential energy distribution along the field line is presented and discussed. For precipitation of magnetospherically trapped hot plasma, this model yields moment calculations that vary by up to a factor of two for various potential energy structures with the same total potential drop. The differences are much greater for the high-latitude outflow scenario, giving order of magnitude variations depending on the shape of the potential energy distribution.
The composition and structure of neutral exospheres imbedded in moving plasmas can be determined by measurements of the velocity distributions of their pickup ion progeny. In turn, the velocity distributions are dependent on the spatial structure of the neutral source gases. Since Titan's neutral exosphere extends into the Saturn's magnetosphere (or solar wind) and well above its ionopause, it serves as a good place to analyze such characteristics. They are analyzed using pickup ion measurements made by the Cassini Plasma Spectrometer (CAPS) at Titan [e.g., Hartle et al., 2006] and an ion kinetic model. An early version of the model [Hartle and Sittler, 2007] is an expression describing the phase space density of pickup ions, which is derived from the Vlasov equation with an ion source that explicitly accounts for the velocity and spatial variation of the exosphere source gases. The current version used here includes exosphere source gases in three dimensions. A fundamental parameter of the phase space densities is the ratio of the gyroradius to the neutral scale height alpha, = r(sub g)/H. Titan's exosphere structure yields pickup ions whose phase space distributions are beam-like when alpha >> 1 and fluid-like when alpha << 1. Downstream from the source peak, the light pickup ions, with alpha << 1, are easily observed because their phase space densities are almost uniform over the orbit phases. On the other hand, the phase space distributions of the heavier ions, with alpha >> 1, peak over narrow velocity and spatial ranges. This beam-like nature makes it considerably more difficult to observe heavy ions because their downstream positions and viewing directions are narrowly constrained. Examples of these extremes will be discussed.
Composition and structure of neutral constituents in the lunar exosphere can be determined through measurements of phase space distributions of pickup ions borne from the exosphere [1]. An essential point made in an early study [ 1 ] and inferred by recent pickup ion measurements [2, 3] is that much lower neutral exosphere densities can be derived from ion mass spectrometer measurements of pickup ions than can be determined by conventional neutral mass spectrometers or remote sensing instruments. One approach for deriving properties of neutral exospheric source gasses is to first compare observed ion spectra with pickup ion model phase space distributions. Neutral exosphere properties are then inferred by adjusting exosphere model parameters to obtain the best fit between the resulting model pickup ion distributions and the observed ion spectra. Adopting this path, we obtain ion distributions from a new general pickup ion model, an extension of a simpler analytic description obtained from the Vlasov equation with an ion source [4]. In turn, the ion source is formed from a three-dimensional exospheric density distribution, which can range from the classical Chamberlain type distribution to one with variable exobase temperatures and nonthermal constituents as well as those empirically derived. The initial stage of this approach uses the Moon's known neutral He and Na exospheres to deriv e He+ and Na+ pickup ion exospheres, including their phase space distributions, densities and fluxes. The neutral exospheres used are those based on existing models and remote sensing studies. As mentioned, future ion measurements can be used to constrain the pickup ion model and subsequently improve the neutral exosphere descriptions. The pickup ion model is also used to estimate the exosphere sources of recently observed pickup ions on KAGUYA [3]. Future missions carrying ion spectrometers (e.g., ARTEMIS) will be able to study the lunar neutral exosphere with great sensitivity, yielding the necessary ion velocity spectra needed to further analysis of parent neutral exosphere properties.
Pickup ions formed from ionized neutral exospheres in flowing plasmas have phase space distributions that reflect their source's spatial distributions. Phase space distributions of the ions are derived from the Vlasov equation with a delta function source using three.dimensional neutral exospheres. The ExB drift produced by plasma motion picks up the ions while the effects of magnetic field draping, mass loading, wave particle scattering, and Coulomb collisions near a planetary body are ignored. Previously, one.dimensional exospheres were treated, resulting in closed form pickup ion distributions that explicitly depend on the ratio rg/H, where rg is the ion gyroradius and H is the neutral scale height at the exobase. In general, the pickup ion distributions, based on three.dimensional neutral exospheres, cannot be written in closed form, but can be computed numerically. They continue to reflect their source's spatial distributions in an implicit way. These ion distributions and their moments are applied to several bodies, including He(+) and Na(+) at the Moon, H(+2) and CH(+4) at Titan, and H+ at Venus. The best places to use these distributions are upstream of the Moon's surface, the ionopause of Titan, and the bow shock of Venus.
We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.
Approximate integrodifferential equations determining ensemble mean and covariance of particle distribution function of Vlasov plasma, discussing incompressible Navier-Stokes turbulence
Derivation of Vlasov shallow shell equations and their application to non-shallow shells
The solution of filtered Vlasov-Maxwell equations, rather than the Vlasov-Maxwell equations themselves, is shown to ameliorate the velocity space filamentation problem in collisionless plasma models. Exact field solutions and filtered velocity distribution functions are obtained without introducing errors. Proper selection of the filter width is demonstrated to inhibit development of velocity space filamentation and, it is conjectured, position space filamentation. The results of sample calculations of both filtered and nonfiltered field solutions illustrate a high degree of agreement between both solutions, with significant savings in computational time and memory requirements with the filtered solutions.
The Hamiltonian formulations for the perturbed Vlasov–Maxwell equations and the perturbed ideal magnetohydrodynamics (MHD) equations are expressed in terms of the perturbation derivative ∂ F / ∂ ϵ ≡ [ F , S ] of an arbitrary functional F [ ψ ] of the Vlasov–Maxwell fields ψ = ( f , E , B ) or the ideal MHD fields ψ = ( ρ , u , s , B ) , which are assumed to depend continuously on the (dimensionless) perturbation parameter ϵ . In this study, [ , ] denotes the functional Poisson bracket for each set of plasma equations and the perturbation action functional S is said to generate dynamically accessible perturbations of the plasma fields. The new Hamiltonian perturbation formulation introduces a framework for functional perturbation methods in plasma physics and highlights the crucial roles played by polarization and magnetization in Vlasov–Maxwell and ideal MHD perturbation theories. One application considered in this paper is a formulation of plasma stability that guarantees dynamical accessibility and leads to a natural generalization to higher-order perturbation theory.
Collective effects in collisionless plasmas in laboratory in case of Vlasov-Maxwell equations with Landau damping
Parareal is a well-known parallel-in-time algorithm that combines a coarse and fine propagator within a parallel iteration. It allows for large-scale parallelism that leads to significantly reduced computational time compared to serial time-stepping methods. However, like many parallel-in-time methods it can fail to converge when applied to non-diffusive equations such as hyperbolic systems or dispersive nonlinear wave equations. Here, this paper explores the use of exponential integrators within the Parareal iteration. Exponential integrators are particularly interesting candidates for Parareal because of their ability to resolve fast-moving waves, even at the large stepsizes used by coarse propagators. This work begins with an introduction to exponential Parareal integrators followed by several motivating numerical experiments involving the nonlinear Schrödinger equation. These experiments are then analyzed using linear analysis that approximates the stability and convergence properties of the exponential Parareal iteration on nonlinear problems. The paper concludes with two additional numerical experiments involving the dispersive Kadomtsev-Petviashvili equation and the hyperbolic Vlasov-Poisson equation. These experiments demonstrate that exponential Parareal methods offer improved time-to-solution compared to serial exponential integrators when solving certain non-diffusive equations.
Solution of initial-value problem for linearized boltzmann equation for longitudinal and transversal plasma oscillations