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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The effects of dilution on turbulence and transport in C-Mod ohmic plasmas and comparisons with gyrokinetic simulations

Main ion dilution has been predicted by gyrokinetic simulations to have a significant effect on ion thermal transport in C-Mod ohmic plasmas. This effect was verified experimentally with a specific set of experiments on C-Mod in which ohmic deuterium plasmas across the linear ohmic confinement (LOC) through the saturated ohmic confinement (SOC) regimes were diluted by seeding with nitrogen gas (Z=7) injection. The seeding was observed to increase the normalized ion temperature gradients (ITGs) by up to 30% without a corresponding increase in the gyrobohm normalized ion energy flux, indicating a change in either the stiffness or the critical ion temperature gradient associated with ITG turbulence. The seeding also reversed the direction of the intrinsic toroidal rotation in plasmas slightly above the normal intrinsic rotation reversal critical density. GYRO simulations of the seeded and unseeded plasmas show that the seeding affected both the critical gradient and the stiffness. For plasmas in the LOC regime, the dilution primarily increased the critical gradient, while for plasmas in the SOC regime the dilution primarily decreased the stiffness. At r/a = 0.8, where the experimental fluxes were above marginal stability, local GYRO predicted and experimental energy fluxes agreed, except for Q i in the SOC regime where GYRO under-predicted the experimental energy flux. At r/a = 0.6, where the experimental fluxes were close to marginally stable, local GYRO predicted ITG modes to be strongly unstable and are responsible for both Q i and Q e (with Q i > Q e ), as opposed to the experiment where Q i < Q e . In contrast, global GYRO in this region predicted the ITG modes to be closer to marginal stability, and accurately predict the experimental Q i when the T i profile is modified within experimental uncertainties. Further, the fact that Q e is always less than Q i in the r/a = 0.6 simulations with k θ ρ s ≤1 indicates that high-k electron temperature gradient driven (ETG) modes must be included in future simulations and may be responsible for the electron energy transport in this case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Observations and Simulations of Formation of Broad Plasma Depletions Through Merging Process

Broad plasma depletions in the equatorial ionosphere near dawn are region in which the plasma density is reduced by 1-3 orders of magnitude over thousands of kilometers in longitude. This phenomenon is observed repeatedly by the Communication/Navigation Outage Forecasting System (C/NOFS) satellite during deep solar minimum. The plasma flow inside the depletion region can be strongly upward. The possible causal mechanism for the formation of broad plasma depletions is that the broad depletions result from merging of multiple equatorial plasma bubbles. The purpose of this study is to demonstrate the feasibility of the merging mechanism with new observations and simulations. We present C/NOFS observations for two cases. A series of plasma bubbles is first detected by C/NOFS over a longitudinal range of 3300-3800 km around midnight. Each of the individual bubbles has a typical width of approx 100 km in longitude, and the upward ion drift velocity inside the bubbles is 200-400 m/s. The plasma bubbles rotate with the Earth to the dawn sector and become broad plasma depletions. The observations clearly show the evolution from multiple plasma bubbles to broad depletions. Large upward plasma flow occurs inside the depletion region over 3800 km in longitude and exists for approx 5 h. We also present the numerical simulations of bubble merging with the physics-based low-latitude ionospheric model. It is found that two separate plasma bubbles join together and form a single, wider bubble. The simulations show that the merging process of plasma bubbles can indeed occur in incompressible ionospheric plasma. The simulation results support the merging mechanism for the formation of broad plasma depletions.

plasma bubbles↗
DOE OSTI2015

Trapping of intense light in hollow shell

A small hollow shell for trapping laser light is proposed. Two-dimensional particle-in-cell simulation shows that under appropriate laser and plasma conditions a part of the radiation fields of an intense short laser pulse can enter the cavity of a small shell through an over-critical density plasma in an adjacent guide channel and become trapped. The trapped light evolves into a circulating radial wave pattern until its energy is dissipated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification of gyrokinetic particle simulations of neoclassical tearing modes in fusion plasmas

Abstract The ability to simulate neoclassical tearing modes (NTMs) in the gyrokinetic toroidal code (GTC) has been developed and verified, in which ions are treated with a gyrokinetic model and electrons are treated as a resistive fluid. The simulation results demonstrate that the neoclassical bootstrap current effect can destabilize an otherwise stable classical tearing mode. In the cylindrical geometry, GTC simulations in the magnetohydrodynamic limit show quantitative agreement with the modified Rutherford theory, both in terms of the scaling law in the small island limit and in terms of the saturation level and pressure flattening effect in the large island limit. The toroidal effects are slightly destabilizing for the NTM, while the kinetic effects of thermal ions are stabilizing for the NTM and increase its excitation threshold.

Physics↗
DOE OSTI2020

Interaction of ultra relativistic ${e}^{-}{e}^{+}$ fireball beam with plasma

Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron–positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density~1017 cm-3 . The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Developing the science and technology for the Material Plasma Exposure eXperiment

Linear plasma generators are cost effective facilities to simulate divertor plasma conditions of present and future fusion reactors. They are used to address important R&D gaps in the science of plasma material interactions and towards viable plasma facing components for fusion reactors. Next generation plasma generators have to be able to access the plasma conditions expected on the divertor targets in ITER and future devices. The steady-state linear plasma device MPEX will address this regime with electron temperatures of 1–10 eV and electron densities of $10^{21}{\text{}}\!-\!10^{20}$ ${\rm m}^{-3}$. The resulting heat fluxes are about 10 MW ${\rm m}^{-2}$ . MPEX is designed to deliver those plasma conditions with a novel Radio Frequency plasma source able to produce high density plasmas and heat electron and ions separately with electron Bernstein wave (EBW) heating and ion cyclotron resonance heating with a total installed power of 800 kW. The linear device Proto-MPEX, forerunner of MPEX consisting of 12 water-cooled copper coils, has been operational since May 2014. Its helicon antenna (100 kW, 13.56 MHz) and EC heating systems (200 kW, 28 GHz) have been commissioned and 14 MW ${\rm m}^{-2}$ was delivered on target. Furthermore, electron temperatures of about 20 eV have been achieved in combined helicon and ECH heating schemes at low electron densities. Overdense heating with EBW was achieved at low heating powers. The operational space of the density production by the helicon antenna was pushed up to $1.1 \times 10^{20}$ ${\rm m}^{-3}$ at high magnetic fields of 1.0 T at the target. Finally, the experimental results from Proto-MPEX will be used for code validation to enable predictions of the source and heating performance for MPEX. MPEX, in its last phase, will be capable to expose neutron-irradiated samples. In this concept, targets will be irradiated in ORNL's High Flux Isotope Reactor and then subsequently exposed to fusion reactor relevant plasmas in MPEX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Real time animation of space plasma phenomena

In pursuit of real time animation of computer simulated space plasma phenomena, the code was rewritten for the Massively Parallel Processor (MPP). The program creates a dynamic representation of the global bowshock which is based on actual spacecraft data and designed for three dimensional graphic output. This output consists of time slice sequences which make up the frames of the animation. With the MPP, 16384, 512 or 4 frames can be calculated simultaneously depending upon which characteristic is being computed. The run time was greatly reduced which promotes the rapid sequence of images and makes real time animation a foreseeable goal. The addition of more complex phenomenology in the constructed computer images is now possible and work proceeds to generate these images.

Jordan, K. F.↗
DOE OSTI2015

WarpVisit

WarpVisit is an insitu simulation application that integrates the Warp laser plasma accelerator simulation framework with Visit a parallel visualization application. WarpVisit is written in python and supports interactive or live mode where user can connect to Warp with the Visit GUI and batch mode for batch for non-interactive use on high-performance computing resources.

Loring, Burlen↗
DOE OSTI2014

ADVECTIVE TRANSPORT OF INTERSTELLAR PLASMA INTO THE HELIOSPHERE ACROSS THE RECONNECTING HELIOPAUSE

We discuss results of magnetohydrodynamical model simulations of plasma dynamics in the proximity of the heliopause (HP). The model is shown to fit details of the magnetic field variations observed by the Voyager 1 spacecraft during the transition from the heliosphere to the local interstellar medium (LISM). We propose an interpretation of magnetic field structures observed by Voyager 1 in terms of fine-scale physical processes. Our simulations reveal an effective transport mechanism of relatively dense LISM plasma across the reconnecting HP into the heliosphere. The mechanism is associated with annihilation of magnetic sectors in the heliospheric plasma near the HP.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Reproducing continuous radio blackout using glow discharge plasma

A novel plasma generator is described that offers large-scale, continuous, non-magnetized plasma with a 30-cm-diameter hollow structure, which provides a path for an electromagnetic wave. The plasma is excited by a low-pressure glow discharge, with varying electron densities ranging from 10{sup 9} to 2.5 × 10{sup 11} cm{sup −3}. An electromagnetic wave propagation experiment reproduced a continuous radio blackout in UHF-, L-, and S-bands. The results are consistent with theoretical expectations. The proposed method is suitable in simulating a plasma sheath, and in researching communications, navigation, electromagnetic mitigations, and antenna compensation in plasma sheaths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Improved fluid simulations of radio-frequency plasmas using energy dependent ion mobilities

Symmetric and asymmetric capacitively coupled radio-frequency plasmas in oxygen at 40 Pa, 300 V voltage amplitude and a discharge gap of 40 mm are investigated by means of one-dimensional numerical semi-kinetic fluid modeling on the basis of a simplified reaction scheme including the dominant positive and negative ions, background gas, and electrons. An improved treatment, by accounting for the dependence of ion mobilities on E/N, is compared to the standard approach, based on using zero-field mobility values only. The charged particle dynamics as a result of direct electron impact ionization of oxygen, secondary electron release from the electrodes, the spatial distribution of all involved particles as well as impact of geometry and model modification on ion energies is analyzed and compared to independent simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Magnetic field generation in core-sheath jets via the kinetic Kelvin-Helmholtz instability

We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas, we find generation of strong large-scale DC currents and magnetic fields that extend over the entire shear surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas, we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates a transverse structure similar to that produced by the Weibel instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Pedestal fueling simulations with a coupled kinetic plasma–kinetic neutral transport code

The simplified neutral transport routine in the guiding center ion–electron-neutral neoclassical PIC code XGC0 has been replaced by one that calls the Monte Carlo neutral transport code DEGAS2. This provides a more realistic treatment of neutral atoms in the tokamak edge plasma and allows detailed atomic physics and plasma–material interaction processes to be incorporated into the simulations. The spatial profile of the neutral particle source used in DEGAS2 is determined from the fluxes of XGC0 ions to the material surfaces, and DEGAS2 is run in a time dependent manner synchronized with XGC0. The kinetic plasma–kinetic neutral transport capability is demonstrated with example pedestal fueling simulations.

Materials Science↗
DOE OSTI2017

Quasi-optical simulation of the electron cyclotron plasma heating in a mirror magnetic trap

The resonance microwave plasma heating in a large-scale open magnetic trap is simulated taking into account all the basic wave effects during the propagation of short-wavelength wave beams (diffraction, dispersion, and aberration) within the framework of the consistent quasi-optical approximation of Maxwell’s equations. The quasi-optical method is generalized to the case of inhomogeneous media with absorption dispersion, a new form of the quasi-optical equation is obtained, the efficient method for numerical integration is found, and simulation results are verified on the GDT facility (Novosibirsk).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Development and Benchmarking of a Hybrid PIC Code For Dense Plasmas and Fast Ignition (Final Report)

Radiation processes play an important role in the study of both fast ignition and other inertial confinement schemes, such as plasma jet driven magneto-inertial fusion, both in their effect on energy balance, and in generating diagnostic signals. In the latter case, warm and hot dense matter may be produced by the convergence of a plasma shell formed by the merging of an assembly of high Mach number plasma jets. This innovative approach has the potential advantage of creating matter of high energy densities in voluminous amount compared with high power lasers or particle beams. An important application of this technology is as a plasma liner for the flux compression of magnetized plasma to create ultra-high magnetic fields and burning plasmas. HyperV Technologies Corp. has been developing plasma jet accelerator technology in both coaxial and linear railgun geometries to produce plasma jets of sufficient mass, density, and velocity to create such imploding plasma liners. An enabling tool for the development of this technology is the ability to model the plasma dynamics, not only in the accelerators themselves, but also in the resulting magnetized target plasma and within the merging/interacting plasma jets during transport to the target. Welch pioneered numerical modeling of such plasmas (including for fast ignition) using the LSP simulation code. Lsp is an electromagnetic, parallelized, plasma simulation code under development since 1995. It has a number of innovative features making it uniquely suitable for modeling high energy density plasmas including a hybrid fluid model for electrons that allows electrons in dense plasmas to be modeled with a kinetic or fluid treatment as appropriate. In addition to in-house use at Voss Scientific, several groups carrying out research in Fast Ignition (LLNL, SNL, UCSD, AWE (UK), and Imperial College (UK)) also use LSP. A collaborative team consisting of HyperV Technologies Corp., Voss Scientific LLC, FAR-TECH, Inc., Prism Computational Sciences, Inc. and Advanced Energy Systems Inc. joined efforts to develop new physics and numerical models for LSP in several key areas to enhance the ability of LSP to model high energy density plasmas (HEDP). This final report details those efforts. Areas addressed in this research effort include: adding radiation transport to LSP, first in 2D and then fully 3D, extending the EMHD model to 3D, implementing more advanced radiation and electrode plasma boundary conditions, and installing more efficient implicit numerical algorithms to speed complex 2-D and 3-D computations. The new capabilities allow modeling of the dominant processes in high energy density plasmas, and further assist the development and optimization of plasma jet accelerators, with particular attention to MHD instabilities and plasma/wall interaction (based on physical models for ion drag friction and ablation/erosion of the electrodes). In the first funding cycle we implemented a solver for the radiation diffusion equation. To solve this equation in 2-D, we used finite-differencing and applied the parallelized sparse-matrix solvers in the PETSc library (Argonne National Laboratory) to the resulting system of equations. A database of the necessary coefficients for materials of interest was assembled using the PROPACEOS and ATBASE codes from Prism. The model was benchmarked against Prism's 1-D radiation hydrodynamics code HELIOS, and against experimental data obtained from HyperV's separately funded plasma jet accelerator development program. Work in the second funding cycle focused on extending the radiation diffusion model to full 3-D, continued development of the EMHD model, optimizing the direct-implicit model to speed up calculations, add in multiply ionized atoms, and improved the way boundary conditions are handled in LSP. These new LSP capabilities were then used, along with analytic calculations and Mach2 runs, to investigate plasma jet merging, plasma detachment and transport, restrike and advanced jet accelerator design. In addition, a strong linkage to diagnostic measurements was made by modeling plasma jet experiments on PLX to support benchmarking of the code. A large number of upgrades and improvements advancing hybrid PIC algorithms were implemented in LSP during the second funding cycle. These include development of fully 3D radiation transport algorithms, new boundary conditions for plasma-electrode interactions, and a charge conserving equation of state that permits multiply ionized high-Z ions. The final funding cycle focused on 1) mitigating the effects of a slow-growing grid instability which is most pronounced in plasma jet frame expansion problems using the two-fluid Eulerian remap algorithm, 2) extension of the Eulerian Smoothing Algorithm to allow EOS/Radiation modeling, 3) simulations of collisionless shocks formed by jet merging, 4) simulations of merging jets using high-Z gases, 5) generation of PROPACEOS EOS/Opacity databases, 6) simulations of plasma jet transport experiments, 7) simulations of plasma jet penetration through transverse magnetic fields, and 8) GPU PIC code development. The tools developed during this project are applicable not only to the study of plasma jets, but also to a wide variety of HEDP plasmas of interest to DOE, including plasmas created in short-pulse laser experiments performed to study fast ignition concepts for inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗