SEARCH · Search NASA

Results for “Deformation Gradient”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1970

Discrete control of linear distributed systems with application to the deformable primary mirror of a large orbiting telescope

The results are presented of a detailed study of the discrete control of linear distributed systems with specific application to the design of a practical controller for a plant representative of a telescope primary mirror for an orbiting astronomical observatory. The problem of controlling the distributed plant is treated by employing modal techniques to represent variations in the optical figure. Distortion of the mirror surface, which arises primarily from thermal gradients, is countered by actuators working against a backing structure to apply a corrective force distribution to the controlled surface. Each displacement actuator is in series with a spring attached to the mirror by means of a pad intentionally introduced to restrict the excitation of high-order modes. Control is exerted over a finite number of the most significant modes.

Creedon, J. F.↗

NASA NTRS1993

Horizontal stresses induced by vertical processes in planetary lithospheres

Understanding the state of stress in the elastic lithosphere is of fundamental importance for planetary geophysics, as it is the link between the observed geologic structures on the surface and the processes which form and modify these structures. As such, it can provide valuable constraints for the difficult problem of determining interior structure and processes. On the Earth, most large scale, organized deformation can be related to lateral tectonics associated with plate dynamics; however, the tectonics on many extraterrestrial bodies (such as the Moon, Mars, and most of the outer-planet satellites) appears to be primarily vertical in nature, and the horizontal stresses induced by vertical motions and loads are expected to dominate the deformation of their lithospheres. The largest stress contributions from vertical loading come from the flexure of the lithosphere, which induces both bending moments and membrane stresses. We are concerned here only with nonflexural changes in the state of stress induced by processes such as sedimentary and volcanic deposition, erosional denudation, and changes in the thermal gradient that induce uplift or subsidence. This analysis is important both for evaluating stresses for specific regions in which the vertical stress history can be estimated, as well as for applying the proper loading conditions to global stress models. It is also of interest for providing a reference state of stress for interpreting stress measurements in the crust of the Earth.

Banerdt, W. B.↗

NASA NTRS1994

Thermocapillary flow with evaporation and condensation and its effect on liquid retention in low-G fluid acquisition devices

The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.

Schmidt, George R.↗

NASA NTRS1989

Jupiter's Great Red Spot as a shallow water system

Voyager cloud-top velocity data for Jupiter's Great Red Spot (GRS) is used to derive the bottom topography up to a constant that depends on the unknown radius of deformation. The bottom topography is inferred from the Bernoulli streamfunction, kinetic energy per unit mass, and absolute vorticity values derived from the velocity data. The results are used to calculate potential vorticity versus latitude far from the vortex. It is found that the deep atmosphere is in differential motion and that the far-field potential vorticity gradient changes sign at several latitudes. Numerical experiments are conducted to study the time-dependent behavior of the shallow water analog of Jupiter's analog.

Dowling, Timothy E.↗

NASA NTRS1990

Thermomechanical deformation testing and modeling in the presence of metallurgical instabilities

A number of viscoplastic constitutive models were developed to describe deformation behavior under complex combinations of thermal and mechanical loading. Questions remain, however, regarding the validity of procedures used to characterize these models for specific structural alloys. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This experimental study is aimed at determining whether viscoplastic constitutive theories characterized using an isothermal data base can adequately model material response under the complex thermomechanical loading conditions typical of power generation service. The approach adopted was to conduct a series of carefully controlled thermomechanical experiments on a nickel-based superalloy, Hastelloy Alloy X. Previous investigations had shown that this material experiences metallurgical instabilities leading to complex hardening behavior, termed dynamic strain aging. Investigating this phenomenon under full thermomechanical conditions leads to a number of challenging experimental difficulties which up to the present work were unresolved. To correct this situation, a number of advances were made in thermomechanical testing techniques. Advanced methods for dynamic temperature gradient control, phasing control and thermal strain compensation were developed and incorporated into real time test control software. These advances allowed the thermomechanical data to be analyzed with minimal experimental uncertainty. The thermomechanical results were evaluated on both a phenomenological and microstructural basis. Phenomenological results revealed that the thermomechanical hardening trends were not bounded by those displayed under isothermal conditions. For the case of Hastelloy Alloy X (and similar dynamic strain aging materials), this strongly suggests that some form of thermomechanical testing is necessary when characterizing a thermoviscoplastic deformation model. Transmission electron microscopy was used to study the microstructural physics, and analyze the unique phenomenological behavior.

Castelli, Michael G.↗

NASA NTRS1989

Large space antenna concepts for ESGP

It is appropriate to note that 1988 marks the 100th anniversary of the birth of the reflector antenna. It was in 1888 that Heinrich Hertz constructed the first one, a parabolic cylinder made of sheet zinc bent to shape and supported by a wooden frame. Hertz demonstrated the existence of the electromagnetic waves that had been predicted theoretically by James Clerk Maxwell some 22 years earlier. In the 100 years since Hertz's pioneering work the field of electromagnetics has grown explosively: one of the technologies is that of remote sensing of planet Earth by means of electromagnetic waves, using both passive and active sensors located on an Earth Science Geostationary Platform (ESEP). For these purposes some exquisitely sensitive instruments were developed, capable of reaching to the fringes of the known universe, and relying on large reflector antennas to collect the minute signals and direct them to appropriate receiving devices. These antennas are electrically large, with diameters of 3000 to 10,000 wavelengths and with gains approaching 80 to 90 dB. Some of the reflector antennas proposed for ESGP are also electrically large. For example, at 220 GHz a 4-meter reflector is nearly 3000 wavelengths in diameter, and is electrically quite comparable with a number of the millimeter wave radiotelescopes that are being built around the world. Its surface must meet stringent requirements on rms smoothness, and ability to resist deformation. Here, however, the environmental forces at work are different. There are no varying forces due to wind and gravity, but inertial forces due to mechanical scanning must be reckoned with. With this form of beam scanning, minimizing momentum transfer to the space platform is a problem that demands an answer. Finally, reflector surface distortion due to thermal gradients caused by the solar flux probably represents the most challenging problem to be solved if these Large Space Antennas are to achieve the gain and resolution required of them.

Love, Allan W.↗

NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗

NASA NTRS2017

Viscous Aerodynamic Shape Optimization with Installed Propulsion Effects

Aerodynamic shape optimization is demonstrated to tailor the under-track pressure signature of a conceptual low-boom supersonic aircraft. Primarily, the optimization reduces nearfield pressure waveforms induced by propulsion integration effects. For computational efficiency, gradient-based optimization is used and coupled to the discrete adjoint formulation of the Reynolds-averaged Navier Stokes equations. The engine outer nacelle, nozzle, and vertical tail fairing are axi-symmetrically parameterized, while the horizontal tail is shaped using a wing-based parameterization. Overall, 48 design variables are coupled to the geometry and used to deform the outer mold line. During the design process, an inequality drag constraint is enforced to avoid major compromise in aerodynamic performance. Linear elastic mesh morphing is used to deform volume grids between design iterations. The optimization is performed at Mach 1.6 cruise, assuming standard day altitude conditions at 51,707-ft. To reduce uncertainty, a coupled thermodynamic engine cycle model is employed that captures installed inlet performance effects on engine operation.

computational fluid dynamics↗

NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗

NASA NTRS1976

Higher-order effects of initial deformation on the vibrations of crystal plates

A system of plate equations for the thickness-shear and flexural vibrations superposed on large initial deflection due to bending is derived; in the stress-strain relations, the terms associated with the fourth-order elastic stiffness coefficients are retained. An explicit formula for the change in the fundamental cut-off thickness shear frequency is obtained, and the effects of the terms associated with the fourth-order constants appear to be significant for large gradients of the rotation angles.

Markenscoff, X.↗

DOE OSTI2022

Micropillar Compression of Additively Manufactured 316L Stainless Steels after 2 MeV Proton Irradiation: A Comparison Study between Planar and Cross-Sectional Micropillars

A micropillar compression study with two different techniques was performed on proton-irradiated additively manufactured (AM) 316L stainless steels. The sample was irradiated at 360 °C using 2 MeV protons to 1.8 average displacement per atom (dpa) in the near-surface region. A comparison study with mechanical test and microstructure characterization was made between planar and cross-sectional pillars prepared from the irradiated surface. While a 2 MeV proton irradiation creates a relatively flat damage zone up to 12 µm, the dpa gradient by a factor of 2 leads to significant dpa uncertainty along the pillar height direction for the conventional planar technique. Cross-sectional pillars can significantly reduce such dpa uncertainty. From one single sample, three cross-sectional pillars were able to show dpa-dependent hardening. Furthermore, post-compression transmission electron microscopy allows the determination of the deformation mechanism of individual micropillars. Cross-sectional micropillar compression can be used to study radiation-induced mechanical property changes with better resolution and less data fluctuation.

transmission electron microscopy↗

NASA NTRS1999

A Novel Shape Parameterization Approach

This paper presents a novel parameterization approach for complex shapes suitable for a multidisciplinary design optimization application. The approach consists of two basic concepts: (1) parameterizing the shape perturbations rather than the geometry itself and (2) performing the shape deformation by means of the soft objects animation algorithms used in computer graphics. Because the formulation presented in this paper is independent of grid topology, we can treat computational fluid dynamics and finite element grids in a similar manner. The proposed approach is simple, compact, and efficient. Also, the analytical sensitivity derivatives are easily computed for use in a gradient-based optimization. This algorithm is suitable for low-fidelity (e.g., linear aerodynamics and equivalent laminated plate structures) and high-fidelity analysis tools (e.g., nonlinear computational fluid dynamics and detailed finite element modeling). This paper contains the implementation details of parameterizing for planform, twist, dihedral, thickness, and camber. The results are presented for a multidisciplinary design optimization application consisting of nonlinear computational fluid dynamics, detailed computational structural mechanics, performance, and a simple propulsion module.

Samareh, Jamshid A.↗

NASA NTRS2000

Multidisciplinary Aerodynamic-Structural Shape Optimization Using Deformation (MASSOUD)

This paper presents a multidisciplinary shape parameterization approach. The approach consists of two basic concepts: (1) parameterizing the shape perturbations rather than the geometry itself and (2) performing the shape deformation by means of the soft object animation algorithms used in computer graphics. Because the formulation presented in this paper is independent of grid topology, we can treat computational fluid dynamics and finite element grids in a similar manner. The proposed approach is simple, compact, and efficient. Also, the analytical sensitivity derivatives are easily computed for use in a gradient-based optimization. This algorithm is suitable for low-fidelity (e.g., linear aerodynamics and equivalent laminated plate structures) and high-fidelity (e.g., nonlinear computational fluid dynamics and detailed finite element modeling analysis tools. This paper contains the implementation details of parameterizing for planform, twist, dihedral, thickness, camber, and free-form surface. Results are presented for a multidisciplinary design optimization application consisting of nonlinear computational fluid dynamics, detailed computational structural mechanics, and a simple performance module.

Samareh, Jamshid A.↗

NASA NTRS2000

Multidisciplinary Aerodynamic-Structural Shape Optimization Using Deformation (MASSOUD)

This paper presents a multidisciplinary shape parameterization approach. The approach consists of two basic concepts: (1) parameterizing the shape perturbations rather than the geometry itself and (2) performing the shape deformation by means of the soft object animation algorithms used in computer graphics. Because the formulation presented in this paper is independent of grid topology, we can treat computational fluid dynamics and finite element grids in the same manner. The proposed approach is simple, compact, and efficient. Also, the analytical sensitivity derivatives are easily computed for use in a gradient-based optimization. This algorithm is suitable for low-fidelity (e.g., linear aerodynamics and equivalent laminate plate structures) and high-fidelity (e.g., nonlinear computational fluid dynamics and detailed finite element modeling) analysis tools. This paper contains the implementation details of parameterizing for planform, twist, dihedral, thickness, camber, and free-form surface. Results are presented for a multidisciplinary application consisting of nonlinear computational fluid dynamics, detailed computational structural mechanics, and a simple performance module.

Samareh, Jamshid A.↗

NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

NASA NTRS1997

Ice in Channels and Ice-Rock Mixtures in Valleys on Mars: Did They Slide on Deformable Rubble Like Antarctic Ice Streams?

Recent studies of ice streams in Antarctica reveal a mechanism of basal motion that may apply to channels and valleys on Mars. The mechanism is sliding of the ice on deformable water-saturated till under high pore pressures. It has been suggested by Lucchitta that ice was present in outflow channels on Mars and gave them their distinctive morphology. This ice may have slid like Antarctic ice streams but on rubbly weathering products rather than till. However, to generate water under high pore pressures, elevated heatflow is needed to melt the base of the ice. Either volcanism or higher heatflow more than 2 b.y. ago could have raised the basal temperature. Regarding valley networks, higher heatflow 3 b.y. ago could have allowed sliding of ice-saturated overburden at a few hundred meters depth. If the original, pristine valleys were somewhat deeper than they are now, they could have formed by the same mechanism. Recent sounding of the seafloor in front of the Ross Ice Shelf in Antarctica reveals large persistent patterns of longitudinal megaflutes and drumlinoid forms, which bear remarkable resemblance to longitudinal grooves and highly elongated streamlined islands found on the floors of martian outflow channels. The flutes are interpreted to have formed at the base of ice streams during the last glacial advance. Additional similarities of Antarctic ice streams with martian outflow channels are apparent. Antarctic ice streams are 30 to 80 km wide and hundreds of kilometers long. Martian outflow channels have similar dimensions. Ice stream beds are below sea level. Carr determined that most common floor elevations of martian outflow channels lie below martian datum, which may have been close to or below past martian sea levels. The Antarctic ice stream bed gradient is flat and locally may go uphill, and surface slopes are exceptionally. Martian channels also have floor gradients that are shallow or go uphill locally and have low surface gradients. The depth to the bed in ice streams is 1 to 1.5 km. At bankful stage, the depth of the fluid in outflow channels was 1 to 2 km, according to the height of bordering scarps. The similarity between Antarctic ice streams and martian outflow channels suggests that ice may have flowed through and shaped the outflow channels, and that perhaps the mechanism of motion of Antarctic ice streams also operated in outflow channels. In addition, sliding on deformable rubble may explain the formation of small valley networks. The large Siple Coast Antarctic ice streams are thought to slide over longitudinally grooved, deforming till, where much of the movement is within the till. The till is saturated with water at high pore pressures that nearly supports all of the weight of the ice. The small differential between overburden pressure and pore pressure at the bed is more important than the volume of water, but water needs to be supplied to the till interface. For pore pressures to remain high, the ice streams have to act as a seal that blocks the flow of water through them, and the rock underneath has to be of low permeability to prevent the water from draining away.

Lucchitta, B. K.↗

DOE OSTI2021

Development and Implementation of Digital Image Correlation for Out-of-Cell Burst Testing

This report details the progress made towards utilizing in-situ methods to observe and quantify unirradiated Zircaloy-4 (Zr4) fuel cladding deformation during out-of-cell design basis accident conditions, specifically during loss-of-coolant accident (LOCA) burst testing in the Severe Accident Test Station at Oak Ridge National Laboratory. Digital image correlation (DIC) was implemented to calculate cladding strain during burst and infrared (IR) thermography was used to attempt a correlation of strain with temperature and to investigate thermal gradients. For both techniques, early experimentation revealed experimental modifications were necessary for implementation. For DIC, this manifested as shortened cladding lengths which could be centered in front of the viewing chamber bored through the IR furnace side, and for IR thermography this meant conducting burst tests in air without a reaction tube and the incorporation of a SiC shell to eliminate direct IR reflections from the tungsten lamps. In the following, the technical issues and subsequent origin of these modifications are discussed, laying the framework for a proof-of-concept DIC/IR characterization of a burst test in the final section.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗