Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA FLUID EQUATIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Observation of dust acoustic shock wave in a strongly coupled dusty plasma

Dust acoustic shock wave is observed in a strongly coupled laboratory dusty plasma. A supersonic flow of charged microparticles is allowed to perturb a stationary dust fluid to excite dust acoustic shock wave. The evolution process beginning with steepening of initial wave front and then formation of a stable shock structure is similar to the numerical results of the Korteweg-de Vries-Burgers equation. The measured Mach number of the observed shock wave agrees with the theoretical results. Reduction of shock amplitude at large distances is also observed due to the dust neutral collision and viscosity effects. The dispersion relation and the spatial damping of a linear dust acoustic wave are also measured and compared with the relevant theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Nonlinear interaction of proton whistler with kinetic Alfvén wave to study solar wind turbulence

This paper presents the nonlinear interaction between small but finite amplitude kinetic Alfvén wave (KAW) and proton whistler wave using two-fluid model in intermediate beta plasma, applicable to solar wind. The nonlinearity is introduced by modification in the background density. This change in density is attributed to the nonlinear ponderomotive force due to KAW. The solutions of the model equations, governing the nonlinear interaction (and its effect on the formation of localized structures), have been obtained using semi-analytical method in solar wind at 1AU. It is concluded that the KAW properties significantly affect the threshold field required for the filament formation and their critical size (for proton whistler). The magnetic and electric field power spectra have been obtained and their relevance with the recent observations of solar wind turbulence by Cluster spacecraft has been pointed out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2013

Closure and transport theory for high-collisionality electron-ion plasmas

Systems of algebraic equations for a high-collisionality electron-ion plasma are constructed from the general moment equations with linearized collision operators [J.-Y. Ji and E. D. Held, Phys. Plasmas 13, 102103 (2006) and J.-Y. Ji and E. D. Held, Phys. Plasmas 15, 102101 (2008)]. A systematic geometric method is invented and applied to solve the system of equations to find closure and transport relations. It is known that some closure coefficients of Braginskii [S. I. Braginskii, Reviews of Plasma Physics (Consultants Bureau, New York, 1965), Vol. 1] are in error up to 65% for some finite values of x (cyclotron frequency Multiplication-Sign electron-ion collision time) and have significant error in the large-x limit [E. M. Epperlein and M. G. Haines, Phys. Fluids 29, 1029 (1986)]. In this work, fitting formulas for electron coefficients are obtained from the 160 moment (Laguerre polynomial) solution, which converges with increasing moments for x{<=}100 and from the asymptotic solution for large x-values. The new fitting formulas are practically exact (less than 1% error) for arbitrary x and Z (the ion charge number, checked up to Z = 100). The ion coefficients for equal electron and ion temperatures are moderately modified by including the ion-electron collision operator. When the ion temperature is higher than the electron temperature, the ion-electron collision and the temperature change terms in the moment equations must be kept. The ion coefficient formulas from 3 moment (Laguerre polynomial) calculations, precise to less than 0.4% error from the convergent values, are explicitly written.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Numerical investigation on two-fluid model (micropolar-Newtonian) for pulsatile flow of blood in a tapered arterial stenosis with radially variable magnetic field and core fluid viscosity

In the present study, an unsteady two-fluid model of blood through a tapered arterial stenosis with variable viscosity in the presence of variable magnetic field has been investigated. In this model, blood in the core region is assumed to be micropolar and plasma in the peripheral layer as Newtonian. Finite difference method is employed in solving the governing equations. The solutions for velocity, flow rate, wall shear stress and flow resistance are computed numerically. A comparison between the velocity profiles obtained by the present study and the experimental data is made and a good agreement between them is found. The model is used to study the effect of parameters such as taper angle, radially variable viscosity, hematocrit, the coupling number, the micropolar parameter, magnetic field and plasma layer thickness on physiologically important parameters such as wall shear stress distribution in the stenotic region and flow resistance. The results in the case of constant magnetic field and variable magnetic field are compared to study the effects of magnetic field on the flow of blood. It is found that the magnitudes of wall shear stress and flow resistance are higher in the case of variable magnetic field. It is important to note that the flow resistance is higher for magneto-micropolar fluid than the micropolar fluid. The wall shear stress decreases with increasing hematocrit parameter whereas flow resistance increases with hematocrit when the applied pressure gradient is held fixed. The model clearly shows the situation of a patient with a tapered arterial stenosis under feverish condition. Due care has been taken to compare the present theoretical results with the existing ones including experimental results and good agreement between them has been observed both qualitatively and quantitatively.

42 ENGINEERING↗
DOE OSTI2012

Wave resonances and the partition of energy in ideal compressible magnetohydrodynamic fluids

Phase mixing and resonant absorption are two processes that have been under scrutiny for some time because of their role in wave damping and in providing a mechanism for heating space and laboratory plasmas. The accumulation or absorption of energy that develops within a resonant layer is usually attributed to a logarithmic singularity, but it will be shown that this build up of energy is inextricably tied to a discontinuity in the fluid displacement at the resonant point. This change in the dynamics of the problem will be examined by establishing a partition of energy that identifies and isolates the individual resonances within the fluid. The partition is based on a variational description of the Fourier transformed equations and is guided by an electrical model of the MHD system that not only illustrates the resonant structure threading the fluid but also exposes the mechanism driving the resonant absorption process. A simplified version of this model is then constructed to help determine the approximate rate of energy absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical modeling of the SNS H{sup −} ion source

Ion source rf antennas that produce H- ions can fail when plasma heating causes ablation of the insulating coating due to small structural defects such as cracks. Reducing antenna failures that reduce the operating capabilities of the Spallation Neutron Source (SNS) accelerator is one of the top priorities of the SNS H- Source Program at ORNL. Numerical modeling of ion sources can provide techniques for optimizing design in order to reduce antenna failures. There are a number of difficulties in developing accurate models of rf inductive plasmas. First, a large range of spatial and temporal scales must be resolved in order to accurately capture the physics of plasma motion, including the Debye length, rf frequencies on the order of tens of MHz, simulation time scales of many hundreds of rf periods, large device sizes on tens of cm, and ion motions that are thousands of times slower than electrons. This results in large simulation domains with many computational cells for solving plasma and electromagnetic equations, short time steps, and long-duration simulations. In order to reduce the computational requirements, one can develop implicit models for both fields and particle motions (e.g. divergence-preserving ADI methods), various electrostatic models, or magnetohydrodynamic models. We have performed simulations using all three of these methods and have found that fluid models have the greatest potential for giving accurate solutions while still being fast enough to perform long timescale simulations in a reasonable amount of time. We have implemented a number of fluid models with electromagnetics using the simulation tool USim and applied them to modeling the SNS H- ion source. We found that a reduced, single-fluid MHD model with an imposed magnetic field due to the rf antenna current and the confining multi-cusp field generated increased bulk plasma velocities of > 200 m/s in the region of the antenna where ablation is often observed in the SNS source. We report here on comparisons of simulated plasma parameters and code performance using more accurate physical models, such as two-temperature extended MHD models, for both a related benchmark system describing a inductively coupled plasma reactor, and for the SNS ion source. We also present results from scaling studies for mesh generation and solvers in the USim simulation code.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Spectral solver for multi-scale plasma physics simulations with dynamically adaptive number of moments

A spectral method for kinetic plasma simulations based on the expansion of the velocity distribution function in a variable number of Hermite polynomials is presented. The method is based on a set of non-linear equations that is solved to determine the coefficients of the Hermite expansion satisfying the Vlasov and Poisson equations. In this paper, we first show that this technique combines the fluid and kinetic approaches into one framework. Second, we present an adaptive strategy to increase and decrease the number of Hermite functions dynamically during the simulation. The technique is applied to the Landau damping and two-stream instability test problems. Performance results show 21% and 47% saving of total simulation time in the Landau and two-stream instability test cases, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Surrogate approximation of the Grad–Shafranov free boundary problem via stochastic collocation on sparse grids

In magnetic confinement fusion devices, the equilibrium configuration of a plasma is determined by the balance between the hydrostatic pressure in the fluid and the magnetic forces generated by an array of external coils and the plasma itself. The location of the plasma is not known a priori and must be obtained as the solution to a free boundary problem. The partial differential equation that determines the behavior of the combined magnetic field depends on a set of physical parameters (location of the coils, intensity of the electric currents going through them, magnetic permeability, etc.) that are subject to uncertainty and variability. The confinement region is in turn a function of these stochastic parameters as well. In this work, we consider variations on the current intensities running through the external coils as the dominant source of uncertainty. This leads to a parameter space of dimension equal to the number of coils in the reactor. With the aid of a surrogate function built on a sparse grid in parameter space, a Monte Carlo strategy is used to explore the effect that stochasticity in the parameters has on important features of the plasma boundary such as the location of the x-point, the strike points, and shaping attributes such as triangularity and elongation. Furthermore, the use of the surrogate function reduces the time required for the Monte Carlo simulations by factors that range between 7 and over 30.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vacuum magnetic fields with exact quasisymmetry near a flux surface. Part 1. Solutions near an axisymmetric surface

While several results have pointed to the existence of exactly quasisymmetric fields on a surface (Garren & Boozer, Phys. Fluids B, vol. 3, 1991, pp. 2805–2821; 2822–2834; Plunk & Helander, J. Plasma Phys. , vol. 84, 2018, 905840205), we have obtained the first such solutions using a vacuum surface expansion formalism. We obtain a single nonlinear parabolic partial differential equation for a function $\eta$ such the field strength satisfies $B = B(\eta )$ . Closed-form solutions are obtained in cylindrical, slab and isodynamic geometries. Numerical solutions of the full nonlinear equations in general axisymmetric toroidal geometry are obtained, resulting in a class of quasihelical local vacuum equilibria near an axisymmetric surface. The analytic models provide additional insight into general features of the nonlinear solutions, such as localization of the surface perturbations on the inboard side. The local solutions thus obtained can be continued globally only for special initial surfaces.

Physics↗
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 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.↗
DOE OSTI2012

Arbitrary amplitude dust ion acoustic solitary waves in a magnetized suprathermal dusty plasma

The linear and nonlinear dust-ion-acoustic (DIA) wave propagating obliquely with respect to an external magnetic field is studied in a magnetized complex plasma which consists of a cold ion fluid, superthermal electrons, and static dust particles. The propagation properties of two possible modes (in the linear regime) are investigated. It is found that the electron suprathermality and the electron population decrease the phase velocities of both modes, while obliqueness leads to increase of separation between two modes. An energy-like equation derived to describe the nonlinear evolution of DIA solitary waves. The influences of electron suprathermality, obliqueness, and electron population on the existence domain of solitary waves and the soliton characteristics are examined. It is shown that the existence domain of the DIA soliton and its profile are significantly depending on the deviation of electrons from thermodynamic equilibrium, electrons population, and obliqueness. It is also found that the suprathermal plasma supports the DIA solitons with larger amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗