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.展开更多
In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase dra...In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the ra...AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.展开更多
A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kin...A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.展开更多
The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows.The entrainment process significantly influences t...The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows.The entrainment process significantly influences the amplitude of the volume fraction.In this paper,we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth,the solid and fluid volume fractions and the bed morphological evolution.The model employs the Mohr–Coulomb plasticity for the solid stress,and the fluid stress is modeled as a Newtonian viscous stress.The interfacial momentum transfer includes viscous drag and buoyancy.A new extended entrainment rate formula that satisfies the boundary momentum jump condition(Iverson and Ouyang,2015)is presented.In this formula,the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models.A finite volume method using Roe’s Riemann approximation is suggested to solve the equations.Three computational cases are conducted and compared with experiments or previous results.The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.展开更多
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.展开更多
The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study o...The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.展开更多
A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the b...A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.展开更多
The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision ter...The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.展开更多
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 work,the static tensile and free vibration of nanorods are studied via both the strain-driven(Strain D)and stress-driven(Stress D)two-phase nonlocal models with a bi-Helmholtz averaging kernel.Merely adjusting...In this work,the static tensile and free vibration of nanorods are studied via both the strain-driven(Strain D)and stress-driven(Stress D)two-phase nonlocal models with a bi-Helmholtz averaging kernel.Merely adjusting the limits of integration,the integral constitutive equation of the Fredholm type is converted to that of the Volterra type and then solved directly via the Laplace transform technique.The unknown constants can be uniquely determined through the standard boundary conditions and two constrained conditions accompanying the Laplace transform process.In the numerical examples,the bi-Helmholtz kernel-based Strain D(or Stress D)two-phase model shows consistently softening(or stiffening)effects on both the tension and the free vibration of nanorods with different boundary edges.The effects of the two nonlocal parameters of the bi-Helmholtz kernel-based two-phase nonlocal models are studied and compared with those of the Helmholtz kernel-based models.展开更多
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.展开更多
A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carr...A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.展开更多
Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turb...Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.展开更多
RH vacuum degasser is a very important secondary refining device in the production of high quality steels. The flow field of molten steel in RH system plays a significant role in determining productivity of the equipm...RH vacuum degasser is a very important secondary refining device in the production of high quality steels. The flow field of molten steel in RH system plays a significant role in determining productivity of the equipment. The homogeneous model and VOF method were often used to predict the flow field in RH system, but these kinds of models simplified the interaction between gas bubbles and molten steel. In the present work, a numerical model of a whole RH system, including vacuum degasser, immersed legs and ladle,was built based on gas-liquid two-fluid model, and it could be used to analyze the interaction between argon bubbles and molten steel, to understand the effect of the bubble size to the flow field.展开更多
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.展开更多
In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction r...In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.展开更多
基金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.
基金funded by Sichuan Natural Science Foundation Outstanding Youth Science Foundation(No.2024NSFJQ0012)Key project of Regional Innovation and Development Joint Fund of National Natural Science Foundation(No.U23A20669)Sichuan Science and Technology Program(2022ZDZX0041).
文摘In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
基金supported by the National Natural Science Foundation of China(Grant Nos.41404101,41174114,41274130,and 41404102)
文摘AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.
基金the Special Funds for Major State Basic Research of China(G-1999-0222-08)the National Natural Science Foundation of China(50376004)Ph.D.Program Foundation,Ministry of Education of China(20030007028)
文摘A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.
基金Financial support from NSFC(Grant No.41572303,4151001059,41101008)Key Projects in the National Science&Technology Pillar Program(2014BAL05B01)CAS"Light of West China"Program
文摘The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows.The entrainment process significantly influences the amplitude of the volume fraction.In this paper,we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth,the solid and fluid volume fractions and the bed morphological evolution.The model employs the Mohr–Coulomb plasticity for the solid stress,and the fluid stress is modeled as a Newtonian viscous stress.The interfacial momentum transfer includes viscous drag and buoyancy.A new extended entrainment rate formula that satisfies the boundary momentum jump condition(Iverson and Ouyang,2015)is presented.In this formula,the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models.A finite volume method using Roe’s Riemann approximation is suggested to solve the equations.Three computational cases are conducted and compared with experiments or previous results.The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.
基金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.
基金supported by the National Natural Science Foundation of China(11722104,11671150)supported by the National Natural Science Foundation of China(11571280,11331005)+3 种基金supported by the National Natural Science Foundation of China(11331005,11771150)by GDUPS(2016)the Fundamental Research Funds for the Central Universities of China(D2172260)FANEDD No.201315
文摘The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.
基金Supported by the Special Funds for Major State Basic Research Projects, PRC(G1999-0222-08) and the National Natural Science Foundation of China(No. 19872039).
文摘A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.
基金The project supported by the National Natural Science Foundation of China (50176022)
文摘The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.
基金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.
基金the National Natural Science Foundation of China(No.12172169)the China Scholarship Council(CSC)(No.202006830038)。
文摘In this work,the static tensile and free vibration of nanorods are studied via both the strain-driven(Strain D)and stress-driven(Stress D)two-phase nonlocal models with a bi-Helmholtz averaging kernel.Merely adjusting the limits of integration,the integral constitutive equation of the Fredholm type is converted to that of the Volterra type and then solved directly via the Laplace transform technique.The unknown constants can be uniquely determined through the standard boundary conditions and two constrained conditions accompanying the Laplace transform process.In the numerical examples,the bi-Helmholtz kernel-based Strain D(or Stress D)two-phase model shows consistently softening(or stiffening)effects on both the tension and the free vibration of nanorods with different boundary edges.The effects of the two nonlocal parameters of the bi-Helmholtz kernel-based two-phase nonlocal models are studied and compared with those of the Helmholtz kernel-based models.
文摘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.
文摘A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.
文摘Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.
文摘RH vacuum degasser is a very important secondary refining device in the production of high quality steels. The flow field of molten steel in RH system plays a significant role in determining productivity of the equipment. The homogeneous model and VOF method were often used to predict the flow field in RH system, but these kinds of models simplified the interaction between gas bubbles and molten steel. In the present work, a numerical model of a whole RH system, including vacuum degasser, immersed legs and ladle,was built based on gas-liquid two-fluid model, and it could be used to analyze the interaction between argon bubbles and molten steel, to understand the effect of the bubble size to the flow field.
基金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.
基金Graduate Research and Innovation Program in Jiangsu Province(KYZZ16_0286)
文摘In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.
