Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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

To grind or not to grind? The influence of mechanical and thermal treatments on the $\mathrm{Cu/Zn}$ disorder in $\mathrm{Cu}$ 2 $\mathrm{ZnSn}$($\mathrm{S}$ x $\mathrm{Se}$ 1-x ) 4 monograins

We report that Kesterite-type based thin film solar cell technologies are mainly based on polycrystalline absorber layers. A promising low cost alternative technology uses Cu 2 ZnSn(S x Se 1-x ) 4 (CZTSSe) monograins (single crystals of 20–100 μm size) which are fixed in a polymer matrix to form a flexible solar cell. The Cu/Zn disorder is discussed as a possible reason for band tailing causing voltage losses limiting the efficiency of CZTSSe-based devices. The experimental determination of the order parameter Q which is a quantitative measure of Cu/Zn disorder, requires a differentiation between the isoelectronic cations Cu + and Zn 2+ . An in-depth analysis of neutron diffraction data allows the determination of type and concentration of intrinsic point defects including a distinction between Cu and Zn. Neutron diffraction requires large sample volumes, thus monograins offer the unique possibility to correlate structural disorder in CZTSSe with device performance parameters. In this study we tackle the influence of grinding the monograins on stoichiometry deviations, the Cu/Zn disorder as well as intrinsic point defects and optoelectronic properties of CZTSSe monograins. Moreover, an easy methodology based on Raman scattering spectroscopy is proposed for the assessment of Cu/Zn disorder in the CZTSSe compounds.

14 SOLAR ENERGY↗
DOE OSTI2020

A Combined Day-ahead and Real-time approach for Scheduling Real and Reactive power dispatch of Battery Energy Storage System

A combined day-ahead (look-ahead) and real-time optimization approach is proposed to schedule BESS real and reactive power dispatch. The look ahead scheduling agent maximizes profit and essentially performs the function of energy arbitrage by looking at forecasted prices, network constraints and BESS operating conditions to generate real power dispatch schedule. While this schedule is followed, the real time optimization agent computes an optimum reactive power set-point (for a particular objective function) using optimal power flow framework. The concepts are implemented by co-simulating legacy tools and demonstrated on a large real world distribution system. In addition to the profits from energy arbitrage, optimum reactive power set point computation results in achieving desirable network attributes like flat voltage profiles and minimizing voltage deviations from desired value.

Bhatti, Bilal A.↗
DOE OSTI2025

Evidence of isospin-symmetry violation in high-energy collisions of atomic nuclei

Strong interactions preserve an approximate isospin symmetry between up (u) and down (d) quarks, part of the more general flavor symmetry. In the case of K meson production, if this isospin symmetry were exact, it would result in equal numbers of charged (K + and K − ) and neutral (K 0 and ${\overline{K}}^{0}$) mesons produced in collisions of isospin-symmetric atomic nuclei. Here, we report results on the relative abundance of charged over neutral K meson production in argon and scandium nuclei collisions at a center-of-mass energy of 11.9 GeV per nucleon pair. We find that the production of K + and K − mesons at mid-rapidity is (18.4 ± 6.1)% higher than that of the neutral K mesons. Although with large uncertainties, earlier data on nucleus-nucleus collisions in the collision center-of-mass energy range $2.6 \, < \, \sqrt{{s}_{NN}} \, < \, 200$ GeV are consistent with the present result. Using well-established models for hadron production, we demonstrate that known isospin-symmetry breaking effects and the initial nuclei containing more neutrons than protons lead only to a small (few percent) deviation of the charged-to-neutral kaon ratio from unity at high energies. Thus, they cannot explain the measurements. The significance of the flavor-symmetry violation beyond the known effects is 4.7σ when the compilation of world data with uncertainties quoted by the experiments is used. New systematic, high-precision measurements and theoretical efforts are needed to establish the origin of the observed large isospin-symmetry breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Accuracy of power spectra in dissipationless cosmological simulations

Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, of scale-free models with a range of spectral indices n, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At 1 per cent level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, ${\sim}2{-}3 k_\Lambda$ at the very latest times we evolve to. At the $5{{\ \rm per\ cent}}$ level, accuracy extends up to wavenumbers of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple n-dependence, allowing a extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale-free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small-scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_\Lambda$ may need revision in light of our convergence analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Probing late-stage hadronic interactions at high baryon density via $K^{⁎0}$ production in the RHIC beam energy scan program

A precision measurement of the $K^{⁎0}$ meson yield is reported in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.6, 19.6, and 27 GeV using the high-statistics data sample collected by the STAR experiment during the Beam Energy Scan II (BES-II) program at RHIC. The transeverse momentum (p T )-integrated yield ratios $\large{(}K^{⁎0} + \overline{K^{⁎0}}\large{)}/(K^+ + K^{-})$ in central collisions show a suppression relative to peripheral collisions at the (1.7–3.6) σ level, while a thermal model without final-stage rescattering overpredicts this ratio with a deviation of (6.9–8.2) σ. These results indicate a loss of the measured $K^{⁎0}$ signal in central collisions due to re-scattering of its hadronic decay products in the hadronic phase. The p T -integrated yield of charged kaons exhibits an approximate scaling with charged-particle multiplicity, independent of collision energy and system size. A similar trend is observed for the short-lived $K^{⁎0}$ resonance, although significant deviations emerge at lower energies. At BES energies, the $K^{⁎0}/K$ ratio shows stronger suppression than at the highest RHIC and LHC energies within a given multiplicity bin, particularly in central and mid-central collisions. This behavior is consistent with changes in the effective hadronic interaction cross section and is supported by transport model calculations, which indicate dominant meson–baryon interactions at lower energies and meson–meson interactions at higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A new coupling of a GPU-resident large-eddy simulation code with a multiphysics wind turbine simulation tool

The development of new wind farm control strategies can benefit from combined analysis of flow dynamics in the farm and the behavior of individual turbines within one simulation environment. In this work, we present such an environment by developing a new coupling between the large-eddy simulation (LES) code GRASP and the multiphysics wind turbine simulation tool OpenFAST via an actuator line model (ALM). In addition, the implementation of the recently proposed filtered actuator line model (FALM) within the coupling is described. The new ALM implementation is cross-verified with results from four other commonly used research LES codes. The results for the blade loads and the near wake obtained with the new coupling are consistent with the other codes. Deviations are observed in the far wake. The results further indicate that the FALM is able to reduce the lift and power overprediction from which the traditional ALM suffers on coarse LES grids. This new simulation environment paves the way for future wind farm simulations under realistic weather conditions by leveraging GRASP's ability to impose data from large-scale meteorological models as boundary conditions.

17 WIND ENERGY↗
DOE OSTI2022

Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection

This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microstructure in matrix elements

We investigate the simple model of Pennington, Shenker, Stanford and Yang for modeling the density matrix of Hawking radiation, but further include dynamics for EOW branes behind the horizon. This allows interactions that scatter one interior state to another, and also allows EOW loops. At strong coupling, we find that EOW states are no longer random; the ensemble has collapsed, and coupling constants encode the microscopic matrix elements of Hawking radiation. This suggests strong interior dynamics are important for understanding evaporating black holes, without any ensemble average. In this concrete model the density matrix of the radiation deviates from the thermal state, small off-diagonal fluctuations encode equivalences between naively orthogonal states, and bound the entropy from above. For almost evaporated black holes the off-diagonal terms become as large as the diagonal ones, eventually giving a pure state. We also find the unique analytic formula for all Renyi entropies.

2D gravity↗
DOE OSTI2021

Effects of magnetic fields on HPGe tracking detectors

Here, we present a study of magnetic fields effects on the position resolution and energy response of hyper-pure germanium detectors. Our results provide realistic estimates of the potential impact on the resolving power of tracking-arrays from (fringe) magnetic fields present when operating together with large spectrometers. By solving the equations of motion for the electron and holes in the presence of both electric and magnetic fields, we analyzed the drift trajectories of the charge carriers to determine the deviations in the positions at the end point of the trajectories, as well as changes in drift lengths affecting the energy resolution and peak shift due to trapping. Our results show that the major effect is in the deviation of the transverse (to the electric field direction) position and suggest that, if no corrective action is taken in the pulse-shape and tracking data analysis procedures, a field strength ≳ 0.1 T will start to impact the intrinsic position resolution of 2 mm (RMS). At fields above ~1 T, the degradation of the energy response becomes observable.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Modeling Henry's law and phase separations of water–NaCl–organic mixtures with solvation and ion-pairing

Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H 2 O–NaCl–MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor–liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid–liquid equilibrium of H 2 O–NaCl–MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H 2 O–NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid–liquid equilibrium (LLE) boundary of the H 2 O–NaCl–MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water–MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid–liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. Finally, this work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational–vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric ( 18 O 16 O 18 O) and asymmetric ( 18 O 18 O 16 O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric ( 16 O 18 O 16 O) and asymmetric ( 16 O 16 O 18 O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The Challenge of Disproportionate Importance of Temporal Features in Predicting HPC Power Consumption

In this work, we demonstrate the challenges in predicting HPC cluster power consumption in the face of significant temporal skew in power consumption behavioral patterns. Predicting large power swings that extend several megawatts has significant operational value for HPC centers, however, prediction is challenging due to the relative rarity of such events and also due to the abrupt or disjoint deviation from the average power consumption levels. To study the impact of this challenge, we have trained a recurrent neural network (RNN) as a reasonably sophisticated model to predict power consumption of the one-year worth of node power consumption data from the Summit supercomputer located in the Oak Ridge Leadership Computing Facility. By studying the prediction results, we have found that although simple usage of RNN models can provide good results on average power consumption levels, it would fail at predicting the power swings that have more operational value. With such results, we discuss potential next steps in addressing such issues aiming towards a robust usage of power prediction techniques in HPC operations.

Li, Chengcheng↗
DOE OSTI2025

Measurement of the associated production of a top-antitop-quark pair and a Higgs boson decaying into a $b\bar{b}$ pair in pp collisions at $\sqrt{s} = 13$ TeV using the ATLAS detector at the LHC

This paper reports the measurement of Higgs boson production in association with a $t\bar{t}$ pair in the $H$ → $b\bar{b}$ decay channel. The analysis uses 140 fb -1 of 13 TeV proton–proton collision data collected with the ATLAS detector at the Large Hadron Collider. The final states with one or two electrons or muons are employed. An excess of events over the expected background is found with an observed (expected) significance of 4.6 (5.4) standard deviations. The $t\bar{t}$ $H$ cross-section is $σ_{t\bar{t} H}$ = $411^{+101}_{-92}$ fb = $411 ± 54$ (stat.) $^{+85}_{-75}$ (syst.) fb for a Higgs boson mass of 125.09 GeV, consistent with the prediction of the Standard Model of $507^{+35}_{-75}$ fb. The cross-section is also measured differentially in bins of the Higgs boson transverse momentum within the simplified template cross-section framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

LES and RANS Spray Combustion Analysis of OME 3-5 and n-Dodecane

Clean-burning oxygenated and synthetic fuels derived from renewable power, so-called e-fuels, are a promising pathway to decarbonize compression–ignition engines. Polyoxymethylene dimethyl ethers (PODEs or OMEs) are one candidate of such fuels with good prospects. Their lack of carbon-to-carbon bonds and high concentration of chemically bound oxygen effectively negate the emergence of polycyclic aromatic hydrocarbons (PAHs) and even their precursors like acetylene (C 2 H 2 ), enabling soot-free combustion without the soot-NO x trade-off common for diesel engines. The differences in the spray combustion process for OMEs and diesel-like reference fuels like n-dodecane and their potential implications on engine applications include discrepancies in the observed ignition delay, the stabilized flame lift-off location, and significant deviations in high-temperature flame morphology. For CFD simulations, the accurate modeling and prediction of these differences between OMEs and n-dodecane proved challenging. This study investigates the spray combustion process of an OME 3 − 5 mixture and n-dodecane with advanced optical diagnostics, Reynolds-Averaged Navier–Stokes (RANS), and Large-Eddy Simulations (LESs) within a constant-volume vessel. Cool-flame and high-temperature combustion were measured simultaneously via high-speed (50 kHz) imaging with formaldehyde (CH 2 O) planar laser-induced fluorescence (PLIF) representing the former and line-of-sight OH* chemiluminescence the latter. Both RANS and LES simulations accurately describe the cool-flame development process with the formation of CH 2 O. However, CH 2 O consumption and the onset of high-temperature reactions, signaled by the rise of OH* levels, show significant deviations between RANS, LES, and experiments as well as between n-dodecane and OME. A focus is set on the quality of the simulated results compared to the experimentally observed spatial distribution of OH*, especially in OME fuel-rich regions. The influence of the turbulence modeling is investigated for the two distinct ambient temperatures of 900 K and 1200 K within the Engine Combustion Network Spray A setup. The capabilities and limitations of the RANS simulations are demonstrated with the initial cool-flame propagation and periodic oscillations of CH 2 O formation/consumption during the quasi-steady combustion period captured by the LES.

CFD↗
DOE OSTI2026

Vector boson scattering and anomalous quartic couplings in final states with ℓνqq or ℓℓqq plus jets using proton-proton collisions at $\sqrt{s}=13$ TeV

A measurement is presented of the electroweak vector boson scattering production of ZV (V = W, Z) boson pairs associated with two jets in proton-proton collisions at a center-of-mass energy of 13 TeV. The data, corresponding to an integrated luminosity of 138 fb −1 , were collected at the CERN LHC with the CMS detector during the 2016–2018 data-taking period. The analysis targets final states with a pair of isolated electrons or muons from Z boson decays and three or four jets, depending on the momentum of the vector boson that decays into quarks. Signal strength is measured for events characterized by a large invariant mass of two forward jets with a wide pseudorapidity gap between them. The electroweak production of ZV in association with two jets is measured with an observed (expected) significance of 1.3 (1.8) standard deviations. A combination of the analyses of ZV channel and the previously published WV channel in the lepton plus jets final state places constraints on effective field theory parameters that describe anomalous electroweak production of WW, WZ, and ZZ boson pairs in association with two jets. Several world best limits are set on anomalous quartic gauge couplings in terms of dimension-8 standard model effective field theory operators.

Hadron-Hadron Scattering↗
DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗