Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface modeling”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of nanoconfinement and surface charge on iron adsorption on mesoporous silica

Here, we present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(II) adsorption onto negatively charged SiO 2 surfaces and Fe(III) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(II) at deprotonated surface sites (SiO - ), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Molecular Models of Atomically Dispersed Uranium at MoS 2 Surfaces Reveal Cooperative Mechanism of Water Reduction

Single atoms of uranium supported on molybdenum sulfide surfaces (U@MoS 2 ) have been recently demonstrated to facilitate the hydrogen evolution reaction (HER) through electrocatalysis. Theoretical calculations have predicted uranium hydroxide moieties bound to edge-sulfur atoms of MoS 2 as a proposed transition state involved in the HER process. However, the isolation of relevant intermediates involved in this process remains a challenge, rendering mechanistic hypotheses unverified. The present work describes the isolation and characterization of a uranium-hydroxide intermediate on molybdenum sulfide surfaces using [(Cp* 3 Mo 3 S 4 )UCp*], a molecular model of a reduced uranium center supported at MoS 2 . Mechanistic investigations highlight the metalloligand cooperativity with uranium involved in the water activation pathway. The corresponding uranium-oxo analogue, [(Cp* 3 Mo 3 S 4 )Cp*U(═O)], was also accessed from the hydroxide cluster via hydrogen atom transfer and from [(Cp* 3 Mo 3 S 4 )UCp*] through an alternative direct oxygen atom transfer. These results provide an atomistic perspective on the reactivity of low-valent uranium at molybdenum sulfide surfaces toward water, modeling key intermediates associated with the HER of U@MoS 2 catalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

A Hybrid Biophysical‐Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy ( LE ) and sensible heat ( H ) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R 2 = 0.81–0.94) and H (R 2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

evapotranspiration↗
DOE OSTI2026

Data for: A hybrid biophysical-machine learning framework for diurnal surface energy flux estimation using proximal sensing

Thermal-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal datasets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for specific surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of an ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81-0.94) and H (R2 = 0.46-0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical – machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

Agricultural Sciences↗
DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗
DOE OSTI2025

Estimating evapotranspiration from deciduous broadleaf forests by considering the effect of foliar litter

Forest litter cover significantly alters the ground surface resistance, which influences the water and energy transfer processes between the land surface and atmosphere. Penman-Monteith (PM) model is widely used to estimate land surface evapotranspiration (ET). However, it was shown to overestimate ET during the period when foliar litter covers the ground surface. Therefore, incorporating the effect of foliar litter on ground surface resistance can potentially improve ET estimates. Here, in this study, we proposed a foliar litter surface resistance model to describe the effect of litter cover on water vapor transport, then incorporated the litter surface resistance into the PM model (noted as PM-EL model). The performance of the PM-EL model was evaluated using observations from 18 deciduous broadleaf forest (DBF) flux sites across global FLUXNET2015 datasets. Results showed that: (1) both the PM model and the PM-EL model are most sensitive to ground surface resistance among resistance parameters; (2) the proposed litter resistance model is capable of describing the seasonal dynamic of litter in DBF ; (3) the PM-EL model, has been shown to improve ET estimates with the coefficient of determination increased by 8.5 %, the mean absolute error and relative root mean square error decreased by 22.1 % and 20.9 %, respectively. The study contributes to our understanding of the effect of litter cover on ET estimation and provides insights into forest ecohydrology and land surface mass and energy transfer processes.

Deciduous Broadleaf Forest↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2024

Quantifying the Impacts of Land-Cover Change on the Hydrologic Response to Hurricane Ida in the Lower Mississippi River Basin

Abstract The Lower Mississippi River basin (LMRB) has experienced significant changes in land cover and is one of the most vulnerable regions to hurricanes in the United States. Here, we study the impacts of land-cover change on the hydrologic response to Hurricane Ida in LMRB. By using an integrated surface–subsurface hydrologic model, Energy Exascale Earth System Model (E3SM) Land Model coupled with the three-dimensional ParFlow subsurface flow model (ELM-ParFlow), we simulate the effects of land-cover change on the flood volume and peak timing induced by rainfall from Hurricane Ida. The results show that land-cover changes from 1850 to 2015, which resulted in a smoother surface and less vegetation, exacerbated both flood peak time and volume induced by Hurricane Ida. The effects of land-cover changes can be decomposed into two mechanisms: a smoother surface routes more water faster to a watershed outlet and less vegetation allows more water to contribute to surface runoff. By comparing scenarios in which the two mechanisms were isolated, we found that changes in soil moisture due to vegetation cover change have more dominant effects on floods in the southern part and changes in Manning’s coefficient have the largest effect on floods in the northern part of the LMRB. The study provides important insights into the complex relationship between land-use, land-cover, and hydrologic processes in coastal regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Computer-aided design file format for additive manufacturing and methods of file generation

A tessellated output file format describing a computer-aided design (CAD) model including a plurality of vertices, a plurality of curves, a plurality of surfaces, and at least one volume, the tessellated output file format comprising base polygon data, a vertex metadata container, a curve metadata container, and a surface metadata container. The base polygon data defines the CAD model for additive manufacturing and includes a plurality of connected polygons, each polygon including a plurality of nodes, a plurality of edges, and a face. The vertex metadata container includes a listing of CAD model vertices and one polygon node associated with each CAD model vertex. The curve metadata container includes a listing of CAD model curves and at least one polygon edge associated with each CAD model curve. The surface metadata container includes a listing of CAD model surfaces and at least one polygon face associated with each CAD model surface.

Vernon, Gregory John↗
DOE OSTI2023

Biogeophysical Effects of Land-Use and Land-Cover Change Not Detectable in Warmest Month

Land-use and land-cover changes (hereafter simply “land use”) alter climates biogeophysically by affecting surface fluxes of energy and water. Yet, near-surface temperature responses to land use across observational versus model-based studies and spatial-temporal scales can be inconsistent. Here we assess the prevalence of the historical land use signal of daily maximum temperatures averaged over the warmest month of the year (t LU ) using regularized optimal fingerprinting for detection and attribution. We use observations from the Climatic Research Unit and Berkeley Earth alongside historical simulations with and without land use from phase 6 of the Coupled Model Intercomparison Project to reconstruct an experiment representing the effects of land use on climate. To assess the signal of land use at spatially resolved continental and global scales, we aggregate all input data across reference regions and continents, respectively. At both scales, land use does not comprise a significantly detectable set of forcings for two of four Earth system models and their multimodel mean. Furthermore, using a principal component analysis, we find that t LU is mostly composed of the nonlocal effects of land use rather than its local effects. These findings show that, at scales relevant for climate attribution, uncertainties in Earth system model representations of land use are too high relative to the effects of internal variability to confidently assess land use.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigating the Global Biogeophysical Impact of Area and Mass Based Wood Harvest in a Vegetation Demography Model

Wood harvesting alters land surface properties and energy redistribution, but there is a lack of studies estimating these changes on a global scale. We coupled a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator, with the E3SM land model to perform offline model simulation to investigate the land biogeophysical responses, including canopy coverage, leaf area index, albedo, surface roughness length, and energy fluxes, to historical wood harvest on the global scale. In this study, we found 50% less harvested carbon (C) when choosing the area-based harvest rate as driving data that has not been spatially harmonized, compared to reharmonized mass-based harvesting. By considering the uncertainty from reconstruction of historical wood harvest time series and the choice of wood harvest approach in the model, continuous wood harvest (1850–2015) results in 5%–10% of canopy coverage loss, contributing 0.5%–1% increase of albedo over disturbed land, which is much stronger than a non-demographic land surface model. Changes in energy flux from the wood harvest are negligible (<1%), but the responses of land surface properties vary (up to 30%) due to differences in model structure between the single canopy, sun-shade leaf model and vegetation demographic model.

Shu, Shijie [Lawrence Berkeley National Laboratory↗
DOE OSTI2024

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A prescription for the asteroseismic surface correction

In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H–R diagram, we parameterize it as a simple function of surface gravity, effective temperature, and metallicity. We determine this function by fitting a wide range of stars. The absolute amount of the surface correction decreases with luminosity, but the ratio between it and νmax increases, suggesting the surface effect is more important for red giants than dwarfs. Applying the prescription can eliminate unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found a reduction of scatter in the model-derived ages for each star in the same cluster. As an important application, we provide a new revision for the Δν scaling relation that, for the first time, accounts for the surface correction. The values of the correction factor, f Δν , are up to 2 percent smaller than those determined without the surface effect considered, suggesting decreases of up to 4 percent in radii and up to 8 percent in masses when using the asteroseismic scaling relations. This revision brings the asteroseismic properties into an agreement with those determined from eclipsing binaries. The new correction factor and the stellar models with the corrected frequencies are available at https://www.github.com/parallelpro/surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2021

AI-enabled Dynamic Finish Machining Optimization for Sustained Surface Integrity

While machining processes are typically leveraged to establish geometric features, many functional characteristics of advanced materials are directly determined by their machining-induced surface integrity (SI). Current modeling approaches struggle to predict surface integrity, and typically neglect the effects of progressive tool-wear, resulting in inefficient ‘static’ process parameters. We present a novel integrated approach based on model-informed artificial intelligence (AI), which optimizes ‘dynamic’ process parameters in real-time. Here, by maximizing the useful life of a cutting tool over which a required set of SI parameters can be maintained, our paradigm will enable significantly more efficient processing of next-generation materials and components.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Crystallization Mechanisms of Poly(vinylidene Fluoride- co -chlorotrifluoroethylene): Nucleation Transitions, Growth Kinetics, and Microstructure Evolution

The crystallization kinetics of FK-800, a commercial semicrystalline copolymer of chlorotrifluoroethylene and vinylidene fluoride (poly(CTFE-co-VDF)), dictate its microstructure and functional performance. Crystallization occurs between the glass transition (T g ≈ 31 °C) and melting point (T m ≈ 110 °C), with grain boundaries playing a key role in applications such as memristors, where they regulate filament growth. This study integrates in situ atomic force microscopy (AFM), grazing-incidence wide-angle X-ray scattering (GIWAXS), and differential scanning calorimetry (DSC) to characterize nucleation, growth, and morphology across the full crystallization window. Hot-stage AFM reveals a transition from homogeneous to heterogeneous nucleation near 45 °C and a shift from reaction-limited to diffusion-limited growth above 60 °C. Kinetics are analyzed using Lauritzen–Hoffman, Turnbull–Fisher, and Avrami models to extract surface free energies and activation barriers. A multiscale modeling framework is developed in which nucleation densities and growth rates measured by AFM are used to reconstruct Avrami kinetics, which are then validated against coverage data, and extrapolated to predict bulk crystallization behavior observed by DSC. GIWAXS confirms a primarily flat-on chain orientation at all temperatures, linking orientation to domain morphology. This integrative approach quantitatively connects nanoscale crystallization dynamics with macroscopic phase evolution, establishing a quantitative framework for modeling crystallization across time, temperature, and length scales.

Crystallization↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗