Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

A transverse Kelvin-Helmholtz instability in a magnetized plasma

An analysis is conducted of the transverse Kelvin-Helmholtz instability in a magnetized plasma for unstable flute modes. The analysis makes use of a two-fluid model. Details regarding the instability calculation are discussed, taking into account the ion continuity and momentum equations, the solution of a zero-order and a first-order component, and the properties of the solution. It is expected that the linear calculation conducted will apply to situations in which the plasma has experienced no more than a few growth periods.

Kintner, P.↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Derivation of a Multi-Species Mix Model

Multi-species flow takes place when two or more gases, plasmas or fluids are flowing and interacting with one another through particle collisions, chemical reactions, turbulence, and thermonuclear processes; this list is not exhaustive. It has a wide range of applications in geological, biological and engineering systems. In what follows species refers to different species of atoms or molecules that may or may not be partially or fully ionized. The interactions between these species takes place at an atomic or molecular level. The species are electrically neutral. These are very important points, because the coupling between species that are interacting on a mesoscale level or are immiscible would be very different than what is being considered in this work. In addition, in section 4, we present a heuristic argument to enhance the collisional mean free path to incorporate the effects of turbulent eddies by use of Prandtl’s assumption that the integral scale length of a turbulent boundary layer is proportional to the width of the turbulent layer. This leads to coupling terms that are non-Newtonian, because they are a function of the dynamics of the flow. It is also assumed that the turbulent stress terms are much smaller than the species partial pressures. This may be a good assumption if the species are hot, dense plasmas that exist in ICF targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Relativistic Alfvén Waves Entering Charge-starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A /j 0 $\ll$ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Particle-mesh techniques

This is an introduction to numerical Particle-Mesh techniques, which are commonly used to model plasmas, gravitational N-body systems, and both compressible and incompressible fluids. The theory behind this approach is presented, and its practical implementation, both for serial and parallel machines, is discussed. This document is based on a four-hour lecture course presented by the author at the NASA Summer School for High Performance Computational Physics, held at Goddard Space Flight Center.

Macneice, Peter↗
DOE OSTI2020

A cancellation problem in hybrid particle-in-cell schemes due to finite particle size

The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the different spatial discretizations of ions as finite-spatial-size particles and electrons as a grid-based fluid can lead to significant numerical wave dispersion errors in the long wavelength limit (kd i «1, where k is the wavenumber and di is the ion skin-depth). The problem occurs when high-order particle-grid interpolations, or grid-based smoothing, spreads the electric field experienced by the ions across multiple spatial cells and leads to inexact cancellation of electric field terms in the total (ion + electron) momentum equation. Practical requirements on the mesh spacing Δx/d i are suggested to bound these errors from above. The accuracy impact of not respecting these resolution constraints is shown for a non-linear shock problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

STORM: Scrape-off layer turbulence in tokamak fusion reactors

The scrape-off layer of a tokamak fusion reactor carries the plasma exhaust from the hot core plasma to the material surfaces of the reactor vessel. The heat loads imposed by the exhaust are a critical limit on the performance of fusion power plants. Turbulent transport of the plasma regulates the width of the scrape-off layer plasma and must be modelled to understand the intensity of these heat loads. STORM is a plasma turbulence code capable of simulating three dimensional turbulence across the full scrape-off layer of a tokamak fusion reactor, using a drift reduced, collisional fluid model. STORM uses mostly finite difference schemes, with a staggered grid in the direction parallel to the magnetic field. We describe the model, geometry and initialisation options used by STORM, as well as the numerical methods, which are implemented using the BOUT++ plasma simulation framework. BOUT++ has been enhanced alongside the development of STORM, providing better support for staggered grid methods. We summarise these enhancements, including a detailed explanation of the parallel derivative methods, which underwent a major update for version 4 of BOUT++.

BOUT++↗
DOE OSTI2020

A Summary of the UCLA HANE-Laser Experiment: 2011-2020

In the early years of this century, there was renewed interest at DTRA in artificial radiation belts, the dynamics of high altitude nuclear explosions that produced them, the development of large-scale kinetic plasma computer models at LLNL and LANL (particularly in the form of “hybrid” {i.e., particle ions, massless fluid electrons} algorithms), and the building of a laser facility at UCLA, under the direction of Prof. Niemann, which was connected to the Large Plasma Device (LAPD), a DOE user facility. These advances sparked the idea for a new laser experimental program to examine early-time HANE issues, reviving the concept from a former program at NRL in the 90’s. The basic motivation for this new effort can be traced back to the first DTRA artificial radiation belt workshop at Stanford in 2009. Subsequent discussions then led to a formal proposal from UCLA that was submitted to DTRA (Grant Jones) in 2010, vigorously reviewed, and finally approved in 2011, with funding begun in 2012. An historical perspective of this development process was presented at a DTRA review last year. In this document we review and summarize the major achievements of the DTRA-sponsored UCLA HANE-laser experiment over the past eight years. Our purpose is to briefly describe the major achievements of this program, as documented in the included extensive list of journal publications [which does not include all the publications nor any of the many invited and contributed presentations of the UCLA group], the role of the national laboratories in this effort, and how this work has impacted (and will continue to improve) our understanding of high altitude nuclear events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Hydrodynamical constraints on cosmic-ray acceleration in relativistic shocks

A two-fluid hydrodynamical model governing the transport of cosmic rays in a relativistically moving background plasma is developed. The equations are used to discuss the time asymptotic structure of a relativistic, plane-parallel shock wave significantly modified by first-order Fermi acceleration of cosmic rays. The model allows for an anisotropic cosmic-ray pressure tensor with pressure components parallel and perpendicular to the shock normal. It is surmised that the well-known energy threshold problem for diffusive shock acceleration that downstream particles have a large enough velocity to overtake the shock and subsequently scatter in the upstream medium before returning to the downstream medium may be expressed hydrodynamically by the condition V less than V(cr), where V is the fluid velocity relative to the shock and V(cr) is the relativistic cosmic-ray sound speed. Astrophysical implications of the results are briefly discussed.

Webb, G. M.↗
DOE OSTI2022

Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++

Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional modeling of a plasma in a strong azimuthal magnetic field

Three-dimensional magnetohydrodynamic (MHD) simulations are able to model the generation of disc-shaped plasma, driven by laser ablation from a current carrying rod in a pulsed-power machine producing azimuthal magnetic fields of 2–3 MG. The plasma at such extreme conditions is unique in that the parameter space for the plasma β and Hall parameter χ transition from below unity to greater than unity at different stages of the plasma generation. In simulations, the formation of the plasma disc in the azimuthal direction, is driven by heat flux from the laser spot and depends on the set of transport coefficients used in simulations. The most recent set of transport coefficients leads to the formation of plasma ejecta at the back end of the rod, that qualitatively matches experiments. Specifically, the cross-gradient Nernst effect, which twists magnetic field, is shown to have a large effect on the shape of the back end ejecta. In the direction along the axis of the rod there is propagation of perturbations from the disc as observed in experiments. In simulations the period of temperature perturbations is in good agreement with experimental results. Furthermore, an instability due to coupling of heat flux and the magnetic field advection provides a possible explanation for perturbation growth along the axis of the rod, and the instability growth rate is consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Heat pump via charge incompressibility in a collisional magnetized multi-ion plasma

Stratification due to ion-ion friction in a magnetized multiple-ion species plasma is shown to be accompanied by a heat pump effect, identified for the first time here, transferring heat from one ion species to another as well as from one region of space to another. Here, the heat pump is produced via newly identified heating mechanisms associated with charge incompressibility and the Ettingshausen effect. Besides their academic interest, these new effects may have useful applications to plasma technologies that involve rotation or compression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Turbulent generation of outward-traveling interplanetary Alfvenic fluctuations

Additional evidence is presented for the dynamic alignment of the fluctuating plasma fluid velocity and the fluctuating magnetic field in two-dimensional MHD turbulence. The evidence is seen as indicating that the interplanetary fluctuations need not represent time-asymptotic states. Rather, it is contended that interplanetary Alfvenic fluctuations are generated in situ in the early stage of dynamic alignment. In addition to providing an explanation for inertial-range alignment of the fluctuating plasma fluid velocity and the fluctuating magnetic field in the solar wind, the model presented here accounts for the fact that the sign of this correlation (relative to the direction of the mean magnetic field) corresponds to the observed 'direction of propagation', that is, away from the sun. Evidence is found to support the contention that MHD fluctuations in the solar wind resemble turbulence arising from local stirring of the medium.

Matthaeus, W. H.↗