Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Analyzing large data sets from XGC1 magnetic fusion simulations using apache spark

Apache Spark is explored as a tool for analyzing large data sets from the magnetic fusion simulation code XGCI. Implementation details of Apache Spark on the NERSC Edison supercomputer are discussed, including binary file reading, and parameter setup. Here, an unsupervised machine learning algorithm, k-means clustering, is applied to XGCI particle distribution function data, showing that highly turbulent spatial regions do not have common coherent structures, but rather broad, ring-like structures in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Dynamic evolutions of multiple toroidal Alfvén eigenmodes with energetic particles

Kinetic simulations based on a reduced model are carried out to study dynamic evolutions of multiple toroidicity-induced shear Alfvén eigenmodes driven by energetic particles. It is found that passing and trapped particles play quite different roles on the mode evolution. For cases without background damping, passing particles lead to a longer linear growth phase and subsequent enhancement of the nonlinear saturation level when resonant regions with two modes are closely overlapped. On the other hand, trapped particles are mainly responsible for the persistent mode growth in the nonlinear phase. By comparing the results from anisotropic and isotropic simulations, it is suggested that passing particles have a negative impact on the mode development in the nonlinear phase. For the near marginal stability cases with large background damping, the upward and downward frequency chirping in multiple branches are significantly affected by the multi-mode interaction due to overlapping of hole and clump in the distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Parallel implementation of three-dimensional molecular dynamic simulation for laser-cluster interaction

The objective of this article is to report the parallel implementation of the 3D molecular dynamic simulation code for laser-cluster interactions. The benchmarking of the code has been done by comparing the simulation results with some of the experiments reported in the literature. Scaling laws for the computational time is established by varying the number of processor cores and number of macroparticles used. The capabilities of the code are highlighted by implementing various diagnostic tools. To study the dynamics of the laser-cluster interactions, the executable version of the code is available from the author.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Energetic protons from an ultraintense laser interacting with a symmetric parabolic concave target

A scheme of a symmetric parabolic concave target irradiated by an ultraintense laser for efficient proton acceleration is proposed and involved problem is studied by using two-dimensional particle-in-cell (PIC) simulations. Results indicate that on one hand, the laser field is focused by the front parabolic concave surface of target and, on the other hand, more energetic hot electrons will traverse to the rear surface of target due to concave shape. The space-charge-separation field, induced by those hot electrons escaping form parabolic concave rear surface of target, can accelerate protons to relatively high energy with narrow energy spread. The dependence of the efficiency of proton acceleration on the target parameters is examined, and the optimal target parameters are obtained. Particle-in-cell simulations show that the proton peak energy and energy spread are greatly enhanced when the target parameters are chosen optimal, for example, a proton bunch with the maximum energy {approx}27.5 MeV and energy spread {approx}7% can be generated. Some implications of our results to experiments and comparisons with the other works are also discussed briefly.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

On the self-organizing process of large scale shear flows

Self organization is invoked as a paradigm to explore the processes governing the evolution of shear flows. By examining the probability density function (PDF) of the local flow gradient (shear), we show that shear flows reach a quasi-equilibrium state as its growth of shear is balanced by shear relaxation. Specifically, the PDFs of the local shear are calculated numerically and analytically in reduced 1D and 0D models, where the PDFs are shown to converge to a bimodal distribution in the case of finite correlated temporal forcing. This bimodal PDF is then shown to be reproduced in nonlinear simulation of 2D hydrodynamic turbulence. Furthermore, the bimodal PDF is demonstrated to result from a self-organizing shear flow with linear profile. Similar bimodal structure and linear profile of the shear flow are observed in gulf stream, suggesting self-organization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Tail-ion transport and Knudsen layer formation in the presence of magnetic fields

Knudsen layer losses of tail fuel ions could reduce significantly the fusion reactivity of highly compressed cylindrical and spherical targets in inertial confinement fusion (ICF). With the class of magnetized ICF targets in mind, the effect of embedded magnetic fields on Knudsen layer formation is investigated for the first time. The modified energy scaling of ion diffusivity in magnetized hot spots is found to suppress the preferential losses of tail-ions perpendicular to the magnetic field lines to a degree that the tail distribution can be at least partially, if not fully, restored. Two simple threshold conditions are identified leading to the restoration of fusion reactivity in magnetized hot spots. A kinetic equation for tail-ion transport in the presence of a magnetic field is derived, and solutions to the equation are obtained numerically in simulations. Numerical results confirm the validity of the threshold conditions for restored reactivity and identify two different asymptotic regimes of the fusion fuel. While Knudsen layer formation is shown to be suppressed entirely in strongly magnetized cylindrical hot spot cavities, uniformly magnetized spherical cavities demonstrate remnant, albeit reduced, levels of tail-ion depletion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Performance and Thermal Characterization of the NASA-300MS 20 kW Hall Effect Thruster

NASA's Space Technology Mission Directorate is sponsoring the development of a high fidelity 15 kW-class long-life high performance Hall thruster for candidate NASA technology demonstration missions. An essential element of the development process is demonstration that incorporation of magnetic shielding on a 20 kW-class Hall thruster will yield significant improvements in the throughput capability of the thruster without any significant reduction in thruster performance. As such, NASA Glenn Research Center and the Jet Propulsion Laboratory collaborated on modifying the NASA-300M 20 kW Hall thruster to improve its propellant throughput capability. JPL and NASA Glenn researchers performed plasma numerical simulations with JPL's Hall2De and a commercially available magnetic modeling code that indicated significant enhancement in the throughput capability of the NASA-300M can be attained by modifying the thruster's magnetic circuit. This led to modifying the NASA-300M magnetic topology to a magnetically shielded topology. This paper presents performance evaluation results of the two NASA-300M magnetically shielded thruster configurations, designated 300MS and 300MS-2. The 300MS and 300MS-2 were operated at power levels between 2.5 and 20 kW at discharge voltages between 200 and 700 V. Discharge channel deposition from back-sputtered facility wall flux, and plasma potential and electron temperature measurements made on the inner and outer discharge channel surfaces confirmed that magnetic shielding was achieved. Peak total thrust efficiency of 64% and total specific impulse of 3,050 sec were demonstrated with the 300MS-2 at 20 kW. Thermal characterization results indicate that the boron nitride discharge chamber walls temperatures are approximately 100 C lower for the 300MS when compared to the NASA- 300M at the same thruster operating discharge power.

High power↗
DOE OSTI2018

Gyrokinetic simulation study of magnetic island effects on neoclassical physics and micro-instabilities in a realistic KSTAR plasma

Here, we perform gyrokinetic simulations to study the effects of a stationary magnetic island on neoclassical flow and micro-instability in a realistic KSTAR plasma condition. Through the simulations, we aim to analyze a recent KSTAR experiment, which was to measure the details of poloidal flow and fluctuation around a stationary (2, 1) magnetic island [M. J. Choi et al., Nucl. Fusion 57, 126058 (2017)]. From the simulations, it is found that the magnetic island can significantly enhance the equilibrium E x B flow. The corresponding flow shearing is strong enough to suppress a substantial portion of ambient micro-instabilities, particularly ∇T e -driven trapped electron modes. This implies that the enhanced E x B flow can sustain a quasi-internal transport barrier for T e in an inner region neighboring the magnetic island. The enhanced E x B flow has a (2, 1) mode structure with a finite phase shift from the mode structure of the magnetic island. It is shown that the flow shear and the fluctuation suppression patterns implied from the simulations are consistent with the observations on the KSTAR experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The transfer between electron bulk kinetic energy and thermal energy in collisionless magnetic reconnection

By performing two-dimensional particle-in-cell simulations, we investigate the transfer between electron bulk kinetic and electron thermal energy in collisionless magnetic reconnection. In the vicinity of the X line, the electron bulk kinetic energy density is much larger than the electron thermal energy density. The evolution of the electron bulk kinetic energy is mainly determined by the work done by the electric field force and electron pressure gradient force. The work done by the electron gradient pressure force in the vicinity of the X line is changed to the electron enthalpy flux. In the magnetic island, the electron enthalpy flux is transferred to the electron thermal energy due to the compressibility of the plasma in the magnetic island. The compression of the plasma in the magnetic island is the consequence of the electromagnetic force acting on the plasma as the magnetic field lines release their tension after being reconnected. Therefore, we can observe that in the magnetic island the electron thermal energy density is much larger than the electron bulk kinetic energy density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Enhance the terahertz Smith-Purcell superradiant radiation by using dielectric loaded grating

A dielectric loaded grating (DLG) for terahertz Smith-Purcell (SP) device is proposed to enhance the radiation intensity. By using the theoretical analysis and particle-in-cell simulations, the dispersion characteristics and SP superradiant radiation are investigated. Compared with the general metal grating, the usage of DLG can improve the magnitude of electric field and, consequently, strengthen the interaction of the evanescent wave with electron beam, which can improve the growth rate, enhance the SP superradiant radiation, and lower the start current for the operation of SP free-electron laser.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Particle-in-cell simulations on spontaneous thermal magnetic field fluctuations

In this paper an electromagnetic particle code is used to investigate the spontaneous thermal emission. Specifically we perform particle-in-cell simulations employing a non-relativistic isotropic Maxwellian particle distribution to show that thermal fluctuations are related to the origin of spontaneous magnetic field fluctuation. These thermal fluctuations can become seed for further amplification mechanisms and thus be considered at the origin of the cosmological magnetic field, at microgauss levels. Our numerical results are in accordance with theoretical results presented in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ion acceleration in shell cylinders irradiated by a short intense laser pulse

The interaction of a short high intensity laser pulse with homo and heterogeneous shell cylinders has been analyzed using particle-in-cell simulations and analytical modeling. We show that the shell cylinder is proficient of accelerating and focusing ions in a narrow region. In the case of shell cylinder, the ion energy exceeds the ion energy for a flat target of the same thickness. The constructed model enables the evaluation of the ion energy and the number of ions in the focusing region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Linear growth rates of resistive tearing modes with sub-Alfvénic streaming flow

The tearing instability with sub-Alfvénic streaming flow along the external magnetic field is investigated using resistive MHD simulation. It is found that the growth rate of the tearing mode instability is larger than that without the streaming flow. With the streaming flow, there exist two Alfvén resonance layers near the central current sheet. The larger perturbation of the magnetic field in two closer Alfvén resonance layers could lead to formation of the observed cone structure and can largely enhance the development of the tearing mode for a narrower streaming flow. For a broader streaming flow, a larger separation of Alfvén resonance layers reduces the magnetic reconnection. The linear growth rate decreases with increase of the streaming flow thickness. The growth rate of the tearing instability also depends on the plasma beta (β). When the streaming flow is embedded in the current sheet, the growth rate increases with β if β < β{sub s}, but decreases if β > β{sub s}. The existence of the specific value β{sub s} can be attributed to competition between the suppressing effect of β and the enhancing effect of the streaming flow on the magnetic reconnection. The critical value β{sub s} increases with increase of the streaming flow strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Avalanche transport of energetic-ions in magnetic confinement plasmas: nonlinear multiple wave-number simulation

Large burst activity, identified as toroidal Alfvén eigenmode (TAE) avalanche, occurs frequently in neutral-beam heated plasmas in National Spherical Torus Experiment (NSTX). Based on the typical experimental observation of TAE avalanche on NSTX, a self-consistent nonlinear multiple wave-number (k ∥ ≃ n/R, where n toroidal mode-number and R major radius) simulation associated with TAE avalanches is performed using the experimental parameters and profiles before the occurrence of TAE avalanche as the M3D-K input. The wave–wave nonlinear coupling among different modes and the resonant interaction between different modes and energetic-ions during TAE avalanches are identified in the nonlinear multiple wave-number simulations. The resonance overlap during the TAE avalanche is clearly observed in the simulation. It is found that the effective wave–wave coupling and a sufficiently strong drive are two important ingredients for the onset of TAE avalanches. TAE avalanche is considered to be a strongly nonlinear process and it is always accompanied by the simultaneous rapid frequency-chirping and large amplitude bursting of multiple modes and significant energetic-ion losses. The experimental phenomenon is observed on NSTX and is qualitatively reproduced by the simulation results in this work. These findings indicate that the onset of avalanche is triggered by nonlinearity of the system, and are also conducive to understanding the underlying mechanism of avalanche transport of energetic particles in the future burning plasmas, such as International Thermonuclear Experiment Reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Quasimonoenergetic proton bunches generation from doped foil targets irradiated by intense lasers

We propose a scheme to generate 10 MeV-level quasimonoenergetic proton bunches using proton-doped heavy-ion targets irradiated by intense lasers via target normal sheath acceleration. The heavy substrate ions provide a long-life quasi-stable sheath field to accelerate the doped protons at the target rear and consequently a quasimonoenergetic proton bunch is produced. The scheme is demonstrated by two-dimensional particle-in-cell simulations. An exemplificative simulation with parameters of targets made by ion-implant technique, a kind of modern doping process, gives a quasimonoenergetic bunch with peak energy {approx}13MeV, energy spread {approx}24%, and {approx}nC charge at the focused laser intensity 10{sup 20}W/cm{sup 2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗