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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Fluid hydrogen at high density - Pressure ionization

The Helmholtz-free-energy model for nonideal mixtures of hydrogen atoms and molecules by Saumon and Chabrier (1991) is extended to describe dissociation and ionization in similar mixtures in chemical equilibrium. A free-energy model is given that describes partial ionization in the pressure and temperature ionization region. The plasma-phase transition predicted by the model is described for hydrogen mixtures including such components as H2, H, H(+), and e(-). The plasma-phase transition has a critical point at Tc = 15,300 K and Pc = 0.614 Mbar, and thermodynamic instability is noted in the pressure-ionization regime. The pressure dissociation and ionization of fluid hydrogen are described well with the model yielding information on the nature of the plasma-phase transition. The model is shown to be valuable for studying dissociation and ionization in astrophysical objects and in high-pressure studies where pressure and temperature effects are significant.

Saumon, Didier↗
DOE OSTI2025

The impact of non-local parallel electron transport on plasma-impurity reaction rates in tokamak scrape-off layer plasmas

Abstract Plasma-impurity reaction rates are a crucial part of modelling tokamak scrape-off layer (SOL) plasmas. To avoid calculating the full set of rates for the large number of important processes involved, a set of effective rates are typically derived which assume Maxwellian electrons. However, non-local parallel electron transport may result in non-Maxwellian electrons, particularly close to divertor targets. Here, the validity of using Maxwellian-averaged rates in this context is investigated by computing the full set of rate equations for a fixed plasma background from kinetic and fluid SOL simulations. We consider the effect of the electron distribution as well as the impact of the electron transport model on plasma profiles. Results are presented for lithium, beryllium, carbon, nitrogen, neon and argon. It is found that electron distributions with enhanced high-energy tails can result in significant modifications to the ionisation balance and radiative power loss rates from excitation, on the order of 50%–75% for the latter. Fluid electron models with Spitzer-Härm or flux-limited Spitzer-Härm thermal conductivity, combined with Maxwellian electrons for rate calculations, can increase or decrease this error, depending on the impurity species and plasma conditions. Based on these results, we also discuss some approaches to experimentally observing non-local electron transport in SOL plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulation of supersonic molecular beam injection fueling into H-mode plasmas on EAST using BOUT++

Transport dynamic 2D simulations for supersonic molecular beam injection (SMBI) fueling into H-mode deuterium plasmas on the Experimental Advanced Superconducting Tokamak (EAST) is first simulated using a seven-field two-fluid model in the BOUT++ framework. The SMB is assumed to be injected into plasma from the midplane at the low field side with a fixed width and constant molecular flux. The different densities and injection velocities of SMBI are investigated within the upper single-null geometry of EAST. The simulations indicate that the SMBI has a self-shielding effect on molecules' inward transport into the plasma, and the deposition of SMBI leads to a large increase in plasma density and decrease in plasma temperature. There is a velocity threshold for SMB penetrating the pedestal and depositing at the top of density pedestal. However, the deposition point would be back toward the plasma boundary after SMB arriving at the deepest penetration position. Comparing the different molecular injection velocities and densities, the outcomes demonstrate that the depth of deposition is closely related to the injection velocity rather than the molecular density. The simulated results show good agreement with the EAST experiments by comparing the electron density profiles obtained by simulation and experiment, respectively. These results will be helpful for guiding future experiments, as well as the design of the SMBI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗
DOE OSTI2023

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2021

Sheath collapse at critical shallow angle due to kinetic effects

The Debye sheath is known to vanish completely in magnetised plasmas for a sufficiently small electron gyroradius and small angle between the magnetic field and the wall. This angle depends on the current onto the wall. When the Debye sheath vanishes, there is still a potential drop between the wall and the plasma across the magnetic presheath. The magnetic field angle corresponding to the predicted sheath collapse is shown to be much smaller than previous estimates, scaling with the electron-ion mass ratio and not with the square root of the mass ratio. This is shown to be a consequence of the kinetic electron and finite ion orbit width effects, which are not captured by fluid models. The wall potential with respect to the bulk plasma at which the Debye sheath vanishes is calculated. Above this wall potential, it is possible that the Debye sheath will invert.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Plasma density enhancement at the Comet Halley diamagnetic cavity boundary

A one-dimensional single-fluid MHD model was used by Cravens (1989) to predict the existence of a narrow (50 km) layer of enhanced plasma density at the boundary of the Comet Halley diamagnetic cavity. The existence of such a layer was confirmed by measurements made by the Giotto ion mass spectrometer (Goldstein et al., 1989). Here, the time-dependent coupled continuity equations are solved for several species including H3O(+), H2O(+), OH(+), O(+), NH4(+), NH3(+), CH4(+), and CH3(+). For several ion species, factors affecting the magnitude of this enhancement and its relationship to the thickness of the transition layer are investigated. For example, ion species with short chemical lifetimes are shown to have smaller density enhancements than species with longer lifetimes. The cometocentric distance of the cavity boundary also strongly affects the magnitude of the density enhancement; the enhancement increases with increasing distance.

Keller, C. N.↗
NASA NTRS2011

Magnetofluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. In this report, the interstellar pickup protons are treated as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. Where possible, the model is compared with Voyager data. Initial results from generalization to a three-fluid model is described elsewhere in this session.

Goldstein, Melvyn L.↗
NASA NTRS1999

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗
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↗