Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA FLUID EQUATIONS”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Amplitude modulation of quantum-ion-acoustic wavepackets in electron-positron-ion plasmas: Modulational instability, envelope modes, extreme waves

A semirelativistic fluid model is employed to describe the nonlinear amplitude modulation of low-frequency (ionic scale) electrostatic waves in an unmagnetized electron-positron-ion plasma. Electrons and positrons are assumed to be degenerated and inertialess, whereas ions are warm and classical. A multiscale perturbation method is used to derive a nonlinear Schrödinger equation for the envelope amplitude, based on which the occurrence of modulational instability is investigated in detail. Various types of localized ion acoustic excitations are shown to exist, in the form of either bright type envelope solitons (envelope pulses) or dark-type envelope solitons (voids, holes). The plasma configurational parameters (namely, the relativistic degeneracy parameter, the positron concentration, and the ionic temperature) are shown to affect the conditions for modulational instability significantly, in fact modifying the associated threshold as well as the instability growth rate. In particular, the relativistic degeneracy parameter leads to an enhancement of the modulational instability mechanism. Furthermore, the effect of different relevant plasma parameters on the characteristics (amplitude, width) of these envelope solitary structures is also presented in detail. Finally, the occurrence of extreme amplitude excitation (rogue waves) is also discussed briefly. Our results aim at elucidating the formation and dynamics of nonlinear electrostatic excitations in superdense astrophysical regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

The relationship between flux coordinates and equilibrium-based frames of reference in fusion theory

The properties of two local reference frames based on the magnetic field and the current density are investigated for magnetized plasmas in toroidal geometry with a symmetric angle. The magnetic field-based local frame of reference has been well-studied, for example, by Dewar and colleagues [Phys. Fluids 27, 1723 (1984)]. An analogous frame based on the current density vector is possible because it is also divergence free and perpendicular to the gradient of the poloidal flux. The concept of straightness of a vector is introduced and used to elucidate the Boozer and Hamada coordinate systems. The relationship of these local frames to the more well-known Frenet frame of reference, which specifies a curve in terms of curvature and torsion, is given. As an example of the usefulness of this formal relationship, we briefly review the ideal MHD theory and its use. We also present a new annihilation operator, useful for eliminating shorter time scales than the time scale of interest, for deriving the inner layer equations of Glasser et al. [Phys. Fluids 18, 875 (1975)]. Compared to the original derivation that is based on the local frame of reference in terms of the magnetic field, the new annihilation operator that is based on the local frame of reference in terms of the current density simplifies the derivation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Gyrokinetic simulations of magnetic reconnection

Fast magnetic reconnection, believed to be a mechanism for rearranging the magnetic topology and creating energetic particles in many astrophysical and laboratory plasmas, is investigated with the nonlinear gyrokinetic code Gene. After some code-code benchmarking, extensive linear studies are presented, covering all relevant parameter dependencies of two-dimensional slab reconnection. The results are used to ascertain the validity of a fluid model and understand for which parameters it fails to describe the physics correctly. The nonlinear phase is studied for two scenarios: decaying and driven turbulence. In the former case, the initially injected energy is cascading towards the largest scales of the system, whereas a fully turbulent, quasi-stationary state develops if the system is driven through a Krook-type term in the gyrokinetic Vlasov equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Self-similar and diffusive expansion of nonextensive plasmas

Exact analytical self-similar solution is presented for free collisionless expansion of a two-component plasma of inertial ions and nonextensive electrons into vacuum, using the generalized nonextensive velocity distribution for electrons. Furthermore, a hydrodynamic model of plasma expansion in the presence of the ambipolar diffusion caused by collisions among the plasma species, such as electrons and ions, is developed and a Fokker-Planck-like generalized diffusion equation for steady-state expansion of a nonextensive electron-ion plasma is derived. For the case of generalized statistics and in the absence of particle diffusion, the density, velocity, electric potential, and field of expansion profiles are exactly obtained and studied in terms of the self-similar parameter. It is found that superthermal electrons lead to an accelerated expansion of plasma compared to that of Maxwellian electrons. It is also revealed that the nonextensivity parameter plays a fundamental role on the density, velocity, electric potential, and field configuration of the expansion. Therefore, one is able to distinguish three different regimes q < 1, q = 1, and q > 1 for expansion corresponding to sub-nonextensive, extensive, and super-nonextensive statistical profiles for electrons, respectively. Current research can provide useful information and suggests techniques for investigation of the involved statistical mechanism on the role of the energetic electron fluid in the expansion of plasma in strong pulsed laser-matter interaction experiments. It is also shown that the particle diffusion expansion mechanism becomes more dominant for relatively large values of the nonextensivity parameter, q.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

The evolution of soft X-ray-emitting flare loops

We have constructed a numerical model for a cooling flare loop in which the complete set of single-fluid equations in a one-dimensional geometry (i.e., parallel to the magnetic field) is solved. Both evaporative and static boundary conditions for the chromosphere-corona interface have been developed. This model is used to investigate the effects of initial temperature and density, loop geometry, and boundary conditions on the form of the plasma evolution and the soft X-ray emission. The results are then compared with Skylab S-056 observations of the 1973 August 9 flare. For this comparison, and under the present assumptions, we conclude that even highly compact flares must have a multiloop structure similar to large flares, and that both radiative and conductive cooling are necessary to explain the observations. The data appear to be consistent with the predicted emission from a combination of evaporative cooling loops.

Antiochos, S. K.↗
NASA NTRS1988

Multistream hydrodynamic modeling of interhemispheric plasma flow

Interhemispheric plasma flow was simulated using one-stream and two-stream hydrodymic models in order to test the suggestion of Banks et al. (1971) and others that the collision of high-speed flows originating from the conjugate hemispheres will cause the formation of a pair of shocks. The single-fluid hydrodynamic equations were modified to include multiple ion streams, allowing for the possibility of counterstreaming flow. It was found that a counterstreaming of ion streams from conjugate hemispheres does occur during the early stages of the refilling of plamaspheric flux tubes, and that a pair of reverse shocks does form. These shocks form away from the equator, and their subsequent motion creates conditions similar to those predicted by the single-stream hydrodynamic models. The findings support the conclusion of earlier studies that the refilling of the plasmasphere occurs from the equatorial region downward.

Rasmussen, C. E.↗
NASA NTRS1993

Nearly incompressible fluids. II - Magnetohydrodynamics, turbulence, and waves

The theory of nearly incompressible (NI) fluid dynamics developed previously for hydrodynamics is extended to magnetohydrodynamics (MHD). Based on a singular expansion technique, modified systems of fluid equations are obtained for which the effects of compressibility are admitted only weakly in terms of the different possible incompressible solutions. NI MHD represents the interface between the compressible and incompressible magnetofluid descriptions in the subsonic regime. It is shown that three distinct NI descriptions exist corresponding to each of the three possible plasma beta regimes. The detailed theory of weakly compressible corrections to the various incompressible MHD descriptions is presented, and the implications for the solar wind are discussed.

Zank, G. P.↗
NASA NTRS1993

Mixing and unmixedness in plasma jets 1: Near-field analysis

The flow characteristics in the near-field of a plasma jet are simulated with a two-fluid model. This model accounts for both gradient-diffusion mixing and uni-directional sifting motion resulting from pressure-gradient-body-force imbalance. This latter mechanism is believed to be responsible for the umixedness observed in plasma jets. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Transport equations are solved for the individual plasma and ambient gas velocities, temperatures and volume fractions. Empirical relations are employed for the interface transfers of mass, momentum and heat. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the torch. By adjusting the inlet condition, the model adequately reproduces the experimental data.

Ilegbusi, Olusegun J.↗
DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of an implicit electromagnetic capability for a hybrid gyrokinetic ion-fluid electron model

Here we report on the development and implementation of a hybrid kinetic ion–fluid electron model for electromagnetic COGENT simulations of edge plasmas. COGENT is a finite-volume gyrokinetic code that employs a locally field-aligned coordinate system combined with a mapped multi-block grid technology to handle strongly anisotropic edge plasma turbulence. The simulation model involves the long-wavelength limit of the ion gyrokinetic equation coupled to the vorticity and Ohm's law equations for the electromagnetic field perturbations. In order to handle the fast Alfvén wave time scales, an implicit-explicit time integration approach with a physics-based preconditioner is used. The model is successfully applied to the simulations of ion-scale resistive-drift ballooning turbulence in a toroidal annulus geometry. Substantial speed-up over a fully explicit time integration approach is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Theoretical Plasma Physics

Lattice Boltzmann algorithms are a mesoscopic method to solve problems in nonlinear physics which are highly parallelized – unlike the direction solution of the original problem. These methods are applied to both fluid and magnetohydrodynamic turbulence. By introducing entropic constraints one can enforce the positive definiteness of the distribution functions and so be able to simulate fluids at high Reynolds numbers without numerical instabilities. By introducing a vector distribution function for the magnetic field one can enforce the divergence free condition on the magnetic field automatically, without the need of divergence cleaning as needed in most direct numerical solutions of the resistive magnetohydrodynamic equations. The principal reason for the high parallelization of lattice Boltzmann codes is that they consist of a kinetic collisional relaxation step (which is purely local) followed by a simple shift of the relaxed data to neighboring lattice sites. In large eddy simulations, the closure schemes are highly nonlocal – the most famous of these schemes is that due to Smagorinsky. Under a lattice Boltzmann representation the Smagorinsky closure is purely local – being simply a particular moment on the perturbed distribution fucntions. After nonlocal fluid moment models were discovered to represent Landau damping, it was found possible to model these fluid models using an appropriate lattice Boltzmann algorithm. The close to ideal parallelization of the lattice Boltzmann codes permitted us to be Gordon Bell finalists on using the Earth Simulation in Japan. We have also been involved in the radio frequency propagation of waves into a tokamak and into a spherical overdense tokamak plasma. Initially we investigated the use of a quasi-optical grill for the launching of lower hybrid waves into a tokamak. It was found that the conducting walls do not prevent the rods from being properly irradiated, the overloading of the quasi-optical grill is not severe with the electric field only being about three times higher than in the ideal case. Moreover, the quasi-optical grill was significantly fewer structural elements that the multijunction grill. Nevertheless there has not been much interest from experimental fusion groups to implementing these structures. Hence we have returned to optimizing the multijunction grill so that the large number of coupling matrix elements can be efficiently evaluated using symmetry arguments. In overdense plasmas, the standard electromagnetic waves cannot propagate into the plasma center, but are reflected at the plasma edge. By optimizing mode conversion processes (in particular, the O-X-B wave propagation of Ordinary Mode converting to an Extraordinary mode which then converts into an electrostatic Bernstein wave) one can excite within the plasma an electrostatic Bernstein wave that does not suffer density cutoffs and is absorbed on the electron cyclotron harmonics. Finally we have started looking at other mesoscopic lattice algorithms that involve unitary collision and streaming steps. Because these algorithms are unitary they can be run on quantum computers when they become available – unlike their computational cousin of lattice Boltzmann which is a purely classical code. These quantum lattice gas algorithms have been tested successfully on exact analytic soliton collision solution. These calculations are hoped to be able to study Bose Einstein condensed atomic gases and their ground states in an optical lattice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Thermal instability in post-flare plasmas

The cooling of post-flare plasmas is discussed and the formation of loop prominences is explained as due to a thermal instability. A one-dimensional model was developed for active loop prominences. Only the motion and heat fluxes parallel to the existing magnetic fields are considered. The relevant size scales and time scales are such that single-fluid MHD equations are valid. The effects of gravity, the geometry of the field and conduction losses to the chromosphere are included. A computer code was constructed to solve the model equations. Basically, the system is treated as an initial value problem (with certain boundary conditions at the chromosphere-corona transition region), and a two-step time differencing scheme is used.

Antiochos, S. K.↗
NASA NTRS2016

Slow Solar Wind: Observations and Modeling

While it is certain that the fast solar wind originates from coronal holes, where and how the slow solar wind (SSW) is formed remains an outstanding question in solar physics even in the post-SOHO era. The quest for the SSW origin forms a major objective for the planned future missions such as the Solar Orbiter and Solar Probe Plus. Nonetheless, results from spacecraft data, combined with theoretical modeling, have helped to investigate many aspects of the SSW. Fundamental physical properties of the coronal plasma have been derived from spectroscopic and imaging remote-sensing data and in situ data, and these results have provided crucial insights for a deeper understanding of the origin and acceleration of the SSW. Advanced models of the SSW in coronal streamers and other structures have been developed using 3D MHD and multi-fluid equations.

Abbo, L.↗
DOE OSTI2020

Anisotropic Electron Fluid Closure Validated by in Situ Spacecraft Observations in the far Exhaust of Guide-field Reconnection

Anisotropic electron heating predictions of the equations of state of Le et al. (2009); https://doi.org/10.1103/PhysRevLett.102.085001 (“the EoS”) are applied to reconnection in the Earth's magnetosheath. The model is applicable to open systems where electrons are streaming along magnetic field lines into the reconnection region, sourced by fixed external reservoirs of plasma ambient to the reconnection region. While spacecraft observations have previously shown the EoS to hold in the region near the X-line, we find that for an event observed far downstream (~100 d i ) from the X-line both inflows and the exhaust follow the predictions of the EoS. Furthermore, the model underlying the EoS is extended to include additional perpendicular heating terms relevant to the considered event undergoing active magnetic compression with local electron trapping. Finally, the agreement between the spacecraft observations and the EoS emphasizes the roles of the ambient plasma sources and the dynamics of trapped and passing electrons in setting and controlling the anisotropic electron heating at large scale in naturally occurring plasma configurations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A MODEL FOR THE ELECTRICALLY CHARGED CURRENT SHEET OF A PULSAR

Global-scale solutions for the magnetosphere of a pulsar consist of a region of low-lying, closed magnetic field near the star, bounded by opposite-polarity regions of open magnetic field along which the pulsar wind flows into space. Separating these open-field regions is a magnetic discontinuity—an electric current sheet—consisting of generally nonneutral plasma. We have developed a self-consistent model for the internal equilibrium structure of the sheet by generalizing the charge-neutral Vlasov/Maxwell equilibria of Harris and Hoh to allow for net electric charge. The resulting equations for the electromagnetic field are solved analytically and numerically. Our results show that the internal thermal pressure needed to establish equilibrium force balance, and the associated effective current-sheet thickness and magnetization, can differ by orders of magnitude from the Harris/Hoh charge-neutral limit. The new model provides a starting point for kinetic or fluid investigations of instabilities that can cause magnetic reconnection and flaring in pulsar magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗