Search NASA⌕ Search

SEARCH · Search NASA

Results for “SURFACE DIFFUSION”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Efficient tandem solar cells with solution-processed perovskite on textured crystalline silicon

Stacking solar cells with decreasing band gaps to form tandems presents the possibility of overcoming the single-junction Shockley-Queisser limit in photovoltaics. The rapid development of solution-processed perovskites has brought perovskite single-junction efficiencies >20%. However, this process has yet to enable monolithic integration with industry-relevant textured crystalline silicon solar cells. We report tandems that combine solution-processed micrometer-thick perovskite top cells with fully textured silicon heterojunction bottom cells. To overcome the charge-collection challenges in micrometer-thick perovskites, we enhanced threefold the depletion width at the bases of silicon pyramids. Moreover, by anchoring a self-limiting passivant (1-butanethiol) on the perovskite surfaces, we enhanced the diffusion length and further suppressed phase segregation. These combined enhancements enabled an independently certified power conversion efficiency of 25.7% for perovskite-silicon tandem solar cells. These devices exhibited negligible performance loss after a 400-hour thermal stability test at 85°C and also after 400 hours under maximum power point tracking at 40°C.

14 SOLAR ENERGY↗
DOE OSTI2024

Gain Stability of Hamamatsu R5912-MOD Photo-Multipliers at Low Temperature

The Photon Detection System is a crucial component of the ICARUS detector in the Short Baseline Neutrino (SBN) Program at Fermi National Accelerator Laboratory (FNAL). It consists of 360 Photo Multiplier Tubes (PMTs) 8” Hamamatsu R5912-MOD and, since June 2020, it has been operated at the liquid Argon cryogenic temperature. During these years, the PMTs have shown gain losses as an aging effect. Using a climatic chamber at the INFN Sezione di Catania, we confirmed a stable gain for a single PMT 8” Hamamatsu 5912-MOD when operated at room temperature, and a persistent gain loss at temperatures around -70°C, although far from the liquid Argon one. No gain recovery was obtained bringing the PMT back to room temperature. We suspect that dynode deterioration due to temperature gradient plays a significant role in this phenomenon. During our test, the PMT was illuminated by a 520 nm laser source and operated in current mode. The laser beam was split into two optical fibers, the first one transmitting light through Neutral Density filters and then diffusing it over the surface of the PMT photocathode. The other optical fiber was used for an independent reference measurement of the injected light via a bolometric photodetector. In this presentation, the main technical characteristics of the measurement system are shown together with preliminary results which need further investigation to elucidate the underlying mechanisms driving the gain loss and to propose new mitigation strategies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Property Measurements of LiF-NaF-KF Molten Salts Doped with Surrogate Fission Products

Measurements of molten salt properties including transition temperatures, phase behavior, heat capacity, density, volumetric thermal expansion, surface tension, viscosity thermal diffusivity, thermal conductivity, and vapor pressure are performed at Argonne. The properties of several fuel and coolant salts including eutectic LiF-NaF-KF (FLiNaK) have been measured. These properties are suitable for use in evaluating reactor performance during startup and early operating conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Property Measurements of LiF-NaF-KF Molten Salts Doped with Corrosion Products and Oxygen

Measurements of thermophysical properties of molten salts are needed for modeling and simulation activities that support the development of molten salt reactor (MSR) technologies. Properties of interest including transition temperatures, phase behavior, heat capacity, density, volumetric thermal expansion, surface tension, viscosity, thermal diffusivity, thermal conductivity, and vapor pressure are being performed at Argonne. Results of these property measurements are suitable for use in evaluating reactor performance during startup and the early operating life of the reactor. Ingressions of oxygen and moisture into the fuel salt are expected to occur at different times during the operating life of the reactor due to system leaks, maintenance, and refueling activities. The presence of these environmental contaminants induces corrosion of structural materials. The introduction of corrosion-derived species, oxygen and moisture is expected to affect the physical and chemical properties of the salt and operation of the reactor. Previous work performed at Argonne evaluated the effects of fission product dopants on the thermal properties of eutectic LiF-NaF-KF (FLiNaK). Properties of FLiNaK are commonly used to represent those of fluoride-bearing fuel salts. Metallic corrosion products such as chromium and nickel ions together with dissolved oxygen are expected to affect system redox differently than the accumulation of fission products. Work summarized in this report was performed to measure the effects of corrosion products and dissolved oxygen on the phase transitions, and specific heat capacity of FLiNaK. Thermophysical property measurements were made using four salts that were prepared by doping aliquots of a eutectic mixture of FLiNaK with surrogate corrosion products. Controlled additions of CrF 3 and NiF 2 were used as surrogates for corrosion product contamination anticipated during extended reactor operations in which fuel salt is in contact with steel reactor components. Controlled additions of UF 4 from two sources containing known amounts of UO 2 at different concentrations were used to represent oxygen contamination. The phase transitions and specific heat capacities of the four salts were measured by using differential scanning calorimetry. Measurements were made at temperatures spanning the range of 500−900 °C, which is the expected operating range of MSRs. Measured property values were compared to values measured previously with eutectic FLiNaK without dopants. Differences between property values measured for the doped and non-doped salts were compared with the uncertainties of the measurements to determine the significance of the effect of corrosion products and oxygen on salt properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Validating a microphysical prognostic stratospheric aerosol implementation in E3SMv2 using observations after the Mount Pinatubo eruption

Abstract. This paper describes the addition of a stratospheric prognostic aerosol (SPA) capability – developed with the goal of accurately simulating sulfate aerosol formation and evolution in the stratosphere – in the Department of Energy (DOE) Energy Exascale Earth System Model, version 2 (E3SMv2). The implementation includes changes to the four-mode Modal Aerosol Module microphysics in the stratosphere to allow for larger particle growth and more accurate stratospheric aerosol lifetime following the Pinatubo eruption. E3SMv2-SPA reasonably reproduces stratospheric aerosol lifetime, burden, aerosol optical depth, and top-of-atmosphere flux when compared to remote sensing observations. E3SMv2-SPA also has close agreement with the interactive chemistry–climate model CESM2-WACCM (Community Earth System Model version 2–Whole Atmosphere Community Climate Model) – which has a more complete chemical treatment – and the observationally constrained, prescribed volcanic aerosol treatment in E3SMv2. Global stratospheric aerosol size distributions identify the nucleation and growth of sulfate aerosol from volcanically injected SO2 from both major and minor volcanic eruptions from 1991 to 1993. The modeled aerosol effective radius is consistently lower than satellite and in situ measurements (max differences of ∼ 30 %). Comparisons with in situ size distribution samples indicate that this simulated underestimation in both E3SMv2-SPA and CESM2-WACCM is due to overly small accumulation and coarse-mode aerosols 6–18 months post-eruption, with E3SMv2-SPA simulating ∼ 50 % of the coarse-mode geometric mean diameters of observations 11 months post-eruption. Effective radii from the models and observations are used to calculate offline scattering and absorption efficiencies to explore the implications of smaller simulated aerosol size for the Pinatubo climate impacts. Scattering efficiencies at wavelengths of peak solar irradiance (∼ 0.5 µm) are 10 %–80 % higher for daily samples in models relative to observations through 1993, suggesting higher diffuse radiation at the surface and a larger cooling effect in the models due to the smaller simulated aerosol; absorption efficiencies at the peak wavelengths of outgoing terrestrial radiation (∼ 10 µm) are 15 %–40 % lower for daily samples in models relative to observations, suggesting an underestimation in stratospheric heating in the models due to the smaller simulated aerosol. These potential biases are based on aerosol size alone and do not take into account differences in the aerosol number. The overall agreement of E3SMv2-SPA with observations and its similar performance to the well-validated CESM2-WACCM makes E3SMv2-SPA a viable alternative to simulating climate impacts from stratospheric sulfate aerosols.

Brown, Hunter York (ORCID:0000000218147874)↗
DOE OSTI2023

Facile etching for single crystal cathode materials

A recycling and synthesis of charge material for secondary batteries generates single-crystal charge materials for producing batteries with greater charge cycle longevity. Charge material particles undergo a heating for fusing or enhancing grain boundaries between polycrystalline particles. The resulting, more well-defined grain boundaries are easily etched by a relatively weak mineral acid solution. The acid solution removes material at the grain boundaries to separate secondary particles into primary particles along the grain boundaries. The resulting single crystal (monocrystalline) charge material particles are washed and filtered, and typically re-sintered to accommodate any needed lithium (lithium carbonate), and result in a charge material with larger surface area, higher lithium diffusivity and lower cation ordering.

Wang, Yan↗
DOE OSTI2020

Strong intracellular signal inactivation produces sharper and more robust signaling from cell membrane to nucleus

For a chemical signal to propagate across a cell, it must navigate a tortuous environment involving a variety of organelle barriers. In this work we study mathematical models for a basic chemical signal, the arrival times at the nuclear membrane of proteins that are activated at the cell membrane and diffuse throughout the cytosol. Organelle surfaces within human B cells are reconstructed from soft X-ray tomographic images, and modeled as reflecting barriers to the molecules’ diffusion. We show that signal inactivation sharpens signals, reducing variability in the arrival time at the nuclear membrane. Inactivation can also compensate for an observed slowdown in signal propagation induced by the presence of organelle barriers, leading to arrival times at the nuclear membrane that are comparable to models in which the cytosol is treated as an open, empty region. In the limit of strong signal inactivation this is achieved by filtering out molecules that traverse non-geodesic paths.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Switching CO 2 Electroreduction Selectivity Between C 1 and C 2 Hydrocarbons on Cu Gas-Diffusion Electrodes

We report that regulating the selectivity toward a target hydrocarbon product is still the focus of CO 2 electroreduction. Here, we discover that the original surface Cu species in Cu gas-diffusion electrodes plays a more important role than the surface roughness, local pH, and facet in governing the selectivity toward C 1 or C 2 hydrocarbons. The selectivity toward C 2 H 4 progressively increases, while CH 4 decreases steadily upon lowering the Cu oxidation species fraction. At a relatively low electrodeposition voltage of 1.5 V, the Cu gas-diffusion electrode with the highest Cu δ+ /Cu 0 ratio favors the pathways of hydrogenation to form CH 4 with maximum Faradaic efficiency of 65.4% and partial current density of 228 mA cm -2 at -0.83 V vs RHE. At 2.0 V, the Cu gas-diffusion electrode with the lowest Cu δ+ /Cu 0 ratio prefers C–C coupling to form C 2+ products with Faradaic efficiency topping 80.1% at -0.75 V vs RHE, where the Faradaic efficiency of C 2 H 4 accounts for 46.4% and the partial current density of C 2 H 4 achieves 279 mA cm -2 . This work demonstrates that the selectivity from CH 4 to C 2 H 4 is switchable by tuning surface Cu species composition of Cu gas-diffusion electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗
DOE OSTI2022

Spatiotemporal Visualization and Chemical Identification of the Metal Diffusion Layer at the Electrochemical Interface

The diffusion layer created by transition metal dissolution is ubiquitous at the electrochemical solid-liquid interface and plays a key role in determining electrochemical performance. Tracking the spatiotemporal dynamics of the diffusion layer has remained an unresolved challenge. With spatially resolved synchrotron X-ray fluorescence microscopy and micro-X-ray absorption spectroscopy, we demonstrate the in situ visualization and chemical identification of the dynamic diffusion layer near the electrode surface under electrochemical operating conditions. Finally, our method allows for direct mapping of the reactive electrochemical interface and provides insights into engineering the diffusion layer for improving electrochemical performance.

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

Quasi-equilibrium predictions of water desorption kinetics from rapidly-heated metal oxide surfaces

Controlling sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for pulsed power applications. Despite the short time scales involved, quasi-equilibrium ideas based on transition state theory (TST) and Arrhenius temperature dependence have been widely applied to fit desorption activation free energies. Here, we apply molecular dynamics (MD) simulations in conjunction with equilibrium potential-of-mean-force (PMF) techniques to directly compute the activation free energies (ΔG*) associated with desorption of intact water molecules from Fe 2 O 3 and Cr 2 O 3 (0001) surfaces. The desorption free energy profiles are diffuse, without maxima, and have substantial dependences on temperature and surface water coverage. Incorporating the predicted ΔG* into an analytical form gives rate equations that are in reasonable agreement with non-equilibrium molecular dynamics desorption simulations. We also show that different ΔG* analytical functional forms which give similar predictions at a particular heating rate can yield desorption times that differ by up to a factor of four or more when the ramp rate is extrapolated by 8 orders of magnitude. This highlights the importance of constructing a physically-motivated ΔG* functional form to predict fast desorption kinetics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Reconstructing the reflectivity of liquid surfaces from grazing incidence X-ray off-specular scattering data

The capillary wave model of a liquid surface predicts both the X-ray specular reflection and the diffuse scattering around it. A quantitative method is presented to obtain the X-ray reflectivity (XRR) from a liquid surface through the diffuse scattering data around the specular reflection measured using a grazing incidence X-ray off-specular scattering (GIXOS) geometry at a fixed horizontal offset angle with respect to the plane of incidence. With this approach the entire Q z -dependent reflectivity profile can be obtained at a single, fixed incident angle. This permits a much faster acquisition of the profile than with conventional reflectometry, where the incident angle must be scanned point by point to obtain a Q z -dependent profile. The XRR derived from the GIXOS-measured diffuse scattering, referred to in this paper as pseudo-reflectivity, provides a larger Q z range compared with the reflectivity measured by conventional reflectometry. Transforming the GIXOS-measured diffuse scattering profile to pseudo-XRR opens up the GIXOS method to widely available specular XRR analysis software tools. Here the GIXOS-derived pseudo-XRR is compared with the XRR measured by specular reflectometry from two simple vapor–liquid interfaces at different surface tension, and from a hexadecyltrimethylammonium bromide monolayer on a water surface. For the simple liquids, excellent agreement (beyond 11 orders of magnitude in signal) is found between the two methods, supporting the approach of using GIXOS-measured diffuse scattering to derive reflectivities. Pseudo-XRR obtained at different horizontal offset angles with respect to the plane of incidence yields indistinguishable results, and this supports the robustness of the GIXOS-XRR approach. The pseudo-XRR method can be extended to soft thin films on a liquid surface, and criteria are established for the applicability of the approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The nanocaterpillar's random walk: diffusion with ligand–receptor contacts

Particles with ligand–receptor contacts bind and unbind fluctuating “legs” to surfaces, whose fluctuations cause the particle to diffuse. Quantifying the diffusion of such “nanoscale caterpillars” is a challenge, since binding events often occur on very short time and length scales. Here, in this work, we derive an analytical formula, validated by simulations, for the long time translational diffusion coefficient of an overdamped nanocaterpillar, under a range of modeling assumptions. We demonstrate that the effective diffusion coefficient, which depends on the microscopic parameters governing the legs, can be orders of magnitude smaller than the background diffusion coefficient. Furthermore it varies rapidly with temperature, and reproduces the striking variations seen in existing data and our own measurements of the diffusion of DNA-coated colloids. Our model gives insight into the mechanism of motion, and allows us to ask: when does a nanocaterpillar prefer to move by sliding, where one leg is always linked to the surface, and when does it prefer to move by hopping, which requires all legs to unbind simultaneously? We compare a range of systems (viruses, molecular motors, white blood cells, protein cargos in the nuclear pore complex, bacteria such as Escherichia coli, and DNA-coated colloids) and present guidelines to control the mode of motion for materials design.

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

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

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

Enhancing Four-Carbon Olefin Production from Acetylene over Copper Nanoparticles in Metal–Organic Frameworks

Four-carbon olefins, such as 1-butene and 1,3-butadiene, are important chemical feedstocks for the production of adhesives and synthetic rubber. These compounds are found in the C 4 fraction of “green oil” products that can arise during the hydrogenation of acetylene. Here, we demonstrate that control of the catalyst structure increases the yield and productivity of these important olefins with a family of catalyst materials comprising Cu nanoparticles (CuNPs) bound within the pores of Zr-based metal–organic frameworks. Using carbon monoxide as a probe molecule, we characterize the surfaces of these catalytic CuNPs with diffuse reflectance infrared Fourier transform spectroscopy, revealing that the electronic structure of the CuNP surfaces is size-dependent. Furthermore, we find that as the CuNP diameter decreases, the selectivity for C 4 products increases and that lowering the stoichiometric ratio of H 2 /acetylene improves the selectivity and productivity of the catalyst.

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

Probing the Surface Chemistry of Lithium Nitridation

Chemical synthesis of Li 3 N through lithium nitridation has potential to advance rechargeable battery and nitrogen fixation technology. However, studies of the conditions for forming Li 3 N on the lithium surface via nitrogen gas exposure report contradictory findings, such as the spontaneous reaction of Li with pure N 2 , the impossibility of forming Li 3 N through pure Li and N 2 interaction, the requirement of trace H 2 O to catalyze the reaction, and evidence to the contrary. In this study, ambient pressure X-ray photoelectron spectroscopy (APXPS) was applied to evaluate the in situ chemical evolution of the lithium metal surface under nitrogen gas up to 800 mTorr. At pressures ≤10 mTorr, no Li 3 N was detected. At higher pressures, surface Li 3 N rapidly reacts with trace CO 2 . Additionally, because metallic lithium is readily oxidized by trace gases, the atomic nitrogen concentration of the lithium surface remains below 2%. When nitridation follows oxidation by O 2 gas, CO 2 gas, or H 2 O vapor, surface Li 3 N formation is inhibited. These results suggest that nitrogen gas can diffuse through the oxidized lithium metal surface to react with subsurface metallic lithium.

Binding energy↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
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↗