Multiphase flow
Multiphase flow is the simultaneous movement of two or more immiscible thermodynamic phases—such as gases, liquids, and solids—within a conduit, domain, or system, where each phase retains its distinct physical and chemical properties while interacting dynamically at their interfaces. This phenomenon often involves phase changes (diabatic flows like boiling or condensation) or no phase change (adiabatic flows), and is characterized by non-equilibrium conditions where phases may differ in velocity, temperature, and concentration, leading to complex interfacial phenomena such as momentum, heat, and mass transfer.[1][2] The study of multiphase flow spans a broad spectrum of scales, from microscopic dispersed particles to macroscopic industrial systems, and is governed by fundamental conservation laws of mass, momentum, and energy applied to each phase or the mixture as a whole. Key variables include the void fraction (the volume fraction occupied by a dispersed phase, such as gas in a liquid), flow quality (the mass fraction of one phase), and velocity ratios between phases, which influence flow behavior and stability. Interfacial instabilities, like Rayleigh-Taylor or Kelvin-Helmholtz types, play a critical role in determining flow patterns and transitions.[1][3] Multiphase flows exhibit distinct regimes based on phase fractions, velocities, fluid properties, and system geometry, including bubbly flow (dispersed bubbles in continuous liquid), slug flow (large bubbles alternating with liquid slugs), annular flow (liquid film around a gas core), and churn flow (chaotic mixing of phases). These regimes are mapped using empirical correlations, such as those developed by Baker, Taitel-Dukler, or Hewitt-Roberts, which rely on extensive experimental data to predict transitions and pressure drops comprising frictional, gravitational, and accelerational components.[1] Applications of multiphase flow are essential in numerous industries and natural processes, including oil and gas pipelines (where gas-liquid mixtures affect transport efficiency), nuclear reactors and boilers (involving steam-water flows for heat transfer), chemical processing equipment like distillation columns and heat exchangers, and environmental systems such as groundwater remediation or atmospheric aerosols. In biomedical contexts, it models blood flow as a suspension of cells in plasma. The field's challenges, including accurate prediction of flow instabilities and critical heat flux, drive ongoing research in modeling approaches ranging from homogeneous equilibrium models to advanced multi-fluid and computational fluid dynamics (CFD) methods with interface-tracking techniques like volume-of-fluid (VOF).[1][3][2]Introduction
Definition and Scope
Multiphase flow refers to the simultaneous flow of two or more thermodynamic phases—typically gas, liquid, and/or solid—within a shared domain, where the phases remain distinct due to immiscibility and are separated by interfaces.[3] This distinguishes it from single-phase flow, which involves only one uniform phase, such as pure water in a pipe, and excludes fully miscible mixtures like alloys or solutions where components mix homogeneously at the molecular level without persistent interfaces.[1] The thermodynamic phases involved are the solid (rigid structure with fixed shape), liquid (incompressible fluid with free surface), and gas (compressible fluid with low density), with interfaces representing the boundaries where phase properties change abruptly, such as the surface tension-driven boundary between a gas bubble and surrounding liquid.[3] The scope of multiphase flow encompasses immiscible combinations, including gas-liquid (e.g., air bubbles in water), liquid-solid (e.g., slurry flows), gas-solid (e.g., pneumatic transport of particles), liquid-liquid (e.g., oil and water), and more complex three- or multi-phase systems (e.g., oil, water, and gas with sand).[1] Configurations within this scope are broadly classified as dispersed, where one phase forms discrete elements like bubbles, droplets, or particles embedded in a continuous carrier phase, or separated, where phases form continuous streams or layers in contact via interfaces, such as stratified liquid layers under a gas flow.[3] These configurations arise due to differences in phase densities, viscosities, and velocities, influencing the overall flow behavior. A key aspect of the scope involves flow regimes, which describe the spatial distribution and morphology of phases, often visualized through regime maps that plot transitions based on parameters like phase velocities and void fractions.[1] For gas-liquid flows in pipes, common regimes include bubbly (discrete small gas bubbles uniformly dispersed in liquid), slug (elongated gas pockets alternating with liquid plugs), and annular (a fast-moving gas core surrounded by a liquid film on the walls, possibly with entrained droplets).[1] These regimes highlight the diversity within multiphase systems, extending to other phase combinations but primarily characterized for gas-liquid due to their prevalence in engineering contexts.Importance and Relevance
Multiphase flow holds significant scientific relevance as a cornerstone of fluid dynamics, thermodynamics, and transport phenomena, where it elucidates complex interactions among immiscible phases such as liquids, gases, and solids.[4] These flows drive advancements in understanding heat and mass transfer across interfaces, influencing fields from chemical engineering to materials science.[5] However, their study is challenged by inherent nonlinearity—stemming from turbulent mixing and phase transitions—and the intricacies of deformable interfaces, which complicate accurate modeling and simulation.[6][7] In environmental contexts, multiphase flows underpin key natural processes, including ocean currents that involve bubble entrainment and sediment-laden water movement, affecting global circulation and marine ecosystems.[8] Atmospheric precipitation relies on multiphase interactions in cloud microphysics, where water vapor condenses into droplets and ice particles, driving weather patterns and hydrological cycles.[9] Furthermore, these flows govern pollutant dispersion, as aerosols and particulates interact with air or water phases to influence air quality and contamination spread in both atmospheric and aquatic environments.[10][11] Economically, multiphase flow is indispensable in resource extraction industries like oil and gas, where it facilitates the transport of hydrocarbons alongside water and gas phases to meet global energy demands.[12] The sector's annual investments and revenues, with revenues in the exploration and production segment reaching approximately $4.0 trillion in 2025, underscore the immense economic stakes tied to efficient multiphase management in production, pipelines, and refining.[13] Despite progress, notable research gaps remain, particularly in nanoscale multiphase flows, where molecular-scale transport mechanisms in confined spaces challenge existing continuum models and require advanced experimental validation.[14] Post-2020 advances in AI-driven predictions have improved flow regime forecasting and phase distribution accuracy, yet integration with traditional simulations for real-time applications is still underdeveloped.[15][16]Historical Development
Early Observations and Theories
Early observations of multiphase flows in nature date back to ancient times, with detailed accounts of volcanic eruptions providing some of the earliest records. In 79 AD, Pliny the Elder and Pliny the Younger documented the eruption of Mount Vesuvius, describing a massive cloud of ash and gas rising from the volcano, part of the pyroclastic activity involving mixtures of hot gases, ash, and rock fragments that surged down the slopes and buried Pompeii and Herculaneum.[17] These eyewitness descriptions captured the dynamic interaction of solid particulates suspended in gaseous and potentially vapor phases, highlighting the destructive power of such natural multiphase phenomena without any theoretical framework. During the Renaissance, Leonardo da Vinci contributed qualitative insights through his sketches and notes on fluid behaviors, including turbulent mixtures in water flows. In the early 1500s, da Vinci illustrated air-water interactions, such as foam formation at waterfalls and the entrainment of air into water jets impacting surfaces, recognizing these as heterogeneous mixtures rather than uniform fluids.[18] His observations of eddies and vortices in rivers also noted sediment-laden flows, where particles were suspended and transported, laying early groundwork for understanding phase interfaces in natural settings.[19] In the 18th and 19th centuries, naturalists and engineers began systematic studies of river sediment transport as a key multiphase process. More quantitatively, in 1861, Andrew Atkinson Humphreys and Henry L. Abbot conducted surveys of the Mississippi River, measuring sediment loads and observing how suspended solids altered flow resistance and channel morphology through empirical sampling. These works emphasized the role of water as a carrier phase for solid particles, influencing early geomorphological concepts.[20] Foundational theories emerged from extensions of single-phase fluid dynamics to mixtures. Daniel Bernoulli's 1738 treatise Hydrodynamica introduced energy conservation principles for fluids, which were later adapted to account for momentum exchange in dilute suspensions, providing a basis for analyzing pressure drops in multiphase systems.[21] Leonhard Euler advanced this in 1757 with his equations describing inviscid flow, offering a continuum framework that influenced early models of fluid-particle interactions by treating mixtures as averaged properties.[22] In the mid-19th century, George Gabriel Stokes developed his 1851 law for the drag force on spherical particles in viscous liquids, enabling predictions of settling velocities in suspensions and quantifying interfacial effects in dilute multiphase flows.[22] Without modern instrumentation like high-speed imaging or rheometers, early researchers faced significant challenges in quantifying interface dynamics and phase interactions. Observations relied on visual inspections and rudimentary measurements, leading to empirical correlations for phenomena such as sediment capacity in rivers, often based on proportional relationships between flow speed and particle size.[20] These limitations fostered a reliance on analogies from single-phase hydrodynamics, delaying precise theories for complex interfacial tensions and phase slip until experimental advancements in the 20th century.[23]Key Milestones in the 20th and 21st Centuries
In the mid-20th century, significant progress was made in understanding pressure drops in two-phase flows, with the Lockhart-Martinelli method emerging as a foundational correlation. Developed in 1949, this approach proposed a dimensionless parameter to relate the pressure gradient in two-phase flow to those in single-phase flows, enabling predictions of frictional losses in isothermal gas-liquid mixtures across various pipe orientations and flow conditions.[24] The method's simplicity and broad applicability made it a staple for engineering design in pipelines and heat exchangers, influencing subsequent models for decades. Building on empirical data, the 1970s saw advancements in flow regime prediction through mechanistic models. In 1976, Taitel and Dukler introduced a comprehensive flow regime map for horizontal and near-horizontal gas-liquid flows, based on physical mechanisms such as wave stability, flooding, and Kelvin-Helmholtz instability.[25] This model predicted transitions between regimes like stratified, intermittent, and annular flows using superficial velocities and fluid properties, providing a predictive tool that reduced reliance on purely empirical correlations and was widely adopted for oil-gas transport systems.[26] The 1980s and 1990s marked a shift toward advanced experimental diagnostics, enhancing visualization and quantification of multiphase regimes. High-speed imaging techniques, utilizing cameras capable of thousands of frames per second, allowed real-time capture of interfacial dynamics in opaque flows, revealing details of bubble coalescence and droplet entrainment previously inaccessible.[27] Concurrently, nuclear methods such as gamma-ray densitometry and neutron radiography emerged for non-intrusive measurements of phase fractions and velocities in industrial-scale pipes, particularly in nuclear and chemical processes, improving accuracy in regime identification under high-pressure conditions.[28] Prominent researchers like Geoffrey F. Hewitt played a pivotal role in these developments, especially in nuclear applications. Throughout the late 20th century, Hewitt's work at the UK Atomic Energy Research Establishment advanced understanding of annular and wispy flows in boiling channels, including experimental demonstrations of critical heat flux mechanisms that informed reactor safety designs.[29] His contributions, spanning over 100 publications, emphasized integrated experimental and modeling approaches for two-phase flows in vertical geometries relevant to reactor cooling. Entering the 21st century, computational fluid dynamics (CFD) revolutionized multiphase flow analysis by integrating multiphase models into simulation frameworks. Post-2000, advancements in Eulerian-Eulerian and volume-of-fluid (VOF) methods enabled detailed predictions of phase interactions, turbulence, and heat transfer in complex geometries, with applications in optimizing oil recovery and reactor simulations far surpassing earlier empirical limits.[30] These tools, validated against experimental data, facilitated virtual prototyping and reduced physical testing needs in industries like energy and chemicals. In the 2020s, artificial intelligence and machine learning have addressed longstanding challenges in regime prediction and flow assurance. Recent models, such as convolutional neural networks trained on high-speed imaging datasets, achieve over 90% accuracy in classifying two-phase regimes in real-time, outperforming traditional maps in handling variable conditions like inclined pipes or non-Newtonian fluids.[31] Explainable AI frameworks further integrate fluid properties and sensor data for predictive maintenance in wet gas pipelines, marking a paradigm shift toward data-driven multiphase engineering.[32]Fundamentals
Basic Principles and Phases Involved
Multiphase flows involve the simultaneous movement of two or more thermodynamic phases, such as gases, liquids, and solids, where each phase exhibits distinct physical properties that govern their interactions. Gases are characterized by low density and high compressibility, allowing significant volume changes under pressure variations, which influences phenomena like bubble expansion in bubbly flows. Liquids, in contrast, possess higher density and are typically treated as incompressible under standard conditions, though minor compressibility arises at elevated pressures; they also feature surface tension that stabilizes interfaces against deformation. Solids, often appearing as particulate matter, are incompressible with high density and exhibit frictional interactions, leading to behaviors like settling or suspension in carrier fluids.[33] Common multiphase combinations include gas-liquid mixtures, such as bubbly or misty flows, and liquid-solid slurries, where solid particles are suspended in a liquid continuum, as seen in sediment transport or industrial pipelines. In slurries, the solid phase introduces additional drag and friction effects due to particle-particle and particle-fluid interactions, altering the overall rheology compared to single-phase liquids. These phase properties dictate the dominant mechanisms, such as buoyancy-driven separation in gas-liquid systems or viscous hindrance in solid-laden flows.[33][34] The foundational principles of multiphase flows derive from the conservation laws of mass, momentum, and energy, adapted to account for interphase interactions and relative motions. For mass conservation, each phase satisfies a continuity equation incorporating mass transfer terms at interfaces, ensuring the total mass balance across phases. Momentum conservation includes body forces, pressure gradients, and interfacial forces like drag, while energy conservation addresses heat and work exchanges between phases, often assuming local thermal equilibrium or nonequilibrium conditions. These laws form the basis for analyzing flow stability and phase distribution.[33] A key concept in these principles is the volume fraction \alpha_i for phase i, defined as the proportion of a control volume occupied by that phase, satisfying \sum_i \alpha_i = 1. Volume fractions quantify phase distribution and are essential for mixture properties, such as effective density \rho = \sum_i \alpha_i \rho_i, where \rho_i is the phase density. They enable the averaging of single-phase equations over multiphase domains, capturing spatial heterogeneity without resolving individual interfaces.[33][35] Interface dynamics in multiphase flows arise from the relative motion between phases, characterized by slip velocity u_{slip} = u_i - u_j, the difference in velocities of phases i and j. At interfaces, boundary conditions range from no-slip (equal velocities, as in fully coupled flows) to slip conditions (velocity discontinuity due to surface tension or viscosity mismatches), influencing energy dissipation and wave propagation. Slip velocities drive phenomena like phase separation and are modulated by factors such as density ratios and Reynolds numbers.[33][36] Flow patterns in multiphase systems are qualitatively classified by the arrangement of phases, distinguishing dispersed flows—where one phase forms discrete particles, droplets, or bubbles within a continuous phase—and separated flows, featuring distinct continuous streams or layers with well-defined interfaces. In dispersed patterns, the continuous phase dominates momentum transfer, while the dispersed phase affects local viscosity and turbulence; examples include bubbly flows (gas dispersed in liquid) and particle-laden flows (solids in gas). These patterns emerge from balances between inertial, viscous, and interfacial forces, setting the stage for regime transitions.[33]Key Parameters and Properties
In multiphase flows, particularly gas-liquid systems, the void fraction \alpha represents the volumetric fraction of the gas phase within the total volume of the mixture, serving as a fundamental parameter for characterizing phase distribution and flow behavior. This scalar quantity, often denoted as \alpha_g for the gas phase, directly influences transport properties and is derived from local or average measurements across the flow cross-section.[37] Accurate determination of \alpha is essential for applications such as nuclear reactor safety analysis and oil-gas transport, where it helps predict heat transfer and pressure gradients.[38] Common measurement techniques for void fraction include the quick-closing valve method, which isolates a pipe segment to directly compute the gas volume via mass balance of the drained phases, providing time-averaged values with uncertainties typically below 5% in steady flows.[37] Non-intrusive alternatives, such as gamma-ray densitometry, utilize attenuation of gamma radiation through the flow to infer \alpha based on the linear relationship between beam intensity and phase densities, enabling real-time profiling in opaque pipelines with resolutions down to 1% void fraction.[38] These techniques, often combined with optical or impedance methods, account for flow unsteadiness and phase interfaces to ensure reliability across varying conditions. Closely related to void fraction is the concept of holdup, defined as the time-averaged volumetric fraction of a specific phase present in a flow conduit over an extended period, reflecting the in-situ accumulation due to phase interactions.[39] For instance, liquid holdup H_L quantifies the fraction of the pipe volume occupied by liquid, which can exceed the input volumetric fraction owing to differential phase velocities, and is typically measured via rapid sampling or radiation-based tomography.[40] The slip ratio, denoted as S = u_g / u_l where u_g and u_l are the actual velocities of the gas and liquid phases, respectively, captures the relative velocity difference between phases, often exceeding unity in vertical flows due to buoyancy effects.[41] This ratio links holdup to input flow rates through relations like H_L = \frac{\beta S}{(1 - \beta) + \beta S}, where \beta is the input liquid fraction, and is pivotal for interpreting phase segregation in pipelines.[42] Pressure drop in multiphase systems arises from three primary components: frictional losses due to shear at walls and interfaces, accelerational effects from momentum changes as phases expand or accelerate, and gravitational contributions from hydrostatic head variations across density gradients. The total pressure gradient \frac{dP}{dz} is expressed as \frac{dP}{dz} = \frac{dP_f}{dz} + \frac{dP_a}{dz} + \frac{dP_g}{dz}, where frictional \frac{dP_f}{dz} dominates in horizontal flows and scales with mixture velocity and viscosity, accelerational \frac{dP_a}{dz} becomes significant in expanding gases, and gravitational \frac{dP_g}{dz} = g \sum \alpha_i \rho_i \sin \theta reflects elevation and orientation.[43] Experimental validations in vertical gas-liquid flows show that frictional components can account for up to 70% of the total drop in turbulent regimes, while gravitational terms prevail in low-velocity stratified flows.[44] Beyond these, effective mixture properties like density and viscosity are critical for bulk flow descriptions. The effective density \rho_\text{eff} is computed as the volume-weighted sum \rho_\text{eff} = \sum \alpha_i \rho_i, where \alpha_i and \rho_i are the volume fraction and density of phase i, providing a homogeneous approximation for conservation equations in dispersed flows.[45] For viscosity, multiphase mixtures often exhibit non-Newtonian behaviors, particularly in suspensions where particle interactions lead to shear-thinning or -thickening, with apparent viscosity \mu_\text{eff} increasing nonlinearly with solid volume fraction \phi via models like \mu_\text{eff} = \mu_c (1 - \phi/\phi_m)^{-2.5\phi_m}, where \mu_c is the carrier fluid viscosity and \phi_m the maximum packing.[46] Rheological models for such complex fluids in multiphase contexts, including emulsions and slurries, incorporate microstructural effects like aggregation to predict yield stress and thixotropy, though challenges persist in capturing transient behaviors under high shear.[47]Types of Multiphase Flows
Two-Phase Flows
Two-phase flows involve the simultaneous movement of two distinct phases, such as gas-liquid, liquid-solid, or gas-solid mixtures, where interactions between the phases govern the overall behavior. These flows serve as foundational configurations in multiphase systems, influencing transport efficiency, pressure gradients, and phase distribution in conduits. The classification of two-phase flows relies on the morphology of the interfaces and relative velocities, with transitions determined by physical mechanisms like buoyancy, drag, and interfacial tension.[48] Gas-liquid two-phase flows are among the most studied configurations, exhibiting distinct regimes based on the superficial velocities of each phase—the volume flow rate divided by the cross-sectional area. In horizontal pipes, common regimes include bubbly flow, where discrete gas bubbles are dispersed in a continuous liquid phase; slug flow, characterized by large elongated gas bubbles alternating with liquid slugs; churn flow, a transitional turbulent regime with chaotic mixing; annular flow, featuring a liquid film along the pipe wall surrounding a gas core; and mist flow, where liquid droplets are entrained in a continuous gas phase. These regimes are mapped using criteria such as the Mandhane et al. (1974) flow pattern map, which correlates transitions to superficial gas and liquid velocities for pipes of various diameters. In vertical orientations, regimes differ due to gravitational effects, with bubbly flow persisting at higher gas fractions before transitioning to slug, churn, and annular-mist patterns, as modeled mechanistically by Taitel et al. (1980) based on bubble rise velocity and flooding limits. Transition criteria often involve dimensionless groups like the Froude number, where bubbly-to-slug shifts occur when gas velocity exceeds a critical value for bubble coalescence.[48][49][25] Liquid-solid two-phase flows, such as slurries and suspensions, occur when solid particles are dispersed or suspended in a liquid carrier. In dilute suspensions, particles settle under gravity, with the settling velocity for spherical particles at low Reynolds numbers given by Stokes' law: v_s = \frac{(\rho_s - \rho_l) g d_p^2}{18 \mu_l} where \rho_s and \rho_l are the densities of the solid and liquid, g is gravitational acceleration, d_p is particle diameter, and \mu_l is liquid viscosity; this relation, derived from balancing gravitational and viscous drag forces, applies when inertial effects are negligible (Re < 1).[50] In denser slurries, hindered settling reduces v_s due to particle interactions, leading to homogeneous or heterogeneous flow depending on concentration. These flows are critical in pipelines where maintaining suspension prevents bed formation, with vertical transport aiding suspension via upward flow while horizontal setups require sufficient velocity to counter settling. Gas-solid two-phase flows, exemplified by pneumatic transport, involve solid particles carried by a gas stream and are classified as dilute or dense phase based on solids loading ratio (mass of solids to mass of gas). Dilute phase flows feature suspended particles with gas velocities well above saltation velocity—the threshold where particles deposit on pipe walls due to insufficient drag—typically resulting in uniform distribution and lower pressure drops. Dense phase flows, with higher loadings, exhibit particle clusters or dunes, promoting saltation in horizontal pipelines where particles saltate (bounce) along the bottom, increasing erosion and energy use; this distinction arises from balancing particle inertia and gas drag, as observed in early studies of vertical and horizontal transport. Pipeline-specific aspects of two-phase flows vary significantly with orientation: horizontal flows promote stratification due to gravity, favoring slug or annular regimes in gas-liquid systems and saltation in gas-solid, whereas vertical flows enhance mixing and suspension through buoyancy, reducing holdup in liquid-solid slurries. In non-Newtonian two-phase flows, such as those involving drilling muds—yield-stress fluids used in oil wells—rheological complexity alters regime transitions and pressure losses, with shear-thinning behavior increasing holdup in annular flows compared to Newtonian counterparts, as evidenced by full-scale experiments in vertical wells. Void fraction, the volume fraction of the dispersed phase, provides a key measure of phase distribution in these configurations.Multiphase Flows Involving Three or More Phases
Multiphase flows involving three or more phases introduce significantly greater complexity than two-phase systems due to the increased number of interfaces and interactions among phases, leading to unique flow regimes that cannot be predicted by simple extensions of binary models. In such systems, phases can include combinations like gas-liquid-liquid or gas-liquid-solid, where mutual influences alter distribution, velocity profiles, and stability. For instance, three-phase oil-water-gas flows in pipelines exhibit regimes such as stratified, slug, and dispersed bubble flows, influenced by phase fractions and superficial velocities, as observed in experimental studies using air-water-oil systems in horizontal pipes. These regimes often show transitional behaviors not seen in pairwise interactions, with water tending to form a lower layer while oil and gas occupy upper regions in stratified configurations.[51] The black oil model serves as a foundational approach for simulating three-phase oil-water-gas flows, particularly in reservoir contexts and extended to wellbore and pipeline simulations. This model treats oil, water, and gas as pseudo-components with simplified phase behavior, accounting for gas dissolution in oil and water without full compositional tracking, enabling efficient prediction of flow dynamics. In practice, it incorporates solution gas-oil ratios and formation volume factors to analyze multiphase transport.[52] Extending to four or more phases, gas-liquid-solid flows, such as those in fluidized beds, involve additional challenges like particle suspension amid bubbling gas and liquid circulation. In these systems, regimes transition from packed beds to bubbling and turbulent fluidization as gas and liquid velocities increase, with solids exhibiting slugging or churning behaviors influenced by interphase drag. Phase identification becomes particularly difficult due to overlapping distributions and rapid morphological changes, requiring advanced techniques like deep learning for accurate segmentation of gas voids, liquid films, and solid clusters from imaging data.[53][54] Regime complexities in three-or-more-phase flows often manifest as hierarchical structures, such as gas bubbles embedded within liquid-solid slurries, where bubbles rise through particle-laden liquids, inducing secondary flows and particle segregation. Particle entrainment into gas cores further complicates dynamics, as solids are lifted from wall films into high-velocity gas streams, enhancing mixing but risking erosion in conduits. These interactions create nested morphologies, like entrained particles within bubble wakes in slurries, which demand multidimensional modeling to resolve. Building briefly on two-phase baselines, such regimes add ternary phase effects that shift transition boundaries.[55] Emerging research highlights multiphase flows with three or more phases in microgravity environments, relevant to space applications like propellant management and life support systems post-2010. In reduced gravity, buoyancy absence leads to symmetric phase distributions and prolonged coalescence times, complicating gas-liquid-solid separations in spacecraft tanks, as reviewed in vapor compression heat pump studies aboard the International Space Station. Similarly, biological systems exemplify such flows, with blood modeled as a multiphase mixture of plasma (liquid), red blood cells (solid-like particles), and minor plasma proteins, where cell aggregation and plasma skimming create heterogeneous distributions in microvessels. These areas underscore the need for gravity-independent models to predict phase partitioning in non-terrestrial and physiological contexts.[56][57]Applications
Natural Phenomena
Multiphase flows play a central role in atmospheric processes, particularly in cloud formation, where water vapor condenses into liquid droplets, creating dispersed gas-liquid systems that influence weather patterns. These flows involve the nucleation and growth of droplets on aerosol particles, leading to the formation of clouds as a multiphase mixture of air, vapor, and liquid phases.[4] In rain formation, droplet coalescence within these clouds drives precipitation, as smaller droplets collide and merge under turbulent conditions to form larger raindrops capable of falling through the atmosphere.[58] Aerosol transport further exemplifies multiphase dynamics, with solid or liquid particles suspended in air facilitating chemical reactions and radiative effects that alter cloud properties and precipitation efficiency.[11] In oceanic and hydrological environments, multiphase flows manifest in sediment-laden rivers, where liquid water carries solid particles in a slurry-like mixture, affecting erosion, deposition, and river morphology. These liquid-solid interactions are governed by turbulence and particle settling, influencing sediment transport over large scales. Ocean waves generate foam through gas-liquid multiphase processes, as air bubbles become entrained in breaking waves, enhancing gas exchange across the air-sea interface and contributing to ocean aeration. Volcanic ash plumes represent complex gas-solid-liquid flows, with hot gases carrying solid ash particles and sometimes liquid water droplets, driving plume rise, dispersion, and atmospheric impacts during eruptions.[59][60][61] Geological settings feature multiphase flows in lava flows, where molten silicate liquid contains dispersed gas bubbles that alter viscosity and flow behavior, leading to varied eruption styles and landform development. In groundwater systems, multi-liquid multiphase flows occur when non-aqueous phase liquids, such as organic contaminants, infiltrate aquifers alongside water, creating immiscible phases that complicate contaminant migration and persistence in porous media.[62][63] Ecological impacts of multiphase flows include nutrient transport in ecosystems, where hyporheic exchange between surface water and groundwater facilitates the cycling of dissolved nutrients through interactions in riverbeds, supporting biodiversity and food webs. Recent studies in the 2020s highlight how climate change intensifies multiphase dynamics, such as accelerated ice melt in polar regions, where solid ice transitions to liquid water and interacts with atmospheric gases, altering freshwater inputs to oceans and ecosystems. These shifts, driven by warming, enhance meltwater flows and influence global nutrient distribution and habitat stability.[64][65]Industrial and Engineering Contexts
In the energy sector, multiphase flows are critical for the transportation of oil and gas mixtures through pipelines, where two-phase gas-liquid flows dominate and influence pressure drops, slugging, and overall efficiency. For instance, in subsea oil-gas pipelines, the management of these flows prevents blockages from hydrates and ensures steady production from remote wells.[66] In nuclear reactors, boiling crises—such as departure from nucleate boiling—pose severe risks during multiphase steam-water flows in cooling systems, potentially leading to fuel rod damage if critical heat flux is exceeded.[67] Carbon capture technologies rely on gas-liquid absorption processes, where multiphase flows in packed columns enhance CO2 dissolution into amine solvents, optimizing capture efficiency in post-combustion systems.[68] Chemical processing extensively utilizes multiphase flows, particularly gas-solid systems in fluidized beds for catalytic reactions, where uniform particle suspension improves reaction rates and heat transfer in processes like fluid catalytic cracking.[69] Vapor-liquid flows in distillation columns govern separation efficiency, with countercurrent multiphase interactions determining tray hydraulics and mass transfer rates in petrochemical refining.[70] In other industries, liquid-solid multiphase flows are essential for pharmaceutical mixing, where agitators suspend powders in liquids to achieve homogeneous formulations, addressing challenges like agglomeration in drug production.[71] Food processing employs emulsions as multiphase systems, stabilizing oil-in-water dispersions for products like mayonnaise through shear-induced droplet breakup and coalescence control.[72] In biomedical applications, blood flow is modeled as a multiphase suspension of plasma, red blood cells, and platelets, influencing rheology and thrombosis risks in vascular devices.[73] Key challenges in these contexts include flow assurance in subsea pipelines, where multiphase interactions cause wax deposition and erosion, necessitating advanced monitoring and mitigation strategies to maintain production integrity.[74] Recent advances extend to renewable energy, such as multiphase two-phase flows in geothermal wells, where steam-liquid mixtures drive enhanced heat extraction.[75] Similarly, hydrogen transport in pipelines involves multiphase considerations, supporting safe integration into renewable grids.[76]Modeling Approaches
Empirical and Semi-Empirical Methods
Empirical and semi-empirical methods in multiphase flow provide practical tools for predicting flow behavior by fitting experimental data to correlations, particularly for pressure drops, void fractions, and regime transitions in two-phase systems. These approaches originated from early experimental studies and remain foundational in engineering applications where computational complexity is prohibitive. They typically assume simplified flow structures, such as homogeneous mixtures or separated phases, and are validated against laboratory data under specific conditions like isothermal gas-liquid flows in pipes. A seminal empirical correlation for frictional pressure drop in two-phase gas-liquid flows is the Lockhart-Martinelli parameter, defined asX = \sqrt{ \frac{ (dP/dz)_L }{ (dP/dz)_G } },
where (dP/dz)_L and (dP/dz)_G represent the single-phase frictional pressure gradients for the liquid and gas, respectively, calculated using their individual flow rates and properties. This parameter relates the two-phase pressure drop to single-phase values through multipliers \phi_L^2 = \frac{ (dP/dz)_{TP} }{ (dP/dz)_L } and \phi_G^2 = \frac{ (dP/dz)_{TP} }{ (dP/dz)_G }, where \phi_L = 1 + \frac{C}{X} + \frac{1}{X^2} and \phi_G = 1 + C X + X^2, with C depending on flow regimes (e.g., 20 for turbulent-turbulent). The correlation distinguishes between homogeneous flow models, which assume equal velocities for both phases and treat the mixture as a single fluid with averaged properties, and separated flow models, which account for distinct phase velocities and interfacial effects via the parameter. Homogeneous models simplify calculations but overpredict holdup in stratified flows, while separated models like Lockhart-Martinelli better capture pressure drops in horizontal pipes by incorporating void fraction data.[77] Flow regime maps empirically delineate transitions between patterns such as bubbly, slug, annular, or stratified flows based on superficial velocities or mass fluxes. The Baker chart, developed for horizontal gas-liquid flows, plots dimensionless parameters \lambda = \frac{\rho_L}{\rho_G} \left( \frac{\mu_G}{\mu_L} \right)^{0.1} and mass velocity ratio \frac{G_G}{G_L}, dividing regimes into dispersed, homogeneous, separated, and intermittent zones using boundaries derived from oil-gas pipeline experiments. These maps aid in selecting appropriate pressure drop correlations but are limited to specific fluid pairs and orientations. Mechanistic semi-empirical models extend this by predicting void fraction \alpha, as in the drift-flux approach, where the mixture velocity j = \sum \alpha_i u_i relates to superficial velocities via the phase velocity u_g = C_0 j + V_{gj}, with void fraction \alpha = \frac{j_g}{u_g}, C_0 as the distribution coefficient (typically 1.2 for bubbly flows) and V_{gj} as the drift velocity. This model improves void fraction predictions in vertical flows by incorporating relative velocity effects from experiments. Semi-empirical tools like the Chisholm parameter refine two-phase multipliers for non-circular geometries or bends, defining K_1 = \sum \left( \frac{\Delta P_i}{\Delta P_{lo}} \right)^{1/2} - 1 for individual phase contributions to total pressure drop, where \Delta P_{lo} is the liquid-only drop. The two-phase multiplier is then \phi_{lo}^2 = 1 + \frac{K_1^2}{X^2} + \frac{1}{X^2}, calibrated against evaporating mixtures in tubes. This approach enhances accuracy for frictional losses but struggles with extrapolation to three-phase or high-viscosity flows due to unaccounted interactions. Recent advancements in the 2020s integrate machine learning for regime classification, using supervised algorithms like random forests or neural networks on pressure or conductance signals to achieve over 90% accuracy in identifying patterns, surpassing traditional maps in dynamic or opaque systems.[78][79]