Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity estimation”

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

Supernova electron-neutrino interactions with xenon in the nEXO detector

Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-𝛽 decay detector (∼5 metric ton, 90% 136 Xe, 10% 134 Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the Model of Argon Reaction Low Energy Yields (MARLEY) event generator. We find good agreement between MARLEY’s predictions and existing theoretical calculations of the inclusive cross sections at supernova neutrino energies. The interactions modeled by MARLEY were simulated within the nEXO simulation framework and were run through an example reconstruction algorithm to determine the detector’s efficiency for reconstructing these events. The simulated data, incorporating the detector response, were used to study the ability of nEXO to reconstruct the incident electron-neutrino spectrum and these results were extended to a larger xenon detector of the same isotope enrichment. We estimate that nEXO will be able to observe electron-neutrino interactions with xenon from supernovae as far as 5–8 kpc from Earth, while the ability to reconstruct incident electron-neutrino spectrum parameters from observed interactions in nEXO is limited to closer supernovae.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Study of Key 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu Resonances to Understand 44 Ti Nucleosynthesis in Supernovae

An important validation of nucleosynthesis models of core-collapse supernovae is the comparison of radioisotope predictions to abundances inferred from observations of 𝛾 rays emitted in remnants. One such isotope, 44 Ti, is especially sensitive to the 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu reaction rate. Despite this importance, no experimentally constrained rate exists for this reaction. It is thus crucial to determine this rate. GRETINA ORRUBA: Dual Detectors for Experimental Structure Studies (GODDESS) was used along with the first use of the Enge split-pole spectrograph at Notre Dame to measure structure properties of 58 Cu via the 58 Ni ⁢( 3 He,𝑡)⁢ 58 Cu reaction. Here, the combined analysis of these complementary data sets allows for precise determination of 58 Cu level energies, and spins are constrained from analysis of the 𝛾-ray decay. 44 Ti yields were found to change by over 25% compared to previous estimates in model calculations using the new experimentally constrained reaction rate.

energy levels↗
DOE OSTI2026

Cross-correlation of SPT-3G D1 CMB lensing and DES Y3 galaxy lensing

Measurements of the weak lensing of galaxies and of the cosmic microwave background (CMB) provide direct probes of the cosmic matter density field, but the two observables are sensitive to different spatial scales, redshift ranges, and survey systematics. Their cross-correlation thus enables consistency checks of the theoretical model and of potential systematics in either dataset. We present measurements of the cross-correlation between CMB lensing and cosmic shear over $\sim$1,300 deg$^2$ of the sky using the SPT-3G D1 CMB lensing maps and the Dark Energy Survey Year 3 (DES Y3) shear catalogs. For the first time, we measure this cross-correlation at high significance ($\sim 14σ$) when using a polarization-only CMB lensing reconstruction that is expected to be robust against biases induced by extragalactic foregrounds. We test a variety of other CMB lensing estimators that include temperature information and exhibit different tradeoffs between foreground biases and noise, as well as a shear sample that consists of blue, star-forming galaxies and has been shown to be less impacted by galaxy intrinsic alignments. Assuming $Λ$CDM and marginalizing over uncertainties in intrinsic alignments, baryonic feedback, and various nuisance parameters, we obtain a constraint on the amplitude of matter clustering $S_8 \equiv σ_8 \sqrt{Ω_m / 0.3} = 0.833^{+0.047}_{-0.061}$, consistent with both the primary CMB results from Planck and shear-only results from DES Y3. By combining our measurement with Planck, we find mild constraints on the astrophysical processes that impact the cross-correlation. We obtain a constraint on the intrinsic alignment amplitude of the DES sample that is competitive with that from shear-only analyses, and we find a lower limit on the strength of baryonic feedback.

Ouellette, A. [Illinois U., Urbana] (ORCID:0000000↗
DOE OSTI2024

Soil application of high-lignin fermentation byproduct to increase the sustainability of liquid biofuel production from crop residues

When digestates from anaerobic digestion of crop residues are added to soil, a considerable body of information indicates that soil organic carbon (SOC) levels are comparable to those when crop residues are left in the field. This occurs although the amount of digestate added to soil is diminished by digestion and implies that digestion increases the proportion of carbon inputs stabilized as SOC. Here we examine the likelihood and implications of these features being manifested for soil application of high lignin-fermentation byproduct (HLFB) from liquid biofuel production. We show that steady-state SOC levels are much less sensitive to crop residue removal with HLFB return than without it, and provide an example supporting the feasibility of foregoing process energy and coproduct revenue when HLFB is returned to the soil. Informed by this review and analysis, we expect with moderate confidence that long-term SOC levels for soils amended with HLFB from some liquid cellulosic biofuel processes will not be substantially lower than those occurring when crop residues are left in the field. We have high confidence that the economically optimum rate of fertilizer nitrogen (N) application and N 2 O emissions will be lower at most sites for HLFB return to the soil than if crop residues were left in the field. We estimate that the per hectare N demand for processing crop residues to liquid biofuels is about a third of the per hectare demand for crop production, giving rise to an opportunity to use N twice and thereby realize cost savings and environmental benefits. These observations support but do not prove the hypothesis that a ‘win-win’ is possible wherein large amounts of liquid biofuel feedstock can be obtained from crop residues while improving the economics and sustainability of food and feed production. A research agenda aimed at exploring and testing this hypothesis is offered.

09 BIOMASS FUELS↗
DOE OSTI2025

Quantifying the Impact of LSST u -band Survey Strategy on Photometric Redshift Estimation and the Detection of Lyman-break Galaxies

The Vera C. Rubin Observatory will conduct the Legacy Survey of Space and Time (LSST), promising to discover billions of galaxies out to redshift 7, using six photometric bands (ugrizy) spanning the near-ultraviolet to the near-infrared. The exact number of and quality of information about these galaxies will depend on survey depth in these six bands, which in turn depends on the LSST survey strategy, i.e., how often and how long to expose in each band. u-band depth is especially important for photometric redshift (photo-z) estimation and for detection of high-redshift Lyman-break galaxies (LBGs). In this paper, we use a simulated galaxy catalog and an analytic model for the LBG population to study how recent updates and proposed changes to Rubin’s u-band throughput and LSST survey strategy impact photo-z accuracy and LBG detection. We find that proposed variations in u-band strategy have a small impact on photo-z accuracy for z < 1.5 galaxies, but the outlier fraction, scatter, and bias for higher-redshift galaxies vary by up to 50%, depending on the survey strategy considered. The number of u-band dropout LBGs at z ∼ 3 is also highly sensitive to the u-band depth, varying by up to 500%, while the number of griz-band dropouts is only modestly affected. Under the new u-band strategy recommended by the Rubin Survey Cadence Optimization Committee, we predict u-band dropout number densities of 110 deg −2 (3200 deg −2 ) in year 1 (10) of LSST. We discuss the implications of these results for LSST cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advanced Instrumentation for Metal Additive Manufacturing

Laser powder bed fusion (LPBF) is the most widely used process for metal additive manufacturing (AM), particularly where complex geometries provide performance advantages unattainable with traditional manufacturing techniques. However, LPBF is highly sensitive to innate variability in both the powder spreading and fusion steps, often leading to defects such as pores that are difficult to detect yet significantly impair component mechanical properties and fatigue life. This thesis presents a range of novel instruments enabling both precise assessment of powder layer characteristics and in-situ thermal metrology of metal AM to advance the quality control of LPBF. First, leveraging a custom X-ray microscope and a radiation-transport model developed through this work, transmission X-ray imaging is used to study spreading of thin metal powder layers. Effective layer depth is directly mapped at a process-relevant size scale, surpassing optical techniques that can only estimate local deposition from layer surface topography. Layer packing density and quality are shown to be influenced by powder flowability and particle size relative to nominal powder layer thickness. Layer quality is additionally connected to the geometry of the spreading implement and its velocity. This technique and its presented findings enable pairing feedstocks with spreading strategies that create layers with consistent packing density and uniformity. Second, a twofold approach is employed to optically interrogate the laser fusion step of LPBF for observing signatures of defect formation. Aperture division multiplexing is conceptualized, providing for simultaneous laser delivery and high-fidelity infrared (IR) process monitoring through a common optic. In-situ microscopy at 50 μm spatial resolution and at mid-wave IR wavelengths is proven readily achievable with the first purpose-built optic of this type. Next, a bespoke imaging spectrometer, along with a temperature-emissivity separation technique, is used to retrieve accurate process temperatures over a 1000 K range. Data from these instruments are correlated to porosity as fine as 4.3 μm in two LPBF test artifacts, as verified using computed tomography (CT), establishing the viability of robust optically-based component qualification.

Penny, Ryan↗
DOE OSTI2026

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Optimized Viewing Techniques to Minimize Radiation Damage From X-ray Imaging Systems

X-ray inspection of ball grid arrays (BGAs) is typically performed at one or more viewing angles to examine adhesion sites for errors such as voids, joint cracking, or head-in-pillow. During this inspection process, the cir cuit board assembly is subject to ionizing radiation exposure, which can cause trapped charge within oxide layers of semiconductor devices. Some x-ray machines allow for programmable inspection routines, which could be used to optimize radiation exposure to semiconductor components. Using Monte Carlo methods, x-ray inspection of a BGA was simulated to determine a range of acceptable viewing angles. Dose rates to circuit board components were estimated at each inspection angle to determine the view resulting in optimized radiation exposure. Results showed that for each BGA, the maximum unobstructed viewing times without exceeding a 5 Gy dose limit to a single part ranged from 82 to 94 minutes. Using a radiation cost function method, optimized viewing across all components was found. Here, it was observed that for a consistent dose limit applied to silicon-based components, performing inspection with BGAs facing the x-ray source was optimal. A third method was applied, assigning individual dose limits based on empirical data from the NASA Goddard Space Flight Center radiation database. This method showed that optimized viewing maximizes the distance between the radiation source and highly sensitive components. It was also observed that cumulative effects from viewing two BGAs will influence viewing angles, causing the optimal view of one BGA to exist nearly 180° from the other.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Electromagnetic Radiation from Cosmic-Ray Scatterings on Relic Neutrinos

Cosmic-ray scatterings on the cosmic neutrino background induce a flux of gamma-rays and X-rays from boosted meson decays and charged lepton processes. Here we present the first estimate of this flux and its cumulative cosmological contribution. Confronting expectations with Fermi-LAT diffuse gamma-ray data, we find a limit on the cosmic neutrino background overdensity at the level of $η\lesssim 2 \times 10^{4}$ for a lightest neutrino mass of $m_ν\gtrsim0.1$ eV, orders of magnitude stronger than current direct laboratory probes, and comparable to constraints on the cosmic neutrino background obtained with IceCube. We further show that X-ray synchrotron emission from cascade electron-positron pairs in intergalactic magnetic fields provides a complementary, albeit weaker, constraint. We discuss how anisotropic signatures and future gamma-ray data from CTA could further improve bounds on the relic neutrino overdensity, approaching in sensitivity the $Λ$CDM expected density.

Herrera, Gonzalo [MIT, MKI; Harvard U.] (ORCID:000↗
DOE OSTI2025

Discerning structure sensitivity of Ni-core Pd-shell nanoparticles for enhanced nitrite reduction

Denitrification using Pd-based catalysts converts N-oxyanions into harmless products. However, the metal cost and scarcity are barriers to application. In this work, we synthesize water-suspended, structure-controlled Ni-core/Pd-shell nanoparticles (Ni@Pd NPs) that show Pd-catalyzed nitrite (NO 2 − ) reduction catalysis promoted by Ni. The NPs immobilized on Al 2 O 3 (0.47 wt % Pd and 0.74 wt % Ni) have higher NH 4 + selectivity and higher catalytic activity than wet-impregnated Pd/Al 2 O 3 (1.0 wt % Pd), i.e., 315.4 vs. 132.1 Lg surface Pd −1 min −1 . X-ray photoelectron spectroscopy (XPS) and CO-diffuse reflectance infrared Fourier transform spectroscopy (CO-DRIFTS) analyses show that Ni increases Pd electron density, which is proposed to lower the bond dissociation energy barrier for nitric oxide surface intermediate. Low-coordinated Pd atoms correlate to high activity and low NH 4 + selectivity. Using CatCost software, we estimate Ni@Pd/Al 2 O 3 to have a ∼40% lower manufacturing cost compared to wet-impregnated Pd/Al 2 O 3 . This rationally designed catalyst illustrates how nickel promotes Pd catalysis and how Pd usage can be lowered for denitrification and other hydrogenation reactions.

CO-DRIFTS↗
DOE OSTI2026

Surface Observations From Atmospheric Radiation Measurement Sites Constrain the Anthropogenic Contribution to Cloud Droplet Number

Uncertainty in anthropogenic forcing driven by aerosol-cloud interactions (aci) limits our ability to infer the sensitivity of the Earth system to forcing from historical records. The driver of aci is the change in cloud droplet number concentration (N d ) due to changes in aerosol serving as cloud condensation nuclei (CCN). Here, we combine a perturbed parameter ensemble run in a global Earth system model with observations of CCN and single-layer-cloud N d at surface sites in the Azores, the Southern Great Plains, and Ascension Island to provide a constraint on the anthropogenic contribution to present-day N d . These observational lines of evidence constrain the preindustrial to present-day change in N d to be between 11 and 43 cm −3 . This is consistent with the upper end of some previous estimates but has a higher minimum perturbation, pointing to a stronger historical aerosol cooling.

ARM↗
DOE OSTI2025

Transformation rate maps of dissolved organic carbon in the contiguous US

Riverine dissolved organic carbon (DOC) plays a vital role in regional and global carbon cycles. However, the processes of DOC conversion from soil organic carbon (SOC) and leaching into rivers are insufficiently understood, inconsistently represented, and poorly parameterized, particularly in land surface and Earth system models. As a first attempt to fill this gap, we propose a generic formula that directly connects SOC concentration with DOC concentration in headwater streams, where a single parameter, the transformation rate from SOC in the soil to DOC leaching flux (P r ), accounts for the overall processes governing SOC conversion to DOC and leaching from soils (along with runoff) into headwater streams. We then derive high-resolution P r maps over the contiguous US (CONUS) using SOC data from two different sources: the Harmonized World Soil Database v1.2 (HWSD) and SoilGrids 2.0. Both maps are developed following the same five major steps: (1) selecting independent catchments where observed riverine DOC data are available with reasonable quality; (2) estimating catchment-average SOC for the independent catchments; (3) estimating the P r values for these catchments based on the generic formula and catchment-average SOC; (4) developing a predictive model of P r with machine learning (ML) techniques and catchment-scale climate, hydrology, geology, and other attributes; and (5) deriving a national map of P r based on the ML model. For evaluation, we compare the DOC concentration derived using the P r map and the observed DOC concentration values at evaluation catchments. The resulting mean absolute scaled error and coefficient of determination are 0.73 and 0.47 for the HWSD-based model and 0.58 and 0.72 for the SoilGrids-based model, respectively, suggesting the effectiveness of the overall methodology. Efforts to constrain uncertainty and evaluate sensitivity of P r to different factors are discussed. To illustrate the use of such maps, we derive a riverine DOC concentration reanalysis dataset over CONUS. The two P r maps, robustly derived and empirically validated, lay a critical cornerstone for better simulating the terrestrial carbon cycle in land surface and Earth system models. Our findings not only set a foundation for improving our predictive understanding of the terrestrial carbon cycle at the regional and global scales, but also hold promises for informing policy decisions related to decarbonization and climate change mitigation. The data presented in this study are publicly available at https://doi.org/10.5281/zenodo.14563816 (Li et al., 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simulation of wind and solar energy generation over California with E3SM SCREAM regionally refined models at 3.25 km and 800 m resolutions

This study presents wind and solar power generation estimates derived from the US Department of Energy’s Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) Regionally Refined Models (RRM) over California at 3.25 km and 800 m horizontal resolutions, using the Python wrapper for the System Advisor Model (PySAM). The resulting wind and solar generation estimates are compared to monthly capacity factors reported to the Energy Information Administration (EIA), High-Resolution Rapid Refresh (HRRR, 3 km resolution) forecast model, and E3SM North American regionally refined model (NARRM, 25 km resolution). We systematically assess the impacts of generation modeling assumptions, meteorological models, and horizontal resolution. Results show that resolution plays a dominant role for wind energy: increasing from 25 to 3.25 km brings qualitative and quantitative improvements, most notably by resolving the phase error in the seasonal cycle found in coarser simulations. However, further refinement to 800 m offers minimal gains. SCREAM performs better than HRRR for solar power generation in single- and dual-axis tracking systems, likely due to more accurate surface radiation. The sensitivity of PySAM to system configuration, particularly for axis-tracking modeling in photovoltaics, is also highlighted. Overall, SCREAM-RRM shows strong potential for high-resolution energy assessments, with future progress depending on more in situ observations and clearer quantification of uncertainties in generation modeling.

Geosciences↗
DOE OSTI2026

Feynman diagrams for matter wave interferometry

We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.

Glick, Jonah [Northwestern U.; Fermilab] (ORCID:00↗
DOE OSTI2025

Assessing the impact of climate change on rainfall-triggered landslides: a case study in California

Landslides are widespread natural hazards which take a heavy toll on lives, property, and infrastructure each year. In January 2023 and 2024, widespread landslides occurred in California as a result of extreme winter precipitation. Estimating the regional evolution of landslide hazard in a changing climate is essential for adaptation planning and risk mitigation efforts. The complex effects of climate change on landslide hazard, however, are poorly understood. Here, we use climate and landslide observations to develop a novel data-driven approach for landslide susceptibility assessment under historical and future climate conditions. We apply this framework to investigate the evolution of landslide susceptibility in California and its impacts on the roadway system based on downscaled climate projections under moderate and high greenhouse gas emission scenarios. On average, the results indicate an increase in landslide susceptibility through 2100 in regions where historical susceptibility is moderate to very high. Although the spatio-temporal variations in landslide susceptibility are sensitive to climate model uncertainties in predicting extreme rainfall, the direst impacts are largely circumvented in lower emission scenarios.

Climate Change↗
DOE OSTI2026

HDSense: An efficient method for ranking observable sensitivity

Identifying which observables most effectively constrain model parameters can be computationally prohibitive when considering full likelihoods of many correlated observables. This is especially important for, e.g., hadronization models, where high precision is required to interpret the results of collider experiments. We introduce the High-Dimensional Sensitivity (HDSense) score, a computationally efficient metric for ranking observable sets using only one-dimensional histograms. Derived by profiling over unknown correlations in the Fisher information framework, the score balances total information content against redundancy between observables. We apply HDSense to rank a set observables in terms of their constraining power with respect to five parameters of the Lund string model of hadronization implemented in Pythia using simulated leptonic collider events at the $Z$ pole. Validation against machine-learning--based full-likelihood approximations demonstrates that HDSense successfully identifies near-optimal observable subsets. The framework naturally handles data from multiple experiments with different acceptances and incorporates detector effects. While demonstrated on hadronization models, the methodology applies broadly to generic parameter estimation problems where correlations are unknown or difficult to model.

Assi, Benoît [Cincinnati U.] (ORCID:00000003092433↗
DOE OSTI2025

Characterization of proton-induced damage in thick, p-channel skipper-CCDs

Here, in this work, we characterize the radiation-induced damage in two thick, p-channel skipper-CCDs irradiated unbiased and at room temperature with 217-MeV protons. We evaluate the overall performance of the sensors and demonstrate their single-electron/single-photon sensitivity after receiving a fluence on the order of 10 10 protons/cm 2 . Using the pocket-pumping technique, we quantify and characterize the proton-induced defects from displacement damage. We report an overall trap density of 0.134 traps/pixel for a displacement damage dose of 2.3 × 10 7 MeV/g. Three main proton-induced trap species were identified, V 2 , C i O i and V n O m , and their characteristic trap energies and cross sections were extracted. We found that while divacancies are the most common proton-induced defects, C i O i defects have a greater impact on charge integrity at typical operating temperatures because their emission-time constants are comparable or larger than typical readout times. To estimate ionization damage, we measure the characteristic output transistor curves. We found no threshold voltage shifts after irradiation. Our results highlight the potential of skipper-CCDs for applications requiring high-radiation tolerance and can be used to find the operating conditions in which effects of radiation-induced damage are mitigated.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Compression–tension cell with sample manipulator for in situ X‐ray nanotomography experiments

In situ X-ray nanotomography experiments where tensile or compressive force is applied on the sample require specialized equipment. A compression-tension device with fluid flow-through capability has been designed for X-ray nanotomography beamlines. The compression-tension cell is equipped with a triaxial stage for sample alignment and a high sensitivity loadcell for measurement of applied force. To handle the <100 µm samples used for X-ray nanotomography imaging and for loading samples on the compression-tension cell a sample manipulator has been built. The sample manipulator is capable of selecting a single <100 µm particle for nanotomography scanning while viewing multiple samples under an optical microscope. To test the functionality of these two devices an initial compression experiment involving two glass beads was performed. To demonstrate instrument stability two spherical glass beads were compressed from a no load condition until one of the beads fractured. Nanotomography data were collected at each step of increasing compressive force. The experimentally observed contact area of the spherical glass beads was compared with the theoretical estimate using the Hertz analysis. To demonstrate the fluid flow capability, two calcite grains were compressed against each other under a calcite saturated solution. Surface topological changes were observed for the stressed grain contact area.

X-ray tomography↗