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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Addressing key physics problems in high-energy-density plasmas with a novel kinetic simulation capability

Many important physical processes in inertial confinement fusion (ICF) and dense Z-pinch (DZP) experiments require a kinetic (velocity-space-dependent) description. Conventional particle-in-cell (PIC) methods are poorly suited for high-energy-density (HED) plasmas, due to restrictive time-step constraints and the inability to conserve energy. In a previous LDRD (21-FS-048), we demonstrated that a fully implicit PIC formulation overcomes these limitations: it conserves energy even when coupled with Coulomb collision models and can be solved efficiently with large grid cells and large time steps. Thus, it is feasible to use this method to study kinetic effects in ICF and DZP plasmas on hydro-like time and spatial scales. In this follow-on LDRD, we advanced this methodology into a high-fidelity tool for production-scale simulations and used it to answer key questions relevant to ICF and DZP experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗
DOE OSTI2018

Wavelet methods for studying the onset of strong plasma turbulence

Recent simulations have demonstrated that coherent current sheets dominate the kinetic-scale energy dissipation in strong turbulence of magnetized plasma. Wavelet basis functions are a natural tool for analyzing turbulent flows containing localized coherent structures of different spatial scales. Here, wavelets are used to study the onset and subsequent transition to fully developed turbulence from a laminar state. Originally applied to neutral fluid turbulence, an iterative wavelet technique decomposes the field into coherent and incoherent contributions. In contrast to Fourier power spectra, finite time Lyapunov exponents, and simple measures of intermittency such as non-Gaussian statistics of field increments, the wavelet technique is found to provide a quantitative measure for the onset of turbulence and to track the transition to fully developed turbulence. The wavelet method makes no assumptions about the structure of the coherent current sheets or the underlying plasma model. Temporal evolution of the coherent and incoherent wavelet fluctuations is found to be highly correlated (a Pearson correlation coefficient of >0.9) with the magnetic field energy and plasma thermal energy, respectively. The onset of turbulence is identified with the rapid growth of a background of incoherent fluctuations spreading across a range of scales and a corresponding drop in the coherent components. This is suggestive of the interpretation of the coherent and incoherent wavelet fluctuations as measures of coherent structures (e.g., current sheets) and dissipation, respectively. The ratio of the incoherent to coherent fluctuations Ric is found to be fairly uniform in the turbulent state across different plasma models and provides an empirical threshold of ~0.1 for turbulence onset. The utility of this technique is illustrated through examples. First, it is applied to the Kelvin–Helmholtz instability from different simulation models including fully kinetic, hybrid (kinetic ion/fluid electron), and Hall MHD simulations. Second, the wavelet diagnostic is applied to the development of turbulence downstream of the bowshock in a global magnetosphere simulation. Finally, the wavelet technique is also shown to be useful as a de-noising method for particle simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Improving the homogeneity of alternating current-drive atmospheric pressure dielectric barrier discharges in helium with an additional low-amplitude radio frequency power source: A numerical study

It was proposed in this paper that the homogeneity of the atmospheric pressure discharge driven by an ac power source could be improved by applying an auxiliary low-amplitude rf power source. To verify the idea, a two-dimensional fluid model then was applied to study the atmospheric discharges in helium driven by ac power, low-amplitude rf power, and combined ac and low-amplitude rf power, respectively. Simulation results confirmed that an auxiliary rf power could improve the homogeneity of a discharge driven by an ac power source. It was further found that there existed a threshold voltage of the rf power source leading to the transition from inhomogeneous to homogeneous discharge. As the frequency of the rf power source increased from 2 to 22 MHz, the magnitude of the threshold voltage dropped first rapidly and then to a constant value. When the frequency was over 13.56 MHz, the magnitude of the threshold voltage was smaller than one-sixth of the ac voltage amplitude under the simulated discharge parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear wavenumber shift of large amplitude Langmuir waves

Nonlinear particle-in-cell simulation is carried out to investigate the nonlinear behavior of the Langmuir wave launched with a fixed frequency in a uniform plasma. It is found that in the strong driving case, the launched wave propagates in a phase velocity larger than that predicted by the linear theory; there appears a nonlinear down-shift of wavenumber. The phase velocity of the nonlinear wave and the down-shift of the wavenumber are demonstrated to be determined by the velocity of nonlinearly accelerated resonant electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2013

On gas desorption from the tokamak first wall during edge localized modes

The effect of gas desorption from the tokamak first wall on the pedestal recovery in the H-mode after an edge-localized-mode burst is considered. Results of FACE code simulations of hydrogen desorption from a beryllium wall are presented. It is found that the wall has a significant effect on plasma processes only at sufficiently low temperatures (of about 400 K), which agrees with qualitative estimates obtained earlier in the zero-dimensional approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Nonlinear fishbone dynamics in spherical tokamaks

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. The results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. A substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Prediction of Breakdown in Air and Solid Dielectrics: A Complete Plasma Model from Discharge Initiation to Flashover (Final Technical Report)

A major challenge for most of the technologies of today, and the technologies of the future, is the compact and reliable design of electrical insulation systems for their electrical components. Although controlled low-temperature plasmas (LTPs) or non-thermal plasmas in air at atmospheric pressure are of interest for a wide range of technologies, the appearance of such discharges is considered highly undesirable in insulation systems because they are precursors of a breakdown or flashover, where complete insulation failure takes place. For non-thermal plasmas, often only one mechanism, or the transition from one discharge mechanism to another, has typically been studied to date. On the other hand, in the case of thermal plasma applications, the conditions following the occurrence of an electric arc or flashover have always been considered. Here, there is a big technical gap: there is no paper describing a complete plasma model simulating all discharge mechanisms and their transitions from discharge initiation to electrical arc and flashover. To address this gap, a complete plasma model for both thermal and non-thermal plasmas should be developed using a hydrodynamic approach. The project aimed to build the model. The work has the potential to address industrial applications and is of great importance to higher voltage systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Gyrokinetic studies of trapped electron mode turbulence in the Helically Symmetric eXperiment stellarator

Gyrokinetic simulations of plasma microturbulence in the Helically Symmetric eXperiment are presented. Using plasma profiles relevant to experimental operation, four dominant drift wave regimes are observed in the ion wavenumber range, which are identified as different flavors of density-gradient-driven trapped electron modes. For the most part, the heat transport exhibits properties associated with turbulence driven by these types of modes. Additionally, long-wavelength, radially localized, nonlinearly excited coherent structures near the resonant central flux surface, not predicted by linear simulations, can further enhance flux levels. Integrated heat fluxes are compatible with experimental observations in the corresponding density gradient range. Despite low shearing rates, zonal flows are observed to regulate turbulence but can be overwhelmed at higher density gradients by the long-wavelength coherent structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Monte-Carlo fluid approaches to detached plasmas in non-axisymmetric divertor configurations

Fluid transport modeling in three-dimensional boundaries of toroidal confinement devices is reviewed with the emphasis on a Monte-Carlo approach to simulate detached plasmas. The loss of axisymmetry in such configurations presents a major challenge for numerical implementation of the standard fluid model widely applied to fusion experimental devices. A large-scale effort has been made to address this problem under complementary aspects including different magnetic topologies and numerical techniques. Here in this paper, we give a brief review of the different strategies pioneered and the challenges involved. A more detailed description is provided for the Monte-Carlo code—EMC3-Eirene, where the physics model and the basic idea behind the applied Monte-Carlo method are presented. The focus is put on its applications to detachment studies for stellarators and tokamaks. Here, major achievements and difficulties encountered are described. Model limitations and further development plans are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Benchmarking gyrokinetic simulations in a toroidal flux-tube

A flux-tube model is implemented in the global turbulence code GEM [Y. Chen and S. E. Parker, J. Comput. Phys. 220, 839 (2007)] in order to facilitate benchmarking with Eulerian codes. The global GEM assumes the magnetic equilibrium to be completely given. The initial flux-tube implementation simply selects a radial location as the center of the flux-tube and a radial size of the flux-tube, sets all equilibrium quantities (B, ∇B, etc.) to be equal to the values at the center of the flux-tube, and retains only a linear radial profile of the safety factor needed for boundary conditions. This implementation shows disagreement with Eulerian codes in linear simulations. An alternative flux-tube model based on a complete local equilibrium solution of the Grad-Shafranov equation [J. Candy, Plasma Phys. Controlled Fusion 51, 105009 (2009)] is then implemented. This results in better agreement between Eulerian codes and the particle-in-cell (PIC) method. The PIC algorithm based on the v{sub ||}-formalism [J. Reynders, Ph.D. dissertation, Princeton University, 1992] and the gyrokinetic ion/fluid electron hybrid model with kinetic electron closure [Y. Chan and S. E. Parker, Phys. Plasmas 18, 055703 (2011)] are also implemented in the flux-tube geometry and compared with the direct method for both the ion temperature gradient driven modes and the kinetic ballooning modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Reinvestigation of the charge density distribution in arc discharge fusion system

A continual arc discharge system has been setup and the light intensity of arc discharge has been profiled. The mathematical model of local energy density distribution in arc discharge fusion has been simulated which is in good qualitative agreement with light intensity profile of arc discharge in the experiments. Eventually, the local energy density distribution of arc discharge system is able to be precisely manipulated to act as heat source in the fabrication of fused fiber devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

ALFVÉN WAVES IN SIMULATIONS OF SOLAR PHOTOSPHERIC VORTICES

Using advanced numerical magneto-hydrodynamic simulations of the magnetized solar photosphere, including non-gray radiative transport and a non-ideal equation of state, we analyze plasma motions in photospheric magnetic vortices. We demonstrate that apparent vortex-like motions in photospheric magnetic field concentrations do not exhibit 'tornado'-like behavior or a 'bath-tub' effect. While at each time instance the velocity field lines in the upper layers of the solar photosphere show swirls, the test particles moving with the time-dependent velocity field do not demonstrate such structures. Instead, they move in a wave-like fashion with rapidly changing and oscillating velocity field, determined mainly by magnetic tension in the magnetized intergranular downflows. Using time-distance diagrams, we identify horizontal motions in the magnetic flux tubes as torsional Alfvén perturbations propagating along the nearly vertical magnetic field lines with local Alfvén speed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

ISS And Space Environment Interactions Without Operating Plasma Contactor

The International Space Station (ISS) will be the largest, highest power spacecraft placed in orbit. Because of this the design of the electrical power system diverged markedly from previous systems. The solar arrays will operate at 160 V and the power distribution voltage will be 120 V. The structure is grounded to the negative side of the solar arrays so under the right circumstances it is possible to drive the ISS potential very negative. A plasma contactor has been added to the ISS to provide control of the ISS structure potential relative to the ambient plasma. The ISS requirement is that the ISS structure not be greater than 40 V positive or negative of local plasma. What are the ramifications of operating large structures with such high voltage power systems? The application of a plasma contactor on ISS controls the potential between the structure and the local plasma, preventing degrading effects. It is conceivable that there can be situations where the plasma contactor might be non-functional. This might be due to lack of power, the need to turn it off during some of the build-up sequences, the loss of functionality for both plasma contactors before a replacement can be installed, similar circumstances. A study was undertaken to understand how important it is to have the contactor functioning and how long it might be off before unacceptable degradation to ISS could occur. The details of interaction effects on spacecraft have not been addressed until driven by design. This was true for ISS. If the structure is allowed to float highly negative impinging ions can sputter exposed conductors which can degrade the primary surface and also generate contamination due to the sputtered material. Arcing has been known to occur on solar arrays that float negative of the ambient plasma. This can also generate electromagnetic interference and voltage transients. Much of the ISS structure and pressure module surfaces exposed to space is anodized aluminum. The anodization thickness is very thin to provide the required solar absorptance and emittance. For conditions where ISS structure can charge negative a large percentage of the array voltage, the dielectric strength of this layer is low, and dielectric breakdown (arcing) can occur. The energy stored capacitively in the structure can be delivered to the arc. The mechanisms by which this energy is delivered and how much of the energy is available hasn't been fully quantified. Questions have been raised regarding the possibility of whether a sustained arc might result due to current collected by the solar arrays from local plasma. It was postulated that even if dielectric breakdown didn't occur, impacts due to micrometeoroids and space debris could penetrate thin layers of dielectric on ISS and initiate an arc due to the coupling provided by the dense local plasma produced by the impact. This was proven in experiments conducted jointly by MSFC and Auburn University. A target chamber with a simulated ionospheric plasma and a biased, anodized aluminum plate and a 1-microfarad capacitor was used. The plate was then impacted by 75-micron particles accelerated to orbital velocity. Arc discharges were sustained for higher voltages but a threshold appears below which no discharge was initiated. Most items without an exposed power system will float electrically near the local plasma potential. This is true of the Space Shuttle, an Astronaut on EVA, and similar items. The structure of ISS might be at a large negative voltage. Therefore, capacitively stored energy can be transferred during docking, installing external boxes and equipment and Astronaut contact with ISS structure. The circumstances of when this can happen and the resulting effects are evaluated in this study. Also, a crewmember on EVA might be in the vicinity of an arc. All safety aspects of such an encounter including charging, molten particles from the arc site and EMI have been evaluated. This paper will report on the total results of this study focussed on the 4A configuration, scheduled to be complete in November, 2000. Interactions such as arcing, debris induced arcs, sustained arcs, sputtering, contamination from sputtering and arcing, docking interactions and Astronaut safety issues will all be addressed.

Carruth, M. R., Jr.↗
DOE OSTI2022

Analysis of wave mode content in fully turbulent, moderately collisional plasma laboratory experiment and kinetic simulation. Final Report

A final report of the activities of Bryn Mawr College for award DE-SC0018258. The 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↗