Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid Breakup”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multiple Environmental Influences on the Lightning of Cold-Based Continental Cumulonimbus Clouds. Part I: Description and Validation of Model

In this two-part paper, influences from environmental factors on lightning in a convective storm are assessed with a model. In Part I, an electrical component is described and applied in the Aerosol–Cloud model (AC). AC treats many types of secondary (e.g., breakup in ice–ice collisions, raindrop-freezing fragmentation, rime splintering) and primary (heterogeneous, homogeneous freezing) ice initiation. AC represents lightning flashes with a statistical treatment of branching from a fractal law constrained by video imagery. The storm simulated is from the Severe Thunderstorm Electrification and Precipitation Study (STEPS; 19/20 June 2000). The simulation was validated microphysically [e.g., ice/droplet concentrations and mean sizes, liquid water content (LWC), reflectivity, surface precipitation] and dynamically (e.g., ascent) in our 2017 paper. Predicted ice concentrations (~10 L -1 ) agreed—to within a factor of about 2—with aircraft data at flight levels (-10° to -15°C). Here, electrical statistics of the same simulation are compared with observations. Flash rates (to within a factor of 2), triggering altitudes and polarity of flashes, and electric fields, all agree with the coincident STEPS observations. The “normal” tripole of charge structure observed during an electrical balloon sounding is reproduced by AC. It is related to reversal of polarity of noninductive charging in ice–ice collisions seen in laboratory experiments when temperature or LWC are varied. Positively charged graupel and negatively charged snow at most midlevels, charged away from the fastest updrafts, is predicted to cause the normal tripole. Total charge separated in the simulated storm is dominated by collisions involving secondary ice from fragmentation in graupel–snow collisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2025

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

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

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

SAF Combustion & Soot Processes [Slides]

The atomization, mixing, combustion and emissions characteristics of aviation fuels were measured using a novel approach based on a non-continuous injection scheme called the single-hole atomizer (SHA). High-speed microscopy revealed differences between fuels in terms of evaporation and mixing regimes over conditions relevant to modern and next generation aero-engine combustors. Measurements of liquid and vapor penetration, mixing fields, combustion and emissions metrics (ignition delay, lift-off length, PAH formation, soot mass) highlighted the effects of fuels and combustor conditions. The experimental results are being leveraged to adjust and validate chemical and CFD models. Detailed analysis of sampled soot showed subtle differences in soot morphology between fuels. The results revealed the presence of contaminants potentially affecting surface chemistry and the nucleation propensity of water droplets on particles. Chemical mechanisms for NJFCP A-2, C-1 and C-4 showed good performance over a large parameter space. Spray breakup at relight conditions is vastly different from the atomization observed at high pressure. CFD simulations of the SHA target conditions confirmed the good behavior of the C-1 kinetic mechanism. The simulations support the strong relationship between low and high temperature reactions. New altitude chamber facility to enable detailed characterization of the heterogeneous nucleation process of water on aerosol particles.

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

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗
DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Measurement of Sub-GeV $ν_μ$ Charged-Current Coherent Pion Production on Argon in MicroBooNE

We report a measurement of the charged-current coherent pion production cross section on argon using the MicroBooNE liquid argon time projection chamber exposed to the Booster Neutrino Beam at Fermilab. The measurement uses the MicroBooNE data set corresponding to $1.26 \times 10^{21}$ protons on target with a mean neutrino energy of $0.8$~GeV. The flux-averaged cross section is measured to be $(9.1 \pm 1.2_{\text{stat}} \pm 1.2_\text{syst}) \times 10^{-40}\,\text{cm}^2/\text{Ar}$. This result represents the first measurement of charged-current coherent pion production on argon at sub-GeV neutrino energies. Due to its clean two-body kinematics, where the neutrino interacts coherently with the entire nucleus producing a forward muon and pion with no nuclear breakup, this process provides a useful tool for constraining neutrino flux uncertainties in current and future oscillation experiments such as DUNE.

Abratenko, P. [Tufts U.] (ORCID:0000000169455941)↗
DOE OSTI2022

High-repetition-rate, multi-MeV deuteron acceleration from converging heavy water microjets at laser intensities of 10 21 W/cm 2

Herein we demonstrate high repetition-rate deuteron acceleration by irradiating a continuously flowing, ambient temperature liquid heavy water jet with the high-intensity ALEPH laser. The laser delivered up to 5.5 J (120 TW, 1.2 × 10 21 W/cm 2 ) laser energy on target at 0.5 Hz. A high repetition-rate Thomson parabola spectrometer measured the deuteron beam energy spectra on each shot for 60 sequential shots (two minutes). Peak fluxes of 5 x 10 10 deuterons/sr/pulse, corresponding to an average flux of 1.5 x 10 12 deuterons/sr/min, were demonstrated with deuteron energies reaching up to 4.4 MeV. High shot-to-shot stability is observed up to 40%–50% of the maximum deuteron energy. These deuteron beams are suited for fast neutron production through deuteron breakup in a converter yielding energies similar to deuteron–deuteron (D–D, 2.45 MeV) fusion reactions of importance for material damage studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗
DOE OSTI2022

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗
DOE OSTI2024

Are Marine Low Cloud Droplet Concentrations Buffered by Entrained Aitken‐Mode Aerosol (Final Technical Report)

During the summertime, the high-latitude oceans come to life with green phytoplankton, which gain their energy from sunlight and are food for sea creatures small and large. When the phytoplankton are eaten or die, sulfur-rich gases are released under the ocean surface and mix into the air. Observations suggest that, over the Southern Ocean, frequent storms lift this air high into the atmosphere while raining out particulates like salt. As a result, the sulfur-rich air then spawns high concentrations of small ‘Aitken-mode’ aerosol particles. We hypothesize that these particles work their way down into the marine boundary layer, where they can replenish the supply of cloud-condensation nuclei scavenged by frequent precipitation. Further, this process maintains high concentrations of liquid cloud droplets in austral summer, promoting more sunlight to be reflected to space. We call this ‘Aitken buffering’. The primary objective of this project has been to document and test what role Aitken-mode aerosols play in clouds over the Southern Ocean and elsewhere. This effort has included three main components: 1) developing a computer model that realistically simulates the aerosol processes and the small-scale turbulent air motions that move aerosols around and create the clouds, 2) using that model to interpret and extend these observations for process understanding, by allowing different factors that contribute to the aerosol budget, such as surface wind speed, precipitation, surface gas exchange, etc. to be separated, and 3) studying Aitken-mode aerosol and its variability with a focus over the Southern Ocean and Antarctica, using data from Atmospheric Radiation Measurement (ARM) sites and other available observations. Initial computer studies in more idealized conditions found that elevated concentrations of Aitken-mode aerosols above the clouds could help prevent the breakup of those clouds by acting as cloud-condensation nuclei after they were entrained into the cloudy boundary layer. Simulations of a day during the ACE-ENA field campaign showed that Aitken-mode aerosols could also prevent cloud breakup under more realistic conditions. A new method that extracts information about Aitken-mode aerosols from measurements of aerosols onto which cloud droplets can form finds that Aitken-mode aerosols do vary seasonally over the Southern Ocean, with a peak in summertime, as described above. Other work during this project has focused on understanding how patterns of water vapor, clouds and precipitation are coupled within low-lying clouds over the oceans, and also on how cloud droplets cluster within clouds and how the distribution of cloud droplet sizes change as dry air is mixed into clouds, with the latter studies also using observations from ACE-ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Time-resolved spray characterization via unified optical flow and binarization technique

This work leverages an unsupervised machine learning and advanced image processing techniques to characterize the breakup of fuel sprays in a small-scale combustor under reacting conditions, providing valuable insights into near-nozzle flow phenomenology. The proposed methodology integrates an improved optical flow model on a convolutional neural network to extract flow vectors with a binarization technique to assess droplets’ size and shape across the region of interest. The velocimetry approach demonstrates superior performance compared to a state-of-the-art optical flow model when applied to high-speed X-ray phase contrast spray images, achieving more accurate and reliable flow predictions. Moreover, breakup processes are quantified by breakup length and sphericity in accordance with velocity estimations, allowing a more complete characterization of the flow. This study establishes a robust methodology for analyzing spray morphology and primary breakup in compact combustors, contributing valuable means of understanding and optimizing fuel spray behavior in advanced combustion systems.

42 ENGINEERING↗
DOE OSTI2026

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗
DOE OSTI2022

Cosmogenic production of 37 Ar in the context of the LUX-ZEPLIN experiment

We estimate the amount of 37 Ar produced in natural xenon via cosmic-ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth’s surface. We then calculate the resulting 37 Ar concentration in a 10-tonne payload (similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage, and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea-level production rate of 37 Ar in natural xenon is estimated to be 0.024 atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1 tonne/month, the average 37 Ar activity after 10 tons are purified and transported underground is 0.058 − 0.090 𝜇⁢Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic 37 Ar will appear as a noticeable background in the early science data, while decaying with a 35-day half-life. This newly noticed production mechanism of 37 Ar should be considered when planning for future liquid-xenon-based experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗