Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase flows”

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

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Potential seal bypass and caprock storage produced by deformation-band-to-opening-mode-fracture transition at the reservoir/caprock interface

Here, we examined the potential impact on CO 2 transport of zones of deformation bands in reservoir rock that transition to opening-mode fractures within overlying caprock. Sedimentological and petrophysical measurements were collected along an approximately 5 m × 5 m outcrop of the Slick Rock and Earthy Members of the Entrada Sandstone on the eastern flank of the San Rafael Swell, Utah, USA. Measured deformation band permeability (2 mD) within the reservoir facies is about three orders of magnitude lower than the host sandstone. Average permeability of the caprock facies (0.0005 mD) is about seven orders of magnitude lower than the host sandstone. Aperture-based permeability estimates of the opening-mode caprock fractures are high (3.3 × 10 7 mD). High-resolution CO 2 –H 2 O transport models incorporate these permeability data at the millimeter scale. We then varied fault properties at the reservoir/caprock interface between open fractures and deformation bands as part of a sensitivity study. Numerical modeling results suggest that zones of deformation bands within the reservoir strongly compartmentalize reservoir pressures largely blocking lateral, cross-fault flow of supercritical CO 2 . Significant vertical CO 2 transport into the caprock occurred in some scenarios along opening-mode fractures. The magnitude of this vertical CO 2 transport depends on the small-scale geometry of the contact between the opening-mode fracture and the zone of deformation bands, as well as the degree to which fractures penetrate caprock. Finally, the presence of relatively permeable units within the caprock allows storage of significant volumes of CO 2 , particularly when the fracture network does not extend all the way through the caprock.

58 GEOSCIENCES↗
DOE OSTI2020

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Dynamic reduced-order models of integrated physics-specific systems for carbon sequestration

In this study, CO 2 leakage from a deep storage reservoir into a shallow groundwater aquifer through a wellbore is simulated using physics-specific models conceptualized as nonisothermal two-phase flow in a deep storage reservoir and through the wellbore, and reactive transport in the aquifer. Supercritical CO 2 is injected into the reservoir subsystem with uncertain permeabilities of reservoir and injection strategies (rate and duration as decision variables). The simulated pressure and CO 2 saturation at an uncertain location are connected to the wellbore model as a boundary condition to output CO 2 flux to the reactive transport model in the shallow aquifer. Uncertainties are propagated from the deep reservoir model, to the wellbore leakage model, and eventually to the reactive transport model, thus contributing to system performance. Finally, to minimize the risk associated with the CO 2 leakage, we develop dynamic reduced order models for quantifying and managing the uncertainty propagation.

58 GEOSCIENCES↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Experimental Investigation of the Effect of Atomizing Electrospray Nozzles on the Cooling and Reduced Humidification of Air

Separating sensible and latent cooling have shown the potential to reduce energy consumption in A/C systems for building HVAC applications. Several technologies exist based on vapor-compression refrigeration, enthalpy wheels, chemical adsorption and absorption materials, and mechanical cooling. However, their thermodynamics limits and high energy consumption hinder their deployment in hot and humid climatic regions. This paper focuses on separating sensible and latent cooling by using a new thermodynamics process and in-kind (i.e., non-vapor compression-based) technology. Highly electrically charged water droplets were sprayed in the airflow. These droplets attracted water vapor molecules to their surfaces and promoted condensation. The phenomenon was the result of simultaneous dielectrophoresis and electro-diffusion interactions. Studies in the literature have shown that using multiple capillary electrodes reduced air moisture by up to 5% when using nanometer-size droplets in the spray. Unfortunately, these studies were limited to low airflow rates, and the objective of this paper was to investigate how to scale up this in-kind approach to airflows typical of buildings. In the present paper, droplets of micrometer size were utilized to control the humidity for a 5-cfm flow rate. While this airflow was still low for building applications, it was 100 times fold the airflows in the literature studies. The air was tested at 20°C and at 50 and 80% relative humidity. A two-fluid atomizing nozzle produced the droplets in the spray, and high DC electric potential, up to 25 kV, was used to charge the fine droplets electrically. The air atomizing nozzle with high voltage potential resembled an evaporative cooler process. However, a measurable reduction of the absolute humidity of up to 2% was observed compared to the case of the nozzle with no high voltage potential. The entire device had one small nozzle selected from off-the-shelf components and had less than a 9 cm2 footprint area.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Development of Multiple Solid Components Capability in GOTHIC Droplet Fields

GOTHIC is a versatile and generally applicable software package that solves complex thermal hydraulics problems. GOTHIC solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and/or multi-dimensional geometries. GOTHIC version 8.2 includes a new feature to model one or more solid components in each of the modeled droplet fields. A solid component is a material present in the water which has volume, density, and heat capacity. GOTHIC has supported modeling multiple droplet fields since version 8.0 - each field has its own mass, energy, and momentum equation. Recent work by Lane et. al. allowed multiple solid components to be tracked in the continuous liquid field. With the release of version 8.2 the droplet fields can now include one or more solid components within each droplet field. The addition of solid components to the droplet field allows faithful simulation of one or more solid components throughout the liquid and droplet fields. Preliminary expected applications are for scenarios such as GSI-191 and longer term desired applications are aerosol simulation. (authors)

42 ENGINEERING↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
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 OSTI2024

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗
DOE OSTI2016

Enhanced Performance Assessment Models for Generic Deep Geologic Repositories for High-Level Waste and Spent Nuclear Fuel - 16223

Development of an enhanced performance assessment (PA) capability for geologic disposal of spent nuclear fuel and high-level waste has been ongoing for several years in the U.S. repository program. The new Generic Disposal System Analysis (GDSA) modeling and software framework is intended to be flexible enough to evolve through the various phases of repository activities, beginning with generic PA activities in the current Concept Evaluation phase to site-specific PA modeling in the Repository Development phase. The GDSA Framework utilizes modern software and hardware capabilities by being based on open-source software architecture and being configured to run in a massively parallel, high-performance computing (HPC) environment. It consists of two main components, the open-source Dakota uncertainty sampling and analysis software and the PFLOTRAN reactive multi-phase flow and transport simulator. Reference cases or 'generic repositories' have been, and are being developed, based on typical properties for potential salt, clay, and granite host-rock formations and corresponding engineered design concepts for each medium. Past simulations have focused on a generic repository in bedded-salt host rock, while the most recent research has focused on a reference case for a typical clay/shale host rock. A variety of single-realization (i.e., deterministic) and multi-realization (probabilistic) results for the new clay reference case are presented, including an analysis of the effects of heat generation on repository performance, assuming a 100-year out-of-reactor commercial SNF waste form. Order-of-magnitude differences between predicted radionuclide concentrations in thermal versus isothermal simulations imply that mechanistic, coupled-process modeling in three-dimensional (3-D) domains can be important for building confidence in post-closure performance assessments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗
DOE OSTI2019

Effect of Ladle Changeover Condition on Transient Three-Phase Flow in a Five-Strand Bloom Casting Tundish

The transient three-phase flow behavior in a five-strand bloom casting tundish during the ladle changeover process considering a constant casting speed was investigated by the VOF model, and the effect of ladle throughput rate, shroud immersion depth, turbulence inhibitor structure and ladle changeover time on the turbulent flow, and level fluctuations during the filling process were analyzed to optimize the flow pattern in the tundish during ladle changeover. The results showed that the exposed area of molten steel in the impact zone increases as the ladle throughput rate increases during the filling process and that the exposure of molten steel is eliminated when the ladle throughput rate is 1.5Q{sub steady}. As the immersion depth of the shroud increases, a slight decrease in slag entrainment happens, while the exposed area of molten steel first decreases and then increases. The shroud immersion depth should be 310 mm. When the optimized turbulence inhibitor is used, the wave crests of the steel level are greatly reduced, the height of the level fluctuation is alleviated with the maximum value decreasing from 33 mm to approximately 10 mm, and the exposure of molten steel is eliminated. When the ladle changeover time is between 2 and 3 minutes, the exposed area of molten steel during the stable period of the filling process is in the range from 150 to 170 cm{sup 2}, which is substantially less than the exposed area of 500 cm{sup 2} when the ladle changeover time is 4 minutes. Therefore, the ladle changeover time should be controlled within 3 minutes. The flaw detection disqualification rate caused by ladle changeover operation is decreased from 43.8 to 25.8 pct of the total disqualified steel rails after optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Pore-scale displacement mechanisms as a source of hysteresis for two-phase flow in porous media

The macroscopic description of the hysteretic behavior of two–phase flow in porous media remains a challenge. It is not obvious how to represent the underlying pore–scale processes at the Darcy–scale in a consistent way. Darcy–scale thermodynamic models do not completely eliminate hysteresis and our findings indicate that the shape of displacement fronts is an additional source of hysteresis that has not been considered before. This is a shortcoming because effective process behavior such as trapping efficiency of CO 2 or oil production during water flooding are directly linked to pore–scale displacement mechanisms with very different front shape such as capillary fingering, flat frontal displacement, or cluster growth. Here we introduce fluid topology, expressed by the Euler characteristic of the nonwetting phase (χn), as a shape measure of displacement fronts. Using two high–quality data sets obtained by fast X–ray tomography, we show that χn is hysteretic between drainage and imbibition and characteristic for the underlying displacement pattern. In a more physical sense, the Euler characteristic can be interpreted as a parameter describing local fluid connectedness. It may provide the closing link between a topological characterization and macroscopic formulations of two–phase immiscible displacement in porous rock. Since fast X–ray tomography is currently becoming a mature technique, we expect a significant growth in high–quality data sets of real time fluid displacement processes in the future. As a result, the novel measures of fluid topology presented here have the potential to become standard metrics needed to fully explore them.

42 ENGINEERING↗
DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2015

Experimental study on nonmonotonicity of capillary desaturation curves in a 2-D pore-network

Immiscible displacement in a porous medium is important in many applications such as soil remediation and enhanced oil recovery. When gravitational forces are negligible, two-phase immiscible displacement at the pore level is controlled by capillary and viscous forces whose relative importance is quantified through the dimensionless capillary number Ca and the viscosity ratio M between liquid phases. Depending on the values of Ca and M, capillary fingering, viscous fingering, or stable displacement may be observed resulting in a variety of patterns affecting the phase entrapment. The Capillary Desaturation Curve (CDC), which represents the relationship between the residual oils saturation and Ca, is an important relation to describe the phase entrapment at a given Ca. In the present study, we investigate the CDC as influenced by the viscosity ratio. A comprehensive series of experiments using a high-resolution microscope and state-of-the-art micromodels were conducted. The CDCs were calculated and the effects of Ca and M on phase entrapments were quantified. The results show that CDCs are not necessarily monotonic for all M.

capillary desaturation curve↗