Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid models”

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 505 records · Page 28

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗

The value of the 4-day headdown bedrest model for screening countermeasures

In order to evaluate the benefits of periodic exposure to the +G(z) vector as a countermeasure to the physiological responses to minus 6 degree head down bedrest (HDT), we considered a two-tiered approach: (a) to use 4 days HDT as a quick and inexpensive means of screening countermeasures, (b) to use a 60 day HDT to validate the most promising candidates. The approach and results of a 4 day study are described here. Methods: Nine males were admitted to our Human Research Facility for one ambulatory control day followed by 4 days HDT and were released on the next day after completion of a peak oxygen consumption test (VO(sub 2 peak)). A battery of tests was selected and standardized to evaluate the known early effects of HDT on plasma volume, early bone markers, orthostatic tolerance, physical performance, and fluid and electrolytes and their hormone regulation. Fluid sodium (Na) and potassium (K) intake and output in the urine were monitored throughout. Plasma volume was determined with a modified Evans Blue method and orthostatic tolerance with a 60 degree head-up tilt test for 30 minutes - both of which were determined on the ambulatory control day and on day 4 of HDT. Immediately after completion of the tilt test subjects were returned to the minus 6 degree position until the next morning when a VO(sub 2 peak) (horizontal ergometer) was done. This was compared to a similar control test determined on 2 separate occasions before subject admission. Results: Four hours after going HDT produced significant decreases (p less than 0.05) in the circulating concentration of fluid and electrolyte regulating hormones. Plasma volume, orthostatic tolerance and VO(sub 2 peak) changed significantly after 4 days HDT. There was also the expected natriuresis on day 1 of HDT but no significant diuresis. The consistency of the pre-bedrest VO(sub 2 peak) tilt tests and plasma volumes was remarkable. Conclusions: The 4 day HDT model seems highly promising for screening a variety of countermeasures alone and in combination before validating their benefits in extended bedrest or flight experiments.

Vernikos, J.↗

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

Hermes-3: Multi-component plasma simulations with BOUT++

A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. Hermes-3 is built on plasma simulation framework BOUT++, consolidating earlier SD1D and Hermes models into a single code that can be configured at run-time to solve plasma models in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. Here, we describe the improved numerical algorithms and software design that have been implemented in Hermes-3. To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs). Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The theory of magnetohydrodynamic wave generation by localized sources. I - General asymptotic theory

The magnetohydrodynamic wave emission from several localized, periodic, kinematically specified fluid velocity fields are calculated using Lighthill's method for finding the far-field wave forms. The waves propagate through an isothermal and uniform plasma with a constant B field. General properties of the energy flux are illustrated with models of pulsating flux tubes and convective rolls. Interference theory from geometrical optics is used to find the direction of minimum fast-wave emission from multipole sources and slow-wave emission from discontinuous sources. The distribution of total flux in fast and slow waves varies with the ratios of the source dimensions l to the acoustic and Alfven wavelengths.

Collins, William↗

Saturn's Magnetosphere and Properties of Upstream Flow at Titan: Preliminary Results

Using Cassini Plasma Spectrometer (CAPS) Ion Mass Spectrometer (IMS) measurements, we present the ion fluid properties and its ion composition of the upstream flow for Titan's interaction with Saturn's magnetosphere. A 3D ion moments algorithm is used which is essentially model independent with only requirement is that ion flow is within the CAPS IMS 2(pi) steradian field-of-view (FOV) and that the ion 'velocity distribution function (VDF) be gyrotropic. These results cover the period from TA flyby (2004 day 300) to T22 flyby (2006 363). Cassini's in situ measurements of Saturn's magnetic field show it is stretched out into a magnetodisc configuration for Saturn Local Times (SLT) centered about midnight local time. Under those circumstances the field is confined near the equatorial plane with Titan either above or below the magnetosphere current sheet. Similar to Jupiter's outer magnetosphere where a magnetodisc configuration applies, one expects the heavy ions within Saturn's outer magnetosphere to be confined within a few degrees of the current sheet while at higher magnetic latitudes protons should dominate. We show that when Cassini is between dusk-midnight-dawn local time and spacecraft is not within the current sheet that light ions (H, 142) tend to dominate the ion composition for the upstream flow. If true, one may expect the interaction between Saturn's magnetosphere, locally devoid of heavy ions and Titan's upper atmosphere and exosphere to be significantly different from that for Voyager 1, TA and TB when heavy ions were present in the upstream flow. We also present observational evidence for Saturn's magnetosphere interaction with Titan's extended H and H2 corona which can extend approx. 1 Rs from Titan.

Sittler, E. C., Jr.↗

Modeling of small tungsten dust grains in EAST tokamak with NDS-BOUT ++

In order to investigate the transport of small dusts as well as their evolution property along their trajectories, the NDS module is developed under the BOUT++ framework, a highly desirable C++ code package to perform parallel plasma fluid simulations with an arbitrary number of equations in three-dimensional curvilinear coordinates. Due to the severe dust ablation in fusion plasmas, the dust size would decrease from micrometer to nanometer, resulting in impurities. Small dusts in the simulations here are specified as tungsten spheres with the radii on or below the order of submicrometer. The Rayleigh limit is included in the charging process when the dust is ablated to the droplet phase. The simulation results from the NDS module show that a 200 nm radius spherical tungsten dust originated from upper divertor region of EAST Tokamak is ablated completely due to the intense heating from the incoming plasma inside the core region, well consistent with the CCD footage of EAST shot # 81459. Furthermore it is found that the magnetic field dominates the dust transport when the dust radius is below 100 nm during the ablation along the trajectory. Our simulations predict that a 10 nm radius spherical tungsten dust injected from the inner midplane is well constrained by the magnetic field, and it reaches the inner divertor target with a velocity on the order of km/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental Study of Alfvén Wave Reflection from an Alfvén-speed Gradient Relevant to the Solar Coronal Holes

Abstract We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of the wavelength of the incident wave to the length scale of the gradient. Two-fluid simulations using the Gkeyll code qualitatively agree with and support the experimental findings. Our experimental results support models of wave heating that rely on wave reflection at low heights from a smooth Alfvén-speed gradient to drive turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗

Intercomparison of numerical models of flaring coronal loops

The proposed Benchmark Problem consists of an infinitesimal magnetic flux tube containing a low-beta plasma. The field strength is assumed to be so large that the plasma can move only along the flux tube, whose shape remains invariant with time (i.e., the fluid motion is essentially one-dimensional). The flux tube cross section is taken to be constant over its entire length. In planar view the flux tube has a semi-circular shape, symmetric about its midpoint s = s sub max and intersecting the chromosphere-corona interface (CCI) perpendicularly at each foot point. The arc length from the loop apex to the CCI is 10,000 km. The flux tube extends an additional 2000 km below the CCI to include the chromosphere, which initially has a uniform temperature of 8000 K. The temperature at the top of the loop was fixed initially at 2 X 1 million K. The plasma is assumed to be a perfect gas (gamma = 5/3), consisting of pure hydrogen which is considered to be fully ionized at all temperatures. For simplicity, moreover, the electron and ion temperatures are taken to be everywhere equal at all times (corresponding to an artificially enhanced electron-ion collisional coupling). While there was more-or-less unanimous agreement as to certain global properties of the system behavior (peak temperature reached, thermal-wave time scales, etc.), no two groups could claim satisfactory accord when a more detailed comparison of solutions was attempted.

Kopp, R. A.↗

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗

Hydromagnetic waves, turbulence, and collisionless processes in the interplanetary medium

An extended discussion is conducted concerning the origin and evolution of interplanetary hydromagnetic waves and turbulence, and their influence on the large scale dynamics of the solar wind. The solar wind is at present the preeminent medium for the study of hydromagnetic waves and turbulence, providing an opportunity for advancement of understanding of the most fundamental processes of the astrophysical plasmas. All interplanetary fluctuations whose time scale is observed to be greater than 1 sec can be regarded as hydromagnetic fluctuations. It has been found to be simplest, and generally very satisfactory, to model interplanetary variations as fluctuations in an MHD fluid. Attention is given to the classification of wave modes, geometrical hydromagnetics, Alfven wave pressure, rugged invariants, and the kinetic theory of collisionless processes.

Barnes, A.↗

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak (R = 0.9 m, B t = 5 T, and I p = 2 MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.

Thomas, V. A.↗