SEARCH · Engineering Papers

Results for “velocity steps”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Fiber coupled laser ultrasound system using a single mode hollow core fiber for excitation laser: Theory and demonstration for on-machine thickness gauging

Laser ultrasound (LU) is a technique that uses a pump laser and a probe laser to optically generate and detect elastic waves in a material. Despite its advantages over traditional contact transducer-based ultrasound, industrial adoption has been limited by complex optical setups and the inability of multi-mode fibers to deliver a stable Gaussian profile for the pump laser. Here, we report a fully fiber-coupled thermoelastic LU system that uses an anti-resonant hollow-core single-mode fiber to deliver 1 mJ nanosecond pulses of 1064 nm light, while preserving the fundamental Gaussian mode (pump laser). When combined with a fiber-coupled interferometer (probe laser), a small, flexible, and environmentally robust sensor capable of optically generating and detecting high frequency broadband ultrasound is realized. We demonstrate such an LU system implemented in situ on a four-axis precision lathe. High-resolution thickness gauging is performed, before and after precision cutting, by exciting and measuring a zero-group velocity guided wave mode. The measurements are verified with ex-situ traceable coordinate measuring machine data. Mean absolute deviations of 0.1%, of nominal thickness, before cutting, and 0.2% and 0.3%, after stepped and tapered cuts, respectively, are reported. A theoretical background for thermoelastic ultrasound generation in an elastic waveguide is also presented. Attention is given to the effect of the pump laser profile on wave generation to elucidate the importance of using single-mode laser light. The fiber-coupled system demonstrated is well-suited for use in scientific and engineering sensing applications and facilitates the adoption of LU for industrial non-destructive testing.

Engineering↗

DOE OSTI2026

Angular momentum of vacuum bubbles in a first-order phase transition

The formation of primordial black holes (PBHs) during a first-order phase transition (FOPT) in a dark sector has been of recent interest. A quantity that characterizes a black hole is its spin. We carry out the first step towards determining the spin of such PBHs, by calculating the spin of spherical false vacuum bubbles induced by cosmological perturbations. The angular momentum is given by the product of density and velocity perturbations. We carefully track the evolution of background quantities and calculate the transfer functions during the FOPT. We find that the dimensionless spin parameter s = J /(G N M 2 ) of false vacuum bubbles of mass M and angular momentum J , take a wide range of values from 𝒪(10 -5 ) to 𝒪(10) for FOPTs between 10 keV and 100 GeV and a dark sector that is 0.1 to 0.4 times cooler than the visible sector. We also find a scaling relation between the root-mean-square value of the spin, the FOPT time scale, the bubble wall velocity, and the dark sector-to-visible sector temperature ratio.

cosmological phase transitions↗

DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Theoretical background for a fast flow liquid metal divertor experiment

Tokamak Energy has developed a liquid metal experiment featuring a lithium loop that circulates lithium through the HIDRA stellarator at the University of Illinois at Urbana-Champaign. The loop includes a replaceable divertor module, designed to demonstrate fast, steady, open-surface liquid metal flow. The first-generation divertor module was deliberately kept simple—an open-surface chute—to enable modelling validation and establish a benchmark for future design iterations. This paper presents the theoretical foundations of the experimental design. Specifically, we identify module overflooding as the primary challenge and determine the limits of fluid velocity and magnetic field strength necessary to prevent it. The steady-state flow patterns are expected to exhibit relatively slow surface velocities, which do not fully align with the fast-flow concept. Nevertheless, the experiment is designed to achieve controllable, continuous open-surface flow within an operational fusion device, representing a significant step toward the development of liquid metal divertor technology.

Experiment↗

DOE OSTI2020

Calculation Of The First Moment Of Energy Using D-T Reactivity Formalisms Under The Maxwell-Boltzmann Distribution--Part II

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of a nuclear fusion formalism is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium- Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hal reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell-Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. This report continues the step toward building understanding of nuclear fusion reactivity formalisms. It is part of a series of reports with the same goal, [5-10] and is the continuation of the Part I paper for defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, <$E$>, with the fusion cross-section, fusion reactivity and its derivative with ion-temperature. In Part I, three variants of the first moment <$E$> were analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment from Brysk. In Part II, attention is given to the definition of <$K$>, originally described as a ratio of moments from Brysk, and its relationship to the first moment definitions from Part I. One measure of the progress made in these documents is the identification that Brysk’s ratio of the second moment to the first moment ratio, <$K$>, does not correspond to his provided solution of the first moment of energy. Another measure of (our) progress from this work is the comparison of first moment variants. That comparison includes confirming the importance of cross sections defined in terms of energy. The analytical relationships we have developed among important physics quantities are useful tools in validation and verification (V&V). For example, we can calculate the kinetic energy as a mean or as a first moment, or as a function of the first moment. These analytically-determined values can be compared directly with numerically-determined values, supplied to the authors, representing <$E$>.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Steady Laminar Flame for 1-Step Thermal Decomposition Model of PBX 9501

We have recently developed a new 1-step thermal decomposition model of PBX 9501 high explosive. In this report, we examine the behavior of the model in 1-D steady laminar flame calculations. In particular, the dependencies of the steady laminar flame velocity on initial pressure, thermal conductivies and parameters of the thermal decomposition model are investigated. The results of our numerical calculations are compared to available experimental data.

42 ENGINEERING↗

DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗

DOE OSTI2021

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

CaloTrilogy: Toward a Breakthrough in One-Step, End-to-End, Physics-Guided Shower Generation for Modern Calorimeters

High-precision calorimeter simulation at current and future colliders imposes rapidly growing computational demands, motivating the development of machine-learning surrogates for traditional Monte Carlo tools such as Geant4. Flow matching and diffusion-based generative models have become leading approaches for high-dimensional fast simulation because of their sample quality, but typically require ${\cal O}(100)$ function evaluations at inference and often rely on auxiliary networks to constrain global observables, compromising streamlined end-to-end generation. We introduce a unified framework that improves the balance between speed, shower quality, and physics fidelity. The method combines: (i) an average velocity field integrator that enables sampling in one or a few evaluations; (ii) a learned generative prior in shower space, constructed from data rather than random noise; and (iii) physics-guided loss terms that impose inductive biases on key observables during training. These elements are training time regularizers, preserving end-to-end inference with no additional cost. With only one or a few evaluation steps, the model achieves shower quality competitive with state-of-the-art flow and diffusion approaches, tested on several public high granularity calorimeter datasets. The results demonstrate inter-layer shower structure consistent with the underlying physics, providing a strong candidate for future fast simulation workflows.

Jiang, Cheng [Edinburgh U.]↗

DOE OSTI2022

Title Coherent X-ray Studies of Surface Growth and Patterning Processes

X-ray Photon Correlation Spectroscopy (XPCS) is being developed as a tool to study nanoscale dynamics of fluctuations during thin film growth and surface patterning. XPCS examines the evolution of the X-ray scattering speckle pattern in reciprocal space to reveal dynamics information not accessible through any other means. Unlike low-coherence conventional X-ray scattering, which incoherently averages over different regions of a sample, coherent X-ray scattering is sensitive to the detailed structure of a given sample at that moment in time. In thin film growth processes, heterodyning, which occurs due to coherent mixing of two scattered signals, can be used to investigate the relationship between surface growth velocity and defect propagation. And, for polycrystalline thin film growth, the spatial coherence of the X-ray beam can substitute for the missing spatial coherence of the growth process to track layered growth in detail, even in a growth regime in which there are no conventional growth oscillations of the X-ray intensity. This has opened the door for detailed dynamics studies of individual atomic layers during real-world growth and patterning, not just in perfect single-crystal growth cases, which greatly expands the applicability of in-situ X-ray scattering methods. Step-flow dynamics in mounded polycrystalline growth has been studied separately for two thin film organic semiconductors deposited by thermal deposition in a vacuum environment, C60 and diindenoperylene (DIP). Highly oriented polycrystalline thin films are readily obtained in both systems, where mounds are composed of crystalline monolayer-height steps and terraces in a so-called wedding cake morphology. The formation of mounds is understood to be due to significant Ehrlich-Schwoebel step edge barriers that inhibit molecules from hopping down from one layer to the one below. The mounds exhibit local step flow, which can be monitored using coherent X-ray scattering. This is made possible due to heterodyning between scattering from the average mounds and the moving steps, which becomes visible in XPCS analysis. The effect of desorption, i.e. re-evaporation of deposited molecules is found to be important for understanding these processes. The impact of this work is to enable testing of models of step dynamics, the shape of mounds, and merging of mounds to form continuous thin films. Highly ordered polycrystalline thin films deposited on inexpensive substrates have applications in thin film solar cells and other organic electronic devices. In a separate set of experiments, speckle analysis during self-organized ion beam nanopatterning reveals memory stretching back to the beginning of patterning in the early stages and enables measurement of the velocity of self-organized patterns across surfaces providing the possibility of stringent new tests of the theory of pattern formation and motion.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Experimental analysis of convective drying of paper and board

Conventional multi-cylinder drying of paper and board involves a mixture of conductive drying from steam-heated dryer cylinders and convective drying by the flow of heated air over the surface of the paper web in the pockets. Pocket ventilation is a critical component in assisting heat and mass transfer during the drying process but is the primary contributor towards removing evaporated water from the web. Air temperature, velocity, and humidity are critical parameters involved in the convective drying process. This paper covers an experimental study involving the design and development of a small lab-scale setup for convective drying of various grades of paper and board, monitoring multiple parameters like paper temperature, moisture content, air humidity, temperature, and velocity measured in situ as the drying proceeds with continuous and accurate sampling capabilities for all parameters in the sample and the system. Instantaneous drying rates, heat, and mass transfer coefficients were also deduced for every time step till the paper completely dried. Furthermore, the coefficients obtained were also reported in the form of dimensionless correlations, and the results were compared against traditional correlations used in the modelling of paper drying. Furthermore, this data will be useful in process development, modelling, design, and the paper drying process simulation.

42 ENGINEERING↗

DOE OSTI2025

High-resolution lidar observations of sedimentation-induced size sorting of droplets near a laboratory cloud top

Cloud optical properties and precipitation, which are crucial to weather and climate, are strongly influenced by cloud microphysical properties that are still poorly understood. Here, we develop a high-resolution time-correlated single-photon-counting lidar and apply it to observe cloud microphysical properties at one-centimeter range resolution in a convection chamber under well-controlled conditions. Together with concurrent in-situ measurements and theoretical analysis, our lidar observations indicate that although turbulent mixing tends to homogenize the cloud in the bulk region, entrainment and sedimentation cause inhomogeneities in droplet concentrations near the cloud top. Specifically, the topmost region is directly affected by entrainment, and lidar profiles show clear evidence of entrained air and detrained cloud filament. The transition region below exhibits vertical size sorting of cloud droplets caused by sedimentation. Our results suggest that using a single sedimentation velocity for all cloud droplets, as is done in many atmospheric models, overlooks key physics relevant to the microphysical structure near the cloud top. In conclusion, our conceptual model used to describe these measurements can serve as a step toward improving the current modeling of processes in the cloud top region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Bayesian State-Space Modeling Framework for Understanding and Predicting Golden Eagle Movements Using Telemetry Data

Predicting raptor movements through a wind power plant under given atmospheric and topographical conditions is a crucial first step in the overall goal of quantifying the risk of turbine-related collisions and mortalities. Extracting behavioral traits of golden eagles (Aquila chrysaetos) from telemetry data requires the fusion of noisy and sparse movement data (location, heading, velocity) with a stochastic mathematical representation of the eagles' decision-making processes. In this study, we framed this problem in a Bayesian state-space framework where both observations and decision-making are assumed to be stochastic processes connected through hidden states (mode of flight, intent), and the unknown model parameters are assumed to be random variables that are calibrated using the available telemetry data. This framework allowed for rigorous consideration of underlying uncertainties while allowing for both data and prior biological knowledge to contribute to a probabilistic and predictive agent-based movement model. We implemented and applied the Bayesian framework to understand movement behavior of 23 GPS-tagged golden eagles travelling in the western US for years 2019 and 2020. Our preliminary findings show that the Bayesian state-space framework provides a robust inverse modeling apparatus to decode eagle behavioral characteristics from telemetry data. This study was primarily aimed at verifying and validating the framework with selected golden eagle tracks (both long- and short-ranged), with future research aimed at extending the framework to include multi-mode flight, consideration of atmospheric data and uplift mechanisms, eagle-to-eagle interaction, and eagle-to-turbine interaction.

Bayesian modeling↗

DOE OSTI2023

Measurement of Beam Energy in the Fermilab's Linac taken at the transfer Line

Linac is the first machine in the Accelerator chain at Fermilab where 𝐻− particles are accelerated from 35 keV to 400 MeV and travel to the Booster where they are stripped. of the extra electrons to become protons. The Booster requires stable beam energy at injection. The Linac extraction energy tuning has been performed using two Beam Velocity Monitors (a.k.a Griffin detectors) to maintain stable relative energy between the Booster and the Linac. Our goal is to control the Linac energy parameters using Machine Learning techniques to increase reliability and quality of beam delivered from Linac. The first step is to understand the data from the Linac extraction diagnostics; and develop reliable, accurate methods for energy measurement and control before implementing machine learning. There were two methods of energy measurement implemented, and the results of these are compared. The first method was time of flight measurements using the Phase probes inside B eam Position Monitors (BPM) that provide the beam’s phase of arrival data. The second method used the relation between beam position data from the BPMs and dispersion values from the 400MeV MAD-X model to calculate momentum variation. The results of these two energy measurement methods compare well as the Linac energy is varied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Dark Matter Velocity Distributions for Direct Detection: Astrophysical Uncertainties Are Smaller Than They Appear

The sensitivity of direct detection experiments depends on the phase-space distribution of dark matter near the Sun, which can be modeled theoretically using cosmological hydrodynamical simulations of Milky Way–like galaxies. However, capturing the halo-to-halo variation in the local dark matter speeds—a necessary step for quantifying the astrophysical uncertainties that feed into experimental results—requires a sufficiently large sample of simulated galaxies, which has been a challenge. In this Letter, we quantify this variation with nearly 100 Milky Way–like galaxies from the tng50 simulation, the largest sample to date at this resolution. Moreover, we introduce a novel phase-space scaling procedure that endows every system with a reference frame that accurately reproduces the local standard-of-rest speed of our Galaxy, providing a principled way of extrapolating the simulation results to real-world data. The ensemble of predicted speed distributions is well characterized by the standard halo model, a Maxwell-Boltzmann distribution truncated at the escape speed, though the individual distributions can deviate from it, especially at high speeds. The dark matter–nucleon cross section limits placed by these speed distributions vary by ∼ 60% about the median. This places the 1⁢𝜎 astrophysical uncertainty at or below the level of the systematic uncertainty of current ton-scale detectors, even down to the energy threshold. The predicted uncertainty remains unchanged when subselecting on those TNG 50 galaxies with merger histories similar to the Milky Way. Tabulated speed distributions, as well as Maxwell-Boltzmann fits, are provided for use in computing direct detection bounds or projecting sensitivities.

Milky Way↗

DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗