Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constrained transport”

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 271 records · Page 15

Ion-Specific Effects on PuO 2 Nanoparticle Aggregation and Dissolution in Concentrated Electrolytes

Hydrolytic PuO 2 nanoparticles (NPs) are a dominant aqueous Pu-bearing phase in high ionic strength nuclear wastes, yet their reactivity in nonideal brines remains poorly constrained. We quantify how electrolyte identity and concentration control PuO 2 NP aggregation and ligand-assisted dissolution in acidic, high salinity solutions (NaCl, NaNO 3 , NaClO 4 , Na 2 SO 4 , Na 2 C 2 O 4 up to 5 M). A multitechnique workflow combining liquid scintillation counting (operationally defined aqueous [Pu]), scattering/electrokinetic measurements (aggregate size and zeta potential), and spectroscopy (UV–vis, XPS, Raman) resolves electrolyte-dependent partitioning between colloidal and molecular Pu species. Weakly coordinating anions (ClO 4 – , Cl – , NO 3 – ) largely preserve the (aggregated) nanoparticulate fraction but generate distinct dissolved species at high concentration, i.e., Pu(IV)–nitrato complexes in NaNO 3 and Pu(VI)–chloro complexes in NaCl. In contrast, stronger ligands substantially perturb PuO 2 NP stability: sulfate promotes partial dissolution to Pu(IV)–sulfate complexes at low concentration but reduces aqueous [Pu] at higher sulfate levels via secondary Pu(IV) sulfate formation, whereas oxalate drives strong dissolution to aqueous Pu–oxalate complexes. Aged NPs show similar trends with reduced aqueous fractions and more dominant aggregation mechanisms. In conclusion, these spectroscopically constrained speciation data provide a foundation for incorporating PuO 2 NP reactivity into thermodynamic and reactive transport models for high salinity waste and brine environments.

Aggregation↗

Elastic NaxMoS2-carbon-BASE triple interface direct robust solid-solid interface for all-solid-state Na-S batteries

The promises of all-solid-state (ASS) sodium batteries for the next generation energy storage are widely recognized but their developments have been severely constrained by the difficulties to design favorable solid-solid interfaces for unhindered Na-ion transport. Using the most promising ß?-Al2O3 solid state electrolyte (BASE) as a platform, we demonstrate here a triple nanojunction strategy that provides simultaneous strong Na adhesion and continuous Na-ions diffusion at solid-solid interface. Such triple junctions (NaxMoS2-carbon-BASE) were constructed by adhering ternary composite Na anodes containing dispersed 3 wt% MoS2 and 3 wt% carbon on BASE, and provide nearly complete adhesion of Na on BASE with a much smaller contact angle (~ 45o vs. 120o of pristine Na). The composite Na anodes exhibited ~ 3 times improved elastic property and the synergy of NaxMoS2 and carbon provides the required ionic and electronic diffusion channels at solid-solid interface, which significantly improve Na utilization and resist premature failure due to loss of solid-solid contact as Na shrink during high capacity stripping. As a result, Na metal at the triple junction exhibited more than five time reduced charge transfer resistance and at least 200 hours stable battery cycling at practical current densities. The novel anode architecture also enabled high capacity cycling of prototype ASS sodium sulfur batteries when coupled with advanced sulfur cathodes containing intrinsic Na-ions diffusion channels and redox catalytic mediators, leading to stable cycling with specific capacity of 1110 mAh g-1.

ß”-Al2O3 solid electrolyte, solid state batteries,↗

Impacts of Precipitation Events on Concentrations of Oxygenated Gas- and Particle-Phase Compounds Observed in the Amazon

Removal of gases and particles by precipitation (wet deposition) is a critical process that significantly influences the transport and chemical transformation of atmospheric compounds. However, there are few studies that directly measure or constrain the rates of this process under real-world conditions. This work quantifies the net change in ambient concentrations during precipitation events (removal rates) of gas- and particle-phase organic compounds at a surface site near Manaus, Brazil, during the GoAmazon2014/5 campaign. Removal rates of identified and unknown compounds that have been previously classified into source-based clusters are measured during rain events and categorized based on estimated properties of compounds and clusters. Highly oxygenated gases, such as isoprene oxidation products, are removed during precipitation events with a median removal rate of 0.09 h –1 and the fastest analyte is removed at a rate of 0.22 h –1 . Removal rates of particle-phase compounds are observed at roughly this median rate, while less soluble gases, such as terpenes, exhibit low removal rates. These results are roughly in agreement with prior theoretical estimates of wet deposition rates for comparable compounds, providing an empirical point of comparison while noting that our metric reflects the net influence of precipitation events rather than wet deposition alone.

atmospheric chemistry↗

Equation of State and Spin Crossover of (Al, Fe)–Phase H

The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases δ-AlOOH (delta phase), MgSiO 2 (OH) 2 (phase H), and ε-FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core-mantle boundary. Here we have constrained the equation of state, including the effects of a spin crossover in the Fe 3+ atoms, of (Al, Fe)-phase H: Al 0.84 Fe 3+ 0.07 Mg 0.02 Si 0.06 OOH, using powder X-ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)-phase H and δ-(Al, Fe)OOH. The changes in spin state of Fe 3+ in (Al, Fe)-phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared with δ-(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)-phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)-phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid-lower mantle (~1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)-phase H could reproduce some of the seismic signatures of large, low seismic-velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ-AlOOH)-(MgSiO 2 (OH) 2 )-(ε-FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.

58 GEOSCIENCES↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting

Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.

36 MATERIALS SCIENCE↗

Dynamic Analysis of a Six-Cable Parallel Robot for Automated Panelized Building Retrofits

Cable-Driven Parallel Robots (CDPRs) are highly suitable for automated panelized building retrofits, thanks to their compact footprint and high payload-to-weight ratio. A common CDPR configuration featuring eight cables, where the anchors form a rectangular prism in front of the building facade, offers a large wrench feasible workspace and good control versatility. However, installing upper anchor points requires additional support structures, such as towers or beams, increasing setup complexity and posing logistical challenges in construction settings. To mitigate these challenges, we propose a six-cable CDPR model specifically designed for automated panelized building retrofits. Although the feasible workspace is limited, our analysis shows that the proposed CDPR adequately covers the critical areas required for panel installation. To validate that the six-cable system can effectively transport the end effector to the desired installation pose, we calculated optimal trajectories based on a constrained dynamic model. The simulation results of the six-cable CDPR demonstrate promising potential for automated panelized building retrofits, effectively balancing simplicity, cost-effectiveness, and functionality.

Liu, Yifang [ORNL] (ORCID:0000000190817417)↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

NO(x) Concentrations in the Upper Troposphere as a Result of Lightning

Upper tropospheric NO(x) controls, in part, the distribution of ozone in this greenhouse-sensitive region of the atmosphere. Many factors control NO(x) in this region. As a result it is difficult to assess uncertainties in anthropogenic perturbations to NO from aircraft, for example, without understanding the role of the other major NO(x) sources in the upper troposphere. These include in situ sources (lightning, aircraft), convection from the surface (biomass burning, fossil fuels, soils), stratospheric intrusions, and photochemical recycling from HNO3. This work examines the separate contribution to upper tropospheric "primary" NO(x) from each source category and uses two different chemical transport models (CTMS) to represent a range of possible atmospheric transport. Because aircraft emissions are tied to particular pressure altitudes, it is important to understand whether those emissions are placed in the model stratosphere or troposphere and to assess whether the models can adequately differentiate stratospheric air from tropospheric air. We examine these issues by defining a point-by-point "tracer tropopause" in order to differentiate stratosphere from troposphere in terms of NO(x) perturbations. Both models predict similar zonal average peak enhancements of primary NO(x) due to aircraft (approx. = 10-20 parts per trillion by volume (pptv) in both January and July); however, the placement of this peak is primarily in a region of large stratospheric influence in one model and centered near the level evaluated as the tracer tropopause in the second. Below the tracer tropopause, both models show negligible NO(x) derived directly from the stratospheric source. Also, they predict a typically low background of 1 - 20 pptv NO(x) when tropospheric HNO3 is constrained to be 100 pptv of HNO3. The two models calculate large differences in the total background NO(x) (defined as the source of NO(x) from lightning + stratosphere + surface + HNO3) when using identical loss frequencies for NO(x). This difference is primarily due to differing treatments of vertical transport. An improved diagnosis of this transport that is relevant to NO(x) requires either measurements of a surface-based tracer with a substantially shorter lifetime than Rn-222 or diagnosis and mapping of tracer correlations with different source signatures. Because of differences in transport by the two models we cannot constrain the source of NO(x) from lightning through comparison of average model concentrations with observations of NO(x).

Penner, Joyce E.↗

A study of wing body blending for an advanced supersonic transport

Increases in supersonic cruise lift drag ratio were sought at Mach numbers 2.2 and 2.7 using wing body planform and thickness blending. Constrained twist and camber optimization was performed in the presence of nacelles. Wing and fuselage thickness distributions were optimized for either minimum volume wave drag or minimum total pressure wave drag. The zero leading edge suction lift drag ratios were determined for three wing planforms. The magnitude of the effect of leading edge suction on attainable lift drag ratio was defined on one planform and estimation of available leading edge suction was made.

Goebel, T. P.↗

Collaborative Research: Atlantic Ocean Tropical/Subtropical Processes from Seasonal to Decadal Time Scales: Model/Data, Model/Model Comparison and Model/Data Synthesis Through Assimilation

The effort of this first year of research has been focused on the assimilation of TOPEX/Poseidon altimetric data into a primitive equation model of the Atlantic tropical/subtropical circulation. A reduced-rank, stationary Kalman filter has been constructed to assimilate the altimetric sea surface height anomaly (SHA) into the model. The goal is to assess how the inter-hemispheric transports between the Atlantic subtropics and tropics are affected by the assimilation and how the subsurface thermocline structure , and its variability ,is dynamically constrained by the SHA. The model is a reduced-gravity primitive equation GCM of the upper Atlantic Ocean between 30 S and 30 N. The assimilation scheme is an approximation to the extended Kalman filter in which the error covariances of the state estimates are calculated only in a reduced- dimension subspace. The subspace is defined by the leading empirical orthogonal functions calculated from an unconstrained model calculation. Both an identical twin experiment using simulated SHA observations and assimilation of the real TOPEX data were performed. Results from the twin experiments demonstrate the ability of the method to constrain the ocean circulation and the subsurface temperature structure. The impact on the subsurface temperature structure of TOPEX assimilation was assessed using data from expandable bathythermographs. This showed a substantial improvement in the estimated temperature variability only within 13 degrees in latitude around the equator. The impact of TOPEX SHA assimilation on zonally integrated meridional transport across different latitudes was also estimated. Again within 13 degrees from the equator both the mean amplitude and interannual variability of the surface and subsurface transports were significantly enhanced, while the transports were insensitive to the assimilation in the subtropics.

Malannotte-Rizzoli, Paola↗

Spatiotemporal control of liquid crystal structure and dynamics through activity patterning

Active materials are capable of converting free energy into mechanical work to produce autonomous motion, and exhibit striking collective dynamics that biology relies on for essential functions. Controlling those dynamics and transport in synthetic systems has been particularly challenging. Here, we introduce the concept of spatially structured activity as a means of controlling and manipulating transport in active nematic liquid crystals consisting of actin filaments and light-sensitive myosin motors. In this work, simulations and experiments are used to demonstrate that topological defects can be generated at will and then constrained to move along specified trajectories by inducing local stresses in an otherwise passive material. These results provide a foundation for the design of autonomous and reconfigurable microfluidic systems where transport is controlled by modulating activity with light.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Burning conditions and transportation pathways determine biomass-burning aerosol properties in the Ascension Island marine boundary layer

Abstract. African biomass-burning aerosol (BBA) in the southeast Atlantic Ocean (SEA) marine boundary layer (MBL) is an important contributor to Earth's radiation budget, yet its representation remains poorly constrained in regional and global climate models. Data from the Layered Atlantic Smoke Interactions with Clouds (LASIC) field campaign on Ascension Island (7.95° S, 14.36° W) provide insight into how burning conditions, fuel type, transport pathways, and atmospheric processing affect the chemical, microphysical, and optical properties of BBA between June and September 2017. A total of 10 individual plume events characterize the seasonal evolution of the BBA properties. Early-season inefficient fires, determined by low refractory black carbon to above-background carbon monoxide mixing ratios (rBC : ΔCO), led to enhanced concentrations of organic- and sulfate-rich aerosols. Mid-season efficient fires, determined by higher rBC : ΔCO values, led to rBC-enriched BBA. A mix of efficient and inefficient fires later in the season resulted in conflicting BBA properties. Prolonged transport (∼ 10 d) through the MBL and lower free troposphere (FT) facilitated chemical and aqueous-phase processing, which led to a reduction in organic aerosol mass concentrations. This resulted in lower organic aerosol (OA) to rBC (OA : rBC) mass ratios (2–5) in the MBL compared to higher values (5–15) in the nearby FT. These atmospheric and cloud oxidation processes yield more light-absorbing BBA and explain the notably low single-scattering albedo at 530 nm (SSA530) values (< 0.80) observed in the MBL. This study establishes a robust correlation between SSA530 and OA : rBC across the MBL and FT, underscoring the dependency of optical properties on chemical composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of the Multi-Purpose Transportation System for the Space Launch System (SLS) Core Stage (CS) Flight Article

NASA's Multi-Purpose Transportation System (MPTS) is designed to transport the Space Launch System (SLS) vehicle segments by waterway and roadway. It is tasked with transporting the vehicle from where it is manufactured to its intermediate test location and final launch destination. Its design incorporates mechanisms that release degrees of freedom to prevent excessive loading during transit and ensure a successful delivery of the vehicle to its intended destination. In addition to the Core Stage (CS) flight article, the system will also move three Structural Test Articles (STAs), the Dynamic Demonstration Unit (DDU), and a simulated CS Pathfinder (weight, center of gravity, outer mold line dimensions, and overall length) over road terrain at four NASA centers and on the Pegasus barge. The MPTS independently supports the article at both ends while operating as a combined unit through automated monitoring of its released degree of freedom and corrective responses. This allows the system to constrain its payload in a statically determinate manner while traversing highly variable terrain. Multi-body simulation of the transportation route is useful to predict free-body motion within its range of mobility. The MPTS has completed its design and analysis developmental cycles. This paper describes the design challenges encountered in developing this system of large-scale structure, which incorporates complex mechanisms. The unique techniques and methodologies developed for analytical assessment of the hardware are also discussed.

Sandridge, Sarah↗

Black Carbon Concentration from Worldwide Aerosol Robotic Network (AERONET)

Worldwide black carbon concentration measurements are needed to assess the efficacy of the carbon emissions inventory and transport model output. This requires long-term measurements in many regions, as model success in one region or season does not apply to all regions and seasons. AERONET is an automated network of more than 180 surface radiometers located throughout the world. The sky radiance measurements obtained by AERONET are inverted to provide column-averaged aerosol refractive indices and size distributions for the AERONET database, which we use to derive column-averaged black carbon concentrations and specific absorptions that are constrained by the measured radiation field. This provides a link between AERONET sky radiance measurements and the elemental carbon concentration of transport models without the need for an optics module in the transport model. Knowledge of both the black carbon concentration and aerosol absorption optical depth (i.e., input and output of the optics module) will enable improvements to the transport model optics module.

Schuster, Greg↗

Burning conditions and transportation pathways determine biomass-burning aerosol properties in the Ascension Island marine boundary layer

African biomass-burning aerosol (BBA) in the southeast Atlantic Ocean (SEA) marine boundary layer (MBL) is an important contributor to Earth’s radiation budget yet its representation remains poorly constrained in regional and global climate models. Data from the Layered Atlantic Smoke Interactions with Clouds (LASIC) field campaign on Ascension Island (-7.95° N, -14.36° E) detail how fire source regions (burning conditions and fuel type), transport pathways, and longer-term chemical processing affect the chemical, microphysical, and optical properties of the BBA in the remote MBL between June and September of 2017. Ten individual plume events characterize the seasonal evolution of BBA characteristics. Inefficient burning conditions, determined by the mass ratio of refractory black carbon to above-background carbon monoxide (rBC:ΔCO), enhance organic- and sulfate-rich aerosol concentrations in June–July. In contrast, the heart of the burning season exhibited higher rBC:ΔCO values indicative of efficient burning conditions, correlating with more rBC-enriched BBA. Toward the end of the burning season, a mix of burning conditions results in increased variation of the BBA properties. The BBA transit to Ascension Island was predominantly through slow-moving pathways in the MBL and lower free troposphere (FT), facilitating prolonged chemical transformations through heterogeneous and aqueous phase processes. Heterogeneous oxidation can persist for up to 10 days, resulting in a considerable decrease in organic aerosol (OA) mass. OA to rBC mass ratios (OA:rBC) in the MBL between 2 and 5 contrast to higher values of 5 to 15 observed in the nearby FT. Conversely, early-season aqueous-phase processes primarily contributed to aerosol oxidation and some aerosol production, but not appreciable aerosol removal. These two chemical processes yield more light-absorbing BBA in the MBL than in the FT and explain the notably low scattering albedo at 530 nm (SSA 530 ) values (< 0.80) at Ascension Island. This study establishes a robust correlation between SSA 530 and OA:rBC across both MBL and FT, underscoring the dependency of optical properties on chemical composition. These findings highlight how the interplay between chemical composition and atmospheric processing can be improved in global and regional climate models. Questions remain on the mixing of aerosols with different pathway histories, and on what accounts for the doubling of the mass absorption coefficient in the boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Interannual and Seasonal Variability of Biomass Burning Emissions Constrained by Satellite Observations

We present a methodology for estimating the seasonal and interannual variation of biomass burning designed for use in global chemical transport models. The average seasonal variation is estimated from 4 years of fire-count data from the Along Track Scanning Radiometer (ATSR) and 1-2 years of similar data from the Advanced Very High Resolution Radiometer (AVHRR) World Fire Atlases. We use the Total Ozone Mapping Spectrometer (TOMS) Aerosol Index (AI) data product as a surrogate to estimate interannual variability in biomass burning for six regions: Southeast Asia, Indonesia and Malaysia, Brazil, Central America and Mexico, Canada and Alaska, and Asiatic Russia. The AI data set is available from 1979 to the present with an interruption in satellite observations from mid-1993 to mid-1996; this data gap is filled where possible with estimates of area burned from the literature for different regions. Between August 1996 and July 2000, the ATSR fire-counts are used to provide specific locations of emissions and a record of interannual variability throughout the world. We use our methodology to estimate mean seasonal and interannual variations for emissions of carbon monoxide from biomass burning, and we find that no trend is apparent in these emissions over the last two decades, but that there is significant interannual variability.

Duncan, Bryan N.↗

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗