Engineering Papers⌕ Search

Engineering topics

Antoun, Tarabay

Publications and source records attributed to Antoun, Tarabay.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Permeability scaling relationships of volcanic tuff from core to field scale measurements

A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.

Environmental sciences↗
DOE OSTI2022

Evaluation of subsurface transport processes of delayed gas signatures applicable to underground nuclear explosions

Radioactive gas signatures from underground nuclear explosions (UNEs) result from gas-migration processes occurring in the subsurface. The processes considered in this study either drive or retard upward migration of gases from the detonation cavity. The relative importance of these processes is evaluated by simulating subsurface transport in a dual-permeability medium for the multi-tracer Noble Gas Migration Experiment (NGME) originally intended to study some aspects of transport from a UNE. For this experiment, relevant driving processes include weak two-phase convection driven by the geothermal gradient, over pressuring of the detonation cavity, and barometric pumping while gas sorption, dissolution, radioactive decay, and usually diffusion represent retarding processes. From deterministic simulations we found that over-pressuring of the post-detonation chimney coupled with barometric pumping produced a synergistic effect amplifying the tracer-gas reaching the surface. Bounding simulations indicated that the sorption and dissolution of gases, tending to retard transport, were much smaller than anticipated by earlier laboratory studies. The NGME observations themselves show that differences in gas diffusivity have a larger effect on influencing upward transport than do the combined effects of tracer-gas sorption and dissolution, which is consistent with a Sobol’ sensitivity analysis. Both deterministic simulations and those considering parametric uncertainties of transport-related properties predict that the excess in concentration of SF 6 compared to 127 Xe as might be captured in small volumetric samples should be much smaller than the order-of-magnitude contrast found in the large-volume gas samples taken at the site. While extraction of large-volume subsurface gas samples is shown to be capable of distorting in situ gas compositions, the highly variable injection rate of SF 6 into the detonation cavity relative to that of 127 Xe at the start of the field experiment is the most likely explanation for the large difference in observed concentrations.

58 GEOSCIENCES↗
DOE OSTI2022

Implications of Underground Nuclear Explosion Cavity Evolution for Radioxenon Isotopic Composition

Abstract Isotopic ratios of radioxenons sampled in the atmosphere or subsurface can be used to verify the occurrence of an underground nuclear explosion (UNE). Differences in the half-lives of radioactive xenon precursors and their decay-chain networks produce different time-dependent concentration profiles of xenon isotopes allowing isotopic ratios to be used for tracking UNE histories including estimating the time of detonation. In this study, we explore the potential effects of post-detonation cavity processes: precipitation of iodine precursors, gas seepage, and prompt venting on radioxenon isotopic evolution which influences UNE histories. Simplified analytical models and closed-form solutions yielding a potentially idealized radioactive decay/ingrowth chain in a closed and well-mixed system typically have limited application by not including the partitioning of the radionuclide inventory between a gas phase and rock melt created by the detonation and by ignoring gas transport from the cavity to host rock or ground surface. In reality, either subsurface transport or prompt release that is principally responsible for gas signatures violates the closed-system (or batch-mode) assumption. A closed-form solution representing time-dependent source-term activities is extended by considering the cavity partitioning process, slow seepage, and/or prompt release of gases from the cavity and applied to realistic systems.

Sun, Yunwei (ORCID:0000000194801439)↗
DOE OSTI2015

Analysis of Ground Motion from An Underground Chemical Explosion

Here in this paper we investigate the excitation and propagation of far-field seismic waves from the 905 kg trinitrotoluene equivalent underground chemical explosion SPE-3 recorded during the Source Physics Experiment (SPE) at the Nevada National Security Site. The recorded far-field ground motion at short and long distances is characterized by substantial shear-wave energy, and large azimuthal variations in P-and S-wave amplitudes. The shear waves observed on the transverse component of sensors at epicentral distances <50 m suggests they were generated at or very near the source. The relative amplitude of the shear waves grows as the waves propagate away from the source. We analyze and model the shear-wave excitation during the explosion in the 0.01–10 Hz frequency range, at epicentral distances of up to 1 km. We used two simulation techniques. One is based on the empirical isotropic Mueller–Murphy (MM) (Mueller and Murphy, 1971) nuclear explosion source model, and 3D anelastic wave propagation modeling. The second uses a physics-based approach that couples hydrodynamic modeling of the chemical explosion source with anelastic wave propagation modeling. Comparisons with recorded data show the MM source model overestimates the SPE-3 far-field ground motion by an average factor of 4. The observations show that shear waves with substantial high-frequency energy were generated at the source. However, to match the observations additional shear waves from scattering, including surface topography, and heterogeneous shallow structure contributed to the amplification of far-field shear motion. Comparisons between empirically based isotropic and physics-based anisotropic source models suggest that both wave-scattering effects and near-field nonlinear effects are needed to explain the amplitude and irregular radiation pattern of shear motion observed during the SPE-3 explosion.

58 GEOSCIENCES↗
DOE OSTI2015

On the generation of horizontal shear waves by underground explosions in jointed rocks

This paper describes computational studies of non-spherical ground motions generated by spherical explosions in a heavily jointed granite formation. Various factors affecting the shear wave generation are considered, including joint spacing, orientation, persistence and properties. Simulations are performed both in 2D for a single joint set to elucidate the basic response mechanisms, and in 3D for multiple joint sets to realistically represent in situ conditions in a realistic geologic setting. The joints are modeled explicitly using both contact elements and weakness planes in the material. Simulations are performed both deterministically and stochastically to quantify the effects of geologic uncertainties on near field ground motions. The mechanical properties of the rock and the joints as well as the joint spacing and orientation are taken from experimental test data and geophysical logs corresponding to the Climax Stock granitic outcrop, which is the geologic setting of the Source Physics Experiment (SPE). Agreement between simulation results and near field wave motion data from SPE enables newfound understanding of the origin and extent of non-spherical motions associated with underground explosions in fractured geologic media.

58 GEOSCIENCES↗
DOE OSTI2013

Evolution of an interfacial crack on the concrete-embankment boundary

Failure of a dam can have subtle beginnings. A small crack or dislocation at the interface of the concrete dam and the surrounding embankment soil initiated by, for example, a seismic or an explosive event can lead to a catastrophic failure of the dam. The dam may ‘self-rehabilitate’ if a properly designed granular filter is engineered around the embankment. Currently, the design criteria for such filters have only been based on experimental studies. We demonstrate the numerical prediction of filter effectiveness at the soil grain scale. This joint LLNL-ERDC basic research project, funded by the Department of Homeland Security’s Science and Technology Directorate (DHS S&T), consists of validating advanced high performance computer simulations of soil erosion and transport of grain- and dam-scale models to detailed centrifuge and soil erosion tests. Validated computer predictions highlight that a resilient filter is consistent with the current design specifications for dam filters. These predictive simulations, unlike the design specifications, can be used to assess filter success or failure under different soil or loading conditions and can lead to meaningful estimates of the timing and nature of full-scale dam failure.

58 GEOSCIENCES↗
DOE OSTI2013

GEOS Code Development Road Map

GEOS is a massively parallel computational framework designed to enable HPC-based simulations of subsurface reservoir stimulation activities with the goal of optimizing current operations and evaluating innovative stimulation methods. GEOS will enable coupling of different solvers associated with the various physical processes occurring during reservoir stimulation in unique and sophisticated ways, adapted to various geologic settings, materials and stimulation methods. The overall architecture of the framework includes consistent data structures and will allow incorporation of additional physical and materials models as demanded by future applications. Along with predicting the initiation, propagation and reactivation of fractures, GEOS will also generate a seismic source term that can be linked with seismic wave propagation codes to generate synthetic microseismicity at surface and downhole arrays. Similarly, the output from GEOS can be linked with existing fluid/thermal transport codes. GEOS can also be linked with existing, non-intrusive uncertainty quantification schemes to constrain uncertainty in its predictions and sensitivity to the various parameters describing the reservoir and stimulation operations. We anticipate that an implicit-explicit 3D version of GEOS, including a preliminary seismic source model, will be available for parametric testing and validation against experimental and field data by Oct. 1, 2013.

02 PETROLEUM↗