Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2024

Signatures of enhanced spin-triplet superconductivity induced by interfacial properties

While spin-triplet pairing remains elusive in nature, there is a growing effort to realize proximity-induced equal-spin triplet superconductivity in junctions with magnetic regions or an applied magnetic field and common s-wave superconductors. To enhance such spin-triplet contribution, it is expected that junctions with a weak interfacial barrier and strong spin-orbit coupling are desirable. Intuitively, a weak interfacial barrier enables a robust proximity-induced superconductivity and strong spin-orbit coupling promotes spin mixing, converting spin-singlet into spin-triplet superconductivity. In contrast, we reveal a nonmonotonic spin-triplet contribution with the strength of the interfacial barrier and spin-orbit coupling. This picture is established by considering different signatures in conductance and superconducting correlations, as well as by performing self-consistent calculations. As a result, we identify a strongly enhanced spin-triplet superconductivity, realized for an intermediate strength of interfacial barrier and spin-orbit coupling. In junctions with magnetic regions, an enhanced spin-triplet superconductivity leads to a large magnetoanisotropy of conductance and superconducting correlations. Furthermore, this picture of an enhanced spin-triplet superconductivity is consistent with experiments demonstrating a huge increase in the conductance magnetoanisotropy, which we predict can be further enhanced at a finite bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Coupled Multiphysics Modeling of Lithium-Ion Batteries for Automotive Crashworthiness Applications

Considerable advances have been made in battery safety models, but achieving predictive accuracy across a wide range of conditions continues to be challenging. Interactions between dynamically evolving mechanical, electrical, and thermal state variables make model prediction difficult during mechanical abuse scenarios. In this study, we develop a physics-based modeling approach that allows for choosing between different mechanical and electrochemical models depending on the required level of analysis. We demonstrate the use of this approach to connect cell-level abuse response to electrode-level and particle-level transport phenomena. A pseudo-two-dimensional model and simplified single-particle models are calibrated to electrical-thermal cycling data and applied to mechanically induced short-circuit scenarios to understand how the choice of electrochemical model affects the model prediction under abuse scenarios. These models are implemented using user-defined subroutines on ls-dyna finite element software and can be coupled with existing automotive crash safety models.

analysis and design of components↗
DOE OSTI2022

Stress Profile in Coating Layers of TRISO Fuel Particles in Contact with One Another

This work presents a discussion on a series of finite element analyses that assess stress evolution in the coating layers of tristructural isotropic (TRISO) particles in contact with each other while embedded in a matrix. The initial simulations were of applied uniaxial pressure versus matrix elastic modulus. These simulations predicted increasing stress in the silicon carbide coating layers of the TRISO particles with decreasing matrix elastic modulus. The second set of simulations focused on the effects of heating and cooling and the associated dimensional change on the state of stress in the coating layers. The general finding was that there was no significant difference below the coating layer’s deposition temperature. Although, above the deposition temperature, the contacting particles had higher stress compared with those that were separated. The third set of simulations focused on the effects of irradiation, specifically, creep, dimensional change, and swelling. An interface debonding model was introduced since these potential effects have a significant bearing on predicted stresses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Metal-insulator transition and magnetism of SU(3) fermions in the square lattice

We study the SU(3) symmetric Fermi-Hubbard model (FHM) in the square lattice at 1/3-filling using numerically exact determinant quantum Monte Carlo and numerical linked-cluster expansion techniques. We present the different regimes of the model in the T–U plane, which are characterized by local and short-range correlations, and capture signatures of the metal-insulator transition and magnetic crossovers. These signatures are detected as the temperature scales characterizing the rise of the compressibility, and an interaction-dependent change in the sign of the diagonal spin-spin correlation function. The analysis of the compressibility estimates the location of the metal-insulator quantum critical point at U c /t ~ 6, and provides a temperature scale for observing Mott physics at finite T. Furthermore, from the analysis of the spin-spin correlation function we observe that for U/t ≳ 6 and T ~ J = 4⁢t 2 /U there is a development of a short-range two-sublattice (2SL) antiferromagnetic structure, as well as an emerging three-sublattice (3SL) antiferromagnetic structure as the temperature is lowered below T/J ≲ 0.57. This crossover from 2SL to 3SL magnetic ordering agrees with Heisenberg limit predictions, and has observable effects on the density of on-site pairs. Finally, we describe how the features of the regimes in the T–U plane can be explored with alkaline-earth-like atoms in optical lattices with currently achieved experimental techniques and temperatures. Furthermore, the results discussed in this paper provide a starting point for the exploration of the SU(3) FHM upon doping.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Imaginary phonon modes and phonon-mediated superconductivity in Y 2 C 3

For Y 2 C 3 with a superconducting critical temperature ($T_c$) ~18 K, zone-center imaginary optical phonon modes have been found for the high-symmetry I-43d structure due to C dimer wobbling motion and electronic instability from a flat band near Fermi energy. After lattice distortion to the more stable lowest-symmetry P1 structure, these stabilized low-energy phonon modes with a mixed C and Y character carry a strong electronphonon coupling to give rise to the observed sizable $T_c$. Our work shows that compounds with the calculated dynamical instability should not be simply excluded in the high-throughput search for new phonon-mediated superconductors. Moreover, we have studied the phase stability of the I-43d structure by calculating the enthalpy of different structural motifs of binary compounds containing group-IV elements at the 2:3 composition and also exploring the energy landscapes via ab initio molecular dynamics near and out of the I-43d structure. Here, our results show that the I-43d-type structures with C dimers are preferred in the low to medium pressure range. Because of the wobbling motion of the C dimers, there are many local energy minimums with degenerated energies. Thus, the ensemble average of many I-43d-distorted structures with C dimer wobbling motion at finite temperature still gives an overall I-43d structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Half-Integer Conductance Plateau at the ν = 2 / 3 Fractional Quantum Hall State in a Quantum Point Contact

The v = 2/3 fractional quantum Hall state is the hole-conjugate state to the primary Laughlin v = 1/3 state. Here, we investigate transmission of edge states through quantum point contacts fabricated on a GaAs/AlGaAs heterostructure designed to have a sharp confining potential. When a small but finite bias is applied, we observe an intermediate conductance plateau with G = 0.5⁢(e 2 /ℎ). This plateau is observed in multiple QPCs, and persists over a significant range of magnetic field, gate voltage, and source-drain bias, making it a robust feature. Using a simple model that considers scattering and equilibration between counterflowing charged edge modes, we find this half-integer quantized plateau to be consistent with full reflection of an inner counterpropagating -1/3 edge mode while the outer integer mode is fully transmitted. In a QPC fabricated on a different heterostructure which has a softer confining potential, we instead observe an intermediate conductance plateau at G = (1/3)⁢(e 2 /ℎ). These results provide support for a model at v = 2/3 in which the edge transitions from a structure having an inner upstream -1/3 charge mode and outer downstream integer mode to a structure with two downstream 1/3 charge modes when the confining potential is tuned from sharp to soft and disorder prevails.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Characterization of latent image of electron beam resist via critical-dimension resonant soft X-ray scattering

Metrology plays a crucial role in semiconductor manufacturing by providing accurate and precise measurement and characterization of critical parameters. With the development of high-resolution extreme ultraviolet lithography (EUVL) processes, critical dimensions are shrinking to sub-10 nm. Resist materials encounter the challenge of providing heightened sensitivity and a handle on exacerbating stochastic variations. A comprehensive understanding of the chemical profile of the latent image is pivotal for mitigating stochastic effects and optimizing pattern quality. However, the subtle differences in chemistry between the exposed and unexposed regions of the resists make it extremely challenging to characterize the latent images with sub-nanometer precision. Here, we develop the metrology with critical-dimension resonant soft X-ray scattering (CD-RSoXS) to probe the chemical profiles of latent images stored in resist after exposure. The combination of absorption spectroscopy and enhanced scattering contrast makes it possible to characterize the subtle structural and chemical variations in the latent image. Moreover, the results of the measurements are compared with the simulations with a finite element method–based Maxwell solver to extract a detailed profile of the latent and developed images. We demonstrate that the CD-RSoXS technique can provide valuable insights into the high spatial resolution and local chemical sensitivity simultaneously, which is crucial to understanding the resolution limits and stochastic effects in EUVL processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological Hall effect and magnetic states in the Nowotny chimney ladder compound Cr 11 Ge 19

In this report, we have investigated the magnetic and charge transport properties of single crystals of the Nowotny chimney ladder compound Cr 11 Ge 19 and mapped out a comprehensive phase diagram reflecting the complicated interplay between the Dzyaloshinskii-Moriya (DM) interaction, the dipolar interaction, and the magnetic anisotropy. We have identified a set of interesting magnetic phases and attributed a finite topological Hall effect to the recently discovered biskyrmion phase. These data also suggest the existence of an antiskyrmion state at finite fields for temperatures just below the magnetic ordering temperature T c as indicated by a distinct change in sign of the topological Hall effect. Above T c , we discovered a region of enhanced magnetic response corresponding to a disordered phase likely existing near the ferromagnetic critical point under small magnetic fields. Strong spin chirality fluctuations are demonstrated by the large value of the topological Hall resistivity persisting up to 1 T, which is most likely due to the existence of the DM interaction. We argue that changes to the topological Hall effect correspond to different topological spin textures that are controlled by magnetic dipolar and DM interactions that vary in importance with temperature

36 MATERIALS SCIENCE↗
DOE OSTI2020

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2026

Phase change material-to-refrigerant heat exchangers: Experimental validation and uneven melting analysis

Thermal energy storage (TES) using phase change materials (PCMs) enables load shifting and reduces mismatches between the building's thermal demand and the heat pump (HP) system's thermal capacity. While most TES studies have focused on single-phase heat transfer fluids, PCM-to-refrigerant heat exchangers (PRHX), in which the refrigerant undergoes liquid-vapor phase change, remain largely underexplored. Here, to address this gap, this study systematically investigates a shell-and-tube PRHX operating as the condenser in an HP through experiments and simulations. A two-dimensional enthalpy-based finite-volume PRHX model was developed and validated against experimental data, including refrigerant and PCM temperature profiles, PRHX capacity, and system-level performance. Results show three characteristic stages: Stage I, the condenser temperature rapidly increases with PCM in the solid state absorbing sensible heat; Stage II, stable operation during PCM phase change; and Stage III, performance decline once part of the PCM becomes fully melted. Heat transfer analysis revealed that uneven PCM melting along the condenser length was driven primarily by variations in the refrigerant-side heat transfer coefficient, rather than local approach temperature differences. Further study demonstrated that minimizing refrigerant outlet subcooling led to more uniform PCM melting and extended Stages I and II by up to 89%. This paper provides a validated PRHX model, clarifies the mechanisms of uneven PCM melting, and highlights subcooling control as an effective strategy to improve TES performance in HP systems without requiring secondary loops.

Graphite matrix↗
DOE OSTI2020

Simulating transient heat transfer in graphene at finite Knudsen number via the Boltzmann transport equation and molecular dynamics

The phonon Boltzmann transport equation (BTE) with the relaxation time approximation (RTA) scattering model is used to calculate transient temperature profiles in graphene, and the results are compared to analogous molecular dynamics (MD) simulations. For the BTE calculations, the phonon dispersion relation and frequency-dependent scattering rates are obtained from a combination of MD data and semi-empirical power-law expressions for the normal and Umklapp phonon lifetimes. The dimensions and initial temperature conditions of graphene are varied to study the size and temperature dependence of thermal transport physics at the mesoscopic scale. Good quantitative agreement to within 5% is found between the BTE and MD results, over a wide range of temperatures and lengthscales of the temperature variation in the graphene sheet. Small differences are attributed to the inaccuracy of the RTA as applied to graphene, and to neglecting four-phonon scattering in the BTE simulations. The present results may further understanding in applications such as the transient heating of nanoelectronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A pseudo-kinetic model to simulate phase changes in gas hydrate bearing sediments

Herein, modeling of the phase transitions anticipated in gas hydrate bearing sediments (GHBS) is critical for a proper understanding of time-dependent changes in states and volumes (e.g. the production of methane from this type of soils). We propose a new pseudo-kinetic approach to simulate the typical phase changes anticipated in GHBS, using published experimental results involving gas hydrate dissociation that are the basis of a widely used kinetic model. The proposed pseudo-kinetic model is formulated in the pressure-temperature (P-T) plane and assumes a rate of gas hydrate dissociation (or formation) proportional to the distance between the current state and the phase boundary. The model consists of only one parameter and is simple to implement in numerical simulators. A similar concept is used to model ice formation/thawing phenomena, but based on the ice/liquid-water phase boundary. We implemented the pseudo-kinetic model in a fully coupled thermo-hydro-chemo-mechanical (THCM) finite element code and validated it against experimental results performed on the dissociation of synthetic gas hydrate. We also evaluated the pseudo-kinetic model using synthetic cases covering several scenarios associated with gas hydrate formation/dissociation and ice formation/thawing. The model successfully reproduced the gas production test from a natural GHBS core from Korea (scaled gas venting experiment over 14 h), and also the formation of gas hydrate and ice in permafrost in Alaska (over 2 × 10 6 years). -The analyses show the versatility of the proposed pseudo-kinetic approach by applying it to model the different types of phase transitions typically encounter in GHBS. The simple formulation, easy implementation in numerical simulator, and reduced number of parameters (only one per phase change) make this model an attractive option for simulating phase transformations in problems involving GHBS.

58 GEOSCIENCES↗
DOE OSTI2025

Dynamical coupled-channel models for hadron dynamics

Dynamical coupled-channel (DCC) approaches parametrize the interactions and dynamics of two and more hadrons and their response to different electroweak probes. The inclusion of unitarity, three-body channels, and other properties from scattering theory allows for a reliable extraction of resonance spectra and their properties from data. Here, we review the formalism and application of the ANL-Osaka, the Juelich-Bonn-Washington, and other DCC approaches in the context of light baryon resonances from meson, (virtual) photon, and neutrino-induced reactions, as well as production reactions, strange baryons, light mesons, heavy meson systems, exotics, and baryon-baryon interactions. Finally, we also provide a connection of the formalism to study finite-volume spectra obtained in Lattice QCD, and review applications involving modern statistical and machine learning tools.

Amplitude analysis↗
DOE OSTI2025

Geometry and size dependent microstructure and crack formation in Rene 41 superalloy fabricated by laser powder bed fusion

This study provides a systematic investigation into the size and geometry-dependent microstructural evolution and cracking susceptibility of LPBF-fabricated Rene 41. By coupling experimental microstructural analysis with thermal modeling, this research uniquely identifies the relationship between geometry-specific thermal histories, carbide coarsening, and liquation cracking. Four different geometries with varying thickness were fabricated with the identical process parameters. It was found that the grain size and morphology are not affected by the part size. However, the thinner parts exhibited coarser sub-grain structures compared to the thicker ones. The crack formation was observed for the parts with cross-sections smaller than 1 mm, whereas thicker parts had high density without any defects. Further, the cracks were observed in the interdendritic regions, suggesting that liquation cracking was the active micro-crack formation mechanism. The detailed microstructural analysis combined with a thermal finite element analysis proved that the heat extraction efficiency was lower for thinner parts causing a lower cooling rate and coarser carbides, making them more susceptible to constitutional liquation. Microhardness measurements were conducted for each geometry and correlated with the observed microstructural variations. The findings highlight the critical need for geometry-specific optimization of LPBF process parameters to mitigate cracking and achieve microstructural uniformity, offering valuable insights into the fabrication of complex, high-performance aerospace components.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Increased Defect Resistance and Ordering in MnBi 2 (Se 1– x Te x ) 4 via Accurate Diffusion Monte Carlo

Stabilizing materials and controlling defect formation remain key challenges in materials science, particularly for theory, where small energy differences must be resolved for accurate predictions. Here, we applied state-of-the-art theoretical methods to topological materials, focusing on MnBi 2 Te 4 (MBT), which is a promising intrinsic magnetic topological insulator. Antisite defects in MBT alter its electronic structure and magnetism, degrading topological properties and causing experimental inconsistencies. Using diffusion Monte Carlo and density functional theory, we investigated the thermodynamic stability and defect formation in MBT, MnBi 2 Se 4 (MBS), and MnBi 2 (Se 1–x Te x ) 4 . We found that MnBi 2 Se 2 Te 2 can be stable at finite temperatures, with higher defect formation energies due to stronger Mn–Se bonding and reduced internal strain. Se preferentially substitutes Te near Mn instead of Te in the outer layer, encouraging long-range ordering when incorporated. For MnBi 2 (Se 1–x Te x ) 4 , cluster expansion phase diagrams revealed solid solution behavior when x <0.5 and phase separation for larger x. MBT and MBS are topological insulators; therefore, the MnBi 2 (Se 1–x Te x ) 4 family could offer tunable topological behavior and improved stability.

MnBi2Te4↗
DOE OSTI2024

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗