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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

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.↗

NASA NTRS1991

A new approach to plasmasphere refilling: Anomalous plasma effects

During the last 10 months of the grant, both laminar and anomalous plasma processes occurring during the refilling of the outer plasmasphere after magnetic storms are investigated. Theoretical investigations were based on two types of models: (1) two-stream hydrodynamic model in which plasma flows from the conjugate ionospheres are treated as separate fluids and the ion temperature anisotropies are treated self-consistently; and (2) large-scale particle-in-cell code.

Singh, N.↗

NASA NTRS1995

Turbulence in the solar wind: Kinetic effects

Although a casual look at the fluctuating magnetic and velocity fields in the solar wind may be reminiscent of a chaotic and disordered flow, there is, nonetheless. considerable organization and structure in the temporal and spatial evolution of those fluctuations. Much of that evolution is controlled by processes operating on rather large scales for example, in the inner heliosphere, the fluctuations in magnetic and velocity are highly correlated in the sense of outward propagating Alfven waves. This correlation can be destroyed both in time and distance by the velocity gradients present between fast and slow streams and by other nonlinear processes which stir the medium, producing a turbulent cascade of energy from large to small scales. Many aspects of this turbulent evolution can be described using fluid models; however, at some scale the fluid approximation breaks down and a more detailed paradigm is necessary. The breakdown is evident in the power spectrum of magnetic fluctuations at scales approaching the wavelength of ion cyclotron waves. At those scales, as evident in Mariner 10 and other magnetometer data, the spectrum bends over and the fluctuations damp, possibly heating the ambient plasma. Some evidence for heating of the solar wind is present in the Voyager data. Fluid models can be modified to some extent to incorporate aspects of a kinetic treatment. This is done by modifying the dissipation terms in the fluid equations and by including extra terms, such as the Hall term. As the scale lengths of phenomena shrink further and approach the spatial and temporal scales characteristic of electron phenomena, the fluid description must be abandoned altogether and a fully kinetic treatment is required. One example is the generation of Langmuir solitons produced by the electron beams that generate type 3 solar radio bursts.

Goldstein, M. L.↗

NASA NTRS2002

Influence of Coronal Abundance Variations

The Multispecies loop modeling project addresses the modeling of TRACE and SOHO observations as a plasma rather than a single fluid. In the single-fluid approximation the effects of heavy species are considered in an averaged sense. Further, loop abundances are usually taken to be uniform throughout the loop, in spite of observational evidence for considerable variation in coroner abundances.

DeLuca, E. E.↗

DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

The Circulation of the Plasmasphere Fluid during the Erosion Event on September 8, 2017

A strong solar wind pressure pulse triggered the magnetic storm on September 7, 2017. Near the end of September 7, the z-component of the interplanetary magnetic field (IMF Bz) dropped from 9 to -10 nT in 30 min. The IMF Bz remained at the level of -10 nT for 2 hours and then had another rapid drop to -31 nT in 30 min. The sudden plunge of IMF Bz and the associated strong convection electric field stirred up the storm main phase with Dst falling from 0 to -122 nT from 2200 UT on September 7 to 0200 UT on the 8th. Severe plasmasphere erosion was observed on September 8 by multiple spacecraft, such as the Van Allen Probes and the Arase satellite. In this study, we examine the fate of the eroded plasmasphere particles during this event by model simulation as well as satellite data analysis. The simulation tool we use is the Space Weather Modeling Framework (SWMF)/Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) model coupled with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. One of the distinctive capabilities of the SWMF/BATSRUS-CIMI model is that it treats the cold plasmas in the plasmasphere as a separate fluid in the MHD equations. As a result, the transport and circulation of the plasmasphere fluid in the global magnetosphere can be traced and the impacts of this cold fluid on the global magnetosphere can be evaluated. In this paper, we will show how the drainage plume is formed during the storm and how the plasmasphere fluid is transported to the flank and lobe regions and eventually to the plasma sheet and reenters into the plasmasphere. We will validate our simulation by plasmasphere signatures observed in both the inner and outer magnetosphere.

Mei-Ching Fok↗

DOE OSTI2020

Experimental observation and numerical investigation of filamentary structures in magnetized plasmas

Filamentary structures in low-pressure, low-temperature plasmas are produced when strong magnetic fields are applied parallel to the electric field defined by parallel electrodes. Filamentary structures are regions within the plasma that have distinct properties such as optical brightness and extend along the magnetic field lines. In our experiments, an argon, radio frequency discharge is exposed to a strong background magnetic field in the magnetized dusty plasma experiment at Auburn University. Different forms of filamentary structures emerge in the plasma such as columns, target, or spiral-like structures. To investigate the origin and the characteristics of these patterns, we have developed a three-dimensional fluid model that can reproduce the experimental observations which enables us to investigate the underlying mechanisms of the filamentation process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Self-consistent simulation of resistive kink instabilities with runaway electrons

A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Here, given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1998

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

DOE OSTI2022

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1987

Calculation of the surface tension of liquid metals using a one-component-plasma reference system

The one-component-plasma (OCP) model is used as a reference system instead of the traditional hard-sphere fluid to calculate the liquid-vapor interfacial surface tension of liquid metals within the density functional formalism. The calculated surface tensions of the alkali metals are in excellent agreement with experiment. For the polyvalent metal Al, the result obtained is larger than experimental measurements. It is concluded that the OCP system is not suitable to describe the liquid-vapor phase transition in simple metals which have a nominal plasma parameter larger than the usual freezing value of about 178. The calculated interfacial widths in all cases are narrower than the expected experimental values.

Zeng, X. C.↗

DOE OSTI2022

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2008

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗

DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗