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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Ion heating dynamics in solid buried layer targets irradiated by ultra-short intense laser pulses

We investigate bulk ion heating in solid buried layer targets irradiated by ultra-short laser pulses of relativistic intensities using particle-in-cell simulations. Our study focuses on a CD{sub 2}-Al-CD{sub 2} sandwich target geometry. We find enhanced deuteron ion heating in a layer compressed by the expanding aluminium layer. A pressure gradient created at the Al-CD{sub 2} interface pushes this layer of deuteron ions towards the outer regions of the target. During its passage through the target, deuteron ions are constantly injected into this layer. Our simulations suggest that the directed collective outward motion of the layer is converted into thermal motion inside the layer, leading to deuteron temperatures higher than those found in the rest of the target. This enhanced heating can already be observed at laser pulse durations as low as 100 fs. Thus, detailed experimental surveys at repetition rates of several ten laser shots per minute are in reach at current high-power laser systems, which would allow for probing and optimizing the heating dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Parameter spaces for linear and nonlinear whistler-mode waves

We examine the growth of magnetospheric whistler-mode waves which comprises a linear growth phase followed by a nonlinear growth phase. We construct time-profiles for the wave amplitude that smoothly match at the transition between linear and nonlinear wave growth. This matching procedure can only take place over a limited “matching region” in (N{sub h}/N{sub 0},A{sub T})-space, where A{sub T} is the electron thermal anisotropy, N{sub h} is the hot (energetic) electron number density, and N{sub 0} is the cold (background) electron number density. We construct this matching region and determine how the matching wave amplitude varies throughout the region. Further, we specify a boundary in (N{sub h}/N{sub 0},A{sub T})-space that separates a region where only linear chorus wave growth can occur from the region in which fully nonlinear chorus growth is possible. We expect that this boundary should prove of practical use in performing computationally expensive full-scale particle simulations, and in interpreting experimental wave data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Plasma contactor design for electrodynamic tether applications

The plasma contacting process is described and experiments are discussed that suggest the key role that cold ions play in establishing a low impedance plasma bridge that can conduct current in either direction between a contactor electrode and a dilute plasma. A ring cusp contactor is shown to provide from 1000-mA of electron emission to 500-mA of electron collection as its bias relative to a simulated space plasma is varied through an 80-v range.

Wilbur, Paul J.↗
DOE OSTI2020

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Self-driven current generation in turbulent fusion plasmas

Global gyrokinetic simulations with self-consistent coupling of neoclassical and turbulent dynamics demonstrate that turbulence can greatly affect plasma self-driven mean current generation in tokamaks. The current amplitude, profile and associated phase space structures can all be modified. Turbulence can critically reduce the current generation in the collisionless regime, generate current profile corrugation near the rational magnetic surface and nonlocally drive current in the linearly stable region—all these are expected to have a radical impact on broad tokamak physics. Both electron parallel acceleration and residual stress from turbulence play tremendous roles in turbulence-induced current generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MHD modeling of shattered pellet injection in JET

Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of a synthetic fast ion loss detector model for Wendelstein 7-X

Abstract We present the first validated synthetic diagnostic for fast ion loss detectors (FILDs) in the Wendelstein 7-X (W7-X) stellarator. This model has been developed on, and validated against experimental data from, a FILD provided by the National Institute for Fusion Science (NIFS-FILD), with potential future applicability to the existing Faraday Cup FILD (FC-FILD) on W7-X as well as the scintillating FILD (S-FILD) currently under development. A workflow combining Monte Carlo codes BEAMS3D and ASCOT5 is used to track fast ions produced by neutral beam injection from the moment of ionization until they are thermalized or lost from the last closed flux surface, and from there to a virtual plane which serves as a projection of the entrance aperture to the FILD. Simulations in ASCOT5 are analyzed via a geometric method to determine the probability of transmission through the FILD aperture and onto the detector as a function of normalized momentum, pitch angle, gyrophase, and position at the virtual plane. This probability is then applied to the simulated ions arriving from the plasma, producing a simulated signal from a computationally tractable number of simulated fast ions. Simulated signals are presented for two W7-X experiments with neutral beam injection and quantitatively compared with experimental measurements from the NIFS-FILD diagnostic. An estimate of the frequency of charge-exchange with neutral particles in the edge is performed, and it is found that this process may have a significant impact on the measured signals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Preliminary study of the CRRES magnetospheric barium releases

Preliminary theoretical and computational analyses of the Combined Release and Radiation Effects Satellite (CRRES) magnetospheric barium releases are presented. The focus of the studies is on the evolution of the diamagnetic cavity which is formed by the barium ions as they expand outward, and on the structuring of the density and magnetic field during the expansion phase of the releases. Two sets of simulation studies are discussed. The first set is based upon a 2D ideal MHD code and provides estimates of the time and length scales associated with the formation and collapse of the diamagnetic cavity. The second set uses a nonideal MHD code; specifically, the Hall term is included. This additional term is critical to the dynamics of sub-Alfvenic plasma expansions, such as the CRRES barium releases, because it leads to instability of the expanding plasma. Detailed simulations of the G4 and G10 releases were performed. In both cases the expanding plasma rapidly structured: the G4 release structured at time t less than about 3 s and developed scale sizes of about 1-2 km, while the G10 release structured at time t less than about 22 s and developed scale sizes of about 10-15 km. It is also found that the diamagnetic cavity size is reduced from those obtained from the ideal MHD results because of the structure. On the other hand, the structuring allows the formation of plasma blobs which appear to free stream across the magnetic field; thus, the barium plasma can propagate to larger distances traverse to the magnetic field than the case where no structuring occurs. Finally, a new normal mode of the system was discovered which may be excited at the leading edge of the expanding barium plasma.

Huba, J. D.↗
DOE OSTI2022

Spectral studies of soft x-ray radiation of laser-produced plasma of various target materials in a wide spectral range

Here, the work is devoted to experimental study of soft X-ray radiation spectra of laser produced plasma in a wide spectral range 5–100 Å at the “Kamerton” facility (GPI) with pulse duration 70 ps, pulse energy 1–5 J, wavelength 0.53 microns at which the laser intensity (power flux density) of 7×10 14 – 3.5×10 15 W/cm 2 was achieved. A spectrograph was used, which had transmission diffraction gratings with a ratio of the elementary gap to the period of the structure of 0.25 and 0.41. Detection was performed on both UV-4 X-ray photo film and Fuji TR fluorescent imaging plate. Solid samples from Al, Si, Ti, Cu, Ta and W were used as targets. The ionization states of the plasma corresponding to various electron temperatures were calculated, which made it possible to estimate the electron temperature by comparison these calculation results with the experimentally obtained spectra. The estimated electron temperature, which depends on the laser pulse energy and the target material, varied within the range of 100–450 eV. To verify the correctness of the temperature estimations obtained by such comparison a numerical simulation of plasma radiation was carried out by the use of PrismSPECT computer program. It was found that the results of this simulation are in a good agreement with estimations on the base of experimentally obtained spectra. The analysis of these spectra showed that tantalum, tungsten or titanium targets are the best candidates among the tested ones for the use of laser produced plasma as a radiation source in the "water window" spectral range (23–44 Å) for applications in biology and medicine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modelling the limiter ramp-up of WEST for addressing the future challenges of ITER

This paper presents a joint experimental and numerical investigation into the physics of long limited plasma ramp-up in tokamaks with tungsten (W) first walls, a critical phase for ITER operation. The comparison between the average plasma quantities simulated using the SolEdge-HDG code and the measurements taken during three successive WEST discharges after boronisation shows how challenging it is to predict this phase. While simulations reproduce the general trends at the midplane, with reasonable match in density profiles, they consistently underestimate core temperatures possibly due to too large perpendicular heat conductivity. On the contrary, at the high field side (HFS) limiter, simulations overestimate the measured quantities, and highlights the limitation of using Bohm boundary conditions at grazing magnetic angles. Experimental measurements reveal that the boron layer is rapidly eroded, on a timescale comparable to a single ITER discharge. The subsequent transition from a boron-coated to a tungsten wall increases recycling and significantly degrades the core plasma, reducing the electron temperature by nearly half due to W contamination, despite wall parameters remaining stable. Furthermore, comparisons with Langmuir probes, bolometry, reflectometry, and spectroscopy indicate that the experimental far scrape-Off layer (SOL) is significantly wider than simulated. This wide SOL implies that boron erosion extends along the entire HFS limiter rather than being confined to the contact point. This work highlights some characteristics of the plasma during this phase of the discharge and emphasizes the current modelling issues that need to be resolved in order to obtain reliable predictions concerning the ITER ramp-up.

ITER↗
DOE OSTI2020

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Enhanced transportation of energetic electrons in dual-frequency atmospheric microplasmas

A comparative study of electron kinetics between single-frequency (SF) microplasmas and their equivalent dual-frequency (DF) microplasmas with matching effective frequencies in atmospheric-pressure helium discharges was performed using particle-in-cell simulation with a Monte Carlo collision. The effective-frequency concept helps in analyzing DF microplasmas in a fashion similar to SF microplasmas with effective parameters. In this study, the plasma characteristics such as the plasma potential, density, and electron energy probability functions of the SF microplasma and its DF counterpart were almost the same. However, the oscillating sheath edge was pushed further into the electrode for a substantial fraction of the time and the sheath width decreased in DF microplasmas. As a result, the transportation of the energetic electrons ({epsilon} > 4 eV) usable for tailoring the surface chemistry in atmospheric microplasmas is enhanced in DF microplasmas as compared to SF microplasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A Fokker-Planck-Landau collision equation solver on two-dimensional velocity grid and its application to particle-in-cell simulation

An approximate two-dimensional solver of the nonlinear Fokker-Planck-Landau collision operator has been developed using the assumption that the particle probability distribution function is independent of gyroangle in the limit of strong magnetic field. The isotropic one-dimensional scheme developed for nonlinear Fokker-Planck-Landau equation by Buet and Cordier [J. Comput. Phys. 179, 43 (2002)] and for linear Fokker-Planck-Landau equation by Chang and Cooper [J. Comput. Phys. 6, 1 (1970)] have been modified and extended to two-dimensional nonlinear equation. In addition, a method is suggested to apply the new velocity-grid based collision solver to Lagrangian particle-in-cell simulation by adjusting the weights of marker particles and is applied to a five dimensional particle-in-cell code to calculate the neoclassical ion thermal conductivity in a tokamak plasma. Error verifications show practical aspects of the present scheme for both grid-based and particle-based kinetic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Overview of the preliminary design of the ITER plasma control system

An overview of the Preliminary Design of the ITER Plasma Control System (PCS) is described here, which focusses on the needs for 1st plasma and early plasma operation in hydrogen/helium (H/He) up to a plasma current of 15 MA with moderate auxiliary heating power in low confinement mode (L-mode). Candidate control schemes for basic magnetic control, including divertor operation and kinetic control of the electron density with gas puffing and pellet injection, were developed. Commissioning of the auxiliary heating systems is included as well as support functions for stray field topology and real-time plasma boundary reconstruction. Initial exception handling schemes for faults of essential plant systems and for disruption protection were developed. The PCS architecture was also developed to be capable of handling basic control for early commissioning and the advanced control functions that will be needed for future high performance operation. A plasma control simulator is also being developed to test and validate control schemes. To handle the complexity of the ITER PCS, a systems engineering approach has been adopted with the development of a plasma control database to keep track of all control requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Asymmetric modes decomposition in an overmoded relativistic backward wave oscillator

Most of the investigated overmoded relativistic backward wave oscillators (RBWOs) are azimuthally symmetric; thus, they are designed through two dimensional (2-D) particle-in-cell (PIC) simulations. However, 2-D PIC simulations cannot reveal the effect of asymmetric modes on beam-wave interaction. In order to investigate whether asymmetric mode competition needs to be considered in the design of overmoded RBWOs, a numerical method of determining the composition of both symmetric and asymmetric modes in three dimensional (3-D) PIC simulations is introduced in this paper. The 2-D and 3-D PIC simulation results of an X-band overmoded RBWO are analyzed. Our analysis indicates that the 2-D and 3-D PIC simulation results of our device are quite different due to asymmetric mode competition. In fact, asymmetric surface waves, especially EH{sub 11} mode, can lead to serious mode competition when electron beam propagates near the surface of slow wave structures (SWSs). Therefore, additional method of suppressing asymmetric mode competition, such as adjusting the reflections at both ends of SWSs to decrease the Q-factor of asymmetric modes, needs to be utilized in the design of overmoded RBWOs. Besides, 3-D PIC simulation and modes decomposition are essential for designing overmoded RBWOs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Using our new 3-D relativistic particle-in-cell (PIC) code, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic compression. Behind the bow shock, in the jet shock, strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to go beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations and calculated radiation based on first principles. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the technique to calculate emission from electrons based on simulations with a small system. We obtain spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system may generate a jitter/synchrotron spectrum.

Nishikawa, K.-I.↗
DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗