Search NASA⌕ Search

SEARCH · Search NASA

Results for “flow (dynamics)”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhancing Streamflow Reanalysis Across the Conterminous US Leveraging Multiple Gridded Precipitation Data Sets

Streamflow observations, essential for various water resource applications, are often unavailable at critical locations in need. Although different models have been proposed to enhance streamflow predictability at ungauged locations, the challenge extends beyond model fidelity. Differences in meteorologic forcing data sets, precipitation in particular, can significantly affect the accuracy of hydrologic predictions. This challenge intensifies across regions characterized by diverse hydro-climatological and geographical conditions, such as in the conterminous US (CONUS) where a single precipitation product struggles to consistently replicate observed hydrographs, particularly peak flow dynamics. To enhance streamflow predictions, we utilize a VIC-RAPID hydrologic modeling framework driven by multiple commonly used meteorological forcing data sets, such as Daymet, PRISM, ST4, AORC, and their hybrids and create multiple sets of 40-year (1980–2019) hourly, daily, and monthly streamflow reanalysis, Dayflow Version 2, for 2.7 million river reaches across the CONUS. Most forcings lead to skillful streamflow performance, except for ST4 in the mountainous west, where severe radar blockage adversely affects the accuracy. The evaluation using over 6,000 hourly stream gauges shows that hourly AORC and ST4 lead to improved annual peak flow performance over Daymet—driven streamflow (Dayflow V1), particularly in smaller basins, highlighting the value of high temporal resolution forcings in hydrologic predictions. Compared with other benchmark data sets like National Water Model V3.0, AORC-driven VIC-RAPID exhibits improved regional streamflow performance, with comparable peak flow representation. We envision that multi-forcing streamflow reanalysis data can inform regions in need of forcing data enhancement, diagnose hydrologic model performance, and benefit diverse water resource applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Numerical simulation of a powered-lift landing, tracking flow features using overset grids, and simulation of high lift devices on a fighter-lift-and-control wing

Attached as appendices to this report are documents describing work performed on the simulation of a landing powered-lift delta wing, the tracking of flow features using overset grids, and the simulation of flaps on the Wright Patterson Lab's fighter-lift-and-control (FLAC) wing. Numerical simulation of a powered-lift landing includes the computation of flow about a delta wing at four fixed heights as well as a simulated landing, in which the delta wing descends toward the ground. Comparison of computed and experimental lift coefficients indicates that the simulations capture the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect. Power spectra of temporal variations of pressure indicate computed vortex shedding frequencies close to the jet exit are in the experimentally observed frequency range; the power spectra of pressure also provide insights into the mechanisms of lift oscillations. Also, a method for using overset grids to track dynamic flow features is described and the method is validated by tracking a moving shock and vortices shed behind a circular cylinder. Finally, Chimera gridding strategies were used to develop pressure coefficient contours for the FLAC wing for a Mach no. of 0.18 and Reynolds no. of 2.5 million.

Chawla, Kalpana↗
NASA NTRS2005

Space Laboratory on a Table Top: A Next Generative ECLSS design and diagnostic tool

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale-time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations. Preliminary experimental and modeling work in this area will be presented. This involves testing of a single inlet-exit model with detailed 3-D flow visualization and quantitative diagnostics and computational modeling of the system.

Ramachandran, N.↗
NASA NTRS1990

Theoretical investigation of non-equilibrium chemistry and optical radiation in hypersonic flow fields

Future space vehicles returning from distant missions or high earth orbits may enter the upper regions of the atmosphere and use aerodynamic drag to reduce their velocity before they skip out of the atmosphere and enter low earth orbit. The Aeroassist Flight Experiment (AFE) is designed to explore the special problems encountered in such entries. A computer code was developed to calculate the radiative transport along line-or-sight in the general 3-D flow field about an arbitrary entry vehicle, if the temperatures and species concentrations along the line-of-sight are known. The radiative heating calculation at the stagnation point of the AFE vehicle along the entry trajectory was performed, including a detailed line-by-line accounting of the radiative transport in the vacuum ultraviolet (below 200 nm) by the atomic N and O lines. A method was developed for making measurements of the haze particles in the Titan atmosphere above 200 km altitude. Several other tasks of a continuing nature, to improve the technical ability to calculate the nonequilibrium gas dynamic flow field and radiative heating of entry vehicles, were completed or advanced.

Whiting, Ellis E.↗
NASA NTRSDate not supplied

Experimental Identification of Bistable Flow States on the Space Launch System at Liftoff Conditions

Time-averaged global force and moment measurements in wind tunnel tests can obscure complex flow dynamics, including bistability phenomena. This paper presents a detailed analysis into the potential cause of previously unobserved discrepancies among repeat runs during the 202114- by 22-Foot Subsonic Tunnel test of the Space Launch System. During the test, it was observed that liftoff configuration repeatability was not as strong as expected when compared to previous tests. These discrepancies occurred between wind angle azimuths of approximately 150◦ and 210◦, as well as between 330◦ and 30◦. After ruling out common confounding factors such as instrumentation hysteresis and temperature, it was hypothesized that the lack of repeatability was caused by the Coandă effect producing a bistable flow state at different locations between the cylindrical rocket centerbody and solid rocket boosters during the acquisition of data. Evidence to support this hypothesis is provided in this paper in the form of time-averaged force and moment coefficient data, time-dependent force and moment coefficient data, and still images from smoke flow and tuft flow visualization runs.

SLS↗
NASA NTRS2018

Validation of Heat Transfer Correlations in Line Chill-Down Tests of Cryogenic Fluid in SINDA/FLUINT

Line chill-down is an important process in cryogenic tank propellant management, storage, and usage Complex flow dynamics during these processes: boiling heat transfer (film, transition, and nucleate) Understanding boiling phenomena can lead to efficient line chill-down systems that use less propellant, propellant stored, reducing cost for space missions Line Chill-down heat transfer was modelled using SINDA/FLUINT version 5.8 (SF) Multiple chill-down tests were modelled using: heat transfer correlations readily available in SF using HTN/HTC TIES heat transfer empiricisms developed by the University of Florida (UF) based on a series of liquid nitrogen chill-down tests using SF HTU TIES Chill-down tests modelled: liquid nitrogen tests conducted by the University of Florida horizontal flow, upward flow, and downward flow (Reynolds Numbers ranging 850-231,000)liquid hydrogen tests conducted by NASA Glenn Research Center vertical upward flow (Reynolds Number range of 18,400 - 433,000)The flow rate was measured far downstream of the test section, near the system exit. Where to set the flow rate? SF was highly sensitive, and sometime unstable, setting the test flow rate downstream (the outlet) of the test section model and setting the test pressure upstream (the inlet) of the test section model higher flow rate oscillations at the entrance of the model's test section SF was more stable setting the test flow rate upstream (than the downstream flow rate set case)test pressure was used as an inlet (SF plenum) to set the thermodynamic state (temperature and quality) coming into the system setting the appropriate downstream pressure was the unknown. The pressure drops predicted by SF for the downstream set flow rate boundary condition were much smaller than test section measured pressure drops. The multiphase pressure drop correlations used internally in SF may need to be adjusted. Models with an upstream flow rate set assumed a pressure drop that was small

cryogenic modeling↗
DOE OSTI2025

Nonseeded linewise temperature measurements by resonantly ionized photoemission thermometry in a Mach 4 Ludwieg tube

A one-dimensional (1D) thermometry using oxygen-tagging resonantly ionized photoelectron thermometry (O 2 RIPT) was employed to investigate thermal gradients within a Mach 4 Ludwieg tube. The Ludwieg tube is pulsed with a test duration of approximately 100 ms, providing a cold supersonic flow at Mach 4 ideal for studying aerothermal effects. This study focused on measuring freestream temperatures, capturing shock-induced heating behind a detached bow shock from a blunt cylinder, and resolving sharp temperature variations across a bow shock generated by a cylinder. The O 2 RIPT technique produced strong emission signals extending approximately 4 cm long, demonstrating its capability for precise temperature measurements in high-speed wind tunnel environments. The results confirm that O 2 RIPT is well-suited for applications in large-scale aerodynamic testing facilities, particularly in regions with strong compression effects, enabling the resolution of sharp thermal gradients. This method presents a promising solution for thermometry in dynamic flow conditions relevant to various experimental ground-test facilities.

McCord, Walker (ORCID:0000000179747550)↗
NASA NTRS1994

Application of static and dynamic enclosures for determining dimethyl sulfide and carbonyl sulfide exchange in Sphagnum peatlands: Implications for the magnitude and direction of flux

A static enclosure method was applied to determine the exchange of dimethyl sulfide (DMS) and carbonyl sulfide (OCS) between the surface of Sphagnum peatlands and the atmosphere. Measurements were performed concurrently with dynamic (flow through) enclosure measurements with sulfur-free air used as sweep gas. This latter technique has been used to acquire the majority of available data on the exchange of S gases between the atmosphere and the continental surfaces and has been criticized because it is thought to overestimate the true flux of gases by disrupting natural S gas gradients. DMS emission rates determined by both methods were not statistically different between 4 and greater than 400 nmol/sq m/h, indicating that previous data on emissions of at least DMS are probably valid. However, the increase in DMS in static enclosures was not linear, indicating the potential for a negative feedback of enlosure DMS concentrations on efflux. The dynamic enclosure method measured positive OCS flux rates (emission) at all sites, while data using static enclosures indicated that OCS was consumed from the atmosphere at these same sites at rates of 3.7 to 55 nmol/sq m/h. Measurements using both enclosure techniques at a site devoid of vegetation showed that peat was a source of both DMS and OCS. However, the rate of OCS efflux from decomposing peat was more than counterbalanced by OCS consumption by vegetation, including Sphagnum mosses, and net OCS uptake occurred at all sites. We propose that all wetlands are net sinks for OCS.

De Mello, William Z.↗
NASA NTRSDate not supplied

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗
DOE OSTI2026

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗
NASA NTRSDate not supplied

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗
NASA NTRS2020

Intercomparison of surface meltwater routing models for the Greenland ice sheet and influence on subglacial effective pressures

Each summer, large volumes of surface meltwater drain off the Greenland ice sheet (GrIS) surface through moulins to the bed, impacting subglacial hydrology and ice flow dynamics. Supraglacial surface routing delays may propagate to englacial and subglacial hydrologic systems, requiring accurate assessment to correctly estimate subglacial effective pressures. We compare hourly supraglacial moulin discharge simulations from three surface meltwater routing models – the synthetic unit hydrograph (SUH), the bare-ice component of surface routing and lake filling (SRLF), and the rescaled width function (RWF) – for four internally drained catchments on the southwestern Greenland ice sheet surface. The routing models are forced identically using surface runoff from the Modèle Atmosphérique Régionale regional climate model (RCM). For each catchment, simulated moulin hydrographs are input to the SHAKTI subglacial hydrologic model to simulate diurnally varying subglacial effective-pressure variations in the vicinity of a single moulin. Overall, all three routing models produce more realistic moulin discharges than simply using RCM runoff outputs without surface routing but produce significant differences in peak moulin discharge and time to peak. In particular, the RWF yields later, smaller peak moulin discharges than the SUH or SRLF due to its representation of slow interfluve flow between supraglacial meltwater channels, and it can readily accommodate the seasonal evolution of supraglacial stream and river networks. Differences among the three routing models are reflected in a series of simple idealized subglacial hydrology simulations that yield different diurnal effective-pressure amplitudes; however, the supraglacial hydrologic system acts as short-term storage for surface meltwater, and the temporal mean effective pressure is relatively consistent across routing models.

Kang Yang↗
DOE OSTI2024

Epoxy Coating of Biofiber: An Effective Modifier of Biofiber Physical and Flow Properties for Improved Tensile Behavior of Biofiber-Reinforced Biocomposite

Biocomposites combine renewable, plant-based fibers with degradable polymers and are an attractive option for sustainable, lightweight, and cost-effective materials with a low carbon footprint, especially for large-scale additive manufacturing. One of the major challenges in the widespread adoption of biocomposites is that their mechanical performance is significantly inferior to that of synthetic composites. Surface treatment is a common and effective technique to improve the mechanical properties of the biofibers used in biocomposites. This study aims to investigate the physical and flow properties of surface-treated biofibers, as well as the tensile properties of their PLA-based biocomposite, to gain insights into how surface treatment changes the fiber‘s characteristics and biocomposite‘s mechanical properties. Surface treatment was created using a two-component epoxy system by reacting poly(bisphenol A-co-epichlorohydrin) glycidyl end-capped (PBG) and dicyandiamide (DICY). The treatment was tested on two different biofibers (loblolly pine and corn stover fibers) with three different PBG/DICY molar ratios (0.25, 0.5, and 2). Results showed that surface-treated fibers improved the tensile strength and Young‘s modulus of the biocomposites. Loblolly pine biocomposites from fibers treated with a PBG/DICY ratio of 0.25 exhibited the best tensile properties. The surface treatment resulted in a more loosely dispersed fiber bulk structure, as evidenced by less fiber agglomeration into smaller particle sizes, higher fiber sphericity, and lower loose bulk density. This can enhance stress distribution and the overall mechanical performance of the biocomposites. Additionally, surface-treated fibers exhibited better dynamic flow properties.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Integrated Analysis of Interferometric SAR, Satellite Altimetry and Hydraulic Modeling to Quantify Louisiana Wetland Dynamics

Wetland loss in Louisiana has been accelerating due primarily to anthropogenic and nature processes, and is being advocated as a problem with national importance. Accurate measurement or modeling of wetland-wide water level changes, its varying extent, its storage and discharge changes resulting in part from sediment loads, erosion and subsidence are fundamental to assessment of hurricane-induced flood hazards and wetland ecology. Here, we use innovative method to integrate interferometric SAR (InSAR) and satellite radar altimetry for measuring absolute or geocentric water level changes and applied the methodology to remote areas of swamp forest in coastal Louisiana. Coherence analysis of InSAR pairs suggested that the HH polarization is preferred for this type of observation, and polarimetric analysis can help to identi:fy double-bonnce backscattering areas in the wetland. Envisat radar altimeter-measured 18- Hz (along-track sampling of 417 m) water level data processed with regional stackfile method have been used to provide vertical references for water bodies separated by levees. The high-resolution (approx.40 m) relative water changes measured from ALOS PALSAR L-band and Radarsat-l C-band InSAR are then integrated with Envisat radar altimetry to obtain absolute water level. The resulting water level time series were validated with in situ gauge observations within the swamp forest. Furthermore, we compare our water elevation changes with 2D flood modeling from LISFLOOD hydrodynamic model. Our study demonstrates that this new technique allows retrospective reconstruction and concurrent monitoring of water conditions and flow dynamics in wetlands, especially those lacking gauge networks.

Lee, Hyongki↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2026

Enabling microbial electrolysis cell scale-up via electrochemistry-, hydrodynamic-, and microbial ecology-informed framework

Microbial electrolysis cells (MECs) can produce green hydrogen while removing organic contaminants from liquid waste streams by leveraging the metabolic activity of electroactive microorganisms. Despite their potential in a sustainable, circular economy, large-scale MECs that can treat relevant volumes of wastewater have failed to deliver performance proportional to their lab-scale counterparts. The reason behind this lower performance at scale remains unclear. Here, in this study, we developed a combined electrochemistry-, hydrodynamic-, and microbial ecology-informed framework to analyze and optimize MEC performance during scale-up, enabling accurate quantification of major limitations and the identification of strategies to overcome them, ultimately facilitating equivalent performance at scale. Applying this framework to the scale-up of a zero-gap MEC from 9 cm 2 electrode area to 100 cm 2 electrode area, resulted in similar maximum current densities in a 100 cm 2 MEC (21.7 ± 1.1 A/m 2 ) compared to a 9 cm 2 system (25.1 ± 2.7 A/m 2 ), as well as equivalent hydrogen production rates of 69.3 L/L-d (100 cm 2 ) and 67.7 ± 2.4 L/L-d (9 cm 2 ). COMSOL flow dynamics simulations were used to scale up the reactor configuration without negatively affecting electrolyte velocity and distribution in the cell, minimizing the increase in internal resistances during scale-up (11.7 ± 0.5 mΩm 2 at 9 cm 2 ; 19.7 ± 1.3 mΩm 2 at 100 cm 2 ). Microbial community structures were assessed at both scales using high-throughput sequencing, highlighting the differences of populations across electrode dimensions and operational parameters. The framework presented here accelerates the development of effective strategies toward the scale-up of MECs by furthering the understanding of how electrochemical, hydrodynamic, and microbial ecology parameters change as the reactor dimension is increased. Ultimately, this approach contributes to advancing electrochemical biotechnology toward practical deployment in energy-efficient wastewater treatment systems.

Flow path↗
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↗
DOE OSTI2025

Elongated particles in flow: commentary on small-angle scattering investigations

Here, this work thoroughly examines several analytical tools, each possessing a different level of mathematical intricacy, for the purpose of characterizing the orientation distribution function of elongated objects under flow. Our investigation places an emphasis on connecting the orientation distribution to the small-angle scattering spectra measured experimentally. The diverse range of mathematical approaches investigated herein provide insights into the flow behavior of elongated particles from different perspectives and serve as powerful tools for elucidating the complex interplay between flow dynamics and the orientation distribution function.

36 MATERIALS SCIENCE↗