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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma

In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – driven by ionization, charge exchange and wall loss – is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Non-thermal electron acceleration in low Mach number collisionless shocks. II. Firehose-mediated Fermi acceleration and its dependence on pre-shock conditions

Electron acceleration to non-thermal energies is known to occur in low Mach number (M{sub s} ≲ 5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M{sub s} = 3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (≳ 20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. I. THEORY AND HYBRID SIMULATIONS

The nearly circular band of energetic neutral atom emission dominating the field of view of the Interplanetary Boundary Explorer ( IBEX ) satellite, is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath and the interstellar space beyond. Several models for the PUI dynamics of this mechanism have been proposed, each requiring either strong or weak scattering of the initial pitch angle. Conventional wisdom states that ring distributions tend to generate waves and scatter onto a shell on timescales too short for charge exchange to occur. We performed a careful study of ring and thin shell proton distribution stability using theoretical tools and hybrid plasma simulations. We show that the kinetic behavior of a freshly injected proton ring is a far more complicated process than previously thought. In the presence of a warm Maxwellian core, narrower rings could be more stable than broader toroidal distributions. The scattered rings possess a fine structure that can only be revealed using very large numbers of macroparticles in a simulation. It is demonstrated that a “stability gap” in ring temperature exists where the protons could retain large gyrating anisotropies for years, and the wave activity could remain below the level of the ambient magnetic fluctuations in interstellar space. In the directions away from the ribbon, however, a partial shell distribution is more likely to be unstable, leading to significant scattering into one hemisphere in velocity space. The process is accompanied by turbulence production, which is puzzling given the very low level of magnetic fluctuations measured in the outer heliosheath by Voyager 1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

FOCUS: A full-orbit CUDA solver for particle simulations in magnetized plasmas

A Full-Orbit CUDA Solver for particle simulations in plasmas, FOCUS, has been developed. The code follows exact particle trajectories by solving Newton’s equations with the Lorentz force. The code can use fields calculated analytically or numerically and it is also possible to use the information provided by equilibrium reconstruction and transport codes. FOCUS has an elastic collisions module which covers the whole particle energy range in magnetic fusion devices. Moreover, an atomic collision module was also included to simulate the interaction of the test particles with neutral or partially ionized species. Regarding performance, the main feature is that the code runs on GPUs, allowing the simulation of a large number of particles using moderate computational resources. Different versions of the code have been used in several papers. Here we present a complete description of its capabilities, the basis physics included and detailed information about the numerical algorithms employed. Finally, some examples are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Progress with the COGENT Edge Kinetic Code: Implementing the Fokker-Plank Collision Operator

Here, COGENT is a continuum gyrokinetic code for edge plasma simulations being developed by the Edge Simulation Laboratory collaboration. The code is distinguished by application of a fourth-order finite-volume (conservative) discretization, and mapped multiblock grid technology to handle the geometric complexity of the tokamak edge. The distribution function F is discretized in v∥ – μ (parallel velocity – magnetic moment) velocity coordinates, and the code presently solves an axisymmetric full-f gyro-kinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. COGENT capabilities are extended by implementing the fully nonlinear Fokker-Plank operator to model Coulomb collisions in magnetized edge plasmas. The corresponding Rosenbluth potentials are computed by making use of a finite-difference scheme and multipole-expansion boundary conditions. Details of the numerical algorithms and results of the initial verification studies are discussed. (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Full-wave feasibility study of magnetic diagnostic based on O-X mode conversion and oblique reflectometry imaging

An innovative millimeter wave diagnostic is proposed to measure the local magnetic field and the edge current as a function of the minor radius in the pedestal region. The idea behind such diagnostic is to localize and characterize a direction of reduced reflectivity at the O-mode cutoff layer. We modeled the wave scattering and mode-conversion processes by means of the finite-element COMSOL Multiphysics code in two dimensions (2D). Sensitivity studies were performed for parameters mocking up DIII-D plasmas. Simulations confirmed the presence of a minimum in reflectivity of an externally injected O-mode beam, and confirmed that this minimum depends on the magnetic field at the cutoff, as expected from the OX mode conversion physics. This study gives confidence in the feasibility of the diagnostic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

EXPANDED IRON UTA SPECTRA-PROBING THE THERMAL STABILITY LIMITS IN AGN CLOUDS

The Fe unresolved transition arrays (UTAs) produce prominent features in the {approx}15-17 A wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in AGNs. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Preliminary design of a tangentially viewing imaging bolometer for NSTX-U

The infrared imaging video bolometer (IRVB) measures plasma radiated power images using a thin metal foil. Two different designs with a tangential view of NSTX-U are made assuming a 640 × 480 (1280 × 1024) pixel, 30 (105) fps, 50 (20) mK, IR camera imaging the 9 cm × 9 cm × 2 μm Pt foil. The foil is divided into 40 × 40 (64 × 64) IRVB channels. This gives a spatial resolution of 3.4 (2.2) cm on the machine mid-plane. The noise equivalent power density of the IRVB is given as 113 (46) μW/cm{sup 2} for a time resolution of 33 (20) ms. Synthetic images derived from Scrape Off Layer Plasma Simulation data using the IRVB geometry show peak signal levels ranging from ∼0.8 to ∼80 (∼0.36 to ∼26) mW/cm{sup 2}.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Plasma-wall interaction in Hall thrusters with magnetic lens configuration

Some recently developed Hall thrusters utilize a magnetic field configuration in which the field lines penetrate the thruster walls at a high incidence angle. This so-called magnetic lens leads to an electric field pointing away from the walls, which is expected to reduce ion losses and improve thruster efficiency. This configuration also introduces an interesting behavior in the sheath formation. At sufficiently large angles, ions are repelled from the wall, and sheath collapse is expected. We use a plasma simulation code to investigate this phenomenon in detail. We consider the role of the magnetic field incidence angle, secondary electron emission, and a magnetic mirror. Numerical study confirms the theoretical predictions, and at large angles, ions are seen to turn away from the wall. We also consider the role of the magnetic field geometry on ion wall flux and channel erosion, and observe reduction in both quantities as the magnetic field incidence angle is increased.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Coupled radiation transport and plasma physics for ICF simulations [Slides]

Describing high-temperature plasmas requires inclusion of several nonlinear coupled physics, including ion kinetics - for each isotope species present in the plasma, electron kinetics, thermal radiative transfer (TRT), electromagnetic description, and nuclear reactions (fusion burn). Simulating plasmas in ICF hohlraums aims to solve kinetic ions (particles), fluid electrons, energetic electrons, electromagnetics, kinetic radiation (MG DP solver), and laser-plasma interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Koopman wavefunctions and Clebsch variables in Vlasov–Maxwell kinetic theory

Motivated by recent discussions on the possible role of quantum computation in plasma simulations, here, we present different approaches to Koopman's Hilbert-space formulation of classical mechanics in the context of Vlasov–Maxwell kinetic theory. The celebrated Koopman–von Neumann construction is provided with two different Hamiltonian structures: one is canonical and recovers the usual Clebsch representation of the Vlasov density, the other is non-canonical and appears to overcome certain issues emerging in the canonical formalism. Furthermore, the canonical structure is restored for a variant of the Koopman–von Neumann construction that carries a different phase dynamics. Going back to van Hove's prequantum theory, the corresponding Koopman–van Hove equation provides an alternative Clebsch representation which is then coupled to the electromagnetic fields. Finally, the role of gauge transformations in the new context is discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Indeterminacy and instability in Petschek reconnection

We explain two puzzling aspects of Petschek's model for fast reconnection. One is its failure to occur in plasma simulations with uniform resistivity. The other is its inability to provide anything more than an upper limit for the reconnection rate. We have found that previously published analytical solutions based on Petschek's model are structurally unstable if the electrical resistivity is uniform. The structural instability is associated with the presence of an essential singularity at the X-line that is unphysical. By requiring that such a singularity does not exist, we obtain a formula that predicts a specific rate of reconnection. For uniform resistivity, reconnection can only occur at the slow, Sweet-Parker rate. For nonuniform resistivity, reconnection can occur at a much faster rate provided that the resistivity profile is not too flat near the X-line. If this condition is satisfied, then the scale length of the nonuniformity determines the reconnection rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Fast-ignition transport studies: Realistic electron source, integrated particle-in-cell and hydrodynamic modeling, imposed magnetic fields

Transport modeling of idealized, cone-guided fast ignition targets indicates the severe challenge posed by fast-electron source divergence. The hybrid particle-in-cell (PIC) code Zuma is run in tandem with the radiation-hydrodynamics code Hydra to model fast-electron propagation, fuel heating, and thermonuclear burn. The fast electron source is based on a 3D explicit-PIC laser-plasma simulation with the PSC code. This shows a quasi two-temperature energy spectrum and a divergent angle spectrum (average velocity-space polar angle of 52 Degree-Sign ). Transport simulations with the PIC-based divergence do not ignite for >1 MJ of fast-electron energy, for a modest (70 {mu}m) standoff distance from fast-electron injection to the dense fuel. However, artificially collimating the source gives an ignition energy of 132 kJ. To mitigate the divergence, we consider imposed axial magnetic fields. Uniform fields {approx}50 MG are sufficient to recover the artificially collimated ignition energy. Experiments at the Omega laser facility have generated fields of this magnitude by imploding a capsule in seed fields of 50-100 kG. Such imploded fields will likely be more compressed in the transport region than in the laser absorption region. When fast electrons encounter increasing field strength, magnetic mirroring can reflect a substantial fraction of them and reduce coupling to the fuel. A hollow magnetic pipe, which peaks at a finite radius, is presented as one field configuration which circumvents mirroring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Non-thermal electron acceleration in low Mach number collisionless shocks. I. Particle energy spectra and acceleration mechanism

Electron acceleration to non-thermal energies in low Mach number (M{sub s} ≲ 5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M{sub s} = 3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15% of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p ≅ 2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration (SDA). The accelerated electrons are then reflected back upstream where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock for further energization via SDA. In summary, the self-generated waves allow for repeated cycles of SDA, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

79 ASTRONOMY AND ASTROPHYSICS↗