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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Observing Arctic Sea Ice

Our understanding of Arctic sea ice and its wide-ranging influence is deeply rooted in observation. Advancing technologies have profoundly improved our ability to observe Arctic sea ice, document its processes and properties, and describe atmosphere-ice-ocean interactions with unprecedented detail. Yet, our progress toward better understanding the Arctic sea ice system is mired by the stark disparities between observations that tend to be siloed by method, scientific discipline, and application. This article presents a review and philosophical design for observing sea ice and accelerating our understanding of the Arctic sea ice system. We give a brief history of Arctic sea ice observations and showcase the 2018 melt season within the context of five observational themes: spatial heterogeneity, temporal variability, cross-disciplinary science, scalability, and retrieval uncertainty. We synthesize buoy data, optical imagery, satellite retrievals, and airborne measurements to demonstrate how disparate data sets can be woven together to transcend observational-scale issues. The results show that there are limitations to interpreting any single data set alone. However, many of these limitations can be surmounted by combining observations that cross spatial and temporal scales. We conclude the article with pathways toward coordinating across observational platforms in order to: (1) optimize the scientific, operational, and community return on observational investments, and (2) facilitate a richer understanding of Arctic sea ice and its role in the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Mesoscale monitoring of the soil freeze/thaw boundary from orbital microwave radiometry

The fundamental objectives are to test the feasibility of delineating the lateral boundary between frozen and thawed condition in the surface layer of soil from orbital microwave radiometry and secondly to examine the sensitivity of general circulation models to an explicit parameterization of the boundary condition. Physical models were developed to relate emissivity to scene properties and a simulation package was developed to predict brightness temperature as a function of emissivity and physical temperature in order to address issues of heterogeneity, scaling, and scene dynamics. Radiative transfer models were develped for both bare soil surfaces and those obscured by an intervening layer of vegetation or snow. These models relate the emissivity to the physical properties of the soil and to those of the snow or vegetation cover. A SMMR simulation package was developed to evaluate the adequacy of the emission models and the limiting effects of scaling for realistic scenarios incorporating spatially heterogeneous scenes with dynamic moisture and temperature gradients at the pixel scale.

Ulaby, Fawwaz T.↗
DOE OSTI2024

An uncertain future change in aridity over the tropics

Abstract An ensemble of climate models from phase six of the Coupled Model Intercomparison Project shows that temperature and potential evapotranspiration (PET) are projected to increase globally towards the end of the 21st century. However, climate models show a spatially heterogeneous change in precipitation over the tropics. Consequently, future changes in aridity (a measure of water availability) are complex and location-dependent. We assess future changes in aridity using three climate models and several single-forcing experiments. Near-term (2021–2040) changes in aridity are small, and we focus instead on its long-term (2081–2100) changes. We show that the increase in greenhouse gases (GHG) primarily explains the spatial pattern, magnitude and ensemble spread of the long-term future changes in aridity. On this timescale, the effects of changes in emissions of anthropogenic aerosols are moderate compared to the effects of increases in atmospheric GHG concentrations. Model diversity in the responses to GHG concentration is large over northern Africa and North and South America. We suggest the large uncertainty is due to differences between models in simulating the effects of an increase in GHG concentrations on surface air temperature over the North Atlantic Ocean, on the interhemispheric temperature gradient, and on PET over North and South America.

Environmental Sciences & Ecology↗
DOE OSTI2020

Deep Learning Prediction of Interspecies Interactions from Self-organized Spatiotemporal Patterns of Co-evolving Organisms

Microorganisms colonizing natural habits such as soils co-evolve to form specific spatial patterns through interspecies interactions. These self-organized patterns are a key ecological phenotype, which provides critical information on their interaction mechanisms. However, conventional network inference techniques that analyze species population data in bulk have yet to be extended to account for such spatial heterogeneity. Here we proposed supervised deep learning as a new network inference tool for predicting interspecies interactions from spatiotemporal patterns of microbial evolution. Due to lack of biological imaging data that can be used for training deep learning networks, we used in silico data generated from high-fidelity agent-based models to determine model structure and parameters. Even though networks were trained under simple configurations where interaction coefficients are assumed to be spatially invariant, we demonstrated that the resulting model can be utilized to successfully predict spatial variation of interactions in more complex domains (i.e., configured with a context-dependent mixture of interaction coefficients) as well as in simple domains without further training. In the further test against real biological data obtained through imaging experiments of a binary consortium (Pseudomonas fluorescens and a mutant of Escherichia coli), our model also predicted the dramatic shifts in interactions of the two organisms across different environmental contexts. Through various successful demonstrations in this work, the combined use of the agent-based model and machine learning algorithm provides a means to use new type of data - microscopic images - for extracting microbial interactions, therefore presenting itself as a useful tool for the analysis of more complex microbial community interactions.

Lee, Joon-Yong↗
DOE OSTI2022

Impact of the plastic deformation microstructure in metals on the kinetics of recrystallization: A phase-field study

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. Here in the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Buried interface modulation via $\mathrm{PEDOT:PSS}$ ionic exchange for the Sn-Pb mixed perovskite based solar cells

To apply Sn-Pb mixed perovskite solar cells for highly efficient single- or multi-junction devices, understanding device-specific buried interfaces is necessary. Poly [3,4-ethylenedioxythiophene]:poly[styrene sulfonate] (PEDOT:PSS) is primarily used as a hole transport layer in Sn-Pb mixed perovskite solar cells. However, the spatial heterogeneity of PEDOT:PSS, caused by its PEDOT-rich and PSS-rich domains, induces many defects at the buried interface in PEDOT:PSS/perovskite, which limits device performance. Here, we present ionic exchange (IE) of PEDOT:PSS via a combination of methylamine iodide (MAI) and dimethyl sulfoxide (DMSO). Through surface analyses and density functional theory (DFT) simulations, we confirm that the IE process preferentially form PEDOT-I and MA-PSS and that PSS-rich domains bind to DMSO. Thus, the spatial separation of PEDOT:PSS is solved, and the exchanged MA + and I - ions serve as a bridge between PEDOT:PSS and the perovskite, leading to improved physical, chemical, and electrical properties of the buried interface. The Sn-Pb mixed perovskite solar cells using IE-PEDOT:PSS achieve an improved efficiency of 21.3% with an open-circuit voltage of 0.85 V and show better long-term stability. Additionally, IE-PEDOT:PSS works effectively in 2-terminal all-perovskite tandem devices, resulting in an improved efficiency of 23.5% and high reproducibility.

14 SOLAR ENERGY↗
DOE OSTI2023

Echogenic Segmentation for Ultrasonic Measurements of Spatially Distributed Properties in Solids

An ultrasonic (US) pulse propagating between transducer-receiver responds cumulatively to distributed material properties it encounters, such as temperature distribution. Such convolution makes the unique reconstitution of distributions impossible short of applying tomographic techniques. However, spatial heterogeneity in material properties is readily revealed by segmenting the US propagation path with echogenic features (EF), producing a train of echoes delayed as a function of segmentally localized properties. Unfortunately, in low attenuating waveguides (WGs), primary and trailing echoes (TE) interfere, limiting the number, density, and placement of EFs we can use. We report the design guidelines for creating US-segmented WGs for measuring distributed properties, in which the impact of TEs is minimized or eliminated. Two waveguide designs using these principles have been tested in spatially distributed temperature measurements in high-temperature or restrictive environments.

Walton, Kenneth↗
DOE OSTI2023

Integrating LEO and GEO Observations: Toward Optimal Summertime Satellite Precipitation Retrieval

Abstract Reliable quantitative precipitation estimation with a rich spatiotemporal resolution is vital for understanding the Earth’s hydrological cycle. Precipitation estimation over land and coastal regions is necessary for addressing the high degree of spatial heterogeneity of water availability and demand, and for resolving the extremes that modulate and amplify hazards such as flooding and landslides. Advancements in computation power along with unique high spatiotemporal and spectral resolution data streams from passive meteorological sensors aboard geosynchronous Earth-orbiting (GEO) and low Earth-orbiting (LEO) satellites offer exciting opportunities to retrieve information about surface precipitation phenomena using data-driven machine learning techniques. In this study, the capabilities of U-Net–like architecture are investigated to map instantaneous, summertime surface precipitation intensity at the spatial resolution of 2 km. The calibrated brightness temperature products from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) radiometer are combined with multispectral images (visible, near-infrared, and infrared bands) from the Advanced Baseline Imager (ABI) aboard the GOES-R satellites as main inputs to the U-Net–like precipitation algorithm. Total precipitable water and 2-m temperature from the Global Forecast System (GFS) model are also used as auxiliary inputs to the model. The results show that the U-Net–like algorithm can capture fine-scale patterns and intensity of surface precipitation at high spatial resolution over stratiform and convective precipitation regimes. The evaluations reveal the potential of extracting relevant, high spatial features over complex surface types such as mountainous regions and coastlines. The algorithm allows users to interpret the inputs’ importance and can serve as a starting point for further exploration of precipitation systems within the field of hydrometeorology.

Meteorology & Atmospheric Sciences↗
NASA NTRS1993

Simulated retrieval of land surface parameters using multichannel satellite microwave radiometry

The current global change research emphasis on understanding water and energy fluxes at the land-atmosphere interface provided renewed interest in using passive microwave satellite data for land studies. Radiative transfer models of microwave emission and scattering in the soil-vegetation-atmosphere column, though still in the process of development, can be used to estimate the accuracies with which land surface parameters can be derived from satellite data. These parameters include surface soil moisture, surface temperature, vegetation water content, and atmospheric water content. A microwave radiative transfer model is used to develop linear and nonlinear versions of a multichannel retrieval algorithm to simulate retrievals of the surface and atmospheric parameters. These simulations include the effect of additive noise, and examine the effects of nonlinearities in the radiative transfer models. However, effects of spatial heterogeneity are not considered. Satellite data to which the retrieval algorithms may be applied include the SSMR and SSM/I, and in future will include the Multifrequency Imaging Microwave Radiometer (MIMR) to be launched as part of the Earth Observing System (EOS). The Scanning Multichannel Microwave Radiometer (SMMR) on the Nimbus-7 satellite provided data from 1978 to 1987 at five microwave frequencies between 6.6 and 37 GHz. The series of Special Sensor Microwave/Imager (SSM/I) instruments launched on the DMSP satellites provide data from 1987 to the present at four frequencies between 19.35 and 8.5 GHz. Spatial resolutions of the data range from approximately 12 to 120 km depending on the frequency. Examples using data from these sensors are shown to indicate the results of applying retrieval algorithms based on model simulations to real data.

Njoku, Eni G.↗
DOE OSTI2022

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Uncertainty quantification and propagation in lithium-ion battery electrodes using bayesian convolutional neural networks

The complex nature of manufacturing processes stipulates electrodes to possess high variability with increased heterogeneity during production. X-ray computed tomography imaging has proved to be critical in visualizing the complicated stochastic particle distribution of as-manufactured electrodes in lithium-ion batteries. However, accurate prediction of their electrochemical performance necessitates precise evaluation of kinetic and transport properties from real electrodes. Image segmentation that characterizes voxels to particle/pore phase is often meticulous and fraught with subjectivity owing to a myriad of unconstrained choices and filter algorithms. Here we utilize a Bayesian convolutional neural network to tackle segmentation subjectivity and quantify its pertinent uncertainties. Otsu inter-variance and Blind/Referenceless Imaging Spatial Quality Evaluator are used to assess the relative image quality of grayscale tomograms, thus evaluating the uncertainty in the derived microstructural attributes. We analyze how image uncertainty is correlated with the uncertainties and magnitude of kinetic and transport properties of an electrode, further identifying pathways of uncertainty propagation within microstructural attributes. The coupled effect of spatial heterogeneity and microstructural anisotropy on the uncertainty quantification of transport parameters is also understood. This work demonstrates a novel methodology to extract microstructural descriptors from real electrode images through quantification of associated uncertainties and discerning the relative strength of their propagation, thus facilitating feedback to manufacturing processes from accurate image based electrochemical simulations.

25 ENERGY STORAGE↗
DOE OSTI2023

Relationships between summertime surface albedo and melt pond fraction in the central Arctic Ocean: The aggregate scale of albedo obtained on the MOSAiC floe

As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), the HELiX uncrewed aircraft system (UAS) was deployed over the sea ice in the central Arctic Ocean during summer 2020. Albedo measurements were obtained with stabilized pyranometers, and melt pond fraction was calculated from orthomosaic imagery from a surface-imaging multispectral camera. This study analyzed HELiX flight data to provide insights on the temporal and spatial evolution of albedo and melt pond fraction of the MOSAiC floe during the melt season as it drifted south through Fram Strait. The surface albedo distributions showed peak values changing from high albedo (0.55–0.6) to lower values (0.3) as the season advanced. Inspired by methods developed for satellite data, an algorithm was established to retrieve melt pond fraction from the orthomosaic images. We demonstrate that the near-surface observations of melt pond fraction were highly dependent on sample area, offering insight into the influence of subgrid scale features and spatial heterogeneity in satellite observations. Vertical observations conducted with the HELiX were used to quantify the influence of melt pond scales on observed surface albedo as a function of sensor footprint. These scaling results were used to link surface-based measurements collected during MOSAiC to broader-scale satellite data to investigate the influence of surface features on observed albedo. Albedo values blend underlying features within the sensor footprint, as determined by the melt pond size and concentration. This study framed the downscaling (upscaling) problem related to the airborne (surface) observations of surface albedo across a variety of spatial scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Heterogeneity of the Dominant Causes of Performance Loss in End-of-Life Cathodes and Their Consequences for Direct Recycling

Recycling Li-ion batteries from electric vehicles is critical for reducing costs and supporting the development of a domestic battery supply chain. Direct recycling of cathodes, like LiNixMnyCozO2 (NMC), is attractive due to its low cost, energy use, and emissions compared to traditional recycling techniques. However, a comprehensive understanding of the active material properties at end-of-life is needed to guide direct recycling processes and the performance-dependent reuse applications. Here, NMC material from an end-of-life commercial pouch cell is characterized and bench-marked against pristine non-cycled counterparts with respect to capacity, impedance, crystallography, morphology, and microstructure to identify major degradation modes and understand variability in the end-of-life material. The spatial heterogeneity of each property throughout the cell is also quantified. While the degraded material demonstrated similar capacity as the pristine, its impedance and rate capability are severely diminished. Furthermore, samples from the periphery of the electrode layers showed more severe performance loss compared to samples extracted from central regions. The dominant culprit of performance loss is the material microstructure, where the magnitude of particle cracking showed the strongest correlation to the impedance components that are most unfavorably impacted. This work suggests severe cracks in cathode active materials are the primary challenge that direct recycling methods must overcome.

25 ENERGY STORAGE↗
DOE OSTI2026

Selective Bidentate Coordination Reconstructs Residual PbI 2 to Homogenize Interfacial Energetics in Perovskite Solar Cells

Spatially heterogeneous interfacial energetics, often originating from residual lead iodide (PbI 2 ), represent a fundamental bottleneck to both the efficiency and operational stability of perovskite photovoltaics. Conventional PbI 2 passivation strategies based on monodentate ligands or highly polar solvents either interact weakly with PbI 2 or undesirably perturb the underlying three-dimensional perovskite lattice. Here, we report a diammonium bidentate strategy that selectively reconstructs residual PbI 2 into corner-sharing PbI 6 octahedra while preserving the bulk perovskite. Enabled by dual-site coordination, a newly designed thiophene-based bidentate ligand (MeXT) stabilizes PbI 6 units during passivation and establishes spatially homogeneous interfacial energetics. Perovskite solar cells incorporating MeXT achieve 26.19% power conversion efficiency and retain over 80% of their initial efficiency after 1000 h of continuous 1-sun illumination at 75 °C. This dual-anchor coordination passivation strategy establishes a general design principle for selectively treating residual PbI 2 and creating electronically coherent and operationally stable interfaces in perovskite photovoltaics.

ligands↗
DOE OSTI2026

Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone

Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.

58 GEOSCIENCES↗