Numerical simulation of a methane-oxygen rotating detonation rocket engine
Not provided.
Energy & Fuels↗
SEARCH · Search NASA
Results for “numerical mixing”
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.
Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling
Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field
Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.
09 BIOMASS FUELS↗
Stabilization of Low Valent Zirconium Nitrides in Titanium Nitride via Plasma-Enhanced Atomic Layer Deposition and Assessment of Electrochemical Properties
The increasing awareness of metal nitrides as feasible electrocatalysts in numerous energy-relevant transformations motivates greater synthetic control of mixed metal nitride compositions and their valence states. We report plasma-enhanced atomic layer deposition of titanium zirconium nitride, with an emphasis on the experimental parameters that influence the valence states of the two metal ions. We further examined the consequence of various ratios of Ti and Zr on the electrochemically active surface areas (ECSAs). Films with a Zr-to-Ti ratio of 1:2 or 1:1 yielded the highest ECSA; deviation from these values led to a large fraction of electrically insulating or semiconducting titanium or zirconium oxynitrides/oxides. The systematic examination establishes a fundamental structure-property relationship critical to assessing the electrochemical/electrocatalytic properties of metal nitrides.
36 MATERIALS SCIENCE↗
Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity
A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f(vy) or f(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
AWARE: The Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment
The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) performed comprehensive meteorological and aerosol measurements and ground-based atmospheric remote sensing at two Antarctic stations using the most advanced instrumentation available. A suite of cloud research radars, lidars, spectral and broadband radiometers, aerosol chemical and microphysical sampling equipment, and meteorological instrumentation was deployed at McMurdo Station on Ross Island from December 2015 through December 2016. A smaller suite of radiometers and meteorological equipment, including radiosondes optimized for surface energy budget measurement, was deployed on the West Antarctic Ice Sheet between 4 December 2015 and 17 January 2016. AWARE provided Antarctic atmospheric data comparable to several well-instrumented high Arctic sites that have operated for many years and that reveal numerous contrasts with the Arctic in aerosol and cloud microphysical properties. These include persistent differences in liquid cloud occurrence, cloud height, and cloud thickness. Antarctic aerosol properties are also quite different from the Arctic in both seasonal cycle and composition, due to the continent’s isolation from lower latitudes by Southern Ocean storm tracks. Antarctic aerosol number and mass concentrations are not only non-negligible but perhaps play a more important role than previously recognized because of the higher sensitivities of clouds at the very low concentrations caused by the large-scale dynamical isolation. Antarctic aerosol chemical composition, particularly organic components, has implications for local cloud microphysics. The AWARE dataset, fully available online in the ARM Program data archive, offers numerous case studies for unique and rigorous evaluation of mixed-phase cloud parameterization in climate models.
54 ENVIRONMENTAL SCIENCES↗
Factors that Impact the Redox Kinetics Associated with Sr1-xAxFeO3 (A = Ca, Ba) in Air Separations Technologies
Perovskite oxides, including SrFeO3-δ, have garnered significant interest for their impressive oxygen storage thermodynamics and kinetics. While traditional oxygen production is performed through cryogenic separation methods at the industrial scale, oxygen carrier-based separations have offered an attractive alternative at the modular scale. Perovskites are considered ideal oxygen carriers, with their structural flexibility minimizing both the thermodynamic barrier and the structural strain related to the oxygen uptake/release cycle. Our latest work has identified a few major factors that directly impact the oxygen release thermodynamics and kinetics for Sr1-xAxFeO3-δ (A = Ca, Ba) materials including elemental composition, crystallographic symmetry, surface area, and pretreatment conditions. This presentation will discuss these factors using Sr1¬-xCaxFeO3 and Ba1-xSrxFeO3 to illustrate their individual and collective impacts on the design of high-performance oxygen carriers when operating at 300-700 °C. Additionally, we will discuss how each of these factors can be leveraged for the numerous catalytic applications also dependent on the redox properties of mixed metal oxides.
Popczun, Eric↗
Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning
Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.
54 ENVIRONMENTAL SCIENCES↗
The Scaling and Units of the Elastic Response Term for Rayleigh Waves that is Output by Computer Programs in Seismology (CPS)
We report on the scaling and units of the Rayleigh wave elastic response function A R (ω) that is output from the widely used Computer Programs in Seismology (CPS) to aid in the modeling of ground motion sourced by atmospheric explosions. The program uses mixed units (km, second, km/s, gm/cc) to keep A R (ω) near 10 0 and prevent any numerical underflow or overflow. We compare two models for the response of an elastic half space to the output from CPS. Our application inputs the recommended, mixed unit geological models to determine how researchers must scale this output to obtain physical units for A R (ω) that represents the amplitude scaling for the minimum group velocity (Airy phase) contribution to Rayleigh waves. We determine that a CPS user must scale the output for A R (ω) by 10 -12 (m/km) 2 (g/cc/m 3 /kg) to obtain MKS (meter, kg, second) units and then must multiply this result by the vertical component eigenfunction squared, that is evaluated at the free surface.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Investigation of finite-volume methods to capture shocks and turbulence spectra in compressible flows
The aim of the present paper is to provide a comparison between several finite-volume methods of different numerical accuracy: the second-order Godunov method with PPM interpolation and the high-order finite-volume WENO method. In this work, the results show that while on a smooth problem the high-order method performs better than the second-order one, when the solution contains a shock all the methods collapse to first-order accuracy. In the context of the decay of compressible homogeneous isotropic turbulence with shocklets, the actual overall order of accuracy of the methods reduces to second-order, despite the use of fifth-order reconstruction schemes at cell interfaces. Most important, results in terms of turbulent spectra are similar regardless of the numerical methods employed, except that the PPM method fails to provide an accurate representation in the high-frequency range of the spectra. It is found that this specific issue comes from the slope-limiting procedure and a novel hybrid PPM/WENO method is developed that has the ability to capture the turbulent spectra with the accuracy of a high-order method, but at the cost of the second-order Godunov method. Overall, it is shown that virtually the same physical solution can be obtained much faster by refining a simulation with the second-order method and carefully chosen numerical procedures, rather than running a coarse high-order simulation. Our results demonstrate the importance of evaluating the accuracy of a numerical method in terms of its actual spectral dissipation and dispersion properties on mixed smooth/shock cases, rather than by the theoretical formal order of convergence rate.
97 MATHEMATICS AND COMPUTING↗
Sensitivity of Global Ocean Deoxygenation to Vertical and Isopycnal Mixing in an Ocean Biogeochemistry Model
Large-scale loss of oxygen under global warming is termed “ocean deoxygenation” and is caused by the imbalance between physical supply and biological consumption of oxygen in the ocean interior. Significant progress has been made in the theoretical understanding of ocean deoxygenation; however, many questions remain unresolved. The oxygen change in the tropical thermocline is poorly understood, with diverging projections among different models. Physical oxygen supply is controlled by a suite of processes that transport oxygen-rich surface waters into the interior ocean, which is expected to weaken due to increasing stratification under global warming. Using a numerical model and a series of sensitivity experiments, the role of ocean mixing is examined in terms of effects on the mean state and the response to a transient warming. Both vertical and horizontal (isopycnal) mixing coefficients are systematically varied over a wide range, and the resulting oxygen distributions in equilibrated and transient simulations are examined. The spatial patterns of oxygen loss are sensitive to both vertical and isopycnal mixing, and the sign of tropical oxygen trend under climate warming can reverse depending on the choice of mixing parameters. An elevated level of isopycnal mixing disrupts the vertical advective-diffusive balance of the tropical thermocline, increasing the mean state oxygen as well as the magnitude of the transient oxygen decline. These results provide first-order explanations for the diverging behaviors of simulated tropical oxygen with respect to ocean mixing parameters.
58 GEOSCIENCES↗
Disjunctive linear separation conditions and mixed-integer formulations for aircraft conflict resolution
In this paper, we address the aircraft conflict resolution problem in air traffic control. We introduce new mixed-integer programming formulations for aircraft conflict resolution with speed, heading and altitude control which are based on disjunctive linear separation conditions. We first examine the two-dimensional aircraft conflict resolution problem with speed and heading control represented as continuous decision variables. We show that the proposed disjunctive linear separation conditions are equivalent to the classical nonlinear conditions for aircraft separation. Further, we characterise conflict-free trajectories based on aircraft velocity bounds and propose a simple pre-processing algorithm to identify aircraft pairs which are either always conflict-free, or which cannot be separated using speed and heading control only. We then incorporate altitude control and propose a lexicographic optimisation formulation that aims to minimise the number of flight level changes before resolving outstanding conflicts via two-dimensional velocity control. The proposed mixed-integer programming formulations are nonconvex, and we propose convex relaxations, decomposition methods and constraint generation algorithms to solve the two-dimensional and lexicographic optimisation formulations to guaranteed optimality. Numerical experiments on four types of conflict resolution benchmarking instances are conducted to test the performance of the proposed mixed-integer formulations. Further, the proposed method is compared against two benchmarks based on state-of-the-art approaches for the aircraft conflict resolution problem. Our numerical results show that the proposed method largely outperforms both benchmarks in terms of runtime and is able to solve significantly more instances to global optimality.
97 MATHEMATICS AND COMPUTING↗
Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment
Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.
54 ENVIRONMENTAL SCIENCES↗
Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit
Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.
42 ENGINEERING↗
A unified hp-HDG framework for Friedrichs' PDE systems
Here, this work proposes a unified hp-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified hp-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified hp-formulations and demonstrate the effectiveness of hp-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, and singularities.
97 MATHEMATICS AND COMPUTING↗
A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations
Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.
54 ENVIRONMENTAL SCIENCES↗
Blending and cavern formation within non-Newtonian fluids in stirred tanks: Application to nuclear waste fluid processing
A numerical approach for predicting the time-accurate fluid flow and mixing properties of non-Newtonian fluids is presented. This approach, which is based on the lattice Boltzmann method, is used to characterize the blending in shear-thinning fluids, as well as cavern formation in yield-stress fluids. Predictions compare favorably to measured data and expectations from first-principles theory. Importantly, across the range of fluids and systems here, the simulations require no reparameterization or retuning between scenarios. This generality is exploited to model a nuclear waste processing unit operation.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Evaluation of several relevant fractionation processes as possible explanation for radioxenon isotopic activity ratios in samples taken near underground nuclear explosions in shafts and tunnels
Gas samples taken from two historic underground nuclear tests done in 1989 at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), were examined to determine how xenon isotopes fractionate because of early-time cavity processes, transport through the rock, or dispersal through tunnels. Xenon isotopes are currently being used to distinguish civilian sources of xenon in the atmosphere from sources associated with underground nuclear explosions (UNEs). The two nuclear tests included (1) BARNWELL, a test conducted in a vertical shaft approximately 600 m below ground surface at Pahute Mesa, and (2) DISKO ELM, a horizontal line-of-sight test done in P-tunnel approximately 261 m below the surface of Aqueduct Mesa. Numerical flow and transport models developed for the two sites had mixed success when attempting to match the observed xenon isotope ratios. At the BARNWELL site, the simulated xenon isotope ratios were consistent with measurements from the chimney and ground surface and appeared to have been affected primarily by fractionation during subsurface transport. At the DISKO ELM site, samples taken from two elevations in the chimney failed to show the degree of fractionation predicted by the models during transport and did not show evidence for significant fractionation due to early-time condensation of refractory xenon-precursor radionuclides into the melt glass. Gas samples taken from the adjacent tunnels in the days following the test showed mixed evidence for early-time separation of xenon isotopes from their iodine precursors.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗