SEARCH · Search NASA

Results for “Plasma codes”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height and width by taking the separatrix temperature and density calculated from SOLPS as inputs to EPED. An empirical proxy function based on data from a subset of these discharges is then used to infer the ratio of density at the pedestal top to that at the separatrix based on the electron temperature at the divertor target. This allows the density at the top of the pedestal to be provided to EPED as well, so that knowledge of this quantity is no longer required a priori, while SOLPS radial transport coefficients are needed instead. The gas puff rate provided to SOLPS is then a primary model input for this forward model coupling (along with radial transport coefficients, which are kept fixed through each scan in this study). Comparisons with measurements from DIII-D density scans through a range of target conditions find significant pedestal pressure degradation as the density is increased towards detachment in three different divertor configurations, consistent with the experimental data in these ballooning-limited pedestal regimes.

core-edge integration↗

DOE OSTI2023

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

ATOMIC Simulations and Experimental Data for Sodium and Copper Mixtures

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

ATOMIC Simulations and Experimental Data for CaCO3 Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1] and are used to develop a statistical characterization of matrix effects. Our overall suite of simulations includes contains several sets of simulations: training and validation sets of simulations for three and four element mixtures of calcium, carbon, oxygen, and nitrogen (included to account for atmosphere) along with simulations of the individual elements. The 4-element simulations include the mixture of all four elements mentioned and for each of the four individual elements. The 3-element simulations include output for the mixture of calcium, carbon, oxygen and for these three individual elements. The training data were produced using a 600-run design, shown in Figure 1, that varies input parameters temperature (T), electron density (Ne), and proportion of the elements calcium, carbon, oxygen, and nitrogen (Ca; C; O; N) for the 4 element output. The 3-element output includes all parameters except for the proportion of nitrogen. The element proportions (all the variables but T and Ne) sum to one and are unused in the single-element simulations. The validation data was produced with a 80-run design shown in Figure 2. The training and validation simulation outputs for the 4-element simulations for the mixture and for the single element calcium are shown as sample simulations in Figures 3 and 4 respectively. The simulations produce spectra over a range of 190nm - 950nm that roughly mimics the range collected by the SciAps Z-300 LIBS instrument that was used for the experimental data. The measured spectra for a CaCO3 (which may include contribution from Earth's atmosphere) in the experiment is shown in in Figure 5. All files are kept in directories whose names indicate the elemental composition (CaCO3, Ca, C, O, or N), number of elements (3 or 4), and purpose (training, which is not labeled in the file name, or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. [1] J Colgan, EJ Judge, DP Kilcrease, and JE Barefield II. Ab-initio modeling of an iron laser-induced plasma: Comparison between theoretical and experimental atomic emission spectra. Spectrochimica Acta Part B: Atomic Spectroscopy, 97:65{73}, 2014.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗

NASA NTRS1987

Real time animation of space plasma phenomena

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

Jordan, K. F.↗

NASA NTRS1993

Wave generation by contaminant ions

Investigations dealing with the generation of waves by contaminant ions are reported. The studies included the properties of the velocity distribution function of such ions. It was found that it takes about one ion-cyclotron period for the distribution function to transform from a beam to a ring distribution. A linear instability analysis was performed to examine the possibility of wave excitation by an ion beam streaming perpendicular to the ambient magnetic field in an ionospheric type of plasma. A simulation code was developed to study the nonlinear behavior of the waves excited by beam and ring types of ion distributions. The code treats both electrons and ions as particles kinetically to also factor in the possible involvement of lower-hybrid waves.

Singh, N.↗

NASA NTRS1999

Coronal Physics and the Chandra Emission Line Project

With the launch of the Chandra X-ray Observatory, high resolution X-ray spectroscopy of cosmic sources has begun. Early, deep observations of three stellar coronal sources will provide not only invaluable calibration data, but will also give us benchmarks for plasma spectral modeling codes. These codes are to interpret data from stellar coronae, galaxies and clusters of galaxies. supernova remnants and other astrophysical sources, but they have been called into question in recent years as problems with understanding moderate resolution ASCA and EUVE data have arisen. The Emission Line Project is a collaborative effort to improve the models, with Phase 1 being the comparison of models with observed spectra of Capella, Procyon, and HR, 1099. Goals of these comparisons are (1) to determine and verify accurate and robust diagnostics and (2) to identify and prioritize issues in fundamental spectroscopy which will require further theoretical and/or laboratory work. A critical issue in exploiting the coronal data for these purposes is to understand the extent to which common simplifying assumptions (coronal equilibrium, time-independence, negligible optical depth) apply. We will discuss recent advances in our understanding of stellar coronae in this context.

Brickhouse, Nancy↗

NASA NTRS2010

Recent Simulation Results on Ring Current Dynamics Using the Comprehensive Ring Current Model

Plasma sheet conditions and electromagnetic field configurations are both crucial in determining ring current evolution and connection to the ionosphere. In this presentation, we investigate how different conditions of plasma sheet distribution affect ring current properties. Results include comparative studies in 1) varying the radial distance of the plasma sheet boundary; 2) varying local time distribution of the source population; 3) varying the source spectra. Our results show that a source located farther away leads to a stronger ring current than a source that is closer to the Earth. Local time distribution of the source plays an important role in determining both the radial and azimuthal (local time) location of the ring current peak pressure. We found that post-midnight source locations generally lead to a stronger ring current. This finding is in agreement with Lavraud et al.. However, our results do not exhibit any simple dependence of the local time distribution of the peak ring current (within the lower energy range) on the local time distribution of the source, as suggested by Lavraud et al. [2008]. In addition, we will show how different specifications of the magnetic field in the simulation domain affect ring current dynamics in reference to the 20 November 2007 storm, which include initial results on coupling the CRCM with a three-dimensional (3-D) plasma force balance code to achieve self-consistency in the magnetic field.

Zheng, Yihua↗

NASA NTRS2010

On the Role of Global Magnetic Field Configuration in Affecting Ring Current Dynamics

Plasma and field interaction is one important aspect of inner magnetospheric physics. The magnetic field controls particle motion through gradient, curvature drifts and E cross B drift. In this presentation, we show how the global magnetic field affects dynamics of the ring current through simulations of two moderate geomagnetic storms (20 November 2007 and 8-9 March 2008). Preliminary results of coupling the Comprehensive Ring Current Model (CRCM) with a three-dimensional plasma force balance code (to achieve self-consistency in both E and B fields) indicate that inclusion of self-consistency in B tends to mitigate the intensification of the ring current as other similar coupling efforts have shown. In our approach, self-consistency in the electric field is already an existing capability of the CRCM. The magnetic self-consistency is achieved by computing the three-dimensional magnetic field in force balance with anisotropic ring current ion distributions. We discuss the coupling methodology and its further improvement. In addition, comparative studies by using various magnetic field models will be shown. Simulation results will be put into a global context by analyzing the morphology of the ring current, its anisotropy and characteristics ofthe interconnected region 2 field-aligned currents.

Zheng, Y.↗

DOE OSTI2024

Laser-Plasma Interaction Modeling in LLNL Radiation-Hydrodynamic Codes

Explore the source record for details and available documents.

Physics - Plasma physics↗

DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗

DOE OSTI2022

Verification of neoclassical toroidal viscosity induced by energetic particles

The thermal particles contributed neoclassical toroidal viscosity (NTV) have been successfully developed and explored by many impressive works such as the study by Shaing et al. [Phys. Plasmas 10, 1443 (2003)] and Zhu et al. [Phys. Rev. Lett. 96, 225002 (2006)]. In this work, the scope of the NTV study is extended to explore the contribution of energetic particles (EPs) through both theory and experiments. In theory, the existence of the NTV torque due to the precessional drift resonance of trapped EPs is identified based on the equivalence between the NTV torque and the perturbed drift kinetic energy [J. Park, Phys. Plasmas 18, 110702 (2011)]. Toroidal modeling with the Magneto Resistive Spectrum - drift Kinetic code [Y. Liu, Phys. Plasmas 15, 112503 (2008)], based on this equivalence, indicates that trapped EPs can contribute a significant amount of the NTV torque. Meanwhile, this work also focuses on developing the dedicated DIII-D experiments in the presence of the n = 2 external magnetic perturbation to verify the EP induced NTV (EP-NTV) by measuring the change of the NTV torque while varying the angle and the voltage of the neutral beam injection. However, the developed experiments have been unable to create conditions necessary to clearly demonstrate the presence of EP-NTV. The main challenge is separating the resonant and non-resonant momentum transport responses in the plasma. Finally, the experience, gained from this study, can help the further exploration of EP-NTV in the future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1987

A possible explanation of the electron temperature enhancement in the wake of a satellite

A numerical simulation was performed in an attempt to explain the enhancement of electron temperature in the wake of an ionospheric satellite in terms of counterstreaming plasma expansions. The simulation was conducted using full kinetic treatments (Vlasov-Poisson code) for both the plasma electrons and ions. The calculations indicated that at a distance of a few Debye lengths from the surface of the body the wake/ambient electron-temperature ratio reaches a value of about 2. This result is in good agreement with some in situ and laboratory experimental results.

Singh, Nagendra↗

DOE OSTI2020

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗

NASA NTRS1992

The dynamics of low-beta plasma clouds as simulated by a three-dimensional, electromagnetic particle code

A 3D electromagnetic and relativistic particle simulation code is used to investigate the dynamics of low-beta plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. It is found that, in a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. Sheath clouds are accelerated up to relativistic energies, which indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions as energetic sheath electrons scatter in the upper atmosphere.

Neubert, T.↗