Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential strain”

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 145 records · Page 8

Coupled momentum balance and phase-field solver with fenicsx module

Code solves momentum balance and phase-field equations simultaneously. The differential equations are solved on a discretized domain with appropriate boundary and initial conditions using finite element method. Primary purpose of the code is to simulate brittle fracture under dynamic loading. Constitutive equations are that of linear elasticity with degradation of stress due to fracture. Small strain formulation is used.

Zecevic, Milovan↗

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes

Introduction: Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi. Objectives: This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution. Methods: We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution. Results: We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation. Conclusion: Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

59 BASIC BIOLOGICAL SCIENCES↗

Low-abundance populations distinguish microbiome performance in plant cell wall deconstruction

Abstract Background Plant cell walls are interwoven structures recalcitrant to degradation. Native and adapted microbiomes can be particularly effective at plant cell wall deconstruction. Although most understanding of biological cell wall deconstruction has been obtained from isolates, cultivated microbiomes that break down cell walls have emerged as new sources for biotechnologically relevant microbes and enzymes. These microbiomes provide a unique resource to identify key interacting functional microbial groups and to guide the design of specialized synthetic microbial communities. Results To establish a system assessing comparative microbiome performance, parallel microbiomes were cultivated on sorghum ( Sorghum bicolor L. Moench) from compost inocula. Biomass loss and biochemical assays indicated that these microbiomes diverged in their ability to deconstruct biomass. Network reconstructions from gene expression dynamics identified key groups and potential interactions within the adapted sorghum-degrading communities, including Actinotalea , Filomicrobium , and Gemmatimonadetes populations. Functional analysis demonstrated that the microbiomes proceeded through successive stages that are linked to enzymes that deconstruct plant cell wall polymers. The combination of network and functional analysis highlighted the importance of cellulose-degrading Actinobacteria in differentiating the performance of these microbiomes. Conclusions The two-tier cultivation of compost-derived microbiomes on sorghum led to the establishment of microbiomes for which community structure and performance could be assessed. The work reinforces the observation that subtle differences in community composition and the genomic content of strains may lead to significant differences in community performance.

59 BASIC BIOLOGICAL SCIENCES↗

Development of dual‐inducible duet‐expression vectors for tunable gene expression control and CRISPR interference‐based gene repression in Pseudomonas putida KT2440

Summary The development of P. putida as an industrial host requires a sophisticated molecular toolbox for strain improvement, including vectors for gene expression and repression. To augment existing expression plasmids for metabolic engineering, we developed a series of dual‐inducible duet‐expression vectors for P. putida KT2440. A number of inducible promoters (P lac , P tac , P tetR/tetA and P bad ) were used in different combinations to differentially regulate the expression of individual genes. Protein expression was evaluated by measuring the fluorescence of reporter proteins (GFP and RFP). Our experiments demonstrated the use of compatible plasmids, a useful approach to coexpress multiple genes in P. putida KT2440. These duet vectors were modified to generate a fully inducible CRISPR interference system using two catalytically inactive Cas9 variants from S. pasteurianus (dCas9) and S. pyogenes (spdCas9). The utility of developed CRISPRi system(s) was demonstrated by repressing the expression of nine conditionally essential genes, resulting in growth impairment and prolonged lag phase for P. putida KT2440 growth on glucose. Furthermore, the system was shown to be tightly regulated, tunable and to provide a simple way to identify essential genes with an observable phenotype.

59 BASIC BIOLOGICAL SCIENCES↗

Probing thermodynamics of radiogenic helium and defects in $δ$-plutonium alloys and interactions with adsorbed environmental gases

Differential scanning calorimetry coupled with simultaneous evolved gas analysis (DSC-EGA) on aged δ-Pu samples shows that most radiogenic helium remains trapped within the Pu matrix at temperatures very close to, or slightly above, the melting temperature. Our results indicate that helium release from 50-year-old δ-Pu occurs as a burst just below the melting temperature (>0.994 T m ), with subsequent pressure oscillations as temperature increases. Subordinate quantities of H 2 were also released along with helium. Here, the helium emission tails off and ceases above ~720 – 750°C. In a δ-Pu alloy aged 6 years, the initial helium burst occurs slightly above melting (~1.015 – 1.042 T m ), with a discrete, larger helium spike occurring between 670 and 686°C. The proximity of helium release to the liquidus transition presented challenges in the deconvolution of overlapping process enthalpies, the liquidus endotherm, and the exotherm resulting from bubble collapse, annealing and gas expulsion. Helium’s strong affinity for vacancy binding in a 2He-vac configuration is predicted by Density Functional Theory (DFT) modeling. The measured stored energy associated with the He release events in a 50-year-old δ alloy is on the order of ~ 10–11 J/g, which is significantly higher than stored energies measured in the sub-solidus regimes (~2 J/g) that are related to the solid-state annealing of processing- and radiation-induced defects. This implies that aged δ Pu alloys have a remarkable resilience to accommodate the lattice strain produced by the internal pressure of the helium bubbles and provides further insight into the thermodynamic behavior of aged δ Pu.

36 MATERIALS SCIENCE↗

Large electromechanical strain and unconventional domain switching near phase convergence in a Pb-free ferroelectric

In many ferroelectrics, large electromechanical strains are observed near regions of composition- or temperature- driven phase coexistence. Phenomenologically, this is attributed to easy re-orientation of the polarization vector and/or phase transition, although their effects are highly convoluted and difficult to distinguish experimentally. Here, we used synchrotron X-ray scattering and digital image correlation to differentiate between the microscopic mechanisms leading to large electrostrains in an exemplary Pb-free piezoceramic Sn-doped barium calcium zirconate titanate. Large electrostrains of ~0.2% measured at room-temperature are attributed to an unconventional effect, wherein polarization switching is aided by a reversible phase transition near the tetragonal-orthorhombic phase boundary. Additionally, electrostrains of ~0.1% or more could be maintained from room temperature to 140°C due to a succession of different microscopic mechanisms. In situ X-ray diffraction elucidates that while 90° domain reorientation is pertinent below the Curie temperature ( T C ), isotropic distortion of polar clusters is the dominant mechanism above T C .

36 MATERIALS SCIENCE↗

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗

Transmission of Human Influenza A Virus in Pigs Selects for Adaptive Mutations on the HA Gene

Influenza A viruses (FLUAV) cause respiratory diseases in many host species, including humans and pigs. The spillover of FLUAV between swine and humans has been a concern for both public health and the swine industry. With the emergence of the triple reassortant internal gene (TRIG) constellation, establishment of human-origin FLUAVs in pigs has become more common, leading to increased viral diversity. However, little is known about the adaptation processes that are needed for a human-origin FLUAV to transmit and become established in pigs. Here, we generated a reassortant FLUAV (VIC11pTRIG) containing surface gene segments from a human FLUAV strain and internal gene segments from the 2009 pandemic and TRIG FLUAV lineages and demonstrated that it can replicate and transmit in pigs. Sequencing and variant analysis identified three mutants that emerged during replication in pigs, which were mapped near the receptor binding site of the hemagglutinin (HA). The variants replicated more efficiently in differentiated swine tracheal cells compared to the virus containing the wildtype human-origin HA, and one of them was present in all contact pigs. These results show that variants are selected quickly after replication of human-origin HA in pigs, leading to improved fitness in the swine host, likely contributing to transmission.

60 APPLIED LIFE SCIENCES↗

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗

A newly characterized malaria antigen on erythrocyte and merozoite surfaces induces parasite inhibitory antibodies

We previously identified a Plasmodium falciparum (Pf) protein of unknown function encoded by a single-copy gene, PF3D7_1134300, as a target of antibodies in plasma of Tanzanian children in a whole-proteome differential screen. Here we characterize this protein as a blood-stage antigen that localizes to the surface membranes of both parasitized erythrocytes and merozoites, hence its designation as Pf erythrocyte membrane and merozoite antigen 1 (PfEMMA1). Mouse anti-PfEMMA1 antisera and affinity-purified human anti-PfEMMA1 antibodies inhibited growth of P. falciparum strains by up to 68% in growth inhibition assays. Following challenge with uniformly fatal Plasmodium berghei (Pb) ANKA, up to 40% of mice immunized with recombinant PbEMMA1 self-cured, and median survival of lethally infected mice was up to 2.6-fold longer than controls (21 vs. 8 d, P = 0.005). Furthermore, high levels of naturally acquired human anti-PfEMMA1 antibodies were associated with a 46% decrease in parasitemia over 2.5 yr of follow-up of Tanzanian children. Together, these findings suggest that antibodies to PfEMMA1 mediate protection against malaria.

60 APPLIED LIFE SCIENCES↗

The Influence of Bacteria on the Growth, Lipid Production, and Extracellular Metabolite Accumulation by Phaeodactylum tricornutum (Bacillariophyceae)

Here, to examine the impact of heterotrophic bacteria on microalgal physiology, we co-cultured the diatom Phaeodactylum tricornutum with six bacterial strains to quantify bacteria-mediated differences in algal biomass, total intracellular lipids, and for a subset, extracellular metabolite accumulation. A Marinobacter isolate significantly increased algal cell concentrations, dry biomass, and lipid content compared to axenic algal cultures. Two other bacterial strains from the Bacteroidetes order, of the genera Algoriphagus and Muricauda, significantly lowered P. tricornutum biomass, leading to overall decreased lipid accumulation. These three bacterial co-cultures (one mutualistic, two competitive) were analyzed for extracellular metabolites via untargeted liquid chromatography mass spectrometry to compare against bacteria-free cultures. Over 80% of the extracellular metabolites differentially abundant in at least one treatment were in higher concentrations in the axenic cultures, in agreement with the hypothesis that the co-cultured bacteria incorporated algal-derived organic compounds for growth. Furthermore, the extracellular metabolite profiles of the two growth-inhibiting cultures were more similar to one another than the growth-promoting co-culture, linking metabolite patterns to ecological role. Our results show that algal–bacterial interactions can influence the accumulation of intracellular lipids and extracellular metabolites, and suggest that utilization and accumulation of compounds outside the cell play a role in regulating microbial interactions.

59 BASIC BIOLOGICAL SCIENCES↗

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Initial development of an in-situ, passive material surveillance test article for monitoring high temperature reactor structural components

This report describes the development of a family of passively-actuated surveillance test articles, suitable for use in an in-situ reactor test program monitoring the degradation in structural material properties caused by the reactor environment. The test articles apply creep-fatigue-type cyclic load to a test section driven by changes in the reactor coolant temperature converted to mechanical force by differential thermal expansion. The report describes three types of test articles, all of which apply the load to a 316H test section using Alloy 617 as the driver material. The three types of specimens each serve a unique purpose: (1) initial proof-of-concept, (2) validation of the models used to size the specimens to achieve a given strain range, and (3) demonstration that the articles can be miniaturized for use in an actual operating reactor. In addition to the specimen design, fabrication, and testing, this report describes the outline of two types of in-situ monitoring programs, using the test article. Monitoring programs such as these could be a critical component in gaining regulator acceptance for Molten Salt Reactor (MSR) designs using structural materials for which there is little long-term experimental environmental testing or operating experience available for calibrating the component structural design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Early immune response to Coccidioides is characterized by robust neutrophil and fibrotic macrophage recruitment and differentiation

Coccidioidomycosis, or Valley fever, is an emerging respiratory disease caused by soil-dwelling fungi of the Coccidioides genus that is expected to spread from the southwest into the central U.S. by 2050. While 60% of infections are asymptomatic, the other 40% of patients experience a range of symptoms, from self-limiting pneumonia to life-threatening disseminated disease. The immunological events that underlie the progression to severe disease remain underdefined. Here, we probed the early immune response to Coccidioides using a high dose of an attenuated strain of Coccidioides posadasii in a mouse model of infection coupled with single-cell RNA sequencing. At 24 h post-infection, robust immune infiltration is detected in the lung, marked by high levels of inflammatory PD-L1 + neutrophils and fungal-contact-dependent pro-fibrotic Spp1 + macrophages. These findings elucidate the early dynamics of the host response to Coccidioides and provide a deeper understanding of host-pathogen interactions in the lung.

Coccidioides↗

Harnessing redox proteomics to study metabolic regulation and stress response in lignin-fed Rhodococci

Abstract Background Rhodococci are studied for their bacterial ligninolytic capabilities and proclivity to accumulate lipids. Lignin utilization is a resource intensive process requiring a variety of redox active enzymes and cofactors for degradation as well as defense against the resulting toxic byproducts and oxidative conditions. Studying enzyme expression and regulation between carbon sources will help decode the metabolic rewiring that stymies lignin to lipid conversion in these bacteria. Herein, a redox proteomics approach was applied to investigate a fundamental driver of carbon catabolism and lipid anabolism: redox balance. Results A consortium of Rhodococcus strains was employed in this study given its higher capacity for lignin degradation compared to monocultures. This consortium was grown on glucose vs. lignin under nitrogen limitation to study the importance of redox balance as it relates to nutrient availability. A modified bottom–up proteomics workflow was harnessed to acquire a general relationship between protein abundance and protein redox states. Global proteomics results affirm differential expression of enzymes involved in sugar metabolism vs. those involved in lignin degradation and aromatics metabolism. As reported previously, several enzymes in the lipid biosynthetic pathways were downregulated, whereas many involved in β-oxidation were upregulated. Interestingly, proteins involved in oxidative stress response were also upregulated perhaps in response to lignin degradation and aromatics catabolism, which require oxygen and reactive oxygen species and generate toxic byproducts. Enzymes displaying little-to-no change in abundance but differences in redox state were observed in various pathways for carbon utilization (e.g., β‑ketoadipate pathway), lipid metabolism, as well as nitrogen metabolism (e.g., purine scavenging/synthesis), suggesting potential mechanisms of redox-dependent regulation of metabolism. Conclusions Efficient lipid production requires a steady carbon and energy flux while balancing fundamental requirements for enzyme production and cell maintenance. For lignin, we theorize that this balance is difficult to establish due to resource expenditure for enzyme production and stress response. This is supported by significant changes to protein abundances and protein cysteine oxidation in various metabolic pathways and redox processes.

09 BIOMASS FUELS↗

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.↗