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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Radiation/convection coupling in rocket motor and plume analysis

A method for describing radiation/convection coupling to a flow field analysis was developed for rocket motors and plumes. The three commonly used propellant systems (H2/O2, RP-1/O2, and solid propellants) radiate primarily as: molecular emitters, non-scattering small particles (soot), and scattering larger particles (Al2O3), respectively. For the required solution, the divergence of the radiation heat flux was included in the energy equation, and the local, volume averaged intensity was determined by a solution to the radiative transfer equation. A rigorous solution to this problem is intractable, therefore, solution methods which use the ordinary and improved differential approximation were developed. This radiation model was being incorporated into the FDNS code, a Navier-Stokes flowfield solver for multiphase, turbulent combusting flows.

Saladino, A. J.↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2015

Quantifying hydrate solidification front advancing using method of characteristics

Here, we develop a one-dimensional analytical solution based on the method of characteristics to explore hydrate formation from gas injection into brine-saturated sediments within the hydrate stability zone. Our solution includes fully coupled multiphase and multicomponent flow and the associated advective transport in a homogeneous system. Our solution shows that hydrate saturation is controlled by the initial thermodynamic state of the system and changed by the gas fractional flow. Hydrate saturation in gas-rich systems can be estimated by 1 – cl 0 /cl e when Darcy flow dominates, where cl 0 is the initial mass fraction of salt in brine, and cl e is the mass fraction of salt in brine at three-phase (gas, liquid, and hydrate) equilibrium. Hydrate saturation is constant, gas saturation and gas flux decrease, and liquid saturation and liquid flux increase with the distance from the gas inlet to the hydrate solidification front. The total gas and liquid flux is constant from the gas inlet to the hydrate solidification front and decreases abruptly at the hydrate solidification front due to gas inclusion into the hydrate phase. The advancing velocity of the hydrate solidification front decreases with hydrate saturation at a fixed gas inflow rate. This analytical solution illuminates how hydrate is formed by gas injection (methane, CO 2 , ethane, propane) at both the laboratory and field scales.

58 GEOSCIENCES↗
DOE OSTI2015

Numerical Simulation of an Industrial Fluid Catalytic Cracking Regenerator

Fluid catalytic cracking (FCC) is one of the most important conversion processes in petroleum refineries, and the FCC regenerator is a key part of an FCC unit utilized in the recovery of solid catalyst reactivity by burning off the deposited coke on the catalyst surface. Furthermore, a three-dimensional multiphase, multispecies reacting flow computational fluid dynamics (CFD) model was established to simulate the flow and reactions inside an FCC regenerator. The Euler–Euler approach, where the two phases (gas and solid) are considered to be continuous and fully interpenetrating, is employed. The model includes gas–solid momentum exchange, gas–solid heat exchange, gas–solid mass exchange, and chemical reactions. Chemical reactions incorporated into the model simulate the combustion of coke which is present on the catalyst surface. We validated the simulation results using plant data.

42 ENGINEERING↗
DOE OSTI2017

HPC Institutional Computing Project: W15_lesreactiveflow KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software

KIVA-hpFE is a high performance computer software for solving the physics of multi-species and multiphase turbulent reactive flow in complex geometries having immersed moving parts. The code is written in Fortran 90/95 and can be used on any computer platform with any popular complier. The code is in two versions, a serial version and a parallel version utilizing MPICH2 type Message Passing Interface (MPI or Intel MPI) for solving distributed domains. The parallel version is at least 30x faster than the serial version and much faster than our previous generation of parallel engine modeling software, by many factors. The 5th generation algorithm construction is a Galerkin type Finite Element Method (FEM) solving conservative momentum, species, and energy transport equations along with two-equation turbulent model k-ω Reynolds Averaged Navier-Stokes (RANS) model and a Vreman type dynamic Large Eddy Simulation (LES) method. The LES method is capable modeling transitional flow from laminar to fully turbulent; therefore, this LES method does not require special hybrid or blending to walls. The FEM projection method also uses a Petrov-Galerkin (P-G) stabilization along with pressure stabilization. We employ hierarchical basis sets, constructed on the fly with enrichment in areas associated with relatively larger error as determined by error estimation methods. In addition, when not using the hp-adaptive module, the code employs Lagrangian basis or shape functions. The shape functions are constructed for hexahedral, prismatic and tetrahedral elements. The software is designed to solve many types of reactive flow problems, from burners to internal combustion engines and turbines. In addition, the formulation allows for direct integration of solid bodies (conjugate heat transfer), as in heat transfer through housings, parts, cylinders. It can also easily be extended to stress modeling of solids, used in fluid structure interactions problems, solidification, porous media modeling and magneto hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Developing Enzyme Immobilization with Fibrous Membranes: Longevity and Characterization Considerations

Fibrous membranes offer broad opportunities to deploy immobilized enzymes in new reactor and application designs, including multiphase continuous flow-through reactions. Enzyme immobilization is a technology strategy that simplifies the separation of otherwise soluble catalytic proteins from liquid reaction media and imparts stabilization and performance enhancement. Flexible immobilization matrices made from fibers have versatile physical attributes, such as high surface area, light weight, and controllable porosity, which give them membrane-like characteristics, while simultaneously providing good mechanical properties for creating functional filters, sensors, scaffolds, and other interface-active biocatalytic materials. This review examines immobilization strategies for enzymes on fibrous membrane-like polymeric supports involving all three fundamental mechanisms of post-immobilization, incorporation, and coating. Post-immobilization offers an infinite selection of matrix materials, but may encounter loading and durability issues, while incorporation offers longevity but has more limited material options and may present mass transfer obstacles. Coating techniques on fibrous materials at different geometric scales are a growing trend in making membranes that integrate biocatalytic functionality with versatile physical supports. Biocatalytic performance parameters and characterization techniques for immobilized enzymes are described, including several emerging techniques of special relevance for fibrous immobilized enzymes. Diverse application examples from the literature, focusing on fibrous matrices, are summarized, and biocatalyst longevity is emphasized as a critical performance parameter that needs increased attention to advance concepts from lab scale to broader utilization. This consolidation of fabrication, performance measurement, and characterization techniques, with guiding examples highlighted, is intended to inspire future innovations in enzyme immobilization with fibrous membranes and expand their uses in novel reactors and processes.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1989

Vortex generator for flow control

Fluidics flow control of a multiphase supply using a cylindrical chamber is achieved by introducing the supply flow radially into the chamber. The supply flow exits through a port in the center at the chamber. A control fluid is then introduced tangentially about 90.degree. upstream from the supply port. A second control fluid port may be added about 90.degree. upstream from the first control fluid port, but preferably two sets of supply and control ports are added with like ports diametrically opposite each other. The control fluid flows against the circular wall of the control chamber, which introduces a vortex in the flow of the supply flow that decays into a spiral path to the exit port in the center of the chamber. The control flow rate may thus be used to control the spiral path, and therefore the supply flow rate through the exit port.

Collins, Jr., Earl R.↗
NASA NTRSDate not supplied

Launching Liquid Slugs Through Tube Networks: Investigating Multiphase Liquid Slug Acceleration Behaviors in Spacecraft Propulsion Systems using 3D Multiphase VOF CFD Simulations and High-Speed Flow-Visualization Validation Experiments

Propellant management in spacecraft propulsion systems is complex. In many current and future crewed and robotic spacecraft propulsion systems, a number of unintended or unavoidable scenarios can occasionally lead to liquid “slugs” of propellant accumulating in unwanted locations and thereafter being rapidly accelerated through complex and arbitrary tubing networks otherwise filled with vapor or gas. To develop better understanding of the multiphase fluid dynamics associated with these “slug launch” scenarios, NASA Marshal Space Flight Center’s (MSFC) ER42 branch – the propulsion fluid dynamics research branch – conducted 3D multiphase transient Volume of Fluid (VOF) CFD simulations using an in-house tool called LOCI-Stream VOF on semi-canonical tube networks with straight sections, bends, and bend complexes. In addition, high-speed flow visualization validation experiments using propellants were conducted on similar tube network geometries at the NASA White Sands Test Facility (WSTF). These simulations and flow-visualization experiments found these scenarios lead to a variety of unique multiphase fluid dynamic effects involving the leading and trailing edges of the slug in both straight sections and curved bends as well as the progressive erosion, acceleration, and eventual “punch-through” of the driving gas through the liquid slug. Understanding these multiphase fluid dynamic effects has led to unique and often counter-intuitive observations which have been of critical relevance to important engineering concerns associated with these “slug launch” scenarios.

Multiphase Fluid Dynamics↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

A model and numerical method for compressible flows with capillary effects

A new model for interface problems with capillary effects in compressible fluids is presented together with a specific numerical method to treat capillary flows and pressure waves propagation. This new multiphase model is in agreement with physical principles of conservation and respects the second law of thermodynamics. A new numerical method is also proposed where the global system of equations is split into several submodels. Each submodel is hyperbolic or weakly hyperbolic and can be solved with an adequate numerical method. This method is tested and validated thanks to comparisons with analytical solutions (Laplace law) and with experimental results on droplet breakup induced by a shock wave.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2015

Multiphase lattice Boltzmann simulations for porous media applications

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this work, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

42 ENGINEERING↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗
DOE OSTI2021

Active motion of multiphase oil droplets: emergent dynamics of squirmers with evolving internal structure

Synthetic soft matter systems, when driven beyond equilibrium by active processes, offer the potential to achieve dynamical states and functions of a complexity found in living matter. Emulsions offer the basis of a simple yet versatile system for identification of the physicochemical principles underlying active soft matter, but how multiple internal phases within emulsion droplets (e.g., Janus morphologies) organize to impact emergent dynamics is not understood. Here, we create multiphase oil droplets with ultralow interfacial tensions but distinct viscosities, and drive them into motion in aqueous micellar solutions. Preferential solubilization of select components of the oil both drives the droplet motion and yields a progression of internal phase morphological states with distinct symmetries. In this work, we find the active droplets to exhibit five dynamical states during morphogenesis. By quantifying microscopic flow fields, we show that it is possible to map the diverse droplet behaviors to squirmer models of spherical microswimmers in Stokes flow, thus showing that multiphase droplets offer the basis of a versatile platform with which to study and engineer the hydrodynamics of microswimmers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗