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 433 records · Page 24

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

Some features of inverted-V events as seen from simulated double layers

Results from a numerical simulation of a double layer show some features similar to those of inverted-V events. The strong heating of thermal and precipitating electrons is observed along with extremely low frequency fluctuations found during inverted-V events. It is suggested that after the acceleration of auroral electrons by the double layer, the precipitating free electrons are heated by the nonlinear effects of the electron beam plasma instability. Fluctuations and pulsations of auroral electron fluxes during auroral events are caused by a relaxation type of oscillation. The finite extent of one dimensional plasma is simulated by solving the Vlasov and Poisson equations as an initial and boundary value problem.

Singh, N.↗

Magnetic Flux Compression Concept for Nuclear Pulse Propulsion and Power

The desire for fast, efficient interplanetary transport requires propulsion systems having short acceleration times and very high specific impulse attributes. Unfortunately, most highly efficient propulsion systems which are within the capabilities of present day technologies are either very heavy or yield very low impulse such that the acceleration time to final velocity is too long to be of lasting interest, One exception, the nuclear thermal thruster, could achieve the desired acceleration but it would require inordinately large mass ratios to reach the range of desired final velocities. An alternative approach, among several competing concepts that are beyond our modern technical capabilities, is a pulsed thermonuclear device utilizing microfusion detonations. In this paper, we examine the feasibility of an innovative magnetic flux compression concept for utilizing microfusion detonations, assuming that such low yield nuclear bursts can be realized in practice. In this concept, a magnetic field is compressed between an expanding detonation driven diamagnetic plasma and a stationary structure formed from a high temperature superconductor (HTSC). In general, we are interested in accomplishing two important functions: (1) collimation of a hot diamagnetic plasma for direct thrust production; and (2) pulse power generation for dense plasma ignition. For the purposes of this research, it is assumed that rnicrofusion detonation technology may become available within a few decades, and that this approach could capitalize on recent advances in inertial confinement fusion ICF) technologies including magnetized target concepts and antimatter initiated nuclear detonations. The charged particle expansion velocity in these detonations can be on the order of 10 (exp 6)- 10 (exp 7) meters per second, and, if effectively collimated by a magnetic nozzle, can yield the Isp and the acceleration levels needed for practical interplanetary spaceflight. The ability to ignite pure fusion micro-bursts with reasonable levels of input energy is an equally challenging scientific problem. It remains to be seen, however, whether an effective ignition driver can be developed which meets the requirements for practical spaceflight application (namely high power density, compactness, low weight, and low cost). In this paper, system level performance and design issues are examined including generator performance, magnetic flux compression processes, magnetic diffusion processes, high temperature superconductor (HTSC) material properties, plasmadynamic processes, detonation plasma expansion processes, magnetohydrodynamic instabilities, magnetic nozzle performance, and thrust production performance. Representative generator performance calculations based on a simplified skin layer formulation are presented as well as the results of exploratory small-scale laboratory experiments on magnetic flux diffusion in HTSC materials. In addition, planned follow-on scientific feasibility experiments are described which utilize high explosive detonations and high energy gas discharges to simulate the plasma conditions associated with thermonuclear micro-detonations.

Litchford, Ronald J.↗

Effects of eletron heating on the current driven electrostatic ion cyclotron instability and plasma transport processes along auroral field lines

Fluid simulations of the plasma along auroral field lines in the return current region have been performed. It is shown that the onset of electrostatic ion cyclotron (EIC) related anomalous resistivity and the consequent heating of electrons leads to a transverse ion temperature that is much higher than that produced by the current driven EIC instability (CDICI) alone. Two processes are presented for the enhancement of ion heating by anomalous resistivity. The anomalous resistivity associated with the turbulence is limited by electron heating, so that CDICI saturates at transverse temperature that is substantially higher than in the absence of resistivity. It is suggested that this process demonstrates a positive feedback loop in the interaction between CDICI, anomalous resistivity, and parallel large-scale dynamics in the topside ionosphere.

Ganguli, Supriya B.↗

Characterization of uncertainties in electron-argon collision cross sections

Abstract The predictive capability of a plasma discharge model depends on accurate representations of electron-impact collision cross sections, which determine the corresponding reaction rates and electron transport properties. The values of cross sections can be known only approximately either through experiments or simulations and are thus subject to uncertainties. Quantifying the uncertainties in plasma simulations allows us to assess the reliability of simulations and to provide a basis for interpreting discrepancies between simulations and experiments. For such uncertainty quantification of plasma simulations, it is essential to quantify the uncertainties of the underlying cross sections. Although much effort has been committed to calibrate the cross section values, their uncertainties are not well investigated. We characterize uncertainties in electron-argon atom collision cross sections using a Bayesian framework. Six collision processes—elastic momentum transfer, ionization, and four excitations—are characterized with semi-empirical models, which effectively capture the features important to the macroscopic properties of the plasma. A probability model for the uncertain parameters of these semi-empirical models is developed. Specifically, a Gaussian-process likelihood model is proposed to capture discrepancies among data sets, as well as the model-form inadequacies of the semi-empirical models. Two other likelihood models are compared with the proposed Gaussian-process model, to illustrate the importance of the choice of the likelihood model. The cross section models are calibrated using the electron-beam experiments and ab-inito quantum simulations. The resulting calibrated uncertainties capture well the scattering among the data sets. The calibrated cross section models are further validated against swarm-parameter experiments and zero-dimensional Boltzmann equation simulations of widely used cross section datasets.

Chung, Seung Whan (ORCID:0000000302501549)↗

Review of Leading Approaches for Mitigating Hypersonic Vehicle Communications Blackout and a Method of Ceramic Particulate Injection Via Cathode Spot Arcs for Blackout Mitigation

Vehicles flying at hypersonic velocities within the atmosphere become enveloped in a "plasma sheath" that prevents radio communication, telemetry, and most importantly, GPS signal reception for navigation. This radio "blackout" period has been a problem since the dawn of the manned space program and was an especially significant hindrance during the days of the Apollo missions. An appropriate mitigation method must allow for spacecraft to ground control and ground control to spacecraft communications through the reentry plasma sheath. Many mitigation techniques have been proposed, including but not limited to, aerodynamic shaping, magnetic windows, and liquid injection. The research performed on these mitigation techniques over the years will be reviewed and summarized, along with the advantages and obstacles that each technique will need to overcome to be practically implemented. A unique approach for mitigating the blackout communications problem is presented herein along with research results associated with this method. The novel method involves the injection of ceramic metal-oxide particulate into a simulated reentry plasma to quench the reentry plasma. Injection of the solid ceramic particulates is achieved by entrainment within induced, energetic cathode spot flows.

Gillman, Eric D.↗

Assessment of Pole Erosion in a Magnetically Shielded Hall Thruster

Numerical simulations of a 6-kW laboratory Hall thruster called H6 have been performed to quantify the erosion rate at the inner pole. The assessments have been made in two versions of the thruster, namely the unshielded (H6US) and magnetically shielded (H6MS) configurations. The simulations have been performed with the 2-D axisymmetric code Hall2De which employs a new multi-fluid ion algorithm to capture the presence of low-energy ions in the vicinity of the poles. It is found that the maximum computed erosion rate at the inner pole of the H6MS exceeds the measured rate of back-sputtered deposits by ~4.5 times. This explains only part of the surface roughening that was observed after a 150-h wear test, which covered most of the pole area exposed to the plasma. For the majority of the pole surface the computed erosion rates are found to be below the back-sputter rate and comparable to those in the H6US which exhibited little to no sputtering in previous tests. Possible explanations for the discrepancy are discussed.

plasma simulations↗

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preface: C/NOFS Results and Equatorial Ionospheric Dynamics

The Communication/Navigation Outage Forecasting System (C/NOFS) satellite was launched into orbit in April 2008 as part of an ongoing effort to understand and identify plasma irregularities that adversely impact the propagation of radio waves in the upper atmosphere. Combined with recent improvements in radar, airglow, and ground-based studies, as well as state-of-the-art modeling techniques, the C/NOFS mission has led to new insights into equatorial ionospheric electrodynamics. In order to document these advances, the C/NOFS Results and Equatorial Dynamics Technical Interchange Meeting was held in Albuquerque, New Mexico from 12 to 14 March 2013. The meeting was a great success with 55 talks and 22 posters, and covered topics including the numerical simulations of plasma irregularities, the effects of atmospheric tides, stratospheric phenomena, and magnetic storms on the upper atmosphere, causes and predictions of scintillation-causing ionospheric irregularities, current and future instrumentation efforts in the equatorial region. The talks were broken into the following three topical sessions: A. Ambient Ionosphere and Thermosphere B. Transient Phenomena in the Low-Latitude Ionosphere C. New Missions, New Sensors, New Science and Engineering Issues. The following special issue was planned as a follow-up to the meeting. We would like to thank Mike Pinnock, the editors and staff of Copernicus, and our reviewers for their work in bringing this special issue to the scientific community. Our thanks also go to Patricia Doherty and the meeting organizing committee for arranging the C/NOFS Technical Interchange Meeting.

CNOFS↗

Reconstruction of a Broadband Spectrum of Alfvenic Fluctuations

Alfvenic fluctuations in the solar wind exhibit a high degree of velocities and magnetic field correlations consistent with Alfven waves propagating away and toward the Sun. Two remarkable properties of these fluctuations are the tendencies to have either positive or negative magnetic helicity (-1 less than or equal to sigma(sub m) less than or equal to +1) associated with either left- or right- topological handedness of the fluctuations and to have a constant magnetic field magnitude. This paper provides, for the first time, a theoretical framework for reconstructing both the magnetic and velocity field fluctuations with a divergence-free magnetic field, with any specified power spectral index and normalized magnetic- and cross-helicity spectrum field fluctuations for any plasma species. The spectrum is constructed in the Fourier domain by imposing two conditions-a divergence-free magnetic field and the preservation of the sense of magnetic helicity in both spaces-as well as using Parseval's theorem for the conservation of energy between configuration and Fourier spaces. Applications to the one-dimensional spatial Alfvenic propagation are presented. The theoretical construction is in agreement with typical time series and power spectra properties observed in the solar wind. The theoretical ideas presented in this spectral reconstruction provide a foundation for more realistic simulations of plasma waves, solar wind turbulence, and the propagation of energetic particles in such fluctuating fields.

waves↗

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↗

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↗

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.↗

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↗

Cooper minima of atoms in a plasma environment

Recent efforts have extended the calculations of photon ionization and associated Cooper minima to atoms in a plasma environment simulated by screening the atomic potential. Furthermore, in this paper previously unrecognized features of bound electron photon excitations in screened atomic potentials are examined. First, the Cooper minima location as a function of photon energy is shown to be weakly dependent on the screening, which is explained by stationary perturbation theory. Second, the screening transfers transition probability from high to low photon energies relative to the isolated atom result.

Cooper minima↗

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗