Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic inductance”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simultaneous Insight into Dissolution and Aggregation of Metal Sulfide Nanoparticles through Single-Particle Inductively Coupled Plasma Mass Spectrometry

Nanoparticles (NPs) and their colloidal aggregates are ubiquitous and play important roles in the transport and release of metals. Knowledge of their dissolution rates and aggregation behavior in solution are crucial for better prediction of their fate in biogeochemical cycling, ecotoxicity, and environmental remediation. There are however significant technical challenges to accurately obtain such information as a result of the heterogeneity and highly dynamic transformation exhibited by NPs, particularly at relatively low particle concentrations in aqueous systems. Here, we quantitatively examine the simultaneous dissolution and aggregation behavior of metal sulfide NPs using single-particle inductively coupled plasma mass spectrometry (spICP–MS). We focus on nickel sulfide (NiS), with additional data presented for copper sulfide (CuS) and cobalt sulfide (CoS). The kinetics of metal release (dissolution and disaggregation) of NiS was fastest under strongly oxidizing conditions (from 0.04 to >3 min –1 with H 2 O 2 ) and were slower under near-neutral (HEPES buffer/H 2 O) and acidic (1 mM HNO 3 ) conditions (≤0.006 min –1 ). Metal release kinetics in HNO 3 was not faster than in H 2 O or HEPES, suggesting that the solution pH has an influence over both the dissolution kinetics of individual particles and the NP aggregation states, which in combination affect metal release rates over time. Between different metal sulfides, the measured metal release rates were largely consistent with predictions based on the crystallinity, solubility products, and specific surface areas of the NPs, following an order of CoS > NiS > CuS. Here, the spICP–MS approach described here can be easily applied to the characterization of metal release and aggregation of other NPs at low concentrations (~10 5 particles/mL) typically found in natural environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Reaction engineering for materials processing in space: Reduction of ilmenite by hydrogen and carbon monoxide

Oxygen is a consumable material which needs to be produced continuously in most space missions. Its use for propulsion as well as life support makes oxygen one of the largest volume chemicals to be produced in space. Production of oxygen from lunar materials is of particular interest and is very attractive possibility. The kinetics and mechanism of reduction of ilmenite by carbon monoxide and hydrogen at 800 to 1100 C were investigated. The temporal profiles of conversion for carbon monoxide have a sigmoidal shape and indicate the presence of three different stages (induction, acceleration, and deceleration) during the reduction reaction. The apparent activation energy decreases from 18 kcal/mole at 10 percent conversion to 10 kcal/mole at 50 percent conversion. The reaction is first order with respect to carbon monoxide under the experimental conditions studied. Both SEM and EDX analysis show that the diffusion of Fe product away from the reaction front and through the TiO2 phase, followed by the nucleation and growth of a separate Fe phase are important steps affecting the process kinetics. The results from hydrogen reduction show that the mechanism of ilmenite reduction by hydrogen is similar to that by carbon monoxide. However, the titanium dioxide can be further reduced by hydrogen at 800 to 1000 C. The detailed comparison and theoretical modeling of both reduction processes is presented.

Zhao, Y.↗
NASA NTRS1989

Radiation effects in Caenorhabditis elegans - Mutagenesis by high and low LET ionizing radiation

The nematode C. elegans was used to measure the effectiveness of high-energy ionized particles in the induction of three types of genetic lesions. Recessive lethal mutations in a 40-map unit autosomal region, sterility, and X-chromosome nondisjunction or damage were investigated. Induction rates were measured as a function of linear energy transfer, LET(infinity), for nine ions of atomic nunmber 1-57 accelerated at the BEVALAC accelerator. Linear kinetics were observed for all three types of lesions within the dose/fluence ranges tested and were found to vary strongly as a function of particle LET(infinity). Relative biological effectiveness (RBE) values of up to 4.2 were measured, and action cross sections were calculated and compared to mutagenic responses in other systems.

Nelson, Gregory A.↗
NASA NTRS1991

Detonability of hydrocarbon fuels in air

Studies were conducted of the detonation of gas-phase mixtures of n-hexane and JP-4, with oxidizers as varied as air and pure oxygen, measuring detonation velocities and cell sizes as a function of stoichiometry and diluent concentration. The induction length of a one-dimensional Zeldovich-von Neumann-Doering detonation was calculated on the basis of a theoretical model that employed the reaction kinetics of the hydrocarbon fuels used. Critical energy and critical tube diameter are compared for a relative measure of the heavy hydrocarbon fuels studied; detonation sensitivity appears to increase slightly with increasing carbon number.

Beeson, H. D.↗
NASA NTRS2005

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗
NASA NTRS1991

Production of oxygen from lunar ilmenite

The kinetics and the mechanism of reduction of synthetic ilmenite by hydrogen in the temperature range of 807 to 1014 C were investigated. At temperatures below 876 C, the temporal profiles of conversion have a sigmoidal shape and indicate the presence of three different stages (induction, acceleration, and deceleration) during the reduction reaction. The apparent activation energy for the reaction is 22.3 kcal/mole, whereas the intrinsic activation energy is 16.9 kcal/mole. Scanning electron microscopy and energy dispersive x-ray analyses show that the diffusion of Fe product away from the reaction front and through the TiO2 phase, followed by the nucleation and growth of a separate Fe phase is an important step affecting the process kinetics. X-ray diffraction and wavelength dispersive x-ray results indicate that the TiO2 can be reduced to lower oxides of titanium at temperatures higher than 876 C.

Shadman, F.↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

R-Value Measurements Performed on Uranium Targets Irradiated with Fission Spectrum Neutrons FY 2021 for F2019 Project

The separation and characterization of two irradiated uranium targets, a depleted uranium (DU) and a highly enriched uranium (HEU) target, was conducted in April of 2021. The two targets were assembled at Los Alamos National Laboratory (LANL) and irradiated using the critical assembly at the National Criticality Experiments Research Center (NCERC). Splits of the dissolved targets were received by Pacific Northwest National Laboratory (PNNL) after which the PNNL and LANL teams chemically separated the solutions using independent separation schemes and analyzed the separated fractions for short lived actinides and fission products. Chemical separations at PNNL were traced with stable or radioactive tracers to allow for the determination of chemical yields, analyzing using either inductively coupled plasma optical emission spectroscopy (ICP-OES), inductively coupled plasma mass spectrometry (ICP-MS) or gamma emission analysis (GEA) depending on the nature of the tracer. Several other analytical techniques were used by PNNL including kinetic phosphorescence analysis (KPA) and thermal ionization mass spectrometry (TIMS) depending on the analyte’s need. Comparisons were made between current and historical PNNL and LANL, as well as literature values. Overall, there was agreement between the two laboratories for the bulk of analytes, with some notable exceptions such as 111 Ag, and 141,143,144 Ce. Included in these comparisons were the short-lived actinides 237 U, 239 Np, the fission products 89 Sr, 91 Y, 95/97 Zr, 99 Mo, 111 Ag, 115/115 mCd, 136/137 Cs, 140 Ba, 141/143/144 Ce, 147 Nd, 153 Sm, 156 Eu, and 161 Tb, providing both total atoms as well as the R-values. The data presented in this report represents the sixth NCERC irradiation of HEU and DU (FY13, FY14, FY15, FY17, FY18, FY21), and their subsequent separation and analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2012

Thrust Stand Measurements of a Conical Inductive Pulsed Plasma Thruster

Inductive Pulsed Plasma Thrusters (iPPT) are spacecraft propulsion devices in which electrical energy is capacitively stored and then discharged through an inductive coil. The thruster is electrodeless, with a time-varying current in the coil interacting with a plasma covering the face of the coil to induce a plasma current. Propellant is accelerated and expelled at a high exhaust velocity (O(10 .. 100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, inductive pulsed plasma thrusters can suffer from both high pulse energy requirements imposed by the voltage demands of inductive propellant ionization, and low propellant utilization efficiencies. A conical coil geometry may o er higher propellant utilization efficiency over that of a at inductive coil, however an increase in propellant utilization may be met with a decrease in axial electromagnetic acceleration, and in turn, a decrease in the total axially-directed kinetic energy imparted to the propellant.

Hallock, Ashley K.↗
NASA NTRS2013

Thrust Stand Measurements of a Conical Inductive Pulsed Plasma Thruster

Inductive Pulsed Plasma Thrusters (iPPT) spacecraft propulsion devices in which electrical energy is capacitively stored and then discharged through an inductive coil. The thruster is electrodeless, with a time-varying current in the coil interacting with a plasma covering the face of the coil to induce a plasma current Propellant is accelerated and expelled at a high exhaust velocity (O(10 -- 100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, inductive pulsed plasma thrusters can suffer from both high pulse energy requirements imposed by the voltage demands of inductive propellant ionization, and low propellant utilization efficiencies. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, inductive pulsed plasma thrusters can suffer from both high pulse energy requirements imposed by the voltage demands of inductive propellant ionization, and low propellant utilization efficiencies. A conical coil geometry may offer higher propellant utilization efficiency over that of a at inductive coil, however an increase in propellant utilization may be met with a decrease in axial electromagnetic acceleration, and in turn, a decrease in the total axially-directed kinetic energy imparted to the propellant.

Hallock, Ashley K.↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Effect of Inductive Coil Geometry and Current Sheet Trajectory of a Conical Theta Pinch Pulsed Inductive Plasma Accelerator

Results are presented demonstrating the e ect of inductive coil geometry and current sheet trajectory on the exhaust velocity of propellant in conical theta pinch pulsed induc- tive plasma accelerators. The electromagnetic coupling between the inductive coil of the accelerator and a plasma current sheet is simulated, substituting a conical copper frustum for the plasma. The variation of system inductance as a function of plasma position is obtained by displacing the simulated current sheet from the coil while measuring the total inductance of the coil. Four coils of differing geometries were employed, and the total inductance of each coil was measured as a function of the axial displacement of two sep- arate copper frusta both having the same cone angle and length as the coil but with one compressed to a smaller size relative to the coil. The measured relationship between total coil inductance and current sheet position closes a dynamical circuit model that is used to calculate the resulting current sheet velocity for various coil and current sheet con gura- tions. The results of this model, which neglects the pinching contribution to thrust, radial propellant con nement, and plume divergence, indicate that in a conical theta pinch ge- ometry current sheet pinching is detrimental to thruster performance, reducing the kinetic energy of the exhausting propellant by up to 50% (at the upper bound for the parameter range of the study). The decrease in exhaust velocity was larger for coils and simulated current sheets of smaller half cone angles. An upper bound for the pinching contribution to thrust is estimated for typical operating parameters. Measurements of coil inductance for three di erent current sheet pinching conditions are used to estimate the magnetic pressure as a function of current sheet radial compression. The gas-dynamic contribution to axial acceleration is also estimated and shown to not compensate for the decrease in axial electromagnetic acceleration that accompanies the radial compression of the plasma in conical theta pinches.

Hallock, Ashley K.↗
NASA NTRS2012

Studies of Thermophysical Properties of Metals and Semiconductors by Containerless Processing Under Microgravity

Electromagnetic levitation under microgravity provides unique opportunities for the investigation of liquid metals, alloys and semiconductors, both above and below their melting temperatures, with minimized disturbances of the sample under investigation. The opportunity to perform such experiments will soon be available on the ISS with the EML payload which is currently being integrated. With its high-performance diagnostics systems EML allows to measure various physical properties such as heat capacity, enthalpy of fusion, viscosity, surface tension, thermal expansion coefficient, and electrical conductivity. In studies of nucleation and solidification phenomena the nucleation kinetics, phase selection, and solidification velocity can be determined. Advanced measurement capabilities currently being studied include the measurement and control of the residual oxygen content of the process atmosphere and a complementary inductive technique to measure thermophysical properties.

Seidel, A.↗
DOE OSTI2022

Time‐Dependent Cation Selectivity of Titanium Carbide MXene in Aqueous Solution

Abstract Electrochemical ion separation is a promising technology to recover valuable ionic species from water. Pseudocapacitive materials, especially 2D materials, are up‐and‐coming electrodes for electrochemical ion separation. For implementation, it is essential to understand the interplay of the intrinsic preference of a specific ion (by charge/size), kinetic ion preference (by mobility), and crystal structure changes. Ti 3 C 2 T z MXene is chosen here to investigate its selective behavior toward alkali and alkaline earth cations. Utilizing an online inductively coupled plasma system, it is found that Ti 3 C 2 T z shows a time‐dependent selectivity feature. In the early stage of charging (up to about 50 min), K + is preferred, while ultimately Ca 2+ and Mg 2+ uptake dominate; this unique phenomenon is related to dehydration energy barriers and the ion exchange effect between divalent and monovalent cations. Given the wide variety of MXenes, this work opens the door to a new avenue where selective ion‐separation with MXene can be further engineered and optimized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
NASA NTRS2005

Dusty-Plasma Particle Accelerator

A dusty-plasma apparatus is being investigated as means of accelerating nanometer- and micrometer-sized particles. Applications for the dusty-plasma particle accelerators fall into two classes: Simulation of a variety of rapidly moving dust particles and micrometeoroids in outer-space environments that include micrometeoroid streams, comet tails, planetary rings, and nebulae and Deposition or implantation of nanoparticles on substrates for diverse industrial purposes that could include hardening, increasing thermal insulation, altering optical properties, and/or increasing permittivities of substrate materials. Relative to prior apparatuses used for similar applications, dusty-plasma particle accelerators offer such potential advantages as smaller size, lower cost, less complexity, and increased particle flux densities. A dusty-plasma particle accelerator exploits the fact that an isolated particle immersed in plasma acquires a net electric charge that depends on the relative mobilities of electrons and ions. Typically, a particle that is immersed in a low-temperature, partially ionized gas, wherein the average kinetic energy of electrons exceeds that of ions, causes the particle to become negatively charged. The particle can then be accelerated by applying an appropriate electric field. A dusty-plasma particle accelerator (see figure) includes a plasma source such as a radio-frequency induction discharge apparatus containing (1) a shallow cup with a biasable electrode to hold the particles to be accelerated and (2) a holder for the substrate on which the particles are to impinge. Depending on the specific design, a pair of electrostatic-acceleration grids between the substrate and discharge plasma can be used to both collimate and further accelerate particles exiting the particle holder. Once exposed to the discharge plasma, the particles in the cup quickly acquire a negative charge. Application of a negative voltage pulse to the biasable electrode results in the initiation of a low-current, high-voltage cathode spot. Plasma pressure associated with the cathode spot as well as the large voltage drop at the cathode spot accelerates the charged particles toward the substrate. The ultimate kinetic energy attained by particles exiting the particle holder depends in part on the magnitude of the cathode spot sheath potential difference, which is proportional to the magnitude of the voltage pulse, and the on the electric charge on the dust. The magnitude of the voltage pulse can be controlled directly, whereas the particle s electric charge can be controlled indirectly by controlling the operating parameters of the plasma apparatus.

Foster, John E.↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗