Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA SIMULATION”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Dielectric barrier discharge-based plasma actuator operation in artificial atmospheres for validation of modeling and simulation

We present an experimental case study of time-resolved force production by an aerodynamic plasma actuator immersed in various mixtures of electropositive (N{sub 2}) and electronegative gases (O{sub 2} and SF{sub 6}) at atmospheric pressure using a fixed AC high-voltage input of 16 kV peak amplitude at 200 Hz frequency. We have observed distinct changes in the discharge structures during both negative- and positive-going voltage half-cycles, with corresponding variations in the actuator's force production: a ratio of 4:1 in the impulse produced by the negative-going half-cycle of the discharge among the various gas mixtures we explored, 2:1 in the impulse produced by the positive-going half-cycle, and cases in which the negative-going half-cycle dominates force production (by a ratio of 1.5:1), where the half-cycles produce identical force levels, and where the positive-going half cycle dominates (by a ratio of 1:5). We also present time-resolved experimental evidence for the first time that shows electrons do play a significant role in the momentum coupling to surrounding neutrals during the negative going voltage half-cycle of the N{sub 2} discharge. We show that there is sufficient macroscopic variation in the plasma that the predictions of numerical models at the microscopic level can be validated even though the plasma itself cannot be measured directly on those spatial and temporal scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Hybrid iterative approach for simulation of radio-frequency fields in plasma

In this work, a novel iterative approach for solving discretized linear wave equations in a frequency domain, which combines time evolution with iterative relaxation schemes, is presented. In this hybrid approach, each iteration cycle consists of evolution of electromagnetic (EM) fields in time over a specified number of field periods followed by several iterative relaxations. Provided that there is sufficient dissipation, both the time evolution and the iterative relaxations contribute to the convergence of the EM fields to the solution of the formulated full wave boundary value problem. Time evolution rapidly distributes EM fields, propagating with group velocity, over the simulation domain, while the iterative relaxations smooth the fields, reducing the numerical errors such that iteration cycles converge to a steady state solution, approximating the solution of the formulated problem. This approach is intended for large scale simulations which are beyond the capabilities of direct solvers presently used for solving wave equations in the frequency domain. Lastly, the technique is demonstrated for solving wave equations on a regular grid using a cold plasma dielectric model with collisions for 2D modeling of EM fields in tokamak in an electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The GBS code for tokamak scrape-off layer simulations

We describe a new version of GBS, a 3D global, flux-driven plasma turbulence code to simulate the turbulent dynamics in the tokamak scrape-off layer (SOL), superseding the code presented by Ricci et al. (2012) [14]. The present work is driven by the objective of studying SOL turbulent dynamics in medium size tokamaks and beyond with a high-fidelity physics model. We emphasize an intertwining framework of improved physics models and the computational improvements that allow them. The model extensions include neutral atom physics, finite ion temperature, the addition of a closed field line region, and a non-Boussinesq treatment of the polarization drift. GBS has been completely refactored with the introduction of a 3-D Cartesian communicator and a scalable parallel multigrid solver. We report dramatically enhanced parallel scalability, with the possibility of treating electromagnetic fluctuations very efficiently. The method of manufactured solutions as a verification process has been carried out for this new code version, demonstrating the correct implementation of the physical model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Magnetized Disruption of Inertially Confined Plasma Flows

Herein, the creation and disruption of inertially collimated plasma flows are investigated through experiment, simulation, and analytical modeling. Supersonic plasma jets are generated by laser-irradiated plastic cones and characterized by optical interferometry measurements. Targets are magnetized with a tunable $B$ field with strengths of up to 5 T directed along the axis of jet propagation. These experiments demonstrate a hitherto unobserved phenomenon in the laboratory, the magnetic disruption of inertially confined plasma jets. This occurs due to flux compression on axis during jet formation and can be described using a Lagrangian-cylinder model of plasma evolution implementing finite resistivity. The basic physical mechanisms driving the dynamics of these systems are described by this model and then compared with two-dimensional radiation-magnetohydrodynamic simulations. Experimental, computational, and analytical results discussed herein suggest that contemporary models underestimate the electrical conductivity necessary to drive the amount of flux compression needed to explain observations of jet disruption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

ISS Plasma Interaction: Measurements and Modeling

Ionospheric plasma interaction effects on the International Space Station are discussed in the following paper. The large structure and high voltage arrays of the ISS represent a complex system interacting with LEO plasma. Discharge current measurements made by the Plasma Contactor Units and potential measurements made by the Floating Potential Probe delineate charging and magnetic induction effects on the ISS. Based on theoretical and physical understanding of the interaction phenomena, a model of ISS plasma interaction has been developed. The model includes magnetic induction effects, interaction of the high voltage solar arrays with ionospheric plasma, and accounts for other conductive areas on the ISS. Based on these phenomena, the Plasma Interaction Model has been developed. Limited verification of the model has been performed by comparison of Floating Potential Probe measurement data to simulations. The ISS plasma interaction model will be further tested and verified as measurements from the Floating Potential Measurement Unit become available, and construction of the ISS continues.

Barsamian, H.↗
NASA NTRS2011

Direct Solar Wind Proton Access into Permanently Shadowed Lunar Polar Craters

Recent analyses of Lunar Prospector neutron spectrometer (LPNS) data have suggested that high abundances of hydrogen exist within cold traps at the lunar poles, and it has often been assumed that hydrogen-bearing volatiles sequestered in permanent shadow are topographically shielded from sputtering by solar wind protons. However, recent simulation results are presented showing that solar wind protons clearly access the floor of an idealized, shadowed lunar crater through a combination of thermal and ambipolar processes, in effect creating a plasma "miniwake". These simulations are the first to model the mini-wake environment in two spatial dimensions with a self-consistent lunar surface-plasma interaction. Progress is reported on constraining the nonzero particle fluxes and energies incident on kilometer-scale shadowed topography, such as a small crater embedded within a larger one. The importance of direct solar wind proton bombardment is discussed within the context of understanding the stability and inventory of hydrogen-bearing volatiles in shadow at the lunar poles. The support of the National Lunar Science institute, the DREAM institute, LPROPS, and the NASA Postdoctoral Program at NASA Goddard Space Flight Center administered by ORAU are gratefully acknowledged.

Zimmerman, M. I.↗
NASA NTRS1992

Effects of equatorially trapped ions on refilling of the plasmasphere

The kinetic aspects of shock formation in response to the equatorial heating of ions are reported by means of small-scale kinetic simulations of countersteaming plasma flows along a spatially varying magnetic field having the features of the dipolar magnetic field lines. Attention is given to the trapping of the field-aligned flowing ions in response to the simulated equatorial heating of ions and the associated self-consistent structures in the electric potentials and the f-plasma flow pattern. Trapping is found to lead to the formation of electrostatic shocks, which propagate away from the 'equator' as seen in the hydrodynamic model of Singh and Torr (1990). The shocks form near the mirror points of the 'equatorially' heated ions. The second stage of the refilling with equatorially trapped ions is reported. A variety of structures in density and temperature distributions as the refilling proceeds slowly is found.

Singh, Nagendra↗
DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗
DOE OSTI2022

IC Report for project “Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED” [Slides]

We’ve achieved our technical goals: (1) Perform large-scale kinetic simulations of SRS and CBET a range of plasma and laser conditions; (2) Develop a new $\delta f$-Gaussian-mixture algorithm to represent non-Maxwellian distribution functions from particle trapping and time-dependent plasma response; (3) Develop physics-based preliminary nonlinear LPI models for coupling to LRT/ Mazinisin and AMP; and (4) Begin implementation of LPI sub-grid models into Mazinisin.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Quiet Plasma Source

Synthesis of plasma from separate ion and electron emitters suppresses electromagnetic interference. Source employs separate emitters for electrons and ions. Plasma source used to simulate variety of astrophysical phenomena and space-plasma effects. For example, used in studies of propagation of electromagnetic waves in plasmas. Serves as interference-free charge neutralizer.

Leung, P. L.↗
NASA NTRSDate not supplied

Nanoflare Heating Frequency of an X-ray Bright Point Observed by MaGIXS

Nanoflares have been considered to be one of the most likely candidates for heating the solar corona to multi-million kelvin temperatures. Individual nanoflares are difficult to detect with today's instruments, but their presence may be established by comparing simulated nanoflare-heated plasma emissions to observed emissions. We present a simulation of emission from an X-ray Bright Point (XBP) detected by the MaGIXS, as well as simultaneous observations from SDO/AIA and Hinode/XRT. To simulate the XBP loops, we utilize the HYDRAD code. The length and magnetic field strength of these loops are determined using potential field extrapolation of SDO/HMI's observed photospheric magnetogram. Each loop is considered to be heated by random nanoflares, the amplitude and frequency of which are governed by the length of the loop and the strength of the magnetic field. The simulated outputs are used to estimate the intensity of spectrally pure maps of Fe-18, Fe-17, Ne-9, O-8, O-9, Ne-9, and so on, which is then compared to the intensity determined from MaGIXS observations. In addition, we derived the intensity maps obtained by AIA and XRT and compared them to the observed data. The composite distribution of the delay time of the nanoflares for which the simulated loops morphology and intensities match with observation shows a peak at 200s-500s, indicating that most of the nanoflares have a high/intermediate frequency.

coronal heating↗
DOE OSTI2015

Mixing of the Interstellar and Solar Plasmas at the Heliospheric Interface

From the ideal MHD perspective, the heliopause is a tangential discontinuity that separates the solar wind plasma from the local interstellar medium plasma. There are physical processes, however, that make the heliopause permeable. They can be subdivided into kinetic and MHD categories. Kinetic processes occur on small length and time scales, and cannot be resolved with MHD equations. On the other hand, MHD instabilities of the heliopause have much larger scales and can be easily observed by spacecraft. The heliopause may also be a subject of magnetic reconnection. In this paper, we discuss mechanisms of plasma mixing at the heliopause in the context of Voyager 1 observations. Numerical results are obtained with a Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), which is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. The code can also treat nonthermal ions and turbulence produced by them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Progress towards modeling tokamak boundary plasma turbulence and understanding its role in setting divertor heat flux widths

The heat flux distributions on divertor targets in H-mode plasmas are grave concerns for future devices. We seek to simulate the tokamak boundary plasma turbulence and heat transport in the edge localized mode-suppressed regimes. The improved BOUT++ model shows that not only I p but also the radial electric field E r plays a crucial role on the turbulence behavior and sets the heat flux width. Instead of calculating E r from the pressure gradient term (diamagnetic E r ), it is measured from the plasma transport equations with the sheath potential in the scrape-off layer and the plasma density and temperature profiles inside the separatrix from the experiment. The simulation results with the new E r model have better agreement with the experiment than using the diamagnetic Er model: (1) The electromagnetic turbulence in enhanced D α H-mode shows the characteristics of quasi-coherent modes (QCMs) and broadband turbulence. The mode spectra are in agreement with the phase contrast imaging data and almost has no change in comparison to the cases which use the diamagnetic Er model; (2) the self-consistent boundary E r is needed for the turbulence simulations to get the consistent heat flux width with the experiment; (3) the frequencies of the QCMs are proportional to E r , while the divertor heat flux widths are inversely proportional to Er; and (4) the BOUT++ turbulence simulations yield a similar heat flux width to the experimental Eich scaling law and the prediction from the Goldston heuristic drift model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Energy Deposition Processes in Titan's Upper Atmosphere

Most of Titan's atmospheric organic and nitrogen chemistry, aerosol formation, and atmospheric loss are driven from external energy sources such as Solar UV, Saturn's magnetosphere, solar wind and galactic cosmic rays. The Solar UV tends to dominate the energy input at lower altitudes of approximately 1100 km but which can extend down to approximately 400 km, while the plasma interaction from Saturn's magnetosphere, Saturn's magnetosheath or solar wind are more important at higher altitudes of approximately 1400 km, but the heavy ion plasma [O(+)] of approximately 2 keV and energetic ions [H(+)] of approximately 30 keV or higher from Saturn's magnetosphere can penetrate below 950km. Cosmic rays with energies of greater than 1 GeV can penetrate much deeper into Titan's atmosphere with most of its energy deposited at approximately 100 km altitude. The haze layer tends to dominate between 100 km and 300 km. The induced magnetic field from Titan's interaction with the external plasma can be very complex and will tend to channel the flow of energy into Titan's upper atmosphere. Cassini observations combined with advanced hybrid simulations of the plasma interaction with Titan's upper atmosphere show significant changes in the character of the interaction with Saturn local time at Titan's orbit where the magnetosphere displays large and systematic changes with local time. The external solar wind can also drive sub-storms within the magnetosphere which can then modify the magnetospheric interaction with Titan. Another important parameter is solar zenith angle (SZA) with respect to the co-rotation direction of the magnetospheric flow. Titan's interaction can contribute to atmospheric loss via pickup ion loss, scavenging of Titan's ionospheric plasma, loss of ionospheric plasma down its induced magnetotail via an ionospheric wind, and non-thermal loss of the atmosphere via heating and sputtering induced by the bombardment of magnetospheric keV ions and electrons. This energy input evidently drives the large positive and negative ions observed below approximately 1100 km altitude with ion masses exceeding 10,000 daltons. We refer to these ions as seed particles for the aerosols observed below 300 km altitude. These seed particles can be formed, for example, from the polymerization of acetylene (C2H2) and benzene (C6H6) molecules in Titan's upper atmosphere to form polycyclic aromatic hydrocarbons (PAH) and/or fullerenes (C60). In the case of fullerenes, which are hollow spherical carbon shells, magnetospheric keV [O(+)] ions can become trapped inside the fullerenes and eventually find themselves inside the aerosols as free oxygen. The aerosols are then expected to fall to Titan's surface as polymerized hydrocarbons with trapped free oxygen where unknown surface chemistry can take place.

Sittler, Edward C., Jr.↗
NASA NTRSDate not supplied

Nanoflare heating in coronal X-ray Bright Points

Understanding the heating of the non-flaring solar corona is an active topic of research in heliophysics. It is well accepted that magnetic fields are mainly responsible for coronal heating. The photospheric driver randomly moves the foot-points of the magnetic field lines, and either generates waves or the quasi-static buildup of magnetic energy, depending on the timescale of motion. Heating by the dissipation of the magnetic energy is termed as DC heating while the dissipation of wave is known as AC heating mechanism. Both the AC and DC heating mechanisms can lead to impulsive heating events, termed nanoflares. The magnitude and frequency of these nanoflares determine whether they can adequately satisfy the coronal heating budget. Thus it is of great importance to study the nanoflares properties to validate their role in coronal heating. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred through indirect observational techniques combined with simulated nanoflare-heated plasma emissions. We combined the observed emission properties of coronal X-ray bright points (XBP) in EUV and X-ray wavelengths with the simulated nanoflare heated plasma to understand nanoflare properties and their contribution in coronal heating. Our results indicate that the heating of coronal XBPs could be explained by frequent nanoflares. Further we demonstrate that the sensitive spectroscopic observations in X-ray wavelength are crucial to diagnose the nanoflare properties.

Biswajit Mondal↗
DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗