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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Pedestal Bifurcation and Resonant Field Penetration at the Threshold of Edge-Localized Mode Suppression in the DIII-D Tokamak

Rapid bifurcations in the plasma response to slowly varying n = 2 magnetic fields are observed as the plasma transitions into and out of edge-localized mode (ELM) suppression. The rapid transition to ELM suppression is characterized by an increase in the toroidal rotation and a reduction in the electron pressure gradient at the top of the pedestal that reduces the perpendicular electron flow there to near zero. These events occur simultaneously with an increase in the inner-wall magnetic response. These observations are consistent with strong resonant field penetration of n = 2 fields at the onset of ELM suppression, based on extended MHD simulations using measured plasma profiles. Spontaneous transitions into (and out of) ELM suppression with a static applied n = 2 field indicate competing mechanisms of screening and penetration of resonant fields near threshold conditions. Magnetic measurements reveal evidence for the unlocking and rotation of tearinglike structures as the plasma transitions out of ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ideal plasma response to vacuum magnetic fields with resonant magnetic perturbations in non-axisymmetric tokamaks

Ideal plasma shielding and amplification of resonant magnetic perturbations in non-axisymmetric tokamak is presented by field line tracing simulation with full ideal plasma response, compared to measurements of divertor lobe structures. Magnetic field line tracing simulations in NSTX with toroidal non-axisymmetry indicate the ideal plasma response can significantly shield/amplify and phase shift the vacuum resonant magnetic perturbations. Ideal plasma shielding for n = 3 mode is found to prevent magnetic islands from opening as consistently shown in the field line connection length profile and magnetic footprints on the divertor target. It is also found that the ideal plasma shielding modifies the degree of stochasticity but does not change the overall helical lobe structures of the vacuum field for n = 3. Furthermore, amplification of vacuum fields by the ideal plasma response is predicted for low toroidal mode n = 1, better reproducing measurements of strong striation of the field lines on the divertor plate in NSTX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling of carbon and tungsten transient dust influx in tokamak edge plasma

The paper presents computer simulation studies of burst injection of carbon and tungsten dust particles in DIII-D-like edge plasmas. The injection causes a large transient influx of the low- and high-Z impurities associated with the dust ablation in the plasmas. The dust transport and the effects of the ablated impurities on the edge plasma dynamics in a modern mid-size tokamak geometry are investigated for low- and high-power plasma discharge conditions. The core plasma contamination with dust-ablated impurities and the factors affecting it are evaluated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Electric potential structures and propagation of electron beams injected from a spacecraft into a plasma

One-dimensional Vlasov simulations are used to study the propagation of electron beams injected from a spacecraft into an ambient plasma and the associated potential structures. It is shown that, for a given beam velocity, the propagation velocity and the potential structure depends on the beam density. In the case of moderate beams, a double layer forms near the beam head which propagates into the ambient plasma much more slowly than the initial beam velocity.

Singh, Nagendra↗
DOE OSTI2020

Microturbulence-mediated route for energetic ion transport and Alfvénic mode amplitude oscillations in tokamaks

The new regimes of Alfvén eigenmode (AE) induced fast ion transport in tokamak plasmas are reported, in which microturbulence plays the role of a mediator of fast ion relaxation. Coulomb collisional scattering alone leads to small AE amplitudes and does not reproduce the steady state regimes observed in experiments. We show that in nonlinear regimes the sufficiently large effective pitch angle scattering due to microturbulence can lead to steady state AE amplitude evolution. This indicates a new route for fast ion losses, which is beyond the scenarios described in “Energetic ion transport by microturbulence is insignificant in tokamaks”. As a result, microturbulence can significantly increase the amplitude of AEs in predictive simulations of burning plasma experiments such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scaling of magnetic dissipation and particle acceleration in ABC fields

Using particle-in-cell numerical simulations with electron–positron pair plasma, we study how the efficiencies of magnetic dissipation and particle acceleration scale with the initial coherence length λ 0 in relation to the system size L of the two-dimensional ‘Arnold–Beltrami–Childress’ (ABC) magnetic field configurations. Topological constraints on the distribution of magnetic helicity in two-dimensional systems, identified earlier in relativistic force-free simulations, that prevent the high-(L/λ 0 ) configurations from reaching the Taylor state, limit the magnetic dissipation efficiency to about ϵ diss ≃60% . Here we find that the peak growth time scale of the electric energy τ E,peak scales with the characteristic value of initial Alfvén velocity β A,ini like τ E,peak ∝(λ 0 /L)β$^{–3}_{A,ini}$ . The particle energy change is decomposed into non-thermal and thermal parts, with non-thermal energy gain dominant only for high initial magnetisation. The most robust description of the non-thermal high-energy part of the particle distribution is that the power-law index is a linear function of the initial magnetic energy fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
NASA NTRS2002

Double Layers in Expanding Plasmas and Their Relevance to the Auroral Plasma Processes

When a dense plasma consisting of a cold and a sufficiently warm electron population expands, a rarefaction shock forms. In the expansion of the polar wind in the magnetosphere, it has been previously shown that when a sufficiently warm electron population also exists, in addition to the usual cold ionospheric one, a discontinuity forms in the electrostatic potential distribution along the magnetic field lines. Despite the lack of spatial resolution and the assumption of quasi-neutrality in the polar wind models, such discontinuities have been called double layers (DLs). Recently similar discontinuities have been invoked to partly explain the auroral acceleration of electrons and ions in the upward current region. By means of one-dimensional Vlasov simulations of expanding plasmas, for the first time we make here the connection between (i) the rarefaction shocks, (ii) the discontinuities in the potential distributions, and (iii) DLs. We show that when plasmas expand from opposite directions into a deep density cavity with a potential drop across it and when the plasma on the high-potential side contains two electron populations, the temporal evolution of the potential and the plasma. distribution generates evolving multiple double layers with an extended density cavity between them. One of the DLs is the rarefaction-shock (RFS) and it forms by the reflections of the cold electrons coming from the high-potential side; it supports a part of the potential drop approximately determined by the hot electron temperature.

Singh, Nagendra↗
DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1990

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗
NASA NTRS1993

Low earth simulation and materials characterization

Oxygen plasma ashers and an electron cyclotron resonance (ECR) sources are currently being used for low Earth orbit (LEO) simulation. The suitability of each of these simulation techniques is considered. Thin film coatings are characterized by optical techniques, including variable-angle spectroscopic ellipsometry, optical spectrophotometry, and laser light scatterometry. Atomic force microscopy (AFM) has been used to characterize the surface morphology of thin aluminum films as a function of substrate temperature during deposition. Results on diamondlike carbon (DLC) films show that DLC degrades with simulated atomic oxygen (AO) exposure at a rate comparable to Kapton polyimide. Since DLC is not as susceptible as Kapton to environmental factors such as moisture absorption, it could potentially provide more accurate measurements of AO fluence on short space flights.

Synowicki, R. A.↗
NASA NTRS2003

Double Layers in Expanding Plasmas and Their Relevance to the Auroral Plasma Processes

When a dense plasma consisting of a cold and a sufficiently warm electron population expands, a rarefaction shock forms [Bezzerides et al., 1978]. In the expansion of the polar wind in the magnetosphere, it has been previously shown that when a sufficiently warm electron population also exists, in addition to the usual cold ionospheric one, a discontinuity forms in the electrostatic potential distribution along the magnetic field lines [Barakat and Schunk, 1984]. Despite the lack of spatial resolution and the assumption of quasi-neutrality in the polar wind models, such discontinuities have been called double layers (DLs). Recently similar discontinuities have been invoked to partly explain the auroral acceleration of electrons and ions in the upward current region [Ergun et al., 2000]. By means of one-dimensional Vlasov simulations of expanding plasmas, for the first time we make here the connection between (1) the rarefaction shocks, (2) the discontinuities in the potential distributions, and (3) DLs. We show that when plasmas expand from opposite directions into a deep density cavity with a potential drop across it and when the plasma on the high-potential side contains hot and cold electron populations, the temporal evolution of the potential and the plasma distribution generates evolving multiple double layers with an ,extended density cavity between them. One of the DLs is the rarefaction-shock (RFS) and it forms by the reflections of the cold electrons coming from the high-potential side; it supports a part of the potential drop approximately determined by the hot electron temperature. The other DLs evolve from charge separations arising either from reflection of ions coming from the low-potential side or stemming from plasma instabilities; they support the rest of the potential drop. The instabilities forming these additional double layers involve electron-ion (e-i) Buneman or ion-ion (i-i) two-stream interactions. The electron-electron two-stream interactions on the high-potential side of the RFS generate electron-acoustic waves, which evolve into electron phase-space holes. The ion population originating from the low-potential side and trapped by the RFS is energized by the e-i and i-i instabilities and it eventually precipitates into the high-potential plasma along with an electron beam. Applications of these findings to the auroral plasma physics are discussed.

Singh, Nagendra↗
DOE OSTI2015

Simulations of radiation damage as a function of the temporal pulse profile in femtosecond X-ray protein crystallography

Serial femtosecond X-ray crystallography of protein nanocrystals using ultrashort and intense pulses from an X-ray free-electron laser has proved to be a successful method for structural determination. However, due to significant variations in diffraction pattern quality from pulse to pulse only a fraction of the collected frames can be used. Experimentally, the X-ray temporal pulse profile is not known and can vary with every shot. This simulation study describes how the pulse shape affects the damage dynamics, which ultimately affects the biological interpretation of electron density. The instantaneously detected signal varies during the pulse exposure due to the pulse properties, as well as the structural and electronic changes in the sample. Here ionization and atomic motion are simulated using a radiation transfer plasma code. Pulses with parameters typical for X-ray free-electron lasers are considered: pulse energies ranging from 10 4 to 10 7 J cm -2 with photon energies from 2 to 12 keV, up to 100 fs long. Radiation damage in the form of sample heating that will lead to a loss of crystalline periodicity and changes in scattering factor due to electronic reconfigurations of ionized atoms are considered here. The simulations show differences in the dynamics of the radiation damage processes for different temporal pulse profiles and intensities, where ionization or atomic motion could be predominant. Overall, the different dynamics influence the recorded diffracted signal in any given resolution and will affect the subsequent structure determination.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗