Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial heterogeneity”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

AI-NERD: Elucidation of relaxation dynamics beyond equilibrium through AI-informed X-ray photon correlation spectroscopy

Abstract Understanding and interpreting dynamics of functional materials in situ is a grand challenge in physics and materials science due to the difficulty of experimentally probing materials at varied length and time scales. X-ray photon correlation spectroscopy (XPCS) is uniquely well-suited for characterizing materials dynamics over wide-ranging time scales. However, spatial and temporal heterogeneity in material behavior can make interpretation of experimental XPCS data difficult. In this work, we have developed an unsupervised deep learning (DL) framework for automated classification of relaxation dynamics from experimental data without requiring any prior physical knowledge of the system. We demonstrate how this method can be used to accelerate exploration of large datasets to identify samples of interest, and we apply this approach to directly correlate microscopic dynamics with macroscopic properties of a model system. Importantly, this DL framework is material and process agnostic, marking a concrete step towards autonomous materials discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DONUT: physics-aware machine learning for real-time X-ray nanodiffraction analysis

Coherent X-ray scattering techniques are critical for investigating the fundamental structural properties of materials at the nanoscale. While advancements have made these experiments more accessible, real-time analysis remains a significant bottleneck, often hindered by artifacts and computational demands. In scanning X-ray nanodiffraction microscopy, which is widely used to spatially resolve structural heterogeneities, this challenge is compounded by the convolution of the divergent beam with the sample’s local structure. To address this, we introduce DONUT (Diffraction with Optics for Nanobeam by Unsupervised Training), a physics-aware neural network designed for the rapid and automated analysis of nanobeam diffraction data. By incorporating a differentiable geometric diffraction model directly into its architecture, DONUT learns to predict crystal lattice strain and orientation in real-time. Crucially, this is achieved without reliance on labeled datasets or pre-training, overcoming a fundamental limitation for supervised machine learning in X-ray science. We demonstrate experimentally that DONUT accurately extracts all features within the data over 200 times more efficiently than conventional fitting methods.

Materials science↗
DOE OSTI2025

The variety and origin of materials accreted by Bennu’s parent asteroid

The first bodies to form in the Solar System acquired their materials from stars, the presolar molecular cloud and the protoplanetary disk. Asteroids that have not undergone planetary differentiation retain evidence of these primary accreted materials. However, geologic processes such as hydrothermal alteration can dramatically change their bulk mineralogy, isotopic compositions and chemistry. Here we analyse the elemental and isotopic compositions of samples from asteroid Bennu to uncover the sources and types of material accreted by its parent body. We show that some primary accreted materials escaped the extensive aqueous alteration that occurred on the parent asteroid, including presolar grains from ancient stars, organic matter from the outer Solar System or molecular cloud, refractory solids that formed close to the Sun, and dust enriched in neutron-rich Ti isotopes. We find Bennu to be richer in isotopically anomalous organic matter, anhydrous silicates, and light isotopes of K and Zn than its closest compositional counterparts, asteroid Ryugu and Ivuna-type (CI) carbonaceous chondrite meteorites. We propose that the parent bodies of Bennu, Ryugu and CI chondrites formed from a common but spatially and/or temporally heterogeneous reservoir of materials in the outer protoplanetary disk.

Kuiper belt↗
DOE OSTI2023

Viral but not bacterial community successional patterns reflect extreme turnover shortly after rewetting dry soils

As central members of soil trophic networks, viruses have the potential to drive substantial microbial mortality and nutrient turnover. Pinpointing viral contributions to terrestrial ecosystem processes remains a challenge, as temporal dynamics are difficult to unravel in the spatially and physicochemically heterogeneous soil environment. In Mediterranean grasslands, the first rainfall after seasonal drought provides an ecosystem reset, triggering microbial activity during a tractable window for capturing short-term dynamics. In this report we simulated precipitation in microcosms from four distinct dry grassland soils and generated 144 viromes, 84 metagenomes and 84 16S ribosomal RNA gene amplicon datasets to characterize viral, prokaryotic and relic DNA dynamics over 10 days. Vastly different viral communities in each soil followed remarkably similar successional trajectories. Wet-up triggered a significant increase in viral richness, followed by extensive compositional turnover. Temporal succession in prokaryotic communities was much less pronounced, perhaps suggesting differences in the scales of activity captured by viromes (representing recently produced, ephemeral viral particles) and total DNA. Still, differences in the relative abundances of Actinobacteria (enriched in dry soils) and Proteobacteria (enriched in wetted soils) matched those of their predicted phages, indicating viral predation of dominant bacterial taxa. Rewetting also rapidly depleted relic DNA, which subsequently reaccumulated, indicating substantial new microbial mortality in the days after wet-up, particularly of the taxa putatively under phage predation. Production of abundant, diverse viral particles via microbial host cell lysis appears to be a conserved feature of the early response to soil rewetting, and results suggest the potential for ‘Cull-the-Winner’ dynamics, whereby viruses infect and cull but do not decimate dominant host populations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Controlled electrochemical nutrient delivery to enhance marine primary productivity

Oceanic photosynthesis contributes to approximately half of the Earth’s net annual primary productivity. Marine photosynthetic productivity has a high degree of heterogeneity due to spatial and temporal co-limitations of light, temperature, and/or nutrients. Across coastal, near-shore, and open ocean regions, insufficient concentrations of key nutrients (e.g., N, P, Fe) can limit primary productivity. Although studies have shown a significant increase in primary productivity with the addition of low doses of trace nutrients, a sustainable approach to reliably deliver and maintain low doses of nutrients and ensure their bioavailability, remains a challenge. Chemical nutrient addition has relied on the use of chelating agents to ensure nutrient bioavailability, but synthetic chelators are persistent environmental pollutants. In this study, we demonstrate for the first time the use of a controlled electrochemical nutrient delivery (CEND) approach to accelerate the growth of phytoplankton without the need for such chelators. Our study uses commercial stainless-steel electrodes to deliver low concentrations of iron to enhance growth rates in the microalga Picochlorum celeri TG2. To demonstrate the process control offered by the CEND method, we evaluate iron delivery as a function of pulse time, pulse frequency, and rest time between pulses. Our data show that at the same total Fe dose of 163 ppb, electrochemical iron delivery can achieve 9.54 ± 0.58 mg biomass/µg Fe, which is comparable to 9.14 ± 0.17 mg biomass/µg Fe achieved with chemical iron additions that include the synthetic chelating agent ethylenediaminetetraacetic acid (EDTA). Further, when different total iron doses (163 ppb, 325 ppb, and 650 ppb) were delivered over 72 h using CEND, biomass yield per iron dose was higher at lower doses: 9.54 ± 0.58 mg biomass/µg Fe at 163 ppb vs. 4.32 ± 0.32 mg biomass/µg Fe at 650 ppb. This highlights the benefits of CEND in delivering frequent and low doses of nutrients for improved process efficiency. Preliminary assessments show both lower cost and reduced greenhouse gas emissions from electrochemical over chemical iron additions with EDTA. The CEND approach opens new pathways to enhance marine primary productivity, without the unintended environmental impacts of synthetic chelators.

analytical chemistry↗
DOE OSTI2025

DONUT: Physics-aware Machine Learning for Real-time X-ray Nanodiffraction Analysis

SF-25-088 Coherent X-ray scattering techniques are critical for investigating the fundamental structural properties of materials at the nanoscale. While advancements have made these experiments more accessible, real-time analysis remains a significant bottleneck, often hindered by artifacts and computational demands. In scanning X-ray nanodiffraction microscopy, which is widely used to spatially resolve structural heterogeneities, this challenge is compounded by the convolution of the divergent beam with the sample’s local structure. To address this, we introduce DONUT (Diffraction with Optics for Nanobeam by Unsupervised Training), a physics-aware neural network designed for the rapid and automated analysis of nanobeam diffraction data. By incorporating a differentiable geometric diffraction model directly into its architecture, DONUT learns to predict crystal lattice strain and orientation in real-time. Crucially, this is achieved without reliance on labeled datasets or pre-training, overcoming a fundamental limitation for supervised machine learning in X-ray science. We demonstrate experimentally that DONUT accurately extracts all features within the data over 200 times more efficiently than conventional fitting methods.

Zhou, Tao [Argonne National Laboratory (ANL), Argo↗
DOE OSTI2020

Observational Study of the Thermodynamics and Morphological Characteristics of a Midlatitude Continental Cold Pool Event

Abstract Convectively generated cold pools are important to the Earth system as they exert strong controls on deep convective-storm initiation, intensity, and life cycle. Despite their importance, efforts to introduce such cold pool controls into weather and climate models lack guidance and/or physical constraints from cold pool observations. This work presents a detailed, purely observational analysis of a cold pool event that took place on 23–24 May 2011 in north-central Oklahoma. The characteristics of the cold pool, and the spatiotemporal evolution of the hydrometeors and dynamics in the proximity of the cold pool, are studied with high-resolution observations. The unprecedented dataset used in this work to study cold pool characteristics includes an enhanced network of surface weather stations, a high-temporal-frequency sounding array, and the NEXRAD and Atmospheric Radiation Measurement (ARM) Southern Great Plains radar networks. The potential use of NEXRAD surveillance scans to estimate height and propagation speed of the leading edge of the cold pool (LECP) is presented in this work. Manual identification and tracking of the LECP from NEXRAD imagery shows a spatial and temporal heterogeneity of the LECP properties. Surprisingly, over its detected life cycle, the LECP speed remains almost constant, even though the strength of the cold pool diminishes in time and its height varies. Radar analysis shows that pulses of graupel and hail within downdrafts in the convective system generating the cold pool appeared to be related to temporary increases in the LECP height.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Report on Properties and Microstructure of 3D Printed Inc-718

The report presents the microstructure and mechanical properties of 3D printed Inconel 718 to assess its potential use as a structural material for the Transformation Challenge Reactor (TCR). The structural components near the outlet of the core will experience significant neutron fluxes and outlet coolant temperatures from the hot standby temperature of 300°C to nearly 550°C at the center of the part. These components must support the core in appropriate loading conditions and require structural analysis at relevant temperatures. Strong spatial and chemical heterogeneity was found in as-built (ASB) Inconel 718. Three heat treatments were designed and conducted to simplify the microstructure and determine how each precipitating phase contributed to the overall strength. Baseline mechanical properties were measured from uniaxial tensile tests on subsize SS-J2 specimens at room temperature and at elevated temperatures of 300, 450, and 600°C. Microstructure electron microscopy was performed on ASB Inconel 718 and heat treated to correlate the observed mechanical properties with nanoscale features. Homogenization of the microstructure led to a highly ductile Inconel with lower strength compared with wrought Inconel 718. The tensile properties of additively manufactured 718 using a standard ASTMrecommended heat treatment were consistent with literature and with the ASTM for the properties of this alloy. A higher fraction of the δ phase led to shorter uniform elongation without altering other engineering properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Tradeoffs in forest resilience to satellite-based estimates of water and productivity losses

Climate change is expected to intensify the global hydrological cycle via increased evapotranspiration (ET). Changes in ET have important consequences for water-related ecosystem services as well as for forest productivity. Forests play a critical role in ET: active trees sequester carbon at the expense of losing water, yet if they are not active, water remains in the trees, but productivity stalls. Here, we estimated the capacity of forests to recover (resilience) from extreme water or productivity losses based on modeled ET through remote sensing methods across the conterminous USA. Results show there is a spatial tradeoff in resilience to water or productivity losses. Forests were more resilient to water loss in the drier west, and more resilient to productivity loss in the wetter east. Still, resilience was highly heterogeneous across the continent. Projections under climate change indicate both declines and increases of up to 13% in forest resilience, with broad regional trends but high spatial heterogeneity within regions. This study brings crucial satellite-based insights on how water dynamics, and consequentially productivity and water provision in forests, might be impacted by climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial mapping of plant N-glycosylation cellular heterogeneity inside soybean root nodules provided insights into legume-rhizobia symbiosis

Although ubiquitously present, information on the function of complex N-glycan posttranslational modification in plants is very limited and it is often neglected. In this work, we adopted an enzyme-assisted matrix-assisted laser desorption/ionization mass spectrometry imaging strategy to visualize the distribution and identity of N-glycans in soybean root nodules at a cellular resolution. We additionally performed proteomics analysis to probe the potential correlation to proteome changes during symbiotic rhizobia-legume interactions. Our ion images reveal that intense N-glycosylation occurs in the sclerenchyma layer, and inside the infected cells within the infection zone, while morphological structures such as the cortex, uninfected cells, and cells that form the attachment with the root are less N-glycosylated. Notably, we observed different N-glycan profiles between soybean root nodules infected with wild type rhizobia and those infected with mutant rhizobia incapable of efficiently fixing atmospheric nitrogen. The majority of complex N-glycan structures, particularly those with characteristic Lewis-a epitopes, are overexpressed in the mutant nodules. Our proteomic results revealed that these glycans likely originated from proteins that maintain the redox balance crucial for proper nitrogen fixation, but also from enzymes involved in N-glycan and phenylpropanoid biosynthesis. These findings indicate possible involvement of Lewis-a glycans in these critical pathways during legume-rhizobia symbiosis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for inverse problems: Determining material properties for heterogeneous systems

Here, we determine spatially varying discontinuous material properties using a domain-decomposition based physics-informed neural networks (PINNs) framework named the Adaptive Interface-PINNs or AdaI-PINNs (Roy et al., 2024). We propose the use of distinct neural networks for the field variables and material properties within each material, utilizing adaptive activation functions. While the neural networks across different materials share the same weights and biases, their activation functions are uniquely tailored using a hyperparameter that influences the slope of the activation function. The proposed framework is tested on several one-dimensional and two-dimensional benchmark examples, and its performance is compared with conventional PINNs and existing domain-decomposition PINNs frameworks, namely, the Multi-domain physics-informed neural network (M-PINN), and the eXtended physics-informed neural networks (XPINNs). The results demonstrate that the proposed approach can determine randomly distributed discontinuous material properties with an L 2 error of $\mathscr{O}$ (10 -3 ) for the material property and the root-mean-square error of $\mathscr{O}$ (10 -3 ) for the primary variable while the other approaches yield errors that are approximately two orders of magnitude larger (that is, $\mathscr{O}$ (10 -1 )). Moreover, the spatial distribution of material properties obtained using the proposed framework is in close agreement with the true distribution, whereas the other approaches fare much worse. Additionally, the proposed approach is approximately 40% faster than its competitors, indicating its potential as a robust alternative for solving inverse problems in heterogeneous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reduced ecosystem resilience quantifies fine-scale heterogeneity in tropical forest mortality responses to drought

Sensitivity of forest mortality to drought in carbon-dense tropical forests remains fraught with uncertainty, while extreme droughts are predicted to be more frequent and intense. Here, in this study, the potential of temporal autocorrelation of high-frequency variability in Landsat Enhanced Vegetation Index (EVI), an indicator of ecosystem resilience, to predict spatial and temporal variations of forest biomass mortality is evaluated against in situ census observations for 64 site-year combinations in Costa Rican tropical dry forests during the 2015 ENSO drought. Temporal autocorrelation, within the optimal moving window of 24 months, demonstrated robust predictive power for in situ mortality (leave-one- out cross-validation R 2 = 0.54), which allows for estimates of annual biomass mortality patterns at 30 m resolution. Subsequent spatial analysis showed substantial fine-scale heterogeneity of forest mortality patterns, largely driven by drought intensity and ecosystem properties related to plant water use such as forest deciduousness and topography. Highly deciduous forest patches demonstrated much lower mortality sensitivity to drought stress than less deciduous forest patches after elevation was controlled. Our results highlight the potential of high-resolution remote sensing to “fingerprint” forest mortality and the significant role of ecosystem heterogeneity in forest biomass resistance to drought.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Spatially Resolved Top-Down Proteomics of Tissue Sections Based on a Microfluidic Nanodroplet Sample Preparation Platform

Conventional proteomics measures the averaged signal from mixed cell populations or bulk tissues, leading to the dilution of significant changes in subpopulations of cells that might serve as important biomarkers. Recent developments in bottom-up proteomics have enabled spatial mapping of cellular heterogeneity in tissue microenvironments. However, bottom-up proteomics cannot precisely infer the abundance changes of intact proteins, which are presented as proteoforms. Herein, we described a spatially resolved top-down proteomics (TDP) platform for proteoform identification and quantification directly on thin tissue sections. The spatial TDP platform consisted of a nanoPOTS (nanodroplet Processing in One pot for Trace Samples)-based sample preparation system and an LCM (laser capture microdissection)-based cell isolation system. We improved the nanoPOTS sample preparation by adding benzonase in the extraction buffer to enhance the coverage of nucleus proteins. Using ~200 cultured cells as model samples, the improved approach increased proteoform identifications from 493 to 700; newly identified proteoforms primarily corresponded to nuclear proteins. To demonstrate the spatial TDP platform in tissue samples, we analyzed LCM-isolated tissue voxels from rat brain cortex and hypothalamus regions. We quantified 426 proteoforms by combining identifications from TopPIC and TDPortal with the quantitation from ProMex. Several proteoforms corresponding to the same gene exhibited mixed abundance profiles between two tissue regions, suggesting potential PTM-specific spatial distributions. The spatial TDP workflow has prospects for biomarker discovery at proteoform level from small tissue sections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Tailoring the Reaction Heterogeneity for Robust Li‐Rich Cathodes

The practical application of Li-rich Mn-based layered oxides (LLO) cathode is hindered by severe capacity and voltage degradation resulting from severe oxygen release and irreversible phase transition. In this work, the reaction heterogeneity, which describes the spatially resolved electrochemical divergence within individual cathode particles, is engineered through compositional gradient design to couple Li + transport kinetics and the anion redox activity between particle interiors and surfaces. It is revealed that Co/Mn concentration gradient within particles creates heterogeneous phase content distribution and structural ordering, inducing surface-bulk reaction heterogeneity that significantly impacts the overall electrochemical performance. Specifically, Li 2 MnO 3 -poor and Co-enriched surface effectively mitigates the oxygen loss and enhances electrochemical reaction kinetics, benefited from the reduced surface redox reactivity and induced highly ordered intra-layered cationic arrangement. Meanwhile, the Li 2 MnO 3 -enriched core with slight Li/Ni intermixing provides high reversible capacity and strong mechanical stability. Consequently, the greatly enhanced anion redox reversibility, Li + diffusion dynamics, and structure stability endow LLO with exceptional electrochemical properties, showing a capacity retention of 86.0% and a reduced voltage decay of 0.518 mV per cycle after 500 cycles at 1 C. This work provides a valuable strategy to tailor the redox chemistry and achieve robust LLO.

25 ENERGY STORAGE↗
DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
DOE OSTI2023

Dynamic In‐Plane Heterogeneous and Inverted Response of Graphite to Fast Charging and Discharging Conditions in Lithium‐Ion Pouch Cells

Solutions for improving fast charging of lithium‐ion batteries have largely focused on alleviating through‐plane lithiation gradients while little is understood about in‐plane heterogeneities and how to resolve them. Herein, high‐speed synchrotron X‐ray diffraction (XRD) resolves graphite lithiation spatially and temporally during 6 C charging and 2 C discharging. At every point during operation, considerable differences in the state of lithiation across the pouch cell are present. Some regions are more responsive to operation than others, reaching full lithiation early during charge and full delithiation during discharge. Other regions within the cell never fully delithiate during discharge, despite a prolonged voltage hold at 2.8 V. Using time‐resolved XRD data, the calculated local current density (mA cm −2 ) at the graphite surface shows an unexpected occurrence of local inverted current densities where regions of graphite are observed to delithiate during charging and lithiate during discharging. A pseudo‐3D model is developed for the graphite electrode with spatially varying microstructural tortuosity to show how microstructural heterogeneity could influence spatial charge dynamics. The model could not predict the complex in‐plane charge behavior observed within the cell. Consequently physics‐based charging protocols based on homogeneous electrode assumptions may underestimate the local variations in charge dynamics and occurrence of lithium plating.

25 ENERGY STORAGE↗
DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗