Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Houston”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Determining Neutrino Mass Ordering with NOvA and Upcoming JUNO Measurements

NOvA has reported a significance of mass ordering determination using ten years of data together with external constraints from reactor-based experiments. The JUNO collaboration is poised to provide a more precise reactor-based constraint on $|Δm^2_{32}|$. In this Letter, we explore the potential impact of this anticipated measurement on the determination of the neutrino mass ordering by NOvA. We find that $3σ$ evidence of the normal ordering is achievable over a range of plausible JUNO measurements within the next five years.

Abubakar, S. [Erciyes U.]↗
DOE OSTI2025

Ca X ML: Chemistry‐informed machine learning explains mutual changes between protein conformations and calcium ions in calcium‐binding proteins using structural and topological features

Proteins' flexibility is a feature in communicating changes in cell signaling instigated by binding with secondary messengers, such as calcium ions, associated with the coordination of muscle contraction, neurotransmitter release, and gene expression. When binding with the disordered parts of a protein, calcium ions must balance their charge states with the shape of calcium-binding proteins and their versatile pool of partners depending on the circumstances they transmit. Accurately determining the ionic charges of those ions is essential for understanding their role in such processes. However, it is unclear whether the limited experimental data available can be effectively used to train models to accurately predict the charges of calcium-binding protein variants. Here, we developed a chemistry-informed, machine-learning algorithm that implements a game theoretic approach to explain the output of a machine-learning model without the prerequisite of an excessively large database for high-performance prediction of atomic charges. We used the ab initio electronic structure data representing calcium ions and the structures of the disordered segments of calcium-binding peptides with surrounding water molecules to train several explainable models. Network theory was used to extract the topological features of atomic interactions in the structurally complex data dictated by the coordination chemistry of a calcium ion, a potent indicator of its charge state in protein. Our design created a computational tool of Ca X ML, which provided a framework of explainable machine learning model to annotate ionic charges of calcium ions in calcium-binding proteins in response to the chemical changes in an environment. Our framework will provide new insights into protein design for engineering functionality based on the limited size of scientific data in a genome space.

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

First observation of strange baryon enhancement with effective energy in pp collisions at the LHC

The production of (multi-)strange hadrons is measured at midrapidity in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the local charged-particle multiplicity in the pseudorapidity interval |η| < 0.5 and of the very-forward energy measured by the ALICE Zero-Degree Calorimeters. The latter provides information on the effective energy, i.e. the energy available for particle production in the collision once subtracted from the centre-of-mass energy. The yields of K$^{0}_{s}$, $Λ$ + $\overline{Λ}$, and $Ξ$ – + $\overline{Ξ}$ + per charged-particle increase with the effective energy. In addition, this work exploits a multi-differential approach to decouple the roles of local multiplicity and effective energy in such an enhancement. The results presented in this article provide new insights into the interplay between global properties of the collision, such as the initial available energy in the event, and the locally produced final hadronic state, connected to the charged-particle multiplicity at midrapidity. Notably, a strong increase of strange baryon production with effective energy is observed for fixed charged-particle multiplicity at midrapidity. These results are discussed within the context of existing phenomenological models of hadronisation implemented in different tunes of the PYTHIA 8 event generator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multiplicity-dependent jet modification from di-hadron correlations in pp collisions at $ \sqrt{s} $ = 13 TeV

Short-range correlations between charged particles are studied via two-particle angular correlations in pp collisions at $\sqrt{s}$ = 13 TeV. The correlation functions are measured as a function of the relative azimuthal angle ∆φ and the pseudorapidity separation ∆η for pairs of primary charged particles within the pseudorapidity interval |η| < 0.9 and the transverse-momentum range 1 < p T < 8 GeV/c. Near-side (|∆φ| < 1.3) peak widths are extracted from a generalised Gaussian fitted over the correlations in full pseudorapidity separation (|∆η| < 1.8), while the per-trigger associated near-side yields are extracted for the short-range correlations (|∆η| < 1.3). Both are evaluated as a function of charged-particle multiplicity obtained by two different event activity estimators. The width of the near-side peak decreases with increasing multiplicity, and this trend is reproduced qualitatively by the Monte Carlo event generators PYTHIA 8, AMPT, and EPOS. However, the models overestimate the width in the low transverse-momentum region (p T < 3 GeV/c). The per-trigger associated near-side yield increases with increasing multiplicity. Although this trend is also captured qualitatively by the considered event generators, the yield is mostly overestimated by the models in the considered kinematic range. The measurement of the shape and yield of the short-range correlation peak can help us understand the interplay between jet fragmentation and event activity, quantify the narrowing trend of the near-side peak as a function of transverse momentum and multiplicity selections in pp collisions, and search for final-state jet modification in small collision systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Measurement of ω meson production in pp collisions at $\sqrt{s}$ = 13 TeV

The p T -differential cross section of ω meson production in pp collisions at $\sqrt{s}$ =13 TeV at midrapidity (|y| < 0.5) was measured with the ALICE detector at the LHC, covering an unprecedented transverse-momentum range of 1.6 < p T < 50 GeV/c. The meson is reconstructed via the ω → π + π − π 0 decay channel. The results are compared with various theoretical calculations: PYTHIA8.2 with the Monash 2013 tune overestimates the data by up to 50%, whereas good agreement is observed with Next-to-Leading Order (NLO) calculations incorporating ω fragmentation using a broken SU(3) model. The ω/π 0 ratio is presented and compared with theoretical calculations and the available measurements at lower collision energies. The presented data triples the p T ranges of previously available measurements. A constant ratio of C ω/π0 = 0.578 ± 0.006 (stat.) ± 0.013 (syst.) is found above a transverse momentum of 4 GeV/c, which is in agreement with previous findings at lower collision energies within the systematic and statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Space-time evolution of particle emission in p–Pb collisions at $\sqrt{{s}_{\text{NN}}} = 5.02$ TeV with 3D kaon femtoscopy

The measurement of three-dimensional femtoscopic correlations between identical charged kaons (K ± K ± ) produced in p–Pb collisions at center-of-mass energy per nucleon pair $\sqrt{{s}_{\text{NN}}}=5.02$ TeV with ALICE at the LHC is presented for the first time. This measurement, supplementary to those in pp and Pb–Pb collisions, allows understanding the particle-production mechanisms at different charged-particle multiplicities and provides information on the dynamics of the source of particles created in p–Pb collisions, for which a general consensus does not yet exist. It is shown that the measured source sizes increase with charged-particle multiplicity and decrease with increasing pair transverse momentum. These trends for K ± K ± are similar to the ones observed earlier in identical charged-pion and ${\text{K}}_{\text{s}}^{0}{\text{K}}_{\text{s}}^{0}$ correlations in Pb–Pb collisions at various energies and in π ± π ± correlations in p–Pb collisions at $\sqrt{{s}_{\text{NN}}}=5.02$ TeV. At comparable multiplicity, the source sizes measured in p–Pb collisions agree within uncertainties with those observed in pp collisions, and there is an indication that they are smaller than those observed in Pb–Pb collisions. The obtained results are also compared with predictions from the hadronic interaction model EPOS 3, which tends to underestimate the source size for the most central collisions and agrees with the data for semicentral and peripheral events. Furthermore, the time of maximal emission for kaons is extracted. It turns out to be comparable with the value obtained in highly peripheral Pb–Pb collisions at the same energy, indicating that the kaon emission evolution is similar to that in p–Pb collisions.

Heavy Ion Experiments↗
DOE OSTI2025

J/ ψ -hadron correlations at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV

We report on the measurement of inclusive, non-prompt, and prompt J/ψ-hadron correlations by the ALICE Collaboration at the CERN Large Hadron Collider in pp collisions at a center-of-mass energy of 13 TeV. The correlations are studied at midrapidity (|y| < 0.9) in the transverse momentum ranges p T < 40 GeV/c for the J/ψ and 0.15 < p T < 10 GeV/c and |η| < 0.9 for the associated hadrons. The measurement is based on minimum bias and high multiplicity data samples corresponding to integrated luminosities of L int = 34 nb −1 and L int = 6.9 pb −1 , respectively. In addition, two more data samples are employed, requiring, on top of the minimum bias condition, a threshold on the tower energy of E = 4 and 9 GeV in the ALICE electromagnetic calorimeters, which correspond to integrated luminosities of L int = 0.9 pb −1 and L int = 8.4 pb −1 , respectively. The azimuthally integrated near and away side yields of associated charged hadrons per J/ψ trigger are presented as a function of the J/ψ and associated hadron transverse momentum. The measurements are discussed in comparison to PYTHIA calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multiplicity-dependent inclusive J/ ψ production at forward rapidity in pp collisions at $\sqrt{s}$ = 13 TeV

This paper presents a study of the inclusive forward J/ψ yield as a function of forward charged-particle multiplicity in pp collisions at $\sqrt{s}$ = 13 TeV using data collected by the ALICE experiment at the CERN LHC. The results are presented in terms of relative J/ψ yields and relative charged-particle multiplicities with respect to these quantities obtained in inelastic collisions having at least one charged particle in the pseudorapidity range |η| < 1. The J/ψ mesons are reconstructed via their decay into μ + μ − pairs in the forward rapidity region (2.5 < y < 4). The relative multiplicity is estimated in the forward pseudorapidity range which overlaps with the J/ψ rapidity region. The results show a steeper-than-linear increase of the J/ψ yields versus the multiplicity. They are compared with previous measurements and theoretical model calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Light neutral-meson production in pp collisions at $\sqrt{\text{s}}$ = 13 TeV

The momentum-differential invariant cross sections of π 0 and η mesons are reported for pp collisions at $\sqrt{s}$ = 13 TeV at midrapidity (|y| < 0.8). The measurement is performed in a broad transverse-momentum range of 0.2 < p T < 200 GeV/c and 0.4 < p T < 60 GeV/c for the π 0 and η, respectively, extending the p T coverage of previous measurements. Transverse-mass-scaling violation of up to 60% at low transverse momentum has been observed, agreeing with measurements at lower collision energies. Transverse Bjorken x (x T ) scaling of the π 0 cross sections at LHC energies is fulfilled with a power-law exponent of n = 5.01 ± 0.05, consistent with values obtained for charged pions at similar collision energies. The data are compared to predictions from next-to-leading order perturbative QCD calculations, where the π 0 spectrum is best described using the CT18 parton distribution function and the NNFF1.0 or BDSS fragmentation function. Expectations from PYTHIA8 and EPOS LHC overestimate the spectrum for the π 0 and are not able to describe the shape and magnitude of the η spectrum. The charged-particle multiplicity dependent π 0 and η p T spectra show the expected change of the spectral shape, characterized by a flatter slope with increasing multiplicity. This is demonstrated across a broad transverse-momentum range and up to events with a charged-particle multiplicity exceeding five times the mean value in minimum bias collisions. The η/π 0 ratio depends on the charged-particle multiplicity for p T < 4 GeV/c. PYTHIA8 and EPOS LHC qualitatively explain this behavior with an increasing contribution from the feed-down of heavier particles to the π 0 spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of correlations among net-charge, net-proton, and net-kaon multiplicity distributions in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations among conserved quantum numbers, such as the net-electric charge, the net-baryon, and the net-strangeness in heavy-ion collisions, are crucial for exploring the QCD phase diagram. In this paper, these correlations are investigated using net-proton number (as a proxy for the net-baryon), net-kaon number (for the net-strangeness), and net-charged particle number in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ALICE detector. The observed correlations deviate from the Poissonian baseline, with a more pronounced deviation at LHC energies than at RHIC. Theoretical calculations of the Thermal-FIST hadron resonance gas model, HIJING, and EPOS LHC event generators are compared with experimental results, where a significant impact of resonance decays is observed. Thermal-FIST calculations under the grand canonical and canonical ensembles highlight significant differences, underscoring the role of local charge conservation in explaining the data. Recent lattice QCD studies have demonstrated that the magnetic field generated by spectator protons in heavy-ion collisions affects susceptibility ratios, in particular those related to the net-electric charge and the net-baryon numbers. The experimental findings are in qualitative agreement with the expectations of lattice QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

First measurement of D* + vector meson spin alignment in Pb–Pb collisions at $\sqrt{{s}_{\text{NN}}}={5}.0{2}$ TeV

The first measurement of prompt D *+ -meson spin alignment in ultrarelativistic heavy-ion collisions with respect to the direction orthogonal to the reaction plane is presented. The spin alignment is quantified by measuring the element ρ 00 of the diagonal spin-density matrix for prompt D *+ mesons with 4 < p T < 30 GeV/c in two rapidity intervals, |y| < 0.3 and 0.3 < |y| < 0.8, in central (0–10%) and midcentral (30–50%) Pb–Pb collisions at $\sqrt{{s}_{\text{NN}}}={5}.0{2}$ TeV. Evidence of spin alignment ρ 00 > 1/3 has been found for p T > 15 GeV/c and 0.3 < |y| < 0.8 with a significance of 3.1σ. The measured spin alignment of prompt D *+ mesons is compared with the one of inclusive J/ψ mesons measured at forward rapidity (2.5 < y < 4).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Study of $\langle {p}_{\text{T}}\rangle$ and its higher moments, and extraction of the speed of sound in Pb-Pb collisions with ALICE

Ultrarelativistic heavy-ion collisions produce a state of hot and dense strongly interacting QCD matter called quark-gluon plasma (QGP). On an event-by-event basis, the volume of the QGP in ultracentral collisions is mostly constant, while its total entropy can vary significantly due to quantum fluctuations, leading to variations in the temperature of the system. Exploiting this unique feature of ultracentral collisions allows for the interpretation of the correlation of the mean transverse momentum ($\langle$p T $\rangle$) of produced charged hadrons and the number of charged hadrons as a measure for the speed of sound, c s . This speed is related to the rate at which compression waves travel in the QGP and is determined by fitting the relative increase in $\langle$p T $\rangle$ with respect to the relative change in the average charged-particle density ($\langle$dN ch /dη$\rangle$) measured at mid-rapidity. This study reports the event-average $\langle$p T $\rangle$ of charged particles as well as the variance, skewness, and kurtosis of the event-by-event transverse momentum per charged particle ([p T ]) distribution in ultracentral Pb-Pb collisions at a center-of-mass energy of 5.02 TeV per nucleon pair using the ALICE detector. Different centrality estimators based on charged-particle multiplicity or the transverse energy of the event are used to select ultracentral collisions. By ensuring a pseudorapidity gap between the region used to define the centrality and the region used to perform the measurement, the influence of biases and their potential effects on the rise of the mean transverse momentum is tested. The measured c$^{2}_{s}$ is found to strongly depend on the exploited centrality estimator and ranges between 0.1146±0.0028 (stat.)±0.0065 (syst.) and 0.4374±0.0006 (stat.)±0.0184 (syst.) in natural units. The self-normalized variance shows a steep decrease towards ultracentral collisions, while the self-normalized skewness variables show a maximum, followed by a fast decrease. These non-Gaussian features are understood in terms of the vanishing of the impact-parameter fluctuations contributing to the event-to-event [p T ] distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multiplicity dependence of ${\Xi }_{\text{c}}^{+}$ and ${\Xi }_{\text{c}}^{0}$ production in pp collisions at $\sqrt{s}=13$ TeV

The first measurement at midrapidity (|y| < 0.5) of the production yield of the strange-charm baryons $Ξ$$^{+}_{c}$ and $Ξ$$^{0}_{c}$ as a function of transverse momentum (p T ) in different charged-particle multiplicity classes in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ALICE experiment at the LHC is reported. The $Ξ$$^{+}_{c}$ baryon is reconstructed via the $Ξ$$^{+}_{c}$ → $Ξ$ – π + π + decay channel in the range 4 < p T < 12 GeV/c, while the $Ξ$$^{0}_{c}$ baryon is reconstructed via both the $Ξ$$^{0}_{c}$ → $Ξ$ – π + and $Ξ$$^{0}_{c}$ → $Ξ$ – e + ν e decay channels in the range 2 < p T < 12 GeV/c. The baryon-to-meson ($Ξ$$^{0}_{c}$ + /D 0 ) and the baryon-to-baryon ($Ξ$$^{0}_{c}$ + /Λ$^{+}_{c}$) production yield ratios show no significant dependence on multiplicity. In addition, the observed yield ratios are not described by theoretical predictions that model charm-quark fragmentation based on measurements at e + e − and e − p colliders, indicating differences in the charm-baryon production mechanism in pp collisions. A comparison with different event generators and tunings, including different modelling of the hadronisation process, is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Imaging the initial condition of heavy-ion collisions and nuclear structure across the nuclide chart

High-energy nuclear collisions encompass three key stages: the structure of the colliding nuclei, informed by low-energy nuclear physics, the initial condition , leading to the formation of quark–gluon plasma (QGP), and the hydrodynamic expansion and hadronization of the QGP, leading to final-state hadron distributions that are observed experimentally. Recent advances in both experimental and theoretical methods have ushered in a precision era of heavy-ion collisions, enabling an increasingly accurate understanding of these stages. However, most approaches involve simultaneously determining both QGP properties and initial conditions from a single collision system, creating complexity due to the coupled contributions of these stages to the final-state observables. To avoid this, we propose leveraging established knowledge of low-energy nuclear structures and hydrodynamic observables to independently constrain the QGP’s initial condition. By conducting comparative studies of collisions involving isobar-like nuclei—species with similar mass numbers but different ground-state geometries—we can disentangle the initial condition’s impacts from the QGP properties. This approach not only refines our understanding of the initial stages of the collisions but also turns high-energy nuclear experiments into a precision tool for imaging nuclear structures, offering insights that complement traditional low-energy approaches. Opportunities for carrying out such comparative experiments at the Large Hadron Collider and other facilities could significantly advance both high-energy and low-energy nuclear physics. Additionally, this approach has implications for the future electron-ion collider. While the possibilities are extensive, we focus on selected proposals that could benefit both the high-energy and low-energy nuclear physics communities. Originally prepared as input for the long-range plan of U.S. nuclear physics, this white paper reflects the status as of September 2022, with a brief update on developments since then.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Impact of anion-exchange membrane degradation on long-term CO 2 electrolysis to ethylene

Electrochemical CO 2 reduction (eCO2R) to ethylene offers a unique opportunity to diversify domestic supply chains for chemical manufacturing. Large-scale deployment of eCO2R to ethylene reactors is currently limited by low energy efficiencies and poor durability. Herein, we demonstrate >100 h of continuous electrolysis at an industrially relevant current density of 200 mA cm −2 in a 25 cm 2 geometric area zero-gap reactor with a full-cell voltage <3.1 V and ethylene Faradaic efficiency >30%. Incorporating expanded polytetrafluoroethylene (PTFE)-supported electrodes into zero-gap reactors allows for durable electrode catalysts that are resistant to flooding, and this enables wider ranges of operating parameters not typically accessible in CO 2 electrolysis. Finally, we identify membrane degradation, evidenced by a loss of membrane ion-exchange capacity, as the leading cause of performance loss over 144 h of electrolysis.

Energy - Conversion↗
DOE OSTI2025

MetaHeart: Metasurface enabled biometrics camouflage

Privacy-invading biometrics monitoring is becoming a prominent security threat as modern sensing systems move to higher operating frequencies (mmWave, sub-THz), increasing sensing resolution and accuracy. As such, developing systems that can protect or obfuscate biometrics from adversarial intrusion becomes pivotal to preserving user privacy. In this work, we develop and implement MetaHeart, a real-time biometrics misinformation system based on reflective, programmable metasurfaces and dynamic phase-front manipulation of radar inferences. MetaHeart’s key goal is to prevent the leakage of a legitimate user’s heartbeat biometrics by spoofing fake heartbeat signals at a malicious, radar-equipped, heart rate sensing intruder. Furthermore, we experimentally demonstrate MetaHeart’s ability to fake Alice’s presence when she is not there and to fool Trudy’s inferences even when Alice is present, achieving an overall accuracy above 98%. Finally, we conduct a robustness analysis to determine MetaHeart’s required spatial placement within the intruder’s monitoring area that would allow for effective spoofing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗