Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuous relaxation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tuning the spin dynamics and magnetic phase transitions of the Cantor alloy via composition and sample processing protocols: A muon spin relaxation study

CrMnFeCoNi, also called the Cantor alloy, is a well-known high-entropy alloy whose magnetic properties have recently become a focus of attention. Here, we present a detailed muon spin relaxation study of the influence of chemical composition and sample processing protocols on the magnetic phase transitions and spin dynamics of several different Cantor alloy samples. Specific samples studied include a pristine equiatomic sample, samples with deficient and excess Mn content, and equiatomic samples magnetized in a field of 9 T or plastically deformed in pressures up to 0.5 GPa. The results confirm the sensitive dependence of the transition temperature on composition and demonstrate that post-synthesis pressure treatments cause the transition to become significantly less homogeneous throughout the sample volume. In addition, we observe critical spin dynamics in the vicinity of the transition in all samples, reminiscent of canonical spin glasses and magnetic materials with ideal continuous phase transitions. Application of an external magnetic field suppresses the critical dynamics in the Mn-deficient sample, while the equiatomic and Mn-rich samples show more robust critical dynamics. The spin-flip thermal activation energy in the paramagnetic phase increases with Mn content, ranging from 3.1⁢(3) × 10 -21 J for 0% Mn to 1.2⁢(2) × 10 -20 J for 30% Mn content. These results shed light on critical magnetic behavior in environments of extreme chemical disorder and demonstrate the tunability of spin dynamics in the Cantor alloy via chemical composition and sample processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

D2NO: Efficient handling of heterogeneous input function spaces with distributed deep neural operators

Neural operators have been applied in various scientific fields, such as solving parametric partial differential equations, dynamical systems with control, and inverse problems. However, challenges arise when dealing with input functions that exhibit heterogeneous properties, requiring multiple sensors to handle functions with minimal regularity. To address this issue, discretization-invariant neural operators have been used, allowing the sampling of diverse input functions with different sensor locations. However, existing frameworks still require an equal number of sensors for all functions. We propose a novel distributed approach to further relax the discretization requirements and solve the heterogeneous dataset challenges. Our method involves partitioning the input function space and processing individual input functions using independent and separate neural networks. A centralized neural network is used to handle shared information across all output functions. This distributed methodology reduces the number of gradient descent back-propagation steps, improving efficiency while maintaining accuracy. Here, we demonstrate that the corresponding neural network is a universal approximator of continuous nonlinear operators and present three numerical examples to validate its performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗
DOE OSTI2023

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous- enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗
DOE OSTI2022

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks: Preprint

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous–enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗
DOE OSTI2023

Spin-squeezed Gottesman-Kitaev-Preskill codes for quantum error correction in atomic ensembles

Gottesman-Kitaev-Preskill (GKP) codes encode a qubit in displaced phase-space combs of a continuous-variable (CV) quantum system and are useful for correcting a variety of high-weight photonic errors. Here we propose atomic ensemble analogs of the single-mode CV GKP code by using the quantum central limit theorem to pull back the phase-space structure of a CV system to the compact phase space of a quantum spin system. We study the optimal recovery performance of these codes under error channels described by stochastic relaxation and isotropic ballistic dephasing processes using the diversity combining approach for calculating channel fidelity. Additionally, we find that the spin GKP codes outperform other spin system codes such as cat codes or binomial codes. Our spin GKP codes based on the two-axis countertwisting interaction and superpositions of SU(2) coherent states are direct spin analogs of the finite-energy CV GKP codes, whereas our codes based on one-axis twisting do not yet have well-studied CV analogs. A state preparation scheme for the spin GKP codes is proposed which uses the linear-combination-of-unitaries method, applicable to both the CV and spin GKP settings. Finally, we discuss a fault-tolerant approximate gate set for quantum computing with spin-GKP-encoded qubits, obtained by translating gates from the CV GKP setting using the quantum central limit theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Multigrid Reduction in Time for Chaotic and Hyperbolic Problems (Final Report)

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallel-in-time methods, with the most common method, parareal, diverging in many cases. The current state of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be further researched and improved (including with a new relaxation method inspired by on Least Squares Shadowing (LSS)) and explored for more complicated problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2021

Coupled valence carrier and core-exciton dynamics in WS 2 probed by few-femtosecond extreme ultraviolet transient absorption spectroscopy

Few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy, performed with optical 500-1000 nm supercontinuum and broadband XUV pulses (30-50 eV), simultaneously probes dynamics of photoexcited carriers in WS$_{2}$ at the W O$_3$ edge (37-45 eV) and carrier-induced modifications of core-exciton absorption at the W N$_{6,7}$ edge (32-37 eV). Access to continuous core-to-conduction band absorption features and discrete core-exciton transitions in the same XUV spectral region in a semiconductor provides a novel means to investigate the effect of carrier excitation on core-exciton dynamics. The core-level transient absorption spectra, measured with either pulse arriving first to explore both core-level and valence carrier dynamics, reveal that core-exciton transitions are strongly influenced by the photoexcited carriers. A $1.2\pm0.3$ ps hole-phonon relaxation time and a $3.1\pm0.4$ ps carrier recombination time are extracted from the XUV transient absorption spectra from the core-to-conduction band transitions at the W O$_{3}$ edge. Global fitting of the transient absorption signal at the W N$_{6,7}$ edge yields $\sim 10$ fs coherence lifetimes of core-exciton states and reveals that the photoexcited carriers, which alter the electronic screening and band filling, are the dominant contributor to the spectral modifications of core-excitons and direct field-induced changes play a minor role. This work provides a first look at the modulations of core-exciton states by photoexcited carriers and advances our understanding of carrier dynamics in metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Initial development of a method for correlating indentation test results to damage accumulation in high temperature structural materials

This report describes an experimental approach based on microindentation creep or relaxation tests for assessing the remaining creep life of a sample material. The objective for this experimental program is to complement the development of passively loaded in situ surveillance specimens. These surveillance specimens could be placed in future advanced reactors to assess the degradation of structural properties, primary creep and creep-fatigue life, when the material is exposed to actual reactor operating conditions. The specimens and the corresponding indentation test protocol could form the basis of a structural health monitoring program designed to ensure the safe continued operation of the component given the relatively unknown effect of the reactor coolant chemistry, radiation exposure, and other environmental effects on the material properties of the constituent structural materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Quantum Zeno Control of Superconducting Qubit Coherence

The Quantum Zeno Effect (QZE) dictates how the dynamics of a quantum system can be modified through continuous or discrete measurements [1]. It has gained increasing attention in quantum computing community in recent years, both due to its inevitable implications for qubit readout as well as for its promise for enabling new methods of quantum state control such as reservoir engineering and Zeno-dragging. Here, we investigate Zeno effects implemented via weak measurements for controlling superconducting-qubit coherence during gate and readout operations in a multi-qubit setup. Building on prior observations [2] that measurement backaction can both suppress and enhance qubit relaxation, we predict and measure QZE-altered qubit coherence times by manipulating the spectral overlap of qubit spectrum with background noise sources. We also discuss modifications and opportunities for controllable QZE due to anharmonic effects in multi-level superconducting atoms. [1] S. Greenfield, A. Kamal, J. Dressel, E. Levenson-Falk, arXiv:2506.12679 (2025) [2] Thorbeck, Z. Xiao, L. Govia, A. Kamal, Phys. Rev. Lett. 132, 090602 (2024)

Seidel, Olivia [Texas U., Arlington; Fermilab]↗
DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗
DOE OSTI2020

In-Operando Study of Shape Retention and Microstructure Development in a Hydrolyzing Sol–Gel Ink during 3D-Printing

3D printing of amorphous and crystalline ceramics is of paramount importance for the fabrication of a wide range of devices with applications across different technology fields. Printed ceramics are remarkably enabled by the sol-gel synthesis method in conjunction with continuous filament direct ink writing. During printing, multiple processes contribute to the evolution of inks including shape retention, chemical conversion, solidification, and microstructure formation. Traditionally, depending on the ink composition and printing environment, several mechanisms have been associated with the shape retention and solidification of 3D printed structures: gelation, rapid solvent evaporation, energy-driven phase transformation, and chemical-driven phase transformation. Understanding the fundamental differences between these mechanisms becomes key since they strongly influence the spatiotemporal evolution of the materials, as the out-of-equilibrium processes inherent to the extrusion, relaxation, and solidification of printed materials have significant effects on the materials properties. Herein, we investigate the shape retention mechanism and the hydrolysis-induced material conversion and microstructure formation during the 3D printing of a water reactive sol-gel ink that transforms into titanium dioxide-based ceramic. This study aims at identifying characteristic mechanisms associated with the material transformation, establishing connections between the microstructure development and the timescales associated with solidification under operando 3D-printing conditions. The investigation of this material's out-of-equilibrium pathways under processing conditions is enabled by time-resolved coherent X-ray scattering, providing simultaneous access to temporospatially resolved microstructural and dynamics information. Furthermore, we explore X-ray speckle tracking as a tool to resolve deformations of the microstructure in a printed filament associated with the deposition of consecutive filaments. Through this work, we aim at providing a fundamental understanding of the relationships behind these transformative processes in 3D printing and their timescales as the basis for achieving unprecedented control over printed materials microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetoelastic coupling, negative thermal expansion, and two-dimensional magnetic excitations in FeAs

We present a temperature-dependent investigation of the local structure and magnetic dynamics of the FeAs binary. The magnetic susceptibility χ (T) result shows an anomalous broad feature up to 550 K, with χ continuing to increase above the Néel antiferromagnetic ordering temperature (T N =70K), peaking at ~250 K, then decreasing gently above. It is remarkable that this peak susceptibility temperature corresponds to the onset of anisotropic negative thermal expansion in both the a and c axes, suggesting that magnetic interactions are affecting the structure even well above the Néel point. A systematic investigation into local bonding correlations, from time-of-flight neutron pair distribution function analyses, shows the octahedral volume around each Fe site growing monotonically while adjacent octahedra tilt toward one another before relaxing away past this peak in the anomalous magnetic susceptibility. We use inelastic-neutron scattering to map spin-wave excitations in FeAs at temperatures above and below the T N . We find magnetic excitations near T N to be very different from the excitations in the ground state at 1.5 K. Spin waves measured at 1.5 K are three dimensional (3D), however, in the vicinity of the magnetic transition, the magnetic fluctuations clearly indicate two-dimensional (2D) character in this intrinsically 3D crystal structure. Unlike the undoped 2D parents of iron-arsenide superconductors, where the magnetic correlations are considerably weaker along the c axis than in the ab plane, inelastic neutron scattering here shows that the spin fluctuations in the 3D FeAs binary are nearly 2D in the bc plane at 90 K. These results demonstrate the importance of short-range correlations in understanding the magnetic properties of transition-metal binaries, and suggest how 2D excitations, even in a 3D structure, can potentially become a breeding ground for unconventional superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Large Exciton Diffusion Coefficients in Two Dimensional Covalent Organic Frameworks with Different Domain Sizes Revealed by Ultrafast Exciton Dynamics

Large singlet exciton diffusion lengths are a hallmark of high performance in organic-based devices such as photovoltaics, chemical sensors, and photodetectors. In this study, exciton dynamics of a two-dimensional covalent organic framework, 2D COF-S, is investigated using ultrafast spectroscopic techniques. After photoexcitation, the COF-S exciton decays via three pathways: (1) excimer formation (4 +/- 2 ps), (2) excimer relaxation (160 +/- 40 ps), and (3) excimer decay (>3 ns). Excitation fluence-dependent transient absorption studies suggest that COF-5 has a relatively large diffusion coefficient (0.08 cm(2)/s). Furthermore, exciton-exciton annihilation processes are characterized as a function of COF-S crystallite domain size in four different samples, which reveal domain-size-dependent exciton diffusion kinetics. These results reveal that exciton diffusion in COF-S is constrained by its crystalline domain size. These insights indicate the outstanding promise of delocalized excitonic processes available in 2D COFs, which motivate their continued design and implementation into optoelectronic devices.

Flanders, Nathan C.↗
DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Ginzburg--Landau functionals in the large-graph limit

Ginzburg–Landau (GL) functionals on graphs, which are relaxations of graph-cut functionals on graphs, have yielded a variety of insights in image segmentation and graph clustering. In this paper, we study large-graph limits of GL functionals by taking a functional-analytic view of graphs as nonlocal kernels. For a graph Wn with n nodes, the corresponding graph GL functional GL W n ϵ is an energy for functions on Wn. We minimize GL functionals on sequences of growing graphs that converge to functions called graphons. For such sequences of graphs, we show that the graph GL functional Γ-converges to a continuous and nonlocal functional that we call the graphon GL functional. We investigate the sharp-interface limits of the graph GL and graphon GL functionals, and we relate these limits to a nonlocal total-variation (TV) functional. We express the limiting GL functional in terms of Young measures and thereby obtain a probabilistic interpretation of the minimization problem in the large-graph limit. Finally, to develop intuition about graphon GL functionals, we determine the GL minimizer for several example families of graphons.

Zhang, Edith↗
DOE OSTI2021

Inference of Rock Flow and Mechanical Properties from Injection-Induced Microseismic Events During Geologic CO 2 Storage

Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.

42 ENGINEERING↗