Search NASA⌕ Search

SEARCH · Search NASA

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,027 records · Page 57

Analyzing large data sets from XGC1 magnetic fusion simulations using apache spark

Apache Spark is explored as a tool for analyzing large data sets from the magnetic fusion simulation code XGCI. Implementation details of Apache Spark on the NERSC Edison supercomputer are discussed, including binary file reading, and parameter setup. Here, an unsupervised machine learning algorithm, k-means clustering, is applied to XGCI particle distribution function data, showing that highly turbulent spatial regions do not have common coherent structures, but rather broad, ring-like structures in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance and Thermal Characterization of the NASA-300MS 20 kW Hall Effect Thruster

NASA's Space Technology Mission Directorate is sponsoring the development of a high fidelity 15 kW-class long-life high performance Hall thruster for candidate NASA technology demonstration missions. An essential element of the development process is demonstration that incorporation of magnetic shielding on a 20 kW-class Hall thruster will yield significant improvements in the throughput capability of the thruster without any significant reduction in thruster performance. As such, NASA Glenn Research Center and the Jet Propulsion Laboratory collaborated on modifying the NASA-300M 20 kW Hall thruster to improve its propellant throughput capability. JPL and NASA Glenn researchers performed plasma numerical simulations with JPL's Hall2De and a commercially available magnetic modeling code that indicated significant enhancement in the throughput capability of the NASA-300M can be attained by modifying the thruster's magnetic circuit. This led to modifying the NASA-300M magnetic topology to a magnetically shielded topology. This paper presents performance evaluation results of the two NASA-300M magnetically shielded thruster configurations, designated 300MS and 300MS-2. The 300MS and 300MS-2 were operated at power levels between 2.5 and 20 kW at discharge voltages between 200 and 700 V. Discharge channel deposition from back-sputtered facility wall flux, and plasma potential and electron temperature measurements made on the inner and outer discharge channel surfaces confirmed that magnetic shielding was achieved. Peak total thrust efficiency of 64% and total specific impulse of 3,050 sec were demonstrated with the 300MS-2 at 20 kW. Thermal characterization results indicate that the boron nitride discharge chamber walls temperatures are approximately 100 C lower for the 300MS when compared to the NASA- 300M at the same thruster operating discharge power.

High power↗

Gyrokinetic simulation study of magnetic island effects on neoclassical physics and micro-instabilities in a realistic KSTAR plasma

Here, we perform gyrokinetic simulations to study the effects of a stationary magnetic island on neoclassical flow and micro-instability in a realistic KSTAR plasma condition. Through the simulations, we aim to analyze a recent KSTAR experiment, which was to measure the details of poloidal flow and fluctuation around a stationary (2, 1) magnetic island [M. J. Choi et al., Nucl. Fusion 57, 126058 (2017)]. From the simulations, it is found that the magnetic island can significantly enhance the equilibrium E x B flow. The corresponding flow shearing is strong enough to suppress a substantial portion of ambient micro-instabilities, particularly ∇T e -driven trapped electron modes. This implies that the enhanced E x B flow can sustain a quasi-internal transport barrier for T e in an inner region neighboring the magnetic island. The enhanced E x B flow has a (2, 1) mode structure with a finite phase shift from the mode structure of the magnetic island. It is shown that the flow shear and the fluctuation suppression patterns implied from the simulations are consistent with the observations on the KSTAR experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhance the terahertz Smith-Purcell superradiant radiation by using dielectric loaded grating

A dielectric loaded grating (DLG) for terahertz Smith-Purcell (SP) device is proposed to enhance the radiation intensity. By using the theoretical analysis and particle-in-cell simulations, the dispersion characteristics and SP superradiant radiation are investigated. Compared with the general metal grating, the usage of DLG can improve the magnitude of electric field and, consequently, strengthen the interaction of the evanescent wave with electron beam, which can improve the growth rate, enhance the SP superradiant radiation, and lower the start current for the operation of SP free-electron laser.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion acceleration in shell cylinders irradiated by a short intense laser pulse

The interaction of a short high intensity laser pulse with homo and heterogeneous shell cylinders has been analyzed using particle-in-cell simulations and analytical modeling. We show that the shell cylinder is proficient of accelerating and focusing ions in a narrow region. In the case of shell cylinder, the ion energy exceeds the ion energy for a flat target of the same thickness. The constructed model enables the evaluation of the ion energy and the number of ions in the focusing region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Avalanche transport of energetic-ions in magnetic confinement plasmas: nonlinear multiple wave-number simulation

Large burst activity, identified as toroidal Alfvén eigenmode (TAE) avalanche, occurs frequently in neutral-beam heated plasmas in National Spherical Torus Experiment (NSTX). Based on the typical experimental observation of TAE avalanche on NSTX, a self-consistent nonlinear multiple wave-number (k ∥ ≃ n/R, where n toroidal mode-number and R major radius) simulation associated with TAE avalanches is performed using the experimental parameters and profiles before the occurrence of TAE avalanche as the M3D-K input. The wave–wave nonlinear coupling among different modes and the resonant interaction between different modes and energetic-ions during TAE avalanches are identified in the nonlinear multiple wave-number simulations. The resonance overlap during the TAE avalanche is clearly observed in the simulation. It is found that the effective wave–wave coupling and a sufficiently strong drive are two important ingredients for the onset of TAE avalanches. TAE avalanche is considered to be a strongly nonlinear process and it is always accompanied by the simultaneous rapid frequency-chirping and large amplitude bursting of multiple modes and significant energetic-ion losses. The experimental phenomenon is observed on NSTX and is qualitatively reproduced by the simulation results in this work. These findings indicate that the onset of avalanche is triggered by nonlinearity of the system, and are also conducive to understanding the underlying mechanism of avalanche transport of energetic particles in the future burning plasmas, such as International Thermonuclear Experiment Reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of Solar Array Designs on High Voltage Operations

As power levels of advanced spacecraft climb above 25 kW, higher solar array operating voltages become attractive. Even in today s satellites, operating spacecraft buses at 100 V and above has led to arcing in GEO communications satellites, so the issue of spacecraft charging and solar array arcing remains a design problem. In addition, micrometeoroid impacts on all of these arrays can also lead to arcing if the spacecraft is at an elevated potential. For example, tests on space station hardware disclosed arcing at 75V on anodized A1 structures that were struck with hypervelocity particles in Low Earth Orbit (LEO) plasmas. Thus an understanding of these effects is necessary to design reliable high voltage solar arrays of the future, especially in light of the Vision for Space Exploration of NASA. In the future, large GEO communication satellites, lunar bases, solar electric propulsion missions, high power communication systems around Mars can lead to power levels well above 100 kW. As noted above, it will be essential to increase operating voltages of the solar arrays well above 80 V to keep the mass of cabling needed to carry the high currents to an acceptable level. Thus, the purpose of this paper is to discuss various solar array approaches, to discuss the results of testing them at high voltages, in the presence of simulated space plasma and under hypervelocity impact. Three different types of arrays will be considered. One will be a planar array using thin film cells, the second will use planar single or multijunction cells and the last will use the Stretched Lens Array (SLA - 8-fold concentration). Each of these has different approaches for protection from the space environment. The thin film cell based arrays have minimal covering due to their inherent radiation tolerance, conventional GaAs and multijunction cells have the traditional cerium-doped microsheet glasses (of appropriate thickness) that are usually attached with Dow Corning DC 93-500 silicone adhesive. In practice, these cover glasses and adhesive do not cover the cell edges. Finally, in the SLA, the entire cell and cell edges are fully encapsulated by a cover glass that overhangs the cell perimeter and the silicone adhesive covers the cell edges providing a sealed environment. These three types of blanket technology have been tested at GRC and Auburn. The results of these tests will be described. For example, 15 modules composed of four state-of-the-art 2x4 cm GaAs solar cells with 150 pm cover glasses connected in two-cell series strings were tested at high voltage, in plasma under hypervelocity impact. A picture of one of the modules is shown in figure 1. These were prepared by standard industry practice from a major supplier and had efficiencies above 18%. The test results and other fabrication factors that influenced the tests will be presented. In addition, results for SLA segments tested under the same conditions will be presented. Testing of thin film blankets at GRC will also be presented. Figure 1 : Typical GaAs Solar Cell Module These results will show significant differences in resistance to arcing that are directly related to array design and manufacturing procedures. Finally, the approaches for mitigating the problems uncovered by these tests will be described. These will lay the foundation for future higher voltage array operation, even including voltages above 300-600 V for direct drive SEP applications.

Brandhorst, Henry W., Jr.↗

Retarding potential analyzer for the Pioneer-Venus Orbiter Mission

The retarding potential analyzer on the Pioneer-Venus Orbiter Mission has been designed to measure most of the thermal plasma parameters within and near the Venusian ionosphere. Parameters include total ion concentration, concentrations of the more abundant ions, ion temperatures, ion drift velocity, electron temperature, and low-energy (0-50 eV) electron distribution function. To accomplish these measurements on a spinning vehicle with a small telemetry bit rate, several functions, including decision functions not previously used in RPA's, have been developed and incorporated into this instrument. The more significant functions include automatic electrometer ranging with background current compensation; digital, quadratic retarding potential step generation for the ion and low-energy electron scans; a current sampling interval of 2 ms throughout all scans; digital logic inflection point detection and data selection; and automatic ram direction detection. Extensive numerical simulation and plasma chamber tests have been conducted to verify adequacy of the design for the Pioneer Mission.

Knudsen, W. C.↗

Influence of axial self-magnetic field component on arcing behavior of spiral-shaped contacts

The transverse magnetic field (TMF) contact design is commonly used in vacuum interrupters. When arcing occurs between the TMF contacts, the contact structure can create a self-induced magnetic field that drives the arc to move and rotate on the contact, and thus local overheating and severe erosion can be avoided. However, TMF contacts could also create an axial self-magnetic component, and the influence of this component on the arc behavior has not been considered to date. In this paper, five different types of Cu-Cr spiral-shaped TMF contacts with three different structures are investigated in a demountable vacuum chamber that contains a high-speed charge-coupled device video camera. It was found that the contact structure greatly influenced the arc behavior, especially in terms of arc rotation and the effective contact area, while contacts with the same slot structure but different diameters showed similar arc behavior and arc motion. The magnetic field distribution and the Lorentz force of each of the three different contact structures are simulated, and the axial self-magnetic field was first taken into consideration for investigation of the TMF contact design. It was found that contact designs that have higher axial self-magnetic field components tend to have arc columns with larger diameters and show poorer arc motion and rotation performance in the experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Integrated System and Application Continuous Performance Monitoring and Analysis Capability

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia's Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia's ElectroMagnetic Plasma In Realistic Environments ( EMPIRE ) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia's HPC users.

97 MATHEMATICS AND COMPUTING↗

Integrated System and Application Continuous Performance Monitoring and Analysis Capability (Final)

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia’s Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia’s ElectroMagnetic Plasma In Realistic Environments (EMPIRE) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia’s HPC users.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear evolution of the Kelvin-Helmholtz instability in the double current sheet configuration

The nonlinear evolution of the Kelvin-Helmholtz (KH) instability driven by a radially antisymmetric shear flow in the double current sheet configuration is numerically investigated based on a reduced magnetohydrodynamic model. Simulations reveal different nonlinear fate of the KH instability depending on the amplitude of the shear flow, which restricts the strength of the KH instability. For strong shear flows far above the KH instability threshold, the linear electrostatic-type KH instability saturates and achieves a vortex flow dominated quasi-steady state of the electromagnetic (EM) KH turbulence with large-amplitude zonal flows as well as zonal fields. The magnetic surfaces are twisted significantly due to strong vortices but without the formation of magnetic islands. However, for the shear flow just over the KH instability threshold, a weak EM-type KH instability is saturated and remarkably damped by zonal flows through modifying the equilibrium shear flow. Interestingly, a secondary double tearing mode (DTM) is excited subsequently in highly damped KH turbulence, behaving as a pure DTM in a flowing plasma as described in Mao et al. [Phys. Plasmas 21, 052304 (2014)]. However, the explosive growth phenomenon is replaced by a gradually growing oscillation due to the extremely twisted islands. As a result, the release of the magnetic energy becomes slow and the global magnetic reconnection tends to be gentle. A complex nonlinear interaction between the EM KH turbulence and the DTMs occurs for the medium shear flows above the KH instability threshold, turbulent EM fluctuations experience oscillatory nonlinear growth of the DTMs, finally achieves a quasi-steady state with the interplay of the fluctuations between the DTMs and the EM KH instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the control of electron heating for optimal laser radiation pressure ion acceleration

We study the onset of electron heating in intense laser–solid interactions and its impact on the spectral quality of radiation pressure accelerated ions in both hole boring and light sail regimes. Two- and three-dimensional particle-in-cell (PIC) simulations are performed over a wide range of laser and target parameters and reveal how the pulse duration, profile, polarization and target surface stability control the electron heating, the dominant ion acceleration mechanisms and the ion spectra. We find that the onset of strong electron heating is associated with the growth of the Rayleigh–Taylor-like instability at the front surface and must be controlled to produce high-quality ion beams, even when circularly polarized lasers are employed. We define a threshold condition for the maximum duration of the laser pulse that allows mitigation of electron heating and radiation pressure acceleration of narrow energy spread ion beams. The model is validated by three-dimensional PIC simulations, and the few experimental studies that reported low energy spread radiation pressure accelerated ion beams appear to meet the derived criteria. The understanding provided by our work will be important in guiding future experimental developments, for example for the ultrashort laser pulses becoming available at state-of-the-art laser facilities, for which we predict that proton beams with $\sim$ 150–250 MeV, $\sim$ 30% energy spread, and a total laser-to-proton conversion efficiency of $\sim$ 20% can be produced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fiber-optic bolometers with high-temperature tolerance and reduced time constants for fusion plasma diagnostics

Fiber-optic bolometers (FOBs) designed for plasma radiation measurement in magnetically confined fusion environments have been previously developed and tested at the DIII-D tokamak. These FOBs utilize a silicon Fabry–Perot interferometer at the fiber tip for temperature measurement and a gold disk as a radiation absorber. This paper presents recent advancements in the development of FOBs with reduced cooling time constants and enhanced temperature tolerance. We also demonstrate high-temperature operation and report findings from irradiation tests. By incorporating a copper sensor holder as a thermal sink and integrating two closely packed sensor heads into a compact design, the cooling time constants were reduced to ∼100 ms in vacuum, while minimizing uncertainties caused by temperature gradients. The FOBs were subjected to heating at 400 °C for five hours, with no observed degradation in noise performance post-heating. High-temperature operation tests were conducted at ∼220 °C using a 404 nm square-wave modulated laser to simulate the plasma radiation. In addition, four FOBs constructed from copper-coated pure-silica fibers were fabricated and underwent irradiation testing, enduring a total γ-radiation dose of 15.2 MGy and a neutron fluence of 1.6 × 10 18 n/cm 2 . As a result, the FOBs maintained their physical integrity and the optical signal shows high visibility throughout the test.

42 ENGINEERING↗