Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface Deformation”

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

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗
DOE OSTI2026

Coming Closer to High Frequency Gravitational Wave Detection With MAGO

In the last years, low frequency gravitational waves (GWs) have been consistently measured by the LIGOVirgo collaboration, but little to no attention has been paid to higher frequencies GWs in the range of 10 kHz to 100 MHz, at which confirmation for current theories or even new physics could be hidden. The MAGO 2.0 project aims at filling this gap in the parameter space using superconducting radio-frequency (SRF) cavities. Exploiting the excellent Q-factors of these resonators, we plan to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. We present the results of the first cold tests ran at DESY and FNAL using the cavity prototype built 20 years ago at the end of the MAGO collaboration, characterizing the RF spectrum, Q-factor and surface resistance. Additionally we introduce the mechanical vibration spectrum characterization and the RF response of the cavity with the injection of a “fake GW” signal using piezoelectric actuators.

Marconato, Giovanni [Hamburg U.]↗
DOE OSTI2025

In Situ Formation of Ripplocations in Hybrid Organic–Inorganic MXenes

Abstract Inorganic–organic hybrid MXenes (h‐MXenes) are a family of 2D transition metal carbides and nitrides functionalized with alkylimido and alkylamido surface groups. Using cryogenic and room temperature scanning transmission electron microscopy (STEM) and electron energy‐loss spectroscopy (EELS), it is shown that ripplocations, a form of a fundamental defect in 2D and layered structures, are abundant in this family of materials. Furthermore, detailed studies of electron probe sample interactions, focusing on structural deformations caused by the electron beam are presented. The findings indicate that at cryogenic temperatures (≈100 K) and below a specific dose threshold, the structure of h‐MXenes remains largely intact. However, exceeding this threshold leads to electron beam‐induced deformation through ripplocations. Interestingly, the deformation behavior, required dose, and resultant structure are highly dependent on temperature. At 100 K, it is demonstrated that the electron beam can induce ripplocations in situ with a high degree of precision.

Chemistry↗
DOE OSTI2021

Out of the Lab: 3D Printing on Non-Ideal Surfaces

In this work, coaxial direct ink writing is utilized to print sensitive or encapsulated materials onto heterogeneous and rough substrates. By enclosing the core fluid within a rigid ceramic shell, continuity is maintained even when printing on substrates traditionally considered challenging. This work presents the development of a coaxial ceramic direct ink writing system and investigates co-flow dynamics through microfluidic principles. A specially designed coaxial nozzle enables the co-extrusion of an alumina shell, with indium-tin-oxide inks serving as the core material. This approach allows the fabrication of core-shell ceramic structures suitable for future high-temperature applications on rough surfaces. The colloidal inks are tailored to achieve the necessary rheological properties and sintering performance. Additionally, flow simulations, guided by microfluidic co-flow principles, are employed to define the processing parameters for coaxial printing, ensuring precise control over the core-shell architectures. Physical modeling is further applied to study core deformation and eccentricity. These simulations are experimentally validated, supporting the successful deposition of coaxial ceramic features on heterogeneous, high-temperature ceramic substrates.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Creep and failure at metal-oxide interfaces

Small-scale bicrystal creep experiments were performed on contacts formed via in situ high-temperature diffusion bonding of metal-oxide interfaces including Ag-ZrO 2 , Pd-ZrO 2 , Pt-ZrO 2 , and Ag-high entropy oxide. This work characterizes deformation and failure at metal-oxide interfaces during mechanical loading. Interfacial sliding can be activated easily, while tensile interfacial creep was not observed at any condition of stress or temperature measured. Plastic strain, instead, localizes within the metal under tensile loading. A variety of mechanisms for plastic strain occur in the metal including lattice dislocation-mediated plasticity, twinning, low-angle grain boundary formation, and low-angle grain boundary creep. Surface and low-angle grain boundary diffusion occur under conditions where no metal-oxide tensile creep is observed, highlighting the significant differences in their interfacial mechanical response. High-temperature interfacial failure occurs when the mean curvature at the contact neck is approximately zero and the applied stresses comparable to brittle fracture stresses. The brittle fracture stresses were measured to be σ ƒ = 180 ± 90 MPa at the Ag-ZrO 2 interface at 225 °C, σ ƒ = 460 ± 160 MPa at the Pd-ZrO 2 interface at 680 °C, and σ ƒ = 640 ± 440 MPa at the Pt-ZrO 2 interface at 1010 °C.

Creep↗
DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗
DOE OSTI2020

Atomic imprinting in the absence of an intrinsic length scale

Bulk metallic glasses (BMGs) have successfully been used to replicate molds that are structured at the nano- and even atomic scale through thermoplastic forming (TPF), an ability that was speculated to be rooted in the glass’ featureless atomic structure. These previous demonstrations of atomically precise imprinting, however, were performed under conditions where mold atomic feature dimensions coincided with the unit cell size of constituents in the BMG. In order to evaluate if accurate atomic-scale replication is possible in general, i.e., independent of the accidental presence of favorable constituent size/feature size relationships, we have used Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 to replicate three different crystalline facets of LaAlO 3 single crystals, each exposing distinct atomic step heights. We find that in all cases, the terraced surface termination can be copied with remarkable fidelity, corroborating that BMGs when thermoplastic formed are capable of adapting to any externally imposed confinement with sub-angstrom precision without being limited by factors related to the specifics of their internal structure. This unprecedented capability of quasi-limitless replication fidelity reveals that the deformation mechanism in the supercooled liquid state of BMGs is essentially homogeneous and suggests TPF of BMGs to be a versatile toolbox for atomic and precision nanoscale imprinting.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Windowed Fourier transform and cross-correlation algorithms for molecular tagging velocimetry

In this study, simulated and experimental molecular tagging velocimetry (MTV) images have been analyzed with a technique commonly used to process grid images on surfaces, the windowed Fourier transform with local spectrum analysis (WFT-LSA). A systematic synthetic image study of the modulation transfer function (MTF) and error tendencies of the WFT-LSA was performed and compared with a PIV-style cross-correlation algorithm to see if advanced strategies such as iterative image deformation can improve analysis of gridded images with high noise levels. Testing of single-pass algorithms showed that in typical MTV images, the WFT-LSA yields significantly lower bias errors than cross-correlation (CC) at displacements greater than 1 pixel but slightly higher random error at all displacements and image conditions. Analysis of the MTF shows that CC provided better resolution of spatial fluctuations than the WFT-LSA in many combinations of grid size and interrogation window. Tests of image deformation algorithms showed that the gap in performance between CC and WFT-LSA is maintained even as both methods improve. Additionally, WFT-LSA and CC methods are applied to real MTV experiments in high-speed gas jet flows. A preliminary analysis of phosphorescence lifetime provided by acetone vapor excited at 266 nm is used for assessing the required gas speeds for making MTV application feasible. The application of WFT-LSA to real MTV images demonstrates the ability of the algorithm to handle further real-world effects that could not be considered in the synthetic image analysis, like reduced signal-to-noise ratio and non-uniform intensity of the tagging grid across the image introduced by the actual laser beam energy distribution. With experimental images, CC is more accurate with shear flows but less robust to high noise levels than WFT-LSA, as predicted by the synthetic image analysis.

42 ENGINEERING↗
DOE OSTI2022

Low-Power, Flexible Sensor Arrays with Solderless Board-to-Board Connectors for Monitoring Soil Deformation and Temperature: Supporting Data

This dataset was used to assess the potential of a Soil Deformation and Temperature Monitoring System developed and presented in the article named "Low-Power, Flexible Sensor Arrays with Solderless Board-to-Board Connectors for Monitoring Soil Deformation and Temperature" and published in Sensors. There are 21 comma-delimited data files (.csv). 14 files contain current measurements performed in a lab setting, enabling the electrical evaluation of the sensor probe. These files are generated by the Keithley DMM6500 multimeter, and list the measured supply current (first column) as a function of time (3rd column). Another set of 6 files is also acquired in lab experiments, but contain soil temperature an deformation measurements, enabling an assessment of the developed device's accuracy. In these files, the first column contains the time (UTC), followed by the sensor's battery voltage and temperature and acceleration values (X, Y, Z) in subsequent columns. The measurements were acquired every 5 seconds. Another .csv file contains field data in a similar format (time, temperature, acceleration) collected at the Teller road (mile 27) site near Nome, Alaska with one probe. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗
DOE OSTI2023

Coupling turbulent flow with blade aeroelastics and control modules in large-eddy simulation of utility-scale wind turbines

We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. Furthermore, the findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.

42 ENGINEERING↗
DOE OSTI2024

Multi frame radiography of supersonic water jets interacting with a foil target

Pulsed-power-driven underwater electrical explosion of cylindrical or conical wire arrays produces supersonic water jets that emerge from a bath, propagating through the air above it. Interaction of these jets with solid targets may represent a new platform for attaining materials at high pressure (>1010 Pa) conditions in a university-scale laboratory. However, measurements of the internal structure of such jets and how they interact with targets are difficult optically due to large densities and density contrasts involved. We utilized multi-frame x-ray radiographic imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility to explore the water jet and its interaction with a 50 μm thick copper foil placed a few mm from the surface of water. The jet was generated with a ∼130 kA-amplitude current pulse of ∼450 ns rise time applied to a conical wire array. X-ray imaging revealed a droplet-type structure of the jet with an average density of <400 kg/m3 propagating with a velocity of ∼1400 m/s. Measurements of deformation and subsequent perforation of the target by the jet suggested pressures at the jet–target interface of ∼5 × 109 Pa. The results were compared to hydrodynamic simulations for better understanding of the jet parameters and their interaction with the foil target. These results can be used in future research to optimize the platform, and extend it to larger jet velocities in the case of higher driving currents supplied to the wire array.

Physics↗
DOE OSTI2024

Insights on seasonal solifluction processes in warm permafrost Arctic landscape using a dense monitoring approach across adjacent hillslopes

Solifluction processes in the Arctic are highly complex, introducing uncertainties in estimating current and future soil carbon storage and fluxes, and assessment of hillslope and infrastructure stability. This study aims to enhance our understanding of triggers and drivers of soil movement of permafrost-affected hillslopes in the Arctic. To achieve this, we established an extensive soil deformation and temperature sensor network, covering 48 locations across multiple hillslopes within a 1 km² watershed on the Seward Peninsula, AK. We report depth-resolved measurements down to 1.8 m depth for May to September 2022, a period conducive to soil movement due to deepening thaw layers and frequent rain events. Over this period, surface movements of up to 334 mm were recorded. In general, these movements occur close to the thawing front, and are initiated as thawing reaches depths of 0.4 to 0.75 m. The largest movements were observed at the top of the south-east facing slope, where soil temperatures are cold (mean annual soil temperatures averaging -1.13°C) and slopes are steeper than 15°. Our analysis highlights three primary factors influencing movements: slope angle, soil thermal conditions, and thaw depth. The latter two significantly impact the generation of pore water pressures at the thaw–freeze interface. Specifically, soil thermal conditions govern the liquid water content, while thaw depth influences both the height of the water column and, consequently, the pressure at the thawing front. These factors affect soil properties, such as cohesion and internal friction angle, which are crucial determinants of slope stability. This underscores the significance of a precise understanding of subsurface thermal conditions, including spatial and temporal variability in soil temperature and thaw depth, when assessing and predicting slope instabilities. Based on our observations, we developed a Factor of Safety proxy that consistently falls below the triggering threshold for all probes exhibiting displacements exceeding 50 mm. This study offers novel insights into patterns and triggers of hillslope movements in the Arctic and provides a venue to evaluate their impact on soil redistribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Nanomechanical squeezing: Toward quantum imaging with atomic force microscopy

Leveraging quantum effects in mechanical oscillators promises major advancements in nanometrology, including sensing and imaging. Microscopic cantilevers serve as versatile force sensors with broad applications in nanoscience, nanotechnology, and, increasingly, quantum sensing. Here, we theoretically model and investigate the quantum mechanics of cantilever probes in atomic force microscopy under conditions of reduced thermal noise. Specifically, our model captures the cantilever interaction with a surface via long-range attractive van der Waals forces and short-range repulsive-adhesive interactions described by the Derjaguin-Muller-Toporov model. We find that when the probe resides within the attractive interaction regime, a coherent state emerges, whereas, in the repulsive regime, a squeezed state forms for the cantilever's deformation state. Our calculations explain the functional role of interaction potentials in the quantum dynamics of macroscopic mechanical systems, which are proving useful in quantum sensing and information processing. The presented calculations can be extended to investigate other interaction forces relevant to atomic force microscopy.

Al-Sawai, Wael M. [Lamar Univ., Beaumont, TX (Unit↗
DOE OSTI2025

Geometric transport signatures of strained multi-Weyl semimetals

The minimal coupling of strain to Dirac and Weyl semimetals, and its modeling as a pseudogauge field has been extensively studied, resulting in several proposed topological transport signatures. In this work, we study the effects of strain on higher winding number Weyl semimetals and show that strain is not a pseudogauge field for any winding number larger than one. Here we focus on the double-Weyl semimetal as an illustrative example to show that the application of strain splits the higher winding number Weyl nodes and produces an anisotropic Fermi surface. Specifically, the Fermi surface of the double-Weyl semimetal acquires nematic order. By extending chiral kinetic theory for such nematic fields, we determine the effective gauge fields acting on the system and show how strain induces anisotropy and affects the geometry of the semiclassical phase space of the double-Weyl semimetal. Further, the strain-induced deformation of the Weyl nodes results in transport signatures related to the covariant coupling of the strain tensor to the geometric tensor associated with the Weyl nodes giving rise to strain-dependent dissipative corrections to the longitudinal as well as the Hall conductance. Thus, by extension, we show that in multi-Weyl semimetals, strain produces geometric signatures rather than topological signatures. Further, we highlight that the most general way to view strain is as a symmetry-breaking field rather than a pseudogauge field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-temperature nanoindentation creep studies on castable and sintered nanostructured low-activation ferritic-martensitic alloys

In this article, we present the creep characteristics of two reduced activation ferritic-martensitic steels of identical starting compositions formed by different fabrication routes: a nanostructured ferritic alloy commonly referred to as a castable nanostructured alloy (CNA) and a sintered nanostructured alloy (SNA) variant. Through a series of nanoindentation experiments spanning a temperature range of 25 °C to 650 °C, with a maximum load of 100 mN, we find creep behaviors in the cast and sintered materials to be remarkably similar. The creep stress exponent (n) for CNA and SNA were found to be in the range of 8–35 and the activation volume was ∼14–42b 3 , underscoring a dominance of dislocation-mediated mechanisms in both alloys. Notably, we observed a decline in the creep stress exponent with increasing temperature, attributable to the heightened influence of thermally activated dislocations. This phenomenon suggests a potential transition in the deformation mechanism towards a thermally activated dislocation climb process, significantly impacting the observed creep behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗