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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,801 records · Page 100

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Study of the Interaction of Plasma with Beryllium That Is a Candidate Material for the First Wall of a Fusion Reactor

The experimental conditions for treatment of a candidate material of the first wall, beryllium, by the plasma have been determined. The plasma parameters obtained using the simulator with a plasma-beam setup during the irradiation of samples as well as the results of material studies obtained by scanning electron microscopy and energy-dispersion analysis of the beryllium surface layer have been considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

The influence of intense electric fields on three-dimensional asymmetric magnetic reconnection

A three-dimensional particle-in-cell simulation of magnetic reconnection in an asymmetric configuration without a guide field and with temperature ratio T{sub i}/T{sub e}>1 demonstrates that intense perpendicular electric fields are produced on the low-density side of the current layer where there is a strong gradient in the plasma density. The simulation shows that the 3-D reconnection rate is unaffected by these intense electric fields, that the electron current layer near the X line remains coherent and does not break up, but that localized regions of strong energy dissipation exist along the low-density separatrices. Near the X line the dominant term in the generalized Ohm's law for the reconnection electric field remains the off-diagonal electron pressure gradient ∂P{sub exy}/∂x. On the low-beta separatrix, however, the anomalous drag −<δnδE{sub y}>/ makes an equally important contribution to that of the pressure gradient to the average E{sub y} field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Benchmark studies of the gyro-Landau-fluid code and gyro-kinetic codes on kinetic ballooning modes

A Gyro-Landau-Fluid (GLF) 3 + 1 model has been recently implemented in BOUT++ framework, which contains full Finite-Larmor-Radius effects, Landau damping, and toroidal resonance [Ma et al., Phys. Plasmas 22, 055903 (2015)]. A linear global beta scan has been conducted using the JET-like circular equilibria (cbm18 series), showing that the unstable modes are kinetic ballooning modes (KBMs). In this work, we use the GYRO code, which is a gyrokinetic continuum code widely used for simulation of the plasma microturbulence, to benchmark with GLF 3 + 1 code on KBMs. To verify our code on the KBM case, we first perform the beta scan based on “Cyclone base case parameter set.” We find that the growth rate is almost the same for two codes, and the KBM mode is further destabilized as beta increases. For JET-like global circular equilibria, as the modes localize in peak pressure gradient region, a linear local beta scan using the same set of equilibria has been performed at this position for comparison. With the drift kinetic electron module in the GYRO code by including small electron-electron collision to damp electron modes, GYRO generated mode structures and parity suggest that they are kinetic ballooning modes, and the growth rate is comparable to the GLF results. However, a radial scan of the pedestal for a particular set of cbm18 equilibria, using GYRO code, shows different trends for the low-n and high-n modes. The low-n modes show that the linear growth rate peaks at peak pressure gradient position as GLF results. However, for high-n modes, the growth rate of the most unstable mode shifts outward to the bottom of pedestal and the real frequency of what was originally the KBMs in ion diamagnetic drift direction steadily approaches and crosses over to the electron diamagnetic drift direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Nonlinear evolution of stimulated Raman scattering near the quarter-critical density

Nonlinear evolution of stimulated Raman scattering (SRS) near the quarter-critical density is studied using one-dimensional (1D) and two-dimensional (2D) particle-in-cell simulations in homogeneous plasmas. In 1D configuration, with two-plasmon decay (TPD) instability excluded, the system evolves into two regimes distinguished by whether density cavities have been formed or not. At low temperatures, a cavity regime characterised by high absorption and low reflection, and at high temperatures nonlinear frequency shift regime due to particle trapping, are observed. Furthermore, a competition between SRS and TPD in 2D simulations reveals that the nonlinear SRS does play a significant role near the quarter-critical density, whose influences were mostly neglected before.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Blob Formation and Ejection in Coronal Jets due to the Plasmoid and Kelvin–Helmholtz Instabilities

We perform 2D resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible to the formation of blobs that are ejected in the jet. Our simulation with low plasma β (Case I) shows that magnetic islands form easily and propagate upward in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 Å and 211 Å) and the cool channel (304 Å) in observations by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory . The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s{sup −1}. These parameters are similar to the properties of blobs observed in extreme-ultraviolet (EUV) jets by AIA. The Kelvin–Helmholtz instability develops in our simulation with moderately high plasma β (Case II) and leads to the formation of bright vortex-like blobs above the multiple high magnetosonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magnetosonic Mach number regions disappear. In the lower plasma β case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current-sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin–Helmholtz instability along the jet are both possible causes of the formation of blobs observed at EUV wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Nonlinear plasma wave in magnetized plasmas

Nonlinear axisymmetric cylindrical plasma oscillations in magnetized collisionless plasmas are a model for the electron fluid collapse on the axis behind an ultrashort relativisically intense laser pulse exciting a plasma wake wave. We present an analytical description of the strongly nonlinear oscillations showing that the magnetic field prevents closing of the cavity formed behind the laser pulse. This effect is demonstrated with 3D PIC simulations of the laser-plasma interaction. An analysis of the betatron oscillations of fast electrons in the presence of the magnetic field reveals a characteristic “Four-Ray Star” pattern.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Assessing the impact of alpha particles on thermal confinement in JET D-T plasmas through global GENE-Tango simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1 keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge–despite this being the shot with the highest fusion output–is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization. These results demonstrate that an increased alpha particle concentration can significantly impact transport properties under simulated JET experimental conditions. However, these findings cannot be directly extrapolated to ITER due to the substantial differences in parameters such as plasma size, magnetic field, plasma current, and thermal pressure.

energetic particles↗

DOE OSTI2022

Modeling deformation and failure in AlSi-polyester abradable sealcoating material using microstructure-based finite element simulation

A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Ion beam emittance from an ECRIS

Simulation of ion beam extraction from an Electron Cyclotron Resonance Ion Source (ECRIS) is a fully 3 dimensional problem, even if the extraction geometry has cylindrical symmetry. Because of the strong magnetic flux density, not only the electrons are magnetized but also the Larmor radius of ions is much smaller than the geometrical dimension of the plasma chamber (Ø 64 × 179 mm). If we assume that the influence of collisions is small on the path of particles, we can do particle tracking through the plasma if the initial coordinates of particles are known. We generated starting coordinates of plasma ions by simulation of the plasma electrons, accelerated stochastically by the 14.5 GHz radio frequency power fed to the plasma. With that we were able to investigate the influence of different electron energies on the extracted beam. Using these assumptions, we can reproduce the experimental results obtained 10 years ago, where we monitored the beam profile with the help of viewing targets. Additionally, methods have been developed to investigate arbitrary 2D cuts of the 6D phase space. To this date, we are able to discuss full 4D information. Currently, we extend our analysis tool towards 5D and 6D, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1992

Interplanetary magnetic field changes and condensations in comet Halley's plasma tail

In a time-dependent three dimensional MHD simulation for cometary plasmas, Schmidt-Voigt (1989) could observe the formation of condensations in the plasma tail after a 90 degree change in the interplanetary magnetic field (IMF) sweeping over the comet. We investigated the IMF measurements of the Vega SC in the vicinity of the comet Halley for 90 degree changes in the clock angle and studied the relation between them and optical observations of condensations in the plasma tail. For the time interval 24 Feb. 86 to 14 Mar. 86, we could not find a correlation between such changes and the release of condensations from the cometary head.

Delva, Magda↗

DOE OSTI2014

Fully three-dimensional simulation and modeling of a dense plasma focus

A dense plasma focus (DPF) is a pulsed-power machine that electromagnetically accelerates and cylindrically compresses a shocked plasma in a Z-pinch. The pinch results in a brief (~ 100 ns) pulse of X-rays, and, for some working gases, also a pulse of neutrons. A great deal of experimental research has been done into the physics of DPF reactions, and there exist mathematical models describing its behavior during the different time phases of the reaction. Two of the phases, known as the inverse pinch and the rundown, are approximately governed by magnetohydrodynamics, and there are a number of well-established codes for simulating these phases in two dimensions or in three dimensions under the assumption of axial symmetry. There has been little success, however, in developing fully three-dimensional simulations. In this work we present three-dimensional simulations of DPF reactions and demonstrate that three-dimensional simulations predict qualitatively and quantitatively different behavior than their two-dimensional counterparts. One of the most important quantities to predict is the time duration between the formation of the gas shock and Z-pinch, and the three-dimensional simulations more faithfully represent experimental results for this time duration and are essential for accurate prediction of future experiments.

Dense plasma focus, magnetohydrodynamics, simulati↗

NASA NTRS1981

Numerical simulation of the generation of electrostatic turbulence in the magnetotail

A two-dimensional plasma model is used to investigate the development of electrostatic turbulence in a magnetized plasma from plasma instabilities. The simulation consists of following the motion of 100,000 ions in their self-consistent electrostatic field. The electrons are treated as a constant neutralizing background. The instabilities modeled are driven by a ring-type velocity distribution and by interpenetrating ion beams in a time-variable magnetic field. Instability growth times are of the order of an ion gyroperiod in the case of the ring distribution and of the order of an ion plasma period in the case of the beam simulation. Maximum potential differences generated are of the order of the ion kinetic energies. These simulations demonstrate the cascade of wave energy to long wavelengths, thus showing the E x B turbulence can be generated from plasma microinstabilities. After the free energy feeding, the instabilities are exhausted, and wave energy at wavelengths less than an ion gyrodiameter decays quickly to equilibrium levels, while longer wavelength modes persist for much longer times. In one model with a time dependent, but spatially uniform, magnetic field the electric field energy at long wavelengths appeared to increase as a result of the increase of the magnetic field.

Swift, D. W.↗

DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Simulation and experimental study of proton bunch self-modulation in plasma with linear density gradients

We present numerical simulations and experimental results of the self-modulation of a long proton bunch in a plasma with linear density gradients along the beam path. Simulation results agree with the experimental results reported [F. Braunmller, T. Nechaeva et al. (AWAKE Collaboration), Phys. Rev. Lett. 125, 264801 (2020)]: with negative gradients, the charge of the modulated bunch is lower than with positive gradients. In addition, the bunch modulation frequency varies with gradient. Simulation results show that dephasing of the wakefields with respect to the relativistic protons along the plasma is the main cause for the loss of charge. The study of the modulation frequency reveals details about the evolution of the self-modulation process along the plasma. In particular for negative gradients, the modulation frequency across time-resolved images of the bunch indicates the position along the plasma where protons leave the wakefields. Simulations and experimental results are in excellent agreement.

43 PARTICLE ACCELERATORS↗

DOE OSTI2014

A plasma source driven predator-prey like mechanism as a potential cause of spiraling intermittencies in linear plasma devices

Three-dimensional global drift fluid simulations are carried out to analyze coherent plasma structures appearing in the NAGDIS-II linear device (nagoya divertor plasma Simulator-II). The numerical simulations reproduce several features of the intermittent spiraling structures observed, for instance, statistical properties, rotation frequency, and the frequency of plasma expulsion. The detailed inspection of the three-dimensional plasma dynamics allows to identify the key mechanism behind the formation of these intermittent events. The resistive coupling between electron pressure and parallel electric field in the plasma source region gives rise to a quasilinear predator-prey like dynamics where the axisymmetric mode represents the prey and the spiraling structure with low azimuthal mode number represents the predator. This interpretation is confirmed by a reduced one-dimensional quasilinear model derived on the basis of the findings in the full three-dimensional simulations. The dominant dynamics reveals certain similarities to the classical Lotka-Volterra cycle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Study of self-consistent particle flows in a plasma blob with particle-in-cell simulations

The self-consistent particle flows in a filamentary coherent structure along the magnetic field line in scrape-off layer (SOL) plasma (plasma blob) have been investigated by means of a three-dimensional electrostatic particle-in-cell simulation code. The presence of the spiral current system composed of the diamagnetic and parallel currents in a blob is confirmed by the particle simulation without any assumed sheath boundary models. Furthermore, the observation of the electron and ion parallel velocity distributions in a blob shows that those distributions are far from Maxwellian due to modification with the sheath formation and that the electron temperature on the higher potential side in a blob is higher than that on the lower potential side. Also, it is found that the ions on the higher potential side are accelerated more intensively along the magnetic field line than those on the lower potential side near the edge. This study indicates that particle simulations are able to provide an exact current closure to analysis of blob dynamics and will bring more accurate prediction of plasma transport in the SOL without any empirical assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

The impact of plasma dynamics on the self-magnetic-pinch diode impedance

In this study, the self-magnetic-pinch diode is being developed as an intense electron beam source for pulsed-power-driven x-ray radiography. The basic operation of this diode has long been understood in the context of pinched diodes, including the dynamic effect that the diode impedance decreases during the pulse due to electrode plasma formation and expansion. Experiments being conducted at Sandia National Laboratories' RITS-6 accelerator are helping to characterize these plasmas using time-resolved and time-integrated camera systems in the x-ray and visible. These diagnostics are analyzed in conjunction with particle-in-cell simulations of anode plasma formation and evolution. The results confirm the long-standing theory of critical-current operation with the addition of a time-dependent anode-cathode gap length. The results may suggest that anomalous impedance collapse is driven by increased plasma radial drift, leading to larger-than-average ion v r × B θ acceleration into the gap.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Parametric instability of a monochromatic Alfven wave: Perpendicular decay in low beta plasma

Two-dimensional hybrid simulations are performed to investigate the parametric decay of a monochromatic Alfven wave in low beta plasma. Both the linearly and left-hand polarized pump Alfven waves are considered in the paper. For the linearly polarized pump Alfven wave, either a parallel or obliquely propagating wave can lead to the decay along the perpendicular direction. Initially, the parametric decay takes place along the propagating direction of the pump wave, and then the decay occurs in the perpendicular direction. With the increase of the amplitude and the propagating angle of the pump wave (the angle between the propagating direction of the pump wave and the ambient magnetic field), the spectral range of the excited waves becomes broad in the perpendicular direction. But the effects of the plasma beta on the spectral range of the excited waves in perpendicular direction are negligible. However, for the left-hand polarized pump Alfven wave, when the pump wave propagates along the ambient magnetic field, the parametric decay occurs nearly along the ambient magnetic field, and there is no obvious decay in the perpendicular direction. Significant decay in the perpendicular direction can only be found when the pump wave propagates obliquely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Refine your search to explore more results.