Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphase modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Bioreactor Optimization through Multi-Phase Flow Models (CRADA Final Report)

Chemical manufacturing uses 29% of energy in the United States and produces 925 million metric tons of CO2 annually. Biomanufacturing offers the potential to leverage America’s rich agricultural resources to produce critical chemicals such as lubricants, pharmaceutical precursors, and components of energetic materials that today are sourced extensively from overseas. The bioreactors used in biomanufacturing applications, such as one developed by Capra Biosciences, involve multiphase flow of biofilm-coated solid support particles that are continuously circulated in a fluidized state within the reactor along with a constant supply of oxygen via an aeration mechanism. In this project, Capra Biosciences and LBNL developed a multiscale modeling framework to simulate the multiphase flows of solid particles in a liquid-gas bubble mixture that occurs in the bioreactor using the current MFIX-Exa software, an opensource multiphase flow solver developed and maintained at LBNL and NETL. By leveraging HPC capabilities, this high-fidelity multiscale model was used to inform design decisions for bioreactor architecture to make them operationally efficient.

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

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Model-parallel Fourier neural operators as learned surrogates for large-scale parametric PDEs

Fourier neural operators (FNOs) are a recently introduced neural network architecture for learning solution operators of partial differential equations (PDEs), which have been shown to perform significantly better than comparable deep learning approaches. Once trained, FNOs can achieve speed-ups of multiple orders of magnitude over conventional numerical PDE solvers. However, due to the high dimensionality of their input data and network weights, FNOs have so far only been applied to two-dimensional or small three-dimensional problems. To remove this limited problem-size barrier, we propose a model-parallel version of FNOs based on domain-decomposition of both the input data and network weights. Here, we demonstrate that our model-parallel FNO is able to predict time-varying PDE solutions of over 2.6 billion variables on Perlmutter using up to 512 A100 GPUs and show an example of training a distributed FNO on the Azure cloud for simulating multiphase CO 2 dynamics in the Earth’s subsurface.

58 GEOSCIENCES↗
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
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 OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2021

Robust hierarchical dispatch for residential distribution network management considering home thermal flexibility and model predictive control

In the transactive energy (TE) paradigm, the devices of participative consumers, or prosumers, may be aggregated and employed to drive operational objectives at the network level. Home heating, ventilation and air-conditioning (HVAC) systems in particular are well-suited to modulate their behaviours based on both home thermal flexibility and requests from the utility grid. This paper develops a robust, hierarchical power dispatch scheme in the context of a residential distribution network. The formulation couples a unique, multiphase linear distribution optimal power flow (OPF) at the upper level with model predictive control (MPC)-based HVAC fleet controllers at the lower level. The proposed approach is tested on nearly 2000 homes with a three-phase distribution network in an intraday market setting, where two major applications are explored and analysed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Molten Halide Salt Surface Tension: Methods and Correlations

Here, this paper reviews various methods for studying surface tension and their applicability to fluoride and chloride molten salt systems, including a comparison of benefits and drawbacks. Such a comparison aids in experiment design based on desired factors such as scale, accuracy, and repeatability. A detailed review is presented for existing literature data regarding the surface tension of molten fluoride and chloride salts. These reference data were compiled and analyzed to determine cross-validated correlation equations for several alkali and alkaline earth fluoride and chloride salts as functions of temperature. These correlations are necessary for reliable multiphysics modeling approaches as well as accurate design and analysis of multiphase molten salt phenomena such as gas sparging and bubble formation/transport. This analysis supports the development of the thermophysical arm of the Molten Salt Thermal Properties Database (MSTDB-TP) managed by Oak Ridge National Laboratory.

Chloride↗
DOE OSTI2024

Role of chemical production and depositional losses on formaldehyde in the Community Regional Atmospheric Chemistry Multiphase Mechanism (CRACMM)

Abstract. Formaldehyde (HCHO) is an important air pollutant with direct cancer risk and ozone-forming potential. HCHO sources are complex because HCHO is both directly emitted and produced from oxidation of most gas-phase reactive organic carbon. We update the secondary production of HCHO in the Community Regional Atmospheric Chemistry Multiphase Mechanism (CRACMM) in the Community Multiscale Air Quality (CMAQ) model. Production of HCHO from isoprene and monoterpenes is increased, correcting an underestimate in the current version. Simulated June–August surface HCHO during peak photochemical production (11:00–15:00 LT, local time) increased by 0.6 ppb (32 %) over the southeastern USA and by 0.2 ppb (13 %) over the contiguous USA. The increased HCHO compares more favorably with satellite-based observations from the TROPOspheric Monitoring Instrument (TROPOMI) and from aircraft-based observations. Evaluation against hourly surface observations indicates a missing nighttime sink that can be improved by increased nighttime deposition, which reduces June–August nocturnal (20:00–04:00 LT) surface HCHO by 1.1 ppb (36 %) over the southeastern USA and 0.5 ppb (29 %) over the contiguous USA. The ability of CRACMM to capture peak levels of HCHO at midday is improved, particularly at sites in the northeastern USA, while peak levels at sites in the southeastern USA are improved, although still lower than observed. Using established risk assessment methods, lifetime exposure of the population in the contiguous USA (∼ 320 million) to ambient HCHO levels predicted here may result in 6200 lifetime cancer cases, with 40 % from controllable anthropogenic emissions of nitrogen oxides and reactive organic compounds. Chemistry updates will be available in CRACMM version 2 (CRACMM2) in CMAQv5.5.

Skipper, T. Nash↗
DOE OSTI2021

NRAP-Open-IAM: FutureGen2 Component Models

This report describes the development and testing of three component models for NRAP-Open-IAM, the National Risk Assessment Partnership’s open-source integrated assessment model. The FutureGen2 Lookup Table Reservoir component model is based on interpolation of data from a set of lookup tables. The lookup tables contain pressures and saturations predicted by multiphase flow simulations performed with a reservoir simulator. The FutureGen2 Above Zone Monitoring Interval (AZMI) component is a surrogate model that can be used to estimate the impact that carbon dioxide (CO 2 ) and brine leaks from the CO 2 storage reservoir at the FutureGen 2.0 site might have had on overlying aquifers or monitoring units were a leak to occur. The model estimates the size of “impact plumes” according to five metrics: pH, Total Dissolved Solids (TDS), pressure, dissolved CO 2 and temperature. The FutureGen2 Aquifer component is similar, but is limited to four metrics: pH, Total Dissolved Solids (TDS), pressure, and dissolved CO 2 . The input parameters for each model are the same, but the Aquifer component is applicable to depths between 100 m and 700 m and the AZMI component is applicable from depths between 700 m and 1050 m.

42 ENGINEERING↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2022

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗
DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

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

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach

Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.

42 ENGINEERING↗