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,045 records · Page 58

Applications of the SCENIC code package to the minority ion-cyclotron heating in Wendelstein 7-X plasmas

We present SCENIC simulations of a W7X 4He plasma with 1% H minority and with an antenna model close to the design foreseen for the W7X ICRF antenna [1, 2]. A high mirror and a standard equilibrium are considered. The injected wave frequency is fixed at 33.8 MHz and 39.6MHz respectively and only fundamental minority heating is considered. Included in this calculation is a new realistic model of the antenna, where it is found that the localization of the antenna geometry tends to break the five-fold periodicity of the system. We assess the heat transfer through the toroidal periods via Coulomb collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Behavior of lithium ions in the turbulent near-wall tokamak plasma under heating of ions and electrons of the main plasma

Turbulent dynamics of the near-wall tokamak plasma is simulated by numerically solving the nonlinear reduced Braginskii magnetohydrodynamic equations with allowance for a lithium ion admixture. The effects of turbulence and radiation of the admixture are analyzed in the framework of a self-consistent approach. The radial distributions of the radiative loss power and the density of Li{sup 0} atoms and Li{sup +1} ions are obtained as functions of the electron and ion temperatures of the main plasma in the near-wall layer. The results of numerical simulations show that supply of lithium ions into the low-temperature near-wall plasma substantially depends on whether the additional power is deposited into the electron or ion component of the main plasma. If the electron temperature in the layer increases (ECR heating), then the ion density drops. At the same time, an increase in the temperature of the main ions (ICR heating) leads to an increase in the density of Li{sup +1} ions. The results of numerical simulations are explained by the different influence of the electron and ion temperatures on the atomic processes governing the accumulation and loss of particles in the balance equations for neutral Li{sup 0} atoms and Li{sup +1} ions in the admixture. The radial profile of the electron temperature and the corresponding distribution of the radiative loss power for different densities of neutral Li{sup 0} atoms on the wall are obtained. The calculations show that the presence of Li{sup +1} ions affects turbulent transport of the main ions. In this case, the electron heat flux increases by 20–30% with increasing Li{sup +1} density, whereas the flux of the main ions drops by nearly the same amount. The radial profile of the turbulent flux of lithium ions is obtained. It is demonstrated that the appearance of the pinch effect is related to the positive density gradient of lithium ions across the calculation layer. For the parameters of the T-10 tokamak, the effect of radiative cooling of the near-wall plasma layer becomes appreciable when the near-wall density of neutral lithium atoms exceeds 7 × 10{sup 11} cm{sup −3}. In this case, the density of radiative loss power in the center of the layer is estimated to be about 500–600 kW/m{sup 3}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three dimensional effects on proton acceleration by intense laser solid target interaction

Multi-dimensional effects on ion acceleration by a normally incident linearly polarized intense laser pulse interacting with a thin solid target have been investigated numerically, where the laser has the peak intensity of 1.37×10{sup 20} W/cm{sup 2}, focused spot size of 6 μm, pulse duration of 33 fs, and total pulse energy about 3 J, which are commercially available now. We have checked the effects of simulation geometries by running one, two, and three dimensional (1D, 2D, 3D) particle-in-cell simulations. 3D simulation results show that, in the case of using a relatively thick target (in the opaque regime, i.e., 2 μm) with the so-called target normal sheath field acceleration mechanism, electrons spread almost uniformly along two transverse directions. While in the case of using an ultra-thin target (in the relativistic-induced transparent regime, i.e., 100 nm) with the so-called break-out afterburner mechanism, electrons spread more quickly along the direction orthogonal to the laser polarization direction especially at the early stage. The transverse spreading of electrons strongly decreases the electron density at the rear side of the target. Such an effect causes different estimation of electron temperatures in different simulation geometries. Usually, 1D and 2D simulations overestimate the temperature; and as a result, the maximum proton energy observed in 1D and 2D simulations is, respectively, about 3 and 2 times of that observed in 3D simulation.

43 PARTICLE ACCELERATORS↗

Transversely pumped laser driven particle accelerator

We present a novel acceleration scheme capable of accelerating electrons and ions in an underdense plasma. Transversely Pumped Acceleration (TPA) uses multiple arrays of counter-propagating laser beamlets that focus onto a central acceleration axis. Tuning the injection timing and the spacing between the adjacent beamlets allows for precise control over the position and velocity of the intersection point of the counter-propagating beam arrays. This results in an accelerating structure that propagates orthogonal to the direction of laser propagation. We present the theory that sets the injection timing of the incoming pulses to accelerate electrons and ions with a tunable phase velocity plasma wave. Simulation results are also presented which demonstrate 1.12 GeV proton beams accelerated in 3.6 mm of plasma and electron acceleration gradients on the order of 1 TeV/m in a scheme that circumvents dephasing. This work has potential applications as a compact accelerator for medical physics and high energy physics colliders.

Laser-produced plasmas↗

Gyrokinetic turbulence cascade via predator-prey interactions between different scales

Gyrokinetic simulations in a closed fieldline geometry are presented to explore the physics of nonlinear transfer in plasma turbulence. As spontaneously formed zonal flows and small-scale turbulence demonstrate “predator-prey” dynamics, a particular cascade spectrum emerges. The electrostatic potential and the density spectra appear to be in good agreement with the simple theoretical prediction based on Charney-Hasegawa-Mima equation | ϕ{sup ~}{sub k} |{sup 2}∼| n{sup ~}{sub k} |{sup 2}∝k{sup −3}/(1+k{sup 2}){sup 2}, with the spectra becoming anisotropic at small scales. The results indicate that the disparate scale interactions, in particular, the refraction and shearing of larger scale eddies by the self-consistent zonal flows, dominate over local interactions, and contrary to the common wisdom, the comprehensive scaling relation is created even within the energy injection region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of nonlinear processes in a beam-plasma system

In the present paper we consider the efficiency of the electromagnetic radiation generation due to various nonlinear processes in the beam-plasma system. The beam and plasma parameters were chosen close to the parameters in the experiment on the GOL-3 facility (BINP SB RAS). The model of the collisionless plasma is described by system of the Vlasov-Maxwell equations with periodic boundary conditions. The parallel numerical algorithm is based on the particles-in-cell method (PIC) with mixed Euler-Lagrangian domain decomposition. Various scenarios of nonlinear evolution in the beam-plasma system under the influence of an external magnetic field in case of a low density beam were studied. The energy transfer from one unstable mode to the others modes was observed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A 1D coupled transport/cold plasma wave model for parallel ponderomotive density modification near RF actuators

The plasma density profile close to RF actuators plays a vital role in the efficiency of power coupling from the antenna to the core plasma of fusion devices. Nonlinear interactions can be significant when considering the interaction of applied RF power and the tenuous SOL plasma, and the inclusion of this region in computational analysis is essential to providing robust predictive simulations of plasma response to RF application. Of particular interest is the observed density modification during RF operation. Here we present a self consistent description of the nonlinear ponderomotive force in the antenna near field region, via coupling a 1D field aligned transport model with a 1D cold plasma wave model. Results are presented, which show a parallel electric field of 300 Vm -1 drove a density reduction of ~ 20% over 25 RF periods.

Barnett, Rhea↗

MAGNETOHYDRODYNAMIC SIMULATION OF A SIGMOID ERUPTION OF ACTIVE REGION 11283

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in an MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

79 ASTRONOMY AND ASTROPHYSICS↗

Black Holes: Connecting Theory and Simulation to Observation [Slides]

Agenda: Review the relevant basics of black holes and black hole astrophysics. Discuss the underlying principles and methodology behind simulating infalling plasma, along with recent advances. Show how we generate synthetic X-ray spectra starting from simulation data. Connect these X-ray spectra to real observations. What are the current challenges facing this technique, and how might they be resolved?

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗

What happens to full-f gyrokinetic transport and turbulence in a toroidal wedge simulation?

Here, in order to save the computing time or to fit the simulation size into a limited computing hardware in a gyrokinetic turbulence simulation of a tokamak plasma, a toroidal wedge simulation may be utilized in which only a partial toroidal section is modeled with a periodic boundary condition in the toroidal direction. The most severe restriction in the wedge simulation is expected to be in the longest wavelength turbulence, i.e., ion temperature gradient (ITG) driven turbulence. The global full-f gyrokinetic code XGC1 is used to compare the transport and turbulence properties from a toroidal wedge simulation against the full torus simulation in an ITG unstable plasma in a model toroidal geometry. It is found that (1) the convergence study in the wedge number needs to be conducted all the way down to the full torus in order to avoid a false convergence, (2) a reasonably accurate simulation can be performed if the correct wedge number N can be identified, (3) the validity of a wedge simulation may be checked by performing a wave-number spectral analysis of the turbulence amplitude |δΦ| and assuring that the variation of δΦ between the discrete kθ values is less than 25% compared to the peak |δΦ|, and (4) a frequency spectrum may not be used for the validity check of a wedge simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modeling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collaborative Research: Analysis of wave mode content in fully turbulent, moderately collisional plasma laboratory experiment and kinetic simulation

Our major goal for the project is to make magnetic and Langmuir-probe-based density fluctuations measurements using identical diagnostic and acquisition setups on two different laboratory-based turbulent plasma devices: SSX at Swarthmore College and a new plasma source in development at Bryn Mawr College (BMX). We propose to generate single-point correlation and k-filtering metrics using these experimental measurements. These metrics will be compared directly to kinetic simulations of the experimental environments in order to determine wave mode content of the laboratory turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HiPACE++ v0.1

Hipace++ is the new implementation of the DESY-LBNL code Hipace (not distributed nor licensed) for plasma acceleration, relying on the quasi-static particle-in-cell (PIC) method, enabling fast and accurate simulations of plasma accelerators for problems satisfying the quasi-static approximation. Plasma acceleration is a priority in the ATAP division, and such problems represent a wide range of parameters. Hipace++ is an HPC (high-performance computing) code running on supercomputers and it provides the first efficient implementation of the quasi-static PIC method on Graphic Processing Units (GPU), showing 10x speedup and using much less resources than CPU implementations for production simulations. It currently runs on Nvidia GPUs and will ultimately be portable to more accelerated HPC platforms. The porting strategy is novel, and can be adapted to other flavors of quasi-static PIC. Besides, Hipace++ aims at using mesh refinement for higher efficiency, thus enabling simulations unfeasible today. For this purpose, Hipace++ depends on the open-source library AMReX (Adaptive Mesh Refinement at EXascale) developed at LBNL.

Vay, Jean-Luc↗

Physics design for a lithium vapor box divertor experiment on Magnum-PSI

The lithium vapor box divertor is a potential solution for power exhaust in toroidal confinement devices. The divertor plasma interacts with a localized, dense cloud of lithium vapor, leading to volumetric radiation, cooling, recombination, and detachment. To minimize contamination of the core plasma, lithium vapor is condensed on cool (300°C to 400°C) baffles upstream of the detachment point. Before implementing this in a toroidal plasma device with a slot divertor geometry, we consider an experiment with a scaled baffled-pipe geometry in the high-power linear plasma device Magnum-PSI. Three 15 cm-scale open cylinders joined by 6 cm diameter ‘nozzles’ are positioned on the plasma beam axis upstream of a target. The central box may be loaded with several tens of grams of lithium, which can be evaporated at 650°C to produce a vapor predicted, using a simple plasma-neutral interaction model, to be dense enough to cause volumetric detachment in the plasma. The power delivered to the target and box walls as measured by increases in their temperatures after a 10 s plasma pulse can be compared to determine the effectiveness of the vapor in detaching the plasma. Direct Simulation Monte Carlo simulations are performed to estimate the flow rates of lithium vapor between the boxes and to estimate the trapping of H2 delivered by the plasma in the boxes, which could inadvertently lead to detachment. Details of the geometry, simulations, and possible diagnostic techniques are presented.

divertor↗

NASA Laboratory Astrophysics Workshop 2006 Introductory Remarks

NASA Laboratory Astrophysics Workshop 2006, is the fourth in a series of workshops held at four year intervals, to assess the laboratory needs of NASA's astrophysics missions - past, current and future. Investigators who need laboratory data to interpret their observations from space missions, theorists and modelers, experimentalists who produce the data, and scientists who compile databases have an opportunity to exchange ideas and understand each other's needs and limitations. The multi-wavelength character of these workshops allows cross-fertilization of ideas, raises awareness in the scientific community of the rapid advances in other fields, and the challenges it faces in prioritizing its laboratory needs in a tight budget environment. Currently, we are in the golden age of Space Astronomy, with three of NASA s Great Observatories, Hubble Space Telescope (HST), Chandra X-Ray Observatory (CXO), and Spitzer Space Telescope (SST), in operation and providing astronomers and opportunity to perform synergistic observations. In addition, the Far Ultraviolet Spectroscopic Explorer (FUSE), XMM-Newton, HETE-2, Galaxy Evolution Explorer (GALEX), INTEGRAL and Wilkinson Microwave Anisotropy Probe (WMAP), are operating in an extended phase, while Swift and Suzaku are in their prime phase of operations. The wealth of data from these missions is stretching the Laboratory Astrophysics program to its limits. Missions in the future, which also need such data include the James Webb Space Telescope (JWST), Space Interferometry Mission (SIM), Constellation-X (Con-X), Herschel, and Planck. The interpretation of spectroscopic data from these missions requires knowledge of atomic and molecular parameters such as transition probabilities, f-values, oscillator strengths, excitation cross sections, collision strengths, which have either to be measured in the laboratory by simulating space plasma and interactions therein, or by theoretical calculations and modeling. Once the laboratory data are obtained, a key step to making them available to the observer is the creation and maintenance of critically compiled databases. Other areas of study, that are important for understanding planet formation, and for detection of molecules that are indicators of life, are also supported by the Laboratory Astrophysics program. Some examples are: studies of ices and dust grains in a space environment; nature and evolution of interstellar carbon-rich dust; and polycyclic aromatic hydrocarbons. In addition, the program provides an opportunity for the investigation of novel ideas, such as simulating radiative shock instabilities in plasmas, in order to understand jets observed in space. A snapshot of the currently funded program, mission needs, and relevance of laboratory data to interpreting observations, will be obtained at this workshop through invited and contributed talks and poster papers. These will form the basis for discussions in splinter groups. The Science Organization Committee will integrate the results of the discussions into a coherent White Paper, which will provide guidance to NASA in structuring the Laboratory Astrophysics program in subsequent years, and also to the scientific community in submitting research proposals to NASA for funding.

Hasan, Hashima↗

Strong Coulomb coupling influences ion and neutral temperatures in atmospheric pressure plasmas

Molecular dynamics simulations are used to model ion and neutral temperature evolution in partially-ionized atmospheric pressure plasma at different ionization fractions. Results show that ion–ion interactions are strongly coupled at ionization fractions as low as 10 –5 and that the temperature evolution is influenced by effects associated with the strong coupling. Specifically, disorder-induced heating is found to rapidly heat ions on a timescale of the ion plasma period (~10 s ps) after an ionization pulse. This is followed by the collisional relaxation of ions and neutrals, which cools ions and heats neutrals on a longer (~ns) timescale. Slight heating then occurs over a much longer (~100 s ns) timescale due to ion-neutral three-body recombination. An analytic model of the temperature evolution is developed that agrees with the simulation results. Finally, a conclusion is that strong coupling effects are important in atmospheric pressure plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗