The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） p...The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.展开更多
In this work, we obtain the global existence and uniqueness of classical solu-tions to a viscous liquid-gas two-phase model with mass-dependent viscosity and vacuum in one dimension, where the initial vacuum is allowe...In this work, we obtain the global existence and uniqueness of classical solu-tions to a viscous liquid-gas two-phase model with mass-dependent viscosity and vacuum in one dimension, where the initial vacuum is allowed. We get the upper and lower bounds of gas and liquid masses n and m by the continuity methods which we use to study the compressible Navier-Stokes equations.展开更多
In this paper,the Cauchy problem for a two-phase model with a magnetic field in three dimensions is considered.Based on a new linearized system with respect to(c−c_(∞),P−P_(∞),u,H)for constants c_(∞)≥0 and P_(∞)&...In this paper,the Cauchy problem for a two-phase model with a magnetic field in three dimensions is considered.Based on a new linearized system with respect to(c−c_(∞),P−P_(∞),u,H)for constants c_(∞)≥0 and P_(∞)>0,the existence theory of global strong solution is established when the initial data is close to its equilibrium in three dimensions for the small H^(2) initial data.We improve the existence results obtained by Wen and Zhu in[40]where an additional assumption that the initial perturbations are bounded in L^(1)-norm was needed.The energy method combined with the low-frequency and high-frequency decomposition is used to derive the decay of the solution and hence the global existence.As a by-product,the time decay estimates of the solution and its derivatives in the L^(2)-norm are obtained.展开更多
A two-phase model for the prediction of macrosegregation formed during solidification is presented. This model incorporates the descriptions of heat transfer, melt convection, solute transport, and solid movement on t...A two-phase model for the prediction of macrosegregation formed during solidification is presented. This model incorporates the descriptions of heat transfer, melt convection, solute transport, and solid movement on the system scale with microscopic relations for grain nucleation and growth. Then the model is used to simulate the solidification of a benchmark industrial 3.3-t steel ingot. Simulations are per- formed to investigate the effects of grain motion and pipe shrinkage formation on the final macrosegregation pattern. The model predictions are compared with experimental data and numerical results from literatures. It is demonstrated that the model is able to express the overall macrosegregation patterns in the ingot. Furthermore, the results show that it is essential to consider the motion of equiaxed grains and the formation of pipe shrinkage in modelling. Several issues for future model improvements are identified.展开更多
In this paper, laminar forced convection of CuO nanofluid is numerically investigated in sudden expansion microchannel with isotherm walls and different expansion ratios (ER). An Eulerian two-fluid model is considered...In this paper, laminar forced convection of CuO nanofluid is numerically investigated in sudden expansion microchannel with isotherm walls and different expansion ratios (ER). An Eulerian two-fluid model is considered to simulate the nanofluid flow inside the microchannel and the governing mass, momentum and energy equations for both phases are solved using the finite volume method. Eulerian-Eulerian two-phase model is very efficient because of considering the relative velocity and temperature of the phases and the nanoparticle concentration distribution. In solving the flow equations for both phases, the SIMPLE algorithm is modified for the coupling of the velocity and pressure and the continuity equations for both phases are combined in order to create the pressure correction equations. However, the Eulerian-Eulerian modeling results show higher heat transfer enhancement in comparison to pure water, so that for a 2% copper-water nanofluid, it has been observed a 35% increase of the heat transfer. The heat transfer enhancement increases with increase in Reynolds number and nanoparticle volume concentration, while the pressure drop increases only slightly. An investigation of the expansion ratio of microchannel shows that the average Nusselt number increases with decrease in expansion ratio as well as with increase in Reynolds number. Also, the Bifurcation has been occurred in higher Reynolds number that is different for each expansion ratio of the microchannel.展开更多
A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctua...A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.展开更多
Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishe...Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mec...Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mechanisms of mass-transfer enhancement in TPPBs would enable efficient predictions for further industrial applications.In this study,influences of gradually increasing silicone oil ratio on the TPPB was explored,and a 94.35%reduction of the n-hexane partition coefficient was observed with 0.1 vol.%silicone,which increased to 80.7%along with a 40-fold removal efficiency enhancement in the stabilised removal period.The elimination capacity increased from 1.47 to 148.35 g/(m^(3)·h),i.e.a 101-fold increase compared with that of the single-phase reactors,when 10 vol.%(3 Critical Micelle Concentration)silicone oil was added.The significantly promoted partition coefficient was the main reason for the mass transfer enhancement,which covered the negative influences of the decreased total mass-transfer coefficient with increasing silicone oil volume ratio.The gradually rising stirring rate was benefit to the n-hexane removal,which became negative when the dominant resistance shifted from mass transfer to biodegradation.Moreover,a mass-transfer-reaction kinetic model of the TPPB was constructed based on the balance of n-hexane concentration,dissolved oxygen and biomass.Similar to the mechanism,the partition factor was predicted sensitive to the removal performance,and another five sensitive parameters were found simultaneously.This forecasting method enables the optimisation of TPPB performance and provides theoretical support for hydrophobic VOCs degradation.展开更多
By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using comput...By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using computational fluid dynamics(CFD),the two-phase flow and piston engraving characteristics of a pyrotechnic actuator are investigated.Initially,the current model was utilized to examine the intricate,multi-dimensional flow,and energy conversion characteristics of the propellant grains and combustion gas within the pyrotechnic actuator chamber.It was discovered that the combustion gas on the wall's constant transition from potential to kinetic energy,along with the combined effect of the propellant motion,are what create the pressure oscillation within the chamber.Additionally,a numerical analysis was conducted to determine the impact of various parameters on the pressure oscillation and piston motion,including pyrotechnic charge,pyrotechnic particle size,and chamber structural dimension.The findings show that decreasing the pyrotechnic charge will lower the terminal velocity,while increasing and decreasing the pyrotechnic particle size will reduce the pressure oscillation in the chamber.The pyrotechnic particle size has minimal bearing on the terminal velocity.The results of this investigation offer a trustworthy forecasting instrument for comprehending and creating pyrotechnic actuator designs.展开更多
As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it ...As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it has several limitations when addressing multi-source heat dissipation.Inspired by the transport and heat dissipation system of plants,large trees achieve stable and efficient liquid supply under the influence of two driving forces:capillary force during transpiration in the leaves(pull)and root pressure generated by osmotic pressure in the roots(push).The root pressure provides an effective liquid supply with a driving force exceeding 2 MPa,far greater than the driving force in conventional capillary-pumped two-phase loops.Research has shown that osmotic heat pipes offer a powerful driving force,and combining osmotic pressure with capillary force has significant advantages.Therefore,this paper designs a multi-evaporator,dual-drive two-phase loop,using both osmotic pressure and capillary force to solve the multi-source heat dissipation challenge.First,a transmembrane water flux model for the osmotic pressure-driven device was established to predict the maximum heat transfer capacity of the dual-drive two-phase loop.Then,an experimental setup for a multi-evaporator“osmotic pressure+capillary force”dual-drive two-phase loop was constructed,capable of transferring at least 235 W of power under a reverse gravity condition of 20 m.The study also analyzed the effects of reverse gravity height,heat load distribution among the three evaporators,startup sequence,and varying branch resistances on the performance of the dual-drive two-phase loop.展开更多
The Proton Exchange Membrane(PEM)water electrolyzer is considered one of the promising energy storing means for harnessing variable renewable energy sources to produce hydrogen.Understanding the internal fluid dynamic...The Proton Exchange Membrane(PEM)water electrolyzer is considered one of the promising energy storing means for harnessing variable renewable energy sources to produce hydrogen.Understanding the internal fluid dynamics,which are often challenging to directly observe experimentally,has prompted the use of numerical models to investigate two-phase flow within PEM water electrolyzers.In this study,we provide a comprehensive review of prior research focusing on two-phase modeling of PEM electrolyzers,encompassing both components at mesoscopic scales and the full electrolyzer at the macroscopic level.We delve into the specifics of various modeling approaches for two-phase flow at different scales and summarize and discuss the current state of the art in the field.Presently,two-phase models for the full electrolyzer predominantly employ a macroscopic homogeneous assumption.However,mesoscopic and microscopic models capable of tracking phase interfaces are limited to components.Challenges persist in integrating various modeling scales into a comprehensive electrolyzer model,particularly in coupling two-phase flow between the channels and porous media.Future efforts should focus on developing multi-scale models and simulating two-phase flow under fluctuating input conditions.Additionally,given the structural similarities between PEM water electrolyzers and PEM fuel cells,we compare and discuss differences in two-phase modeling between the two technologies.This work offers the insights for researchers in the field of modeling of PEM water electrolyzers and even fuel cells.展开更多
With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchan...With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.展开更多
Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristic...Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.展开更多
Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was c...Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.展开更多
In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The e...In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The effects of the divergence, straight, and convergence isolators on the rotating detonation wave dynamics and the upstream oblique shock wave propagation mechanism are analyzed. The differences in the rotating detonation wave behaviors between ground and flight operations are clarified.The results indicate that the propagation regimes of the upstream oblique shock wave depend on the isolator configurations and operation conditions. With a divergence isolator, the airflow is accelerated throughout the isolator and divergence section, leading to a maximum Mach number(~1.8) before the normal shock. The total pressure loss reaches the largest, and the detonation pressure drops. The upstream oblique shock wave can be suppressed within the divergence section with the divergence isolator.However, for the straight and convergence isolators, the airflow in the isolator with a larger ψ_(1)(0.3 and0.4) can suffer from the disturbance of the upstream oblique shock wave. The critical incident angle is around 39° at ground operation conditions. The upstream oblique shock wave tends to be suppressed when the engine operates under flight operation conditions. The critical pressure ratio β_(cr0) is found to be able to help in distinguishing the propagation regimes of the upstream oblique shock wave. Slightly below or above the β_(cr0) can obtain different marginal propagation results. The high-speed airflow in the divergence section affects the fuel droplet penetration distance, which deteriorates the reactant mixing and the detonation area. Significant detonation velocity deficits are observed and the maximum velocity deficit reaches 26%. The results indicate the engine channel design should adopt different isolator configurations based on the purpose of total pressure loss or disturbance suppression. This study can provide useful guidance for the channel design of a more complete two-phase rotating detonation engine.展开更多
The exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine m...The exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine milk exosomes.The milk exosome partition behaviors and bovine milk separation were investigated,and the ATPSs and bovine milk whey addition was optimized.The optimal separation conditions were identified as 16%(mass)polyethylene glycol 4000,10%(mass)dipotassium phosphate,and 1%(mass)enzymatic hydrolysis bovine milk whey.During the separation process,bovine milk exosomes were predominantly enriched in the interphase,while protein impurities were primarily found in the bottom phase.The process yielded bovine milk exosomes of 2.0×10^(11)particles per ml whey with high purity(staining rate>90%,7.01×10^(10)particles per mg protein)and high uniformity(polydispersity index<0.03).The isolated exosomes were characterized and identified by transmission electron microscopy,zeta potential and size distribution.The results demonstrated aqueous two-phase extraction possesses a robust capability for the enrichment and separation of exosomes directly from bovine milk whey,presenting a novel approach for the large-scale isolation of exosomes.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffecti...In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffective purging can lead to crystallization of the molten salt,resulting in blockages.To address this issue,understanding the gas-liquid two-phase flow dynamics during high-pressure gas purging is crucial.This study utilizes the Volume of Fluid(VOF)model and adaptive dynamic grids to simulate the gas-liquid two-phase flow during the purging process in a DN50 PN50 conventional molten salt regulating valve.Initially,the reliability of the CFD simulations is validated through comparisons with experimental data and findings from the literature.Subsequently,simulation experiments are conducted to analyze the effects of various factors,including purge flow rates,initial liquid accumulation masses,purge durations,and the profiles of the valve bottom flow channels.The results indicate that the purging process comprises four distinct stages:Initial violent surge stage,liquid discharge stage,liquid partial fallback stage,liquid dissipation stage.For an initial liquid height of 17 mm at the bottom of the valve,the critical purge flow rate lies between 3 and 5 m/s.Notably,the critical purge flow rate is independent of the initial liquid accumulation mass.As the purge gas flow rate increases,the volume of liquid discharged also increases.Beyond the critical purge flow rate,higher purge gas velocities lead to shorter purge durations.Interestingly,the residual liquid mass after purging remains unaffected by the initial liquid accumulation.Additionally,the flow channel profile at the bottom of the valve significantly influences both the critical purge speed and the efficiency of the purging process.展开更多
Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle ...Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle model rely on solving implicit transcendental equations, which require iterative numerical root-finding methods,thereby introducing computational complexity and inefficiency. This paper proposes the high-precision explicit approximate solutions for the two models, directly correlating the geometric parameters with the flow parameters, thus significantly enhancing the efficiency and accuracy of two-phase flow analysis.展开更多
文摘The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.
基金Supported by the National Natural Science Foundation of China(11171340)
文摘In this work, we obtain the global existence and uniqueness of classical solu-tions to a viscous liquid-gas two-phase model with mass-dependent viscosity and vacuum in one dimension, where the initial vacuum is allowed. We get the upper and lower bounds of gas and liquid masses n and m by the continuity methods which we use to study the compressible Navier-Stokes equations.
基金supported by the National Natural Science Foundation of China(11871341 and 12071152).
文摘In this paper,the Cauchy problem for a two-phase model with a magnetic field in three dimensions is considered.Based on a new linearized system with respect to(c−c_(∞),P−P_(∞),u,H)for constants c_(∞)≥0 and P_(∞)>0,the existence theory of global strong solution is established when the initial data is close to its equilibrium in three dimensions for the small H^(2) initial data.We improve the existence results obtained by Wen and Zhu in[40]where an additional assumption that the initial perturbations are bounded in L^(1)-norm was needed.The energy method combined with the low-frequency and high-frequency decomposition is used to derive the decay of the solution and hence the global existence.As a by-product,the time decay estimates of the solution and its derivatives in the L^(2)-norm are obtained.
基金supported by the National Science and Technology Major Project of China (No.2011ZX04014-052)the National Basic Research Priorities Program of China (No.2011CB012900)
文摘A two-phase model for the prediction of macrosegregation formed during solidification is presented. This model incorporates the descriptions of heat transfer, melt convection, solute transport, and solid movement on the system scale with microscopic relations for grain nucleation and growth. Then the model is used to simulate the solidification of a benchmark industrial 3.3-t steel ingot. Simulations are per- formed to investigate the effects of grain motion and pipe shrinkage formation on the final macrosegregation pattern. The model predictions are compared with experimental data and numerical results from literatures. It is demonstrated that the model is able to express the overall macrosegregation patterns in the ingot. Furthermore, the results show that it is essential to consider the motion of equiaxed grains and the formation of pipe shrinkage in modelling. Several issues for future model improvements are identified.
文摘In this paper, laminar forced convection of CuO nanofluid is numerically investigated in sudden expansion microchannel with isotherm walls and different expansion ratios (ER). An Eulerian two-fluid model is considered to simulate the nanofluid flow inside the microchannel and the governing mass, momentum and energy equations for both phases are solved using the finite volume method. Eulerian-Eulerian two-phase model is very efficient because of considering the relative velocity and temperature of the phases and the nanoparticle concentration distribution. In solving the flow equations for both phases, the SIMPLE algorithm is modified for the coupling of the velocity and pressure and the continuity equations for both phases are combined in order to create the pressure correction equations. However, the Eulerian-Eulerian modeling results show higher heat transfer enhancement in comparison to pure water, so that for a 2% copper-water nanofluid, it has been observed a 35% increase of the heat transfer. The heat transfer enhancement increases with increase in Reynolds number and nanoparticle volume concentration, while the pressure drop increases only slightly. An investigation of the expansion ratio of microchannel shows that the average Nusselt number increases with decrease in expansion ratio as well as with increase in Reynolds number. Also, the Bifurcation has been occurred in higher Reynolds number that is different for each expansion ratio of the microchannel.
基金funded by Natural Science Foundation of China(Grants Nos.51279144 and 11432015)
文摘A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.
基金funded by the National Key R&D Program of China,China(Grant No.2023YFB4005500)National Natural Science Foundation of China,China(Grant Nos.52379113 and 52379114).
文摘Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金supported by the National Key Research and Development Program of China(No.2022YFC3702000)the National Natural Science Foundation of China(No.52070169)the Project of Bureau of Science and Technology of Zhoushan,China(No.2022C41013).
文摘Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mechanisms of mass-transfer enhancement in TPPBs would enable efficient predictions for further industrial applications.In this study,influences of gradually increasing silicone oil ratio on the TPPB was explored,and a 94.35%reduction of the n-hexane partition coefficient was observed with 0.1 vol.%silicone,which increased to 80.7%along with a 40-fold removal efficiency enhancement in the stabilised removal period.The elimination capacity increased from 1.47 to 148.35 g/(m^(3)·h),i.e.a 101-fold increase compared with that of the single-phase reactors,when 10 vol.%(3 Critical Micelle Concentration)silicone oil was added.The significantly promoted partition coefficient was the main reason for the mass transfer enhancement,which covered the negative influences of the decreased total mass-transfer coefficient with increasing silicone oil volume ratio.The gradually rising stirring rate was benefit to the n-hexane removal,which became negative when the dominant resistance shifted from mass transfer to biodegradation.Moreover,a mass-transfer-reaction kinetic model of the TPPB was constructed based on the balance of n-hexane concentration,dissolved oxygen and biomass.Similar to the mechanism,the partition factor was predicted sensitive to the removal performance,and another five sensitive parameters were found simultaneously.This forecasting method enables the optimisation of TPPB performance and provides theoretical support for hydrophobic VOCs degradation.
基金supported by the National Natural Science Foundation of China(Grant No.11972194).
文摘By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using computational fluid dynamics(CFD),the two-phase flow and piston engraving characteristics of a pyrotechnic actuator are investigated.Initially,the current model was utilized to examine the intricate,multi-dimensional flow,and energy conversion characteristics of the propellant grains and combustion gas within the pyrotechnic actuator chamber.It was discovered that the combustion gas on the wall's constant transition from potential to kinetic energy,along with the combined effect of the propellant motion,are what create the pressure oscillation within the chamber.Additionally,a numerical analysis was conducted to determine the impact of various parameters on the pressure oscillation and piston motion,including pyrotechnic charge,pyrotechnic particle size,and chamber structural dimension.The findings show that decreasing the pyrotechnic charge will lower the terminal velocity,while increasing and decreasing the pyrotechnic particle size will reduce the pressure oscillation in the chamber.The pyrotechnic particle size has minimal bearing on the terminal velocity.The results of this investigation offer a trustworthy forecasting instrument for comprehending and creating pyrotechnic actuator designs.
基金Science Foundation for Distinguished Young Scholars 2020-JCJQ-ZQ-042 Intelligent and Bionic Spacecraft Thermal Control Technology Inspired by Tree Sap Transport Principle.
文摘As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it has several limitations when addressing multi-source heat dissipation.Inspired by the transport and heat dissipation system of plants,large trees achieve stable and efficient liquid supply under the influence of two driving forces:capillary force during transpiration in the leaves(pull)and root pressure generated by osmotic pressure in the roots(push).The root pressure provides an effective liquid supply with a driving force exceeding 2 MPa,far greater than the driving force in conventional capillary-pumped two-phase loops.Research has shown that osmotic heat pipes offer a powerful driving force,and combining osmotic pressure with capillary force has significant advantages.Therefore,this paper designs a multi-evaporator,dual-drive two-phase loop,using both osmotic pressure and capillary force to solve the multi-source heat dissipation challenge.First,a transmembrane water flux model for the osmotic pressure-driven device was established to predict the maximum heat transfer capacity of the dual-drive two-phase loop.Then,an experimental setup for a multi-evaporator“osmotic pressure+capillary force”dual-drive two-phase loop was constructed,capable of transferring at least 235 W of power under a reverse gravity condition of 20 m.The study also analyzed the effects of reverse gravity height,heat load distribution among the three evaporators,startup sequence,and varying branch resistances on the performance of the dual-drive two-phase loop.
基金the financial supports of National Key Research and Development Program of China(2021YFB4000100)Yang Yang acknowledges the financial support of Chongqing Basic Science and Advanced Technology Research Program(No.CSTB2022NSCQMSX1676).
文摘The Proton Exchange Membrane(PEM)water electrolyzer is considered one of the promising energy storing means for harnessing variable renewable energy sources to produce hydrogen.Understanding the internal fluid dynamics,which are often challenging to directly observe experimentally,has prompted the use of numerical models to investigate two-phase flow within PEM water electrolyzers.In this study,we provide a comprehensive review of prior research focusing on two-phase modeling of PEM electrolyzers,encompassing both components at mesoscopic scales and the full electrolyzer at the macroscopic level.We delve into the specifics of various modeling approaches for two-phase flow at different scales and summarize and discuss the current state of the art in the field.Presently,two-phase models for the full electrolyzer predominantly employ a macroscopic homogeneous assumption.However,mesoscopic and microscopic models capable of tracking phase interfaces are limited to components.Challenges persist in integrating various modeling scales into a comprehensive electrolyzer model,particularly in coupling two-phase flow between the channels and porous media.Future efforts should focus on developing multi-scale models and simulating two-phase flow under fluctuating input conditions.Additionally,given the structural similarities between PEM water electrolyzers and PEM fuel cells,we compare and discuss differences in two-phase modeling between the two technologies.This work offers the insights for researchers in the field of modeling of PEM water electrolyzers and even fuel cells.
基金supported by the Beijing Municipal Science&Technology Commission(Z231100006123010).
文摘With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.
基金the National Natural Science Foundation of China (Nos. 42302143, 42172159)China Geological Survey Project (No. DD20211350)support from the G. Albert Shoemaker endowment
文摘Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.
基金the financial support from the National Natural Science Foundation of China (No.42102127)the Postdoctoral Research Foundation of China (No.2024 M751860)。
文摘Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.
基金supported by the National Natural Science Foundation of China (Grant No. 12202204)the Natural Science Foundation of Jiangsu Province (Grant No. BK20220953)+1 种基金the Fundamental Research Funds for the Central Universitiesthe Science and Technology Association's Young Talent Nurturing Program of Jiangsu Province (Grant No. JSTJ-2024-004)
文摘In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The effects of the divergence, straight, and convergence isolators on the rotating detonation wave dynamics and the upstream oblique shock wave propagation mechanism are analyzed. The differences in the rotating detonation wave behaviors between ground and flight operations are clarified.The results indicate that the propagation regimes of the upstream oblique shock wave depend on the isolator configurations and operation conditions. With a divergence isolator, the airflow is accelerated throughout the isolator and divergence section, leading to a maximum Mach number(~1.8) before the normal shock. The total pressure loss reaches the largest, and the detonation pressure drops. The upstream oblique shock wave can be suppressed within the divergence section with the divergence isolator.However, for the straight and convergence isolators, the airflow in the isolator with a larger ψ_(1)(0.3 and0.4) can suffer from the disturbance of the upstream oblique shock wave. The critical incident angle is around 39° at ground operation conditions. The upstream oblique shock wave tends to be suppressed when the engine operates under flight operation conditions. The critical pressure ratio β_(cr0) is found to be able to help in distinguishing the propagation regimes of the upstream oblique shock wave. Slightly below or above the β_(cr0) can obtain different marginal propagation results. The high-speed airflow in the divergence section affects the fuel droplet penetration distance, which deteriorates the reactant mixing and the detonation area. Significant detonation velocity deficits are observed and the maximum velocity deficit reaches 26%. The results indicate the engine channel design should adopt different isolator configurations based on the purpose of total pressure loss or disturbance suppression. This study can provide useful guidance for the channel design of a more complete two-phase rotating detonation engine.
基金supported by the National Natural Science Foundation of China(22378350).
文摘The exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine milk exosomes.The milk exosome partition behaviors and bovine milk separation were investigated,and the ATPSs and bovine milk whey addition was optimized.The optimal separation conditions were identified as 16%(mass)polyethylene glycol 4000,10%(mass)dipotassium phosphate,and 1%(mass)enzymatic hydrolysis bovine milk whey.During the separation process,bovine milk exosomes were predominantly enriched in the interphase,while protein impurities were primarily found in the bottom phase.The process yielded bovine milk exosomes of 2.0×10^(11)particles per ml whey with high purity(staining rate>90%,7.01×10^(10)particles per mg protein)and high uniformity(polydispersity index<0.03).The isolated exosomes were characterized and identified by transmission electron microscopy,zeta potential and size distribution.The results demonstrated aqueous two-phase extraction possesses a robust capability for the enrichment and separation of exosomes directly from bovine milk whey,presenting a novel approach for the large-scale isolation of exosomes.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
文摘In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffective purging can lead to crystallization of the molten salt,resulting in blockages.To address this issue,understanding the gas-liquid two-phase flow dynamics during high-pressure gas purging is crucial.This study utilizes the Volume of Fluid(VOF)model and adaptive dynamic grids to simulate the gas-liquid two-phase flow during the purging process in a DN50 PN50 conventional molten salt regulating valve.Initially,the reliability of the CFD simulations is validated through comparisons with experimental data and findings from the literature.Subsequently,simulation experiments are conducted to analyze the effects of various factors,including purge flow rates,initial liquid accumulation masses,purge durations,and the profiles of the valve bottom flow channels.The results indicate that the purging process comprises four distinct stages:Initial violent surge stage,liquid discharge stage,liquid partial fallback stage,liquid dissipation stage.For an initial liquid height of 17 mm at the bottom of the valve,the critical purge flow rate lies between 3 and 5 m/s.Notably,the critical purge flow rate is independent of the initial liquid accumulation mass.As the purge gas flow rate increases,the volume of liquid discharged also increases.Beyond the critical purge flow rate,higher purge gas velocities lead to shorter purge durations.Interestingly,the residual liquid mass after purging remains unaffected by the initial liquid accumulation.Additionally,the flow channel profile at the bottom of the valve significantly influences both the critical purge speed and the efficiency of the purging process.
基金supported by the General Research Fund from the Research Grants Council of the Hong Kong Special Administrative Region of China (No. PolyU 15210624)。
文摘Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle model rely on solving implicit transcendental equations, which require iterative numerical root-finding methods,thereby introducing computational complexity and inefficiency. This paper proposes the high-precision explicit approximate solutions for the two models, directly correlating the geometric parameters with the flow parameters, thus significantly enhancing the efficiency and accuracy of two-phase flow analysis.
