SEARCH · Engineering Papers

Results for “time resolution data”

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 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗

DOE OSTI2021

Optimal renormalization of multiscale systems

While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.

Price, Jacob↗

DOE OSTI2023

Measurement of the cross section of the Q=4.4398 MeV 12 C (n, n'γ) reaction from threshold to 16.5 MeV using γ and correlated n–γ detection

The Q=4.4398 MeV 12 C (n, n'γ) cross section was measured using a white incident neutron source through the detection of γ rays only and n–γ coincidences using a segmented liquid scintillator detector array. While the n–γ technique utilized here is more generally applicable to a wide variety of neutron scattering measurements, the γ-only technique was successfully applied to this reaction to exploit the precise time resolution and high efficiency of this detection system to yield results with unprecedented statistical precision and total uncertainties < 2 % from reaction threshold up to 16 MeV incident neutron energy, clearly resolving many features in this reaction that were previously not well known. The γ-only and n–γ results are consistent with each other for the majority of the incident energy range covered in this paper, thereby lending validation to the n–γ technique for future measurements, though significant disagreements are observed between both results and with the ENDF/B-VIII.0 nuclear data evaluation. These differences are particularly noticeable in the recently evaluated energy range below 6.5 MeV, and also near 14 MeV where a “sawtooth”-like feature is observed similar to that in other 12 C+n reaction channels. Finally, both γ-only and n–γ results are presented here with thorough covariance derivations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Shape evolution in neutron-rich odd-even 105–109 Nb isotopes

Background: Neutron-rich nuclei around 𝑍 ≈ 40 are well known for exhibiting multiple shape transitions. Here, this region shows one of the sharpest shape transitions in the nuclear chart, evolving from a spherical vibrator at 𝑁 = 58 to a strongly deformed prolate shape at 𝑁 = 60. The largest deformations are observed for 38 Sr and 40 Zr . This abrupt shape transition disappears at 𝑍 = 36 and below, where a shape transition from spherical to oblate nuclei is predicted. On the other hand, for 𝑍 ≥ 42 and 𝑁 ≥ 60, the shape is known to evolve from axial to triaxial. While the even-𝑍 nuclei in this region have already been extensively studied, new insights can be gained from the studies of odd-𝑍 isotopes for a better understanding of the underlying mechanisms driving these phenomena. Purpose: The 41 Nb nuclei lie at the boundary between axially deformed Zr and triaxially deformed Mo nuclei. This work investigates the nuclear structure of very neutron-rich Nb nuclei up to 𝑁 = 68. The goal is to understand how the nuclear shape evolves as a function of isospin in this isotopic chain and provide new insights into the emergence of triaxial deformation. Methods: The structure of the neutron-rich Nb isotopes was investigated using state-of-the-art high-resolution 𝛾-ray spectroscopy of fission fragments produced via two different fission reactions. The use of 9 Be ⁢( 238 U, 𝑓) inverse kinematics, with a detection system comprising AGATA, EXOGAM, and VAMOS++, enabled the measurement of prompt and delayed 𝛾 rays from isotopically identified fission fragments, and 𝛾−𝛾−𝛾−𝛾 highfold data were obtained from a spontaneous fission source of 252 Cf using the Gammasphere array. Results: The level scheme of 105 Nb has been significantly extended, with the addition of two negative-parity bands observed for the first time. A new level scheme is proposed for 107 Nb, which is not in agreement with an earlier measurement, and new levels and transitions have been added to the very neutron-rich 109 Nb. The degree of triaxiality of the new bands is discussed on the basis of signature splitting analysis. The recently reported level scheme of 99 Nb has been revised. Conclusions: This systematic study on the Nb isotopic chain, compared to Zr and Mo, indicates that while the ground-state band exhibits a triaxial deformation, attributed to a proton hole coupled to a triaxially deformed Mo core, the negative-parity bands, based on isomeric bandheads, display an axially symmetric deformed structure, similar to that observed in the Zr isotopes, revealing the existence of a shape coexistence in the neutron-rich Nb nuclei.

Abushawish, M. [Université Claude-Bernard Lyon 1 (↗

DOE OSTI2021

Feature Analysis, Tracking, and Data Reduction: An Application to Multiphase Reactor Simulation MFiX-Exa for In-Situ Use Case

As we enter the exascale computing regime, powerful supercomputers continue to produce much higher amounts of data than what can be stored for offline data processing. To utilize such high compute capabilities on these machines, much of the data processing needs to happen in situ, when the full high-resolution data is available at the supercomputer memory. In this article, we discuss our MFiX-Exa simulation, which models multiphase flow by tracking a very large number of particles through the simulation domain. In one of the use cases, the carbon particles interact with air to produce carbon dioxide bubbles from the reactor. These bubbles are of primary interest to the domain experts for these simulations. For this particle-based simulation, we propose a streaming technique that can be deployed in situ to efficiently identify the bubbles, track them over time, and use them to down-sample the data with minimal loss in these features.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Are Atmospheric Models Too Cold in the Mountains? The State of Science and Insights from the SAIL Field Campaign

Mountains play an outsized role in water resource availability, and the amount and timing of water they provide depend strongly on temperature. To that end, we ask the question: How well are atmospheric models capturing mountain temperatures? We synthesize results showing that high-resolution, regionally relevant climate models produce 2-m air temperature (T2m) measurements colder than what is observed (a “cold bias”), particularly in snow-covered midlatitude mountain ranges during winter. We find common cold biases in 44 studies across global mountain ranges, including single-model and multimodel ensembles. We explore the factors driving these biases and examine the physical mechanisms, data limitations, and observational uncertainties behind T2m. Our analysis suggests that the biases are genuine and not due to observation sparsity or resolution mismatches. Cold biases occur primarily on mountain peaks and ridges, whereas valleys are often warm biased. Our literature review suggests that increasing model resolution does not clearly mitigate the bias. By analyzing data from the Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign in the Colorado Rocky Mountains, we test various hypotheses related to cold biases and find that local wind circulations, longwave (LW) radiation, and surface-layer parameterizations contribute to the T2m biases in this particular location. We conclude by emphasizing the value of coordinated model evaluation and development efforts in heavily instrumented mountain locations for addressing the root cause(s) of T2m biases and improving predictive understanding of mountain climates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Are Atmospheric Models Too Cold in the Mountains? The State of Science and Insights from the SAIL Field Campaign

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Further improvement and evaluation of nudging in the E3SM Atmosphere Model version 1 (EAMv1): simulations of the mean climate, weather events, and anthropogenic aerosol effects

Abstract. A previous study on the use of nudging in E3SM Atmosphere Model version 1 (EAMv1) had an unresolved issue; i.e., a simulation nudged to EAMv1's own meteorology showed non-negligible deviations from the free-running baseline simulation over some of the subtropical marine stratocumulus and trade cumulus regions. Here, we demonstrate that the deviations can be substantially reduced by (1) changing where the nudging tendency is calculated in the time integration loop of a nudged EAM simulation so as to improve consistency with the free-running baseline and (2) increasing the frequency of the constraining data so as to better capture strong sub-diurnal variations. The fact that modification (2) improves the climate representativeness of the nudged simulations has motivated us to investigate whether the use of newer reanalysis products with higher data frequency can help improve nudged hindcast simulations by better capturing the observed weather events. To answer this question, we present simulations conducted at EAMv1's standard horizontal resolution (approximately 1∘) with nudging towards 6-hourly ERA-Interim reanalysis or 6-hourly, 3-hourly, or hourly ERA5 reanalysis. These simulations are evaluated against the climatology of free-running EAMv1 simulations as well as reanalyses, satellite retrievals, and in situ measurements from the Atmospheric Radiation Measurement user facility. For the 1∘ EAMv1 simulations, we recommend using the relocated nudging tendency calculation and the ERA5 reanalysis at 3-hourly or higher frequency. Simulations used for estimating the anthropogenic aerosol effects often use nudging to help discern signal from noise. The sensitivity of such estimates to the configuration of nudging is investigated in EAMv1, again using the standard 1∘ horizontal resolution. We find that, when estimating the global mean effects, the frequency of constraining data has relatively small impacts, while the choice of nudged variables can change the results substantially. The nudging of air temperature (in addition to horizontal winds) has two non-negligible effects. First, when the constraining data come from reanalysis, the nudging-induced mean bias correction can cause significant changes in the simulated clouds and hence substantially different estimates of the aerosol effects. The impact of the mean bias correction on ice cloud formation has been noted in previous studies and is also seen in EAMv1. For applications like ours, where the preferred configurations of nudging are those capable of providing results consistent with the multi-year free-running simulations, the consequence of the mean bias correction is undesirable. The second important impact of temperature nudging is a significant suppression of adjustments to aerosol forcing, which also causes changes in the estimated aerosol effects. This effect can be seen in simulations nudged to either reanalysis or EAM's own meteorology. These results suggest that nudging horizontal winds but not temperature is a better choice for estimating the anthropogenic aerosol effects.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

915rwpprecipcor

Radar Wind Profiler (RWP) has operational modes: short duration high resolution (referred to as 'high' or 'high mode') and long duration low resolution (referred to as 'low' or 'low mode'). The motivation for this is to create a VAP that quality controls (noise, sig echo filtering), merges, co-grids, and performs initial analyses on these RWP modes. This is done to streamline downstream VAP development on the RWP, as well as enable incorporation of echo properties into the ARSCL chain. The proposed VAP will use the RWPPRECIP b0/b1 level data stream as input. Additional consideration for the incorporation of other ARM data streams to improve this VAP was performed, but not included at this time.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Machine learning–aided real-time detection of keyhole pore generation in laser powder bed fusion

Porosity defects are currently a major factor that hinders the widespread adoption of laser-based metal additive manufacturing technologies. One common porosity occurs when an unstable vapor depression zone (keyhole) forms because of excess laser energy input. With simultaneous high-speed synchrotron x-ray imaging and thermal imaging, coupled with multiphysics simulations, we discovered two types of keyhole oscillation in laser powder bed fusion of Ti-6Al-4V. Amplifying this understanding with machine learning, we developed an approach for detecting the stochastic keyhole porosity generation events with submillisecond temporal resolution and near-perfect prediction rate. Finally, the highly accurate data labeling enabled by operando x-ray imaging allowed us to demonstrate a facile and practical way to adopt our approach in commercial systems.

42 ENGINEERING↗

DOE OSTI2025

Direct reactions with the AT-TPC

Direct reactions are crucial tools for accessing properties of the atomic nucleus. Fundamental and exotic phenomena such as collective modes, pairing, weakbinding effects and evolution of single-particles energies can be investigated in peripheral collisions between a heavy nucleus and a light target. The necessity of using inverse kinematics to reveal how these structural properties change with isospin imbalance renders direct reactions a challenging technique when using the missing mass method. In this scenario, Active Target Time Projection Chambers (AT-TPC) have demonstrated an outstanding performance in enabling these types of reactions even under conditions of very low beam intensities. The AT-TPC of the Facility for Rare Isotope Beams (FRIB) is a next generation multipurpose Active Target. When operated inside a solenoidal magnet, direct reactions benefit from the measurement of the magnetic rigidity that enables particle identification and the determination of the excitation energy with high resolution without the need of auxiliary detectors. Additionally, the AT-TPC can be coupled to a magnetic spectrometer improving even further its spectroscopic investigation capability. In this contribution, we discuss inelastic scattering and transfer reaction data obtained via the AT-TPC and compare them to theory. In particular, we present the results for the 14 C(p,p′) and 12 Be (p,d) 11 Be reactions. For 14 C, we compare the experimental excitation energy of the first 1 – excited state with coupled-cluster calculationsbased on nuclear interactions from chiral effective field theory and with available shell-model predictions. For 12 Be, we determine the theoretical spectroscopic factors of the 12 Be (p,d) 11 Be transfer reaction in the shell modeland compare them to the experimental excitation spectrum from a qualitative standpoint.

active target↗

DOE OSTI2020

Time-Sequenced Flow Field Prediction in an Optima Spark-Ignition Direct-Injection Engine Using Bidirectional Recurrent Neural Network (bi-RNN) with Long Short-Term Memory

To further improve the energy conversion efficiency of internal combustion engine, the transient and complex air flow movement inside the cylinder needs to be better understood and controlled. Although the in-cylinder flow fields are highly stochastic with strong cycle-to-cycle fluctuations, machine learning can still provide an efficient way to learn and regress the complex flow movement process inside the cylinder. In this work, a bidirectional recurrent neural network (bi-RNN) model with long short-term memory was applied to predict the in-cylinder flow fields at different time steps using training data from mull-cycle particle image velocimetry (PIV) measurements. To evaluate the agreement between the true and predicted flow fields, structure and magnitude comparison indices are calculated both globally and locally. The comparison results show that the bi-RNN model can accurately predict the bulk flow and vortex motions from early intake stroke to compression stroke. This work demonstrates that the machine learning model has the potential to predict the underlying dynamics of the interaction between in-cylinder flows and provides a reliable way to improve temporal resolution in PIV flow data to better reveal transient in-cylinder flow features.

Bi-RNN model↗

DOE OSTI2023

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗

DOE OSTI2023

Grand-scale Atmospheric Imaging Apparatus (GAIA) and Wind LiDAR Multi-scale Measurements in the Atmospheric Surface Layer

Understanding the organization and dynamics of turbulence structures in the atmospheric surface layer (ASL) is important for fundamental and applied research in different fields, including weather prediction, snow settling, particle and pollutant transport, and wind energy. The main challenges associated with probing and modeling turbulence in the ASL are: i) the broad range of turbulent scales associated with the different eddies present in high Reynolds-number boundary layers ranging from the viscous scale ($\mathcal{O}$(mm)) up to large energy-containing structures ($\mathcal{O}$(km)); ii) the non-stationarity of the wind conditions and the variability associated with the daily cycle of the atmospheric stability; iii) the interactions among eddies of different sizes populating different layers of the ASL, which contribute to momentum, energy, and scalar turbulent fluxes. Creative and innovative measurement techniques are required to probe near-surface turbulence by generating spatio-temporally-resolved data in the proximity of the ground and, at the same time, covering the entire ASL height with large enough streamwise extent to characterize the dynamics of larger eddies evolving aloft. To this aim, the U.S. National Science Foundation sponsored the development of the Grand-scale Atmospheric Imaging Apparatus (GAIA) enabling super-large snow particle image velocimetry (SLPIV) in the near-surface region of the ASL. This inaugural version of GAIA provides a comprehensive measuring system by coupling SLPIV and two scanning Doppler LiDARs to probe the ASL at an unprecedented resolution. A field campaign performed in 2021–2022 and its preliminary results are presented herein elucidating new research opportunities enabled by the GAIA measuring system.

17 WIND ENERGY↗

DOE OSTI2023

Considerations for the temperature stratification in a pre-burn constant-volume combustion chamber

In recent years, the Engine Combustion Network (ECN) has developed as a worldwide reference for understanding and describing engine combustion processes, successfully bringing together experimental and numerical efforts. Since experiments and numerical simulations both target the same boundary conditions, an accurate characterization of the stratified environment that is inevitably present in experimental facilities is required. The difference between the core-, and pressure-derived bulk-temperature of pre-burn combustion vessels has been addressed in various previous publications. Additionally, thermocouple measurements have provided initial data on the boundary layer close to the injector nozzle, showing a transition to reduced ambient temperatures. The conditions at the start of fuel injection influence physicochemical properties of a fuel spray, including near nozzle mixing, heat release computations, and combustion parameters. To address the temperature stratification in more detail, thermocouple measurements at larger distances from the spray axis have been conducted. Both the temperature field prior to the pre-combustion event that preconditions the high-temperature, high-pressure ambient, as well as the stratification at the moment of fuel injection were studied. To reveal the cold boundary layer near the injector with a better spatial resolution, Rayleigh scattering experiments and thermocouple measurements at various distances close to the nozzle have been carried out. The impact of the boundary layers and temperature stratification are illustrated and quantified using numerical simulations at Spray A conditions. Next to a reference simulation with a uniform temperature field, six different stratified temperature distributions have been generated. These distributions were based on the mean experimental temperature superimposed by a randomized variance, again derived from the experiments. The results showed that an asymmetric flame structure arises in the computed results when the temperature stratification input is used. In these predictions, first-stage ignition is advanced by 24 $\mu$s , while second-stage ignition is delayed by 11 $\mu$s. At the same time a lift-off length difference between the top and the bottom of up to 1.1 mm is observed. Furthermore, the lift-off length is less stable over time. Given the shown dependency, the temperature data is made available along with the vessel geometry data as a recommended basis for future numerical simulations.

42 ENGINEERING↗

DOE OSTI2021

Single Molecule Spectroscopy Studies of Acid–Base Chemical Gradients Using Nile Red as a Probe of Local Surface Acidity

Single molecule spectroscopy studies of local acidity along bifunctional acid-base gradients are reported. Gradients are prepared by directional vapor phase diffusion and subsequent reaction of 3-aminopropyl-trimethoxysilane with a uniform silica film. Gradient formation is confirmed by spectroscopic ellipsometry and by static water contact angle measurements. X-ray photoelectron spectroscopy is used to characterize the nitrogen content and degree of nitrogen protonation along the gradient. Nile Red is employed as the probe dye in single molecule spectroscopy studies of these gradients. While Nile Red is well-known for its solvent sensitivity, it is used here, for the first time, to sense the acid/base properties of the film in two-color wide-field fluorescence imaging experiments. The data reveal broad bimodal distributions of Nile Red emission spectra that vary along the gradient direction. The single-molecule results are consistent with solution-phase ensemble acid/base studies of the dye. The former reveal a gradual transition from a surface dominated by basic aminosilane sites at the high-amine end of the gradient to one dominated by acidic silanol sites at the low-amine end. The sub-diffraction-limited spatial resolution afforded by superlocalization of the single molecules reveals spatial correlations in the acid/base properties of the gradient over ~ 200 nm distances. Furthermore, these studies provide data relevant to the use of aminosilane-modified silica in bifunctional, cooperative chemical catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Inertial Effects During the Process of Supercritical CO 2 Displacing Brine in a Sandstone: Lattice Boltzmann Simulations Based on the Continuum-Surface-Force and Geometrical Wetting Models

Inertial effects during the process of supercritical CO 2 displacing brine in porous media may not be negligible according to recent studies. Capturing the inertial effects of the physical CO 2 -brine system imposes a requirement on the grid resolution and viscosity to surface tension ratio for pore-scale simulations, which some commonly used simulators may not be able to meet. To fulfill the parameter requirement, we combine the continuum-surface-force based color-gradient lattice Boltzmann (LB) multiphase model and the geometrical wetting model and extend the model to 3D under the multiple-relaxation-time framework. We validate the model via simple benchmarks which show significant improvement over the traditional models. Furthermore, we perform 3D drainage simulations in a heterogeneous micromodel where the simulation result agrees well with experimental data, while our previous work fails to reproduce certain displacement patterns in the experiment due to the use of a traditional LB model that cannot fulfill the parameter requirement. Finally, we perform high-fidelity 3D drainage simulations to study the inertial effects in a Bentheimer sandstone sample. These findings show that stronger inertial effects generally help develop more CO 2 flow pathways for the same capillary number which results in higher CO 2 saturation, consistent with the micromodel results. The phenomena can be found in both low and high capillary number cases, indicating that the inertial effects are not dependent on the mean velocity. In addition, the change of the invasion patterns is not proportional to the change of inertial effects, thus exhibiting threshold behavior.

58 GEOSCIENCES↗