Hydrate phase transition may pose risks in pipeline blockage and severe challenges for offshore natural gas hydrate pro-duction.The present work involves the development of a multiphase gas-liquid-solid vertical slug ...Hydrate phase transition may pose risks in pipeline blockage and severe challenges for offshore natural gas hydrate pro-duction.The present work involves the development of a multiphase gas-liquid-solid vertical slug flow hydrodynamic model consi-dering hydrate phase transition kinetics with heat and mass transfer behaviors.The varying gas physical properties due to pressure and temperature variations are also introduced to evaluate vertical slug flow characteristics.The proposed model is used to carry out a series of numerical simulations to examine the interactions between hydrate phase transition and vertical slug flow hydrodynamics.Furthermore,the hydrate volumetric fractions under different pressure and temperature conditions are predicted.The results reveal that hydrate formation and gas expansion cause the mixture superficial velocity,and the gas and liquid fractions,void fraction in liq-uid slug,and unit length tend to decrease.The increase in outlet pressure leads to an increased hydrate formation rate,which not only increases the hydrate volumetric fraction along the pipe but also causes the upward shift of the hydrate phase transition critical point.展开更多
Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions;however,its underlying multiphase reactive flow mechanisms remain poorly understood.We conducted steady-state imbibition relative...Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions;however,its underlying multiphase reactive flow mechanisms remain poorly understood.We conducted steady-state imbibition relative permeability experiments on sandstone from a proposed storage site,comple-mented by in situ X-ray imaging and ex situ analyses using scanning electron microscopy(SEM)and energy-dispersive X-ray spectroscopy(EDS).Despite our use of a brine that was pre-equilibrated with CO_(2),there was a significant reduction in both CO_(2) relative permeability and absolute permeability during multiphase flow due to chemical reactions.This reduction was driven by decreased pore and throat sizes,diminished connectivity,and increased irregularity of pore and throat shapes,as revealed by in situ pore-scale imaging.Mineral dissolution,primarily of feldspar,albite,and calcite,along with precipitation resulting from feldspar-to-kaolinite transformation and fines migration,were identified as contributing factors through SEM-EDS analysis.This work provides a benchmark for storage in mineralogically complex sandstones,for which the impact of chemical reactions on multiphase flow properties has been measured.展开更多
Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid ...Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.展开更多
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.展开更多
This study introduces a Transformer-based multimodal fusion framework for simulating multiphase flow and heat transfer in carbon dioxide(CO_(2))–water enhanced geothermal systems(EGS).The model integrates geological ...This study introduces a Transformer-based multimodal fusion framework for simulating multiphase flow and heat transfer in carbon dioxide(CO_(2))–water enhanced geothermal systems(EGS).The model integrates geological parameters,thermal gradients,and control schedules to enable fast and accurate prediction of complex reservoir dynamics.The main contributions are:(i)development of a workflow that couples physics-based reservoir simulation with a Transformer neural network architecture,(ii)design of physics-guided loss functions to enforce conservation of mass and energy,(iii)application of the surrogate model to closed-loop optimization using a differential evolution(DE)algorithm,and(iv)incorporation of economic performance metrics,such as net present value(NPV),into decision support.The proposed framework achieves root mean square error(RMSE)of 3–5%,mean absolute error(MAE)below 4%,and coefficients of determination greater than 0.95 across multiple prediction targets,including production rates,pressure distributions,and temperature fields.When compared with recurrent neural network(RNN)baselines such as gated recurrent units(GRU)and long short-term memory networks(LSTM),as well as a physics-informed reduced-order model,the Transformer-based approach demonstrates superior accuracy and computational efficiency.Optimization experiments further show a 15–20%improvement in NPV,highlighting the framework’s potential for real-time forecasting,optimization,and decision-making in geothermal reservoir engineering.展开更多
In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer...In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.展开更多
The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied...The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied in many industries.The flow of a liquid containing air is a representative example of a multiphase flow and exhibits complex flow characteristics.In particular,the greater the gas volume fraction（GVF）,the more inhomogeneous the flow becomes.As a result,using a venturi meter to measure the rate of a flow that has a high GVF generates an error.In this study,the cause of the error occurred in measuring the flow rate for the multiphase flow when using the venturi meter for analysis by CFD.To ensure the reliability of this study,the accuracy of the multiphase flow models for numerical analysis was verified through comparison between the calculated results of numerical analysis and the experimental data.As a result,the Grace model,which is a multiphase flow model established by an experiment with water and air,was confirmed to have the highest reliability.Finally,the characteristics of the internal flow Held about the multiphase flow analysis result generated by applying the Grace model were analyzed to find the cause of the uncertainty occurring when measuring the flow rate of the multiphase flow using the venturi meter.A phase separation phenomenon occurred due to a density difference of water and air inside the venturi,and flow inhomogeneity happened according to the flow velocity difference of each phase.It was confirmed that this flow inhomogeneity increased as the GVF increased due to the uncertainty of the flow measurement.展开更多
It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conser...It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conservation equations, were established for annular flow with gas hydrate phase transition during gas kick. The behavior of annular multiphase flow with hydrate phase transition was investigated by analyzing the hydrate-forming region, the gas fraction in the fluid flowing in the annulus, pit gain, bottom hole pressure, and shut-in casing pressure. The simulation shows that it is possible to move the hydrate-forming region away from sea floor by increasing the circulation rate. The decrease in gas volume fraction in the annulus due to hydrate formation reduces pit gain, which can delay the detection of well kick and increase the risk of hydrate plugging in lines. Caution is needed when a well is monitored for gas kick at a relatively low gas production rate, because the possibility of hydrate presence is much greater than that at a relatively high production rate. The shut-in casing pressure cannot reflect the gas kick due to hydrate formation, which increases with time.展开更多
Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （...Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （LCS） defined by the ridges of the finite-time Lyapunov exponent （FTLE） is utilized in this study to elucidate the multiphase interactions in gaseous jets injected into water and time-dependent turbu- lent cavitation under the framework of Navier-Stokes flow computations. For the gaseous jets injected into water, the highlighted phenomena of the jet transportation can be observed by the LCS method, including expansion, bulge, necking/breaking, and back-attack. Besides, the observation of the LCS reveals that the back-attack phenomenon arises from the fact that the injected gas has difficulties to move toward downstream re- gion after the necking/breaking. For the turbulent cavitating flow, the ridge of the FTLE field can form a LCS to capture the front and boundary of the re-entraint jet when the ad- verse pressure gradient is strong enough. It represents a bar- rier between particles trapped inside the circulation region and those moving downstream. The results indicate that the FFLE field has the potential to identify the structures of mul- tiphase flows, and the LCS can capture the interface/barrier or the vortex/circulation region.展开更多
The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore ...The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore network model from digital cores at different confining pressures and evaluated the effect of pressure sensitivity on the multiphase displacement process. In both the pore network model and QEMSCAN scanning, the pore structure was observed to be damaged under a high confining pressure. Due to their different scales, the pores and throats exhibited inhomogeneous changes; further, the throats exhibited a significant variation compared to that exhibited by the pores. Meanwhile, the heterogeneity of the pore structure under the two aforementioned activities was aggravated by the elastic-plastic deformation of the pore structure.The pressure-sensitive effect increased the proportion of mineral particles, such as quartz(the main component of the core skeleton), and reduced the proportion of clay minerals. The clay minerals were originally attached to the pore walls or interspersed in the pores; however, as the pressure increased, the clay minerals accumulated in the pores resulting in blockage of the pores. While simulating the multiphase displacement process, increasing the confining pressure was observed to severely restrict the flowability of oil and water. This study promises to improve the efficiency of reservoir development in terms of oil and gas exploitation.展开更多
Direct numerical simulations have recently emerged as a viable tool to study finite Reynolds number multiphase flows. The approach parallels direct numerical simulations of turbulent flows, but the unsteady motion of ...Direct numerical simulations have recently emerged as a viable tool to study finite Reynolds number multiphase flows. The approach parallels direct numerical simulations of turbulent flows, but the unsteady motion of a deformable phase boundary adds considerable complexity. Here, a front tracking method that has been used to study several multiphase flow problems is described. The Navier Stokes equations are solved by a finite difference/front tracking technique that allows the inclusion of fully deformable interfaces and surface tension, in addition to inertial and viscous effects. A parallel version of the method makes it possible to use large grids and resolve flows containing a few hundred bubbles.展开更多
The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian...The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian points to the Eulerian field.The inter-phase exchange terms in the conservation equations are distributed over the volume encompassing the particle size, as opposed to the Dirac delta function generally used in the point-particle approach.The proposed approach is benchmarked against three different flow configurations in a numerical framework based on large eddy simulation（LES） turbulence closure.First, the flow over a circular cylinder is simulated for a Reynolds number of 3900 at 1 atm pressure.Results show good agreement with experimental data for the mean streamwise velocity and the vortex shedding frequency in the wake region.The calculated flowfield exhibits correct physics, which the conventional point-particle approach fails to capture.The second case deals with diesel jet injection in quiescent environment over a pressure range of 1.1–5.0 MPa.The calculated jet penetration depth closely matches measurements.It decreases with increasing chamber pressure, due to enhanced drag force in a denser fluid environment.Finally, water and acetone jet injection normal to air crossflow is studied at1 atm.The calculated jet penetration and Sauter mean diameter of liquid droplets compare very well with measurements.展开更多
This study involved the analysis and characterization of the multiphase flow phenomenon inside the lower stage cyclone separator used in the clinker burning process.The analysis was performed using both CFD and experi...This study involved the analysis and characterization of the multiphase flow phenomenon inside the lower stage cyclone separator used in the clinker burning process.The analysis was performed using both CFD and experimental research methods.Very few studies are devoted to such types of cyclone separators,which in addition to their basic functions are also responsible for the technological process.Due to the atypical working conditions of these cyclone separators,they are characterized with a complex geometry,which significantly differs from that of the traditional separators.Furthermore,the evaluation of the accuracy and level of reliability of the two models of turbulence closure—k-e RNG and RSM(RANS),and the LES.The results obtained led to the conclusion that for the lower stage cyclone separators,the LES model proved to be the most accurate(both in the case of forecasting the separation efficiency and pressure drop).The performance parameter(in particular the separation efficiency)values obtained for the RSM model were also characterized by high accuracy.The k-e RNG model was characterized by significantly larger deviations.展开更多
In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enabl...In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enables the study of multicomponent flows with different molecular weights, different viscosities and different Schmidt numbers, but also ensures that the distribution function of each component evolves on the same square lattice without invoking ad- ditional interpolations. Furthermore, the Chapman-Enskog analysis shows that the present model results in the correct hydrodynamic equations, and satisfies the indifferentiability principle. The numerical validation exercises further demonstrate that the favorable performance of the present model.展开更多
Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars,but many phenomena concerned are not well understood,especially in the internal jet ...Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars,but many phenomena concerned are not well understood,especially in the internal jet flow of the cutting head at the condition of ultra-high pressure.The multiphase flow in the cutting head is numerically simulated to study the abrasive motion mechanism and wear inside the cutting head at the pressure beyond 300 MPa.Visible predictions of the particles trajectories and wear rate in the cutting head are presented.The influences of the abrasive physical properties,size of the jewel orifice and the operating pressure on the trajectories are discussed.Based on the simulation,a wear experiment is carried out under the corresponding pressures.The simulation and experimental results show that the flow in the mixing chamber is composed of the jet core zone and the disturbance zone,both affect the particles trajectories.The mixing efficiency drops with the increase of the abrasive granularity.The abrasive density determines the response of particles to the effects of different flow zones,the abrasive with medium density gives the best general performance.Increasing the operating pressure or using the jewel with a smaller orifice improves the coherency of p articles trajectories but increases the wear rate of the jewel holder at the same time.Walls of the jewel holder,the entrance of the mixing chamber and the convergence part of the mixing tube are subject to wear out.The computational and experimental results give a qualitative consistency which proves that this numerical method can provide a reliable and visible cognition of the flow characteristics of ultra-high pressure abrasive water jet.The investigation is benefit for improving the machining properties of water jet cutting systems and the optimization design of the cutting head.展开更多
A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a...A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.展开更多
Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronauti...Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronautical and astronautical,transportation,hydraulic and nuclear engineering. In this paper,advances and re-search needs in fundamental studies of dispersed multiphase flows,including the particle/droplet/bubble dynamics,particle-particle,droplet-droplet and bubble-bubble interactions,gas-particle and bubble-liquid turbulence interac-tions,particle-wall interaction,numerical simulation of dispersed multiphase flows,including Reynolds-averaged modeling(RANS modeling),large-eddy simulation(LES) and direct numerical simulation(DNS) are reviewed. The research results obtained by the present author are also included in this review.展开更多
Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information whe...Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.展开更多
Tarbela dam is one of the largest earth filled dam in the world used for power generation and irrigation purposes. Like all reservoirs the sediments inflow in the Tarbela reservoir has resulted in reduction in water s...Tarbela dam is one of the largest earth filled dam in the world used for power generation and irrigation purposes. Like all reservoirs the sediments inflow in the Tarbela reservoir has resulted in reduction in water storage capacity and is also causing damage to the tunnels, power generating units and ultimately to the plant equipment. This numerical study was performed to predict the flow patterns and characteristics in Tarbela dam. Tunnel 3 and 4 inlets;originally on the bed level were raised in the 3-D model and meshed. Analysis was performed using multiphase flow (water and air) for maximum inflow in the reservoir, i.e., considering summer season and discharging water through different locations, i.e., tunnels and spillways. Pressure, velocities, flow rate and free surface height results obtained were found in good agreement with the analytical and existing results where available. Results show uneven discharge through each gate due to maximum velocity near exits and overall stagnant phenomena of water within the reservoir. Maximum velocity was observed along the spillways outlet. Strong vortex motion was observed near the spillways outlet and tunnel inlets. New design of Tunnel 3 and 4 were suggested to WAPDA in order to decrease the sediment inflow and improvements in design of the spillways were suggested.展开更多
A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten st...A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.展开更多
基金supported by the National Natural Science Foundation of China(No.52301355)the Natu-ral Science Foundation of Qingdao Municipality(No.23-2-1-108-zyyd-jch)the China University of Petroleum(East China)Independent Innovation Research Project(Science and Engineering)-Leading Talent Cultivation Fund(No.24CX07001A).
文摘Hydrate phase transition may pose risks in pipeline blockage and severe challenges for offshore natural gas hydrate pro-duction.The present work involves the development of a multiphase gas-liquid-solid vertical slug flow hydrodynamic model consi-dering hydrate phase transition kinetics with heat and mass transfer behaviors.The varying gas physical properties due to pressure and temperature variations are also introduced to evaluate vertical slug flow characteristics.The proposed model is used to carry out a series of numerical simulations to examine the interactions between hydrate phase transition and vertical slug flow hydrodynamics.Furthermore,the hydrate volumetric fractions under different pressure and temperature conditions are predicted.The results reveal that hydrate formation and gas expansion cause the mixture superficial velocity,and the gas and liquid fractions,void fraction in liq-uid slug,and unit length tend to decrease.The increase in outlet pressure leads to an increased hydrate formation rate,which not only increases the hydrate volumetric fraction along the pipe but also causes the upward shift of the hydrate phase transition critical point.
文摘Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions;however,its underlying multiphase reactive flow mechanisms remain poorly understood.We conducted steady-state imbibition relative permeability experiments on sandstone from a proposed storage site,comple-mented by in situ X-ray imaging and ex situ analyses using scanning electron microscopy(SEM)and energy-dispersive X-ray spectroscopy(EDS).Despite our use of a brine that was pre-equilibrated with CO_(2),there was a significant reduction in both CO_(2) relative permeability and absolute permeability during multiphase flow due to chemical reactions.This reduction was driven by decreased pore and throat sizes,diminished connectivity,and increased irregularity of pore and throat shapes,as revealed by in situ pore-scale imaging.Mineral dissolution,primarily of feldspar,albite,and calcite,along with precipitation resulting from feldspar-to-kaolinite transformation and fines migration,were identified as contributing factors through SEM-EDS analysis.This work provides a benchmark for storage in mineralogically complex sandstones,for which the impact of chemical reactions on multiphase flow properties has been measured.
基金Supported by the Guangxi Natural Science Foundation(Grant No.2021GXNSFBA196008)the Guangxi Science and Technology Development Program(Grant No.GuikeAD22035189).
文摘Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.
基金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.
文摘This study introduces a Transformer-based multimodal fusion framework for simulating multiphase flow and heat transfer in carbon dioxide(CO_(2))–water enhanced geothermal systems(EGS).The model integrates geological parameters,thermal gradients,and control schedules to enable fast and accurate prediction of complex reservoir dynamics.The main contributions are:(i)development of a workflow that couples physics-based reservoir simulation with a Transformer neural network architecture,(ii)design of physics-guided loss functions to enforce conservation of mass and energy,(iii)application of the surrogate model to closed-loop optimization using a differential evolution(DE)algorithm,and(iv)incorporation of economic performance metrics,such as net present value(NPV),into decision support.The proposed framework achieves root mean square error(RMSE)of 3–5%,mean absolute error(MAE)below 4%,and coefficients of determination greater than 0.95 across multiple prediction targets,including production rates,pressure distributions,and temperature fields.When compared with recurrent neural network(RNN)baselines such as gated recurrent units(GRU)and long short-term memory networks(LSTM),as well as a physics-informed reduced-order model,the Transformer-based approach demonstrates superior accuracy and computational efficiency.Optimization experiments further show a 15–20%improvement in NPV,highlighting the framework’s potential for real-time forecasting,optimization,and decision-making in geothermal reservoir engineering.
基金Project(07JJ4016) supported by the Natural Science Foundation of Hunan Procvince,China
文摘In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.
基金supported by the Industrial Infrastructure Program through The Korea Institute for Advancement of Technology(KIAT) Grant funded by the Korea government Ministry of Trade,Industry and Energy(Grant N0000502)
文摘The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied in many industries.The flow of a liquid containing air is a representative example of a multiphase flow and exhibits complex flow characteristics.In particular,the greater the gas volume fraction（GVF）,the more inhomogeneous the flow becomes.As a result,using a venturi meter to measure the rate of a flow that has a high GVF generates an error.In this study,the cause of the error occurred in measuring the flow rate for the multiphase flow when using the venturi meter for analysis by CFD.To ensure the reliability of this study,the accuracy of the multiphase flow models for numerical analysis was verified through comparison between the calculated results of numerical analysis and the experimental data.As a result,the Grace model,which is a multiphase flow model established by an experiment with water and air,was confirmed to have the highest reliability.Finally,the characteristics of the internal flow Held about the multiphase flow analysis result generated by applying the Grace model were analyzed to find the cause of the uncertainty occurring when measuring the flow rate of the multiphase flow using the venturi meter.A phase separation phenomenon occurred due to a density difference of water and air inside the venturi,and flow inhomogeneity happened according to the flow velocity difference of each phase.It was confirmed that this flow inhomogeneity increased as the GVF increased due to the uncertainty of the flow measurement.
基金supported by the China National 863 Program (Grant No.2006AA09A106)the Doctoral Program of Higher Education of China (Grant No.20060425502)+1 种基金the National Natural Science Foundation of China (Grant No.50874116)Shandong Province Natural Science Foundation(Grant No.Z2007A01)
文摘It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conservation equations, were established for annular flow with gas hydrate phase transition during gas kick. The behavior of annular multiphase flow with hydrate phase transition was investigated by analyzing the hydrate-forming region, the gas fraction in the fluid flowing in the annulus, pit gain, bottom hole pressure, and shut-in casing pressure. The simulation shows that it is possible to move the hydrate-forming region away from sea floor by increasing the circulation rate. The decrease in gas volume fraction in the annulus due to hydrate formation reduces pit gain, which can delay the detection of well kick and increase the risk of hydrate plugging in lines. Caution is needed when a well is monitored for gas kick at a relatively low gas production rate, because the possibility of hydrate presence is much greater than that at a relatively high production rate. The shut-in casing pressure cannot reflect the gas kick due to hydrate formation, which increases with time.
文摘Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （LCS） defined by the ridges of the finite-time Lyapunov exponent （FTLE） is utilized in this study to elucidate the multiphase interactions in gaseous jets injected into water and time-dependent turbu- lent cavitation under the framework of Navier-Stokes flow computations. For the gaseous jets injected into water, the highlighted phenomena of the jet transportation can be observed by the LCS method, including expansion, bulge, necking/breaking, and back-attack. Besides, the observation of the LCS reveals that the back-attack phenomenon arises from the fact that the injected gas has difficulties to move toward downstream re- gion after the necking/breaking. For the turbulent cavitating flow, the ridge of the FTLE field can form a LCS to capture the front and boundary of the re-entraint jet when the ad- verse pressure gradient is strong enough. It represents a bar- rier between particles trapped inside the circulation region and those moving downstream. The results indicate that the FFLE field has the potential to identify the structures of mul- tiphase flows, and the LCS can capture the interface/barrier or the vortex/circulation region.
文摘The pressure-sensitive effect on the pore structure of sandstone was investigated using X-ray computed micro-tomography and QEMSCAN quantitative mineral analysis. In a physical simulation study, we extracted the pore network model from digital cores at different confining pressures and evaluated the effect of pressure sensitivity on the multiphase displacement process. In both the pore network model and QEMSCAN scanning, the pore structure was observed to be damaged under a high confining pressure. Due to their different scales, the pores and throats exhibited inhomogeneous changes; further, the throats exhibited a significant variation compared to that exhibited by the pores. Meanwhile, the heterogeneity of the pore structure under the two aforementioned activities was aggravated by the elastic-plastic deformation of the pore structure.The pressure-sensitive effect increased the proportion of mineral particles, such as quartz(the main component of the core skeleton), and reduced the proportion of clay minerals. The clay minerals were originally attached to the pore walls or interspersed in the pores; however, as the pressure increased, the clay minerals accumulated in the pores resulting in blockage of the pores. While simulating the multiphase displacement process, increasing the confining pressure was observed to severely restrict the flowability of oil and water. This study promises to improve the efficiency of reservoir development in terms of oil and gas exploitation.
文摘Direct numerical simulations have recently emerged as a viable tool to study finite Reynolds number multiphase flows. The approach parallels direct numerical simulations of turbulent flows, but the unsteady motion of a deformable phase boundary adds considerable complexity. Here, a front tracking method that has been used to study several multiphase flow problems is described. The Navier Stokes equations are solved by a finite difference/front tracking technique that allows the inclusion of fully deformable interfaces and surface tension, in addition to inertial and viscous effects. A parallel version of the method makes it possible to use large grids and resolve flows containing a few hundred bubbles.
基金sponsored by the William R.T.Oakes Endowment of the Georgia Institute of Technology
文摘The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian points to the Eulerian field.The inter-phase exchange terms in the conservation equations are distributed over the volume encompassing the particle size, as opposed to the Dirac delta function generally used in the point-particle approach.The proposed approach is benchmarked against three different flow configurations in a numerical framework based on large eddy simulation（LES） turbulence closure.First, the flow over a circular cylinder is simulated for a Reynolds number of 3900 at 1 atm pressure.Results show good agreement with experimental data for the mean streamwise velocity and the vortex shedding frequency in the wake region.The calculated flowfield exhibits correct physics, which the conventional point-particle approach fails to capture.The second case deals with diesel jet injection in quiescent environment over a pressure range of 1.1–5.0 MPa.The calculated jet penetration depth closely matches measurements.It decreases with increasing chamber pressure, due to enhanced drag force in a denser fluid environment.Finally, water and acetone jet injection normal to air crossflow is studied at1 atm.The calculated jet penetration and Sauter mean diameter of liquid droplets compare very well with measurements.
基金the support of the Interdisciplinary Centre for Mathematical and Computational Modelling (ICM) University of Warsaw under grant no G71-5
文摘This study involved the analysis and characterization of the multiphase flow phenomenon inside the lower stage cyclone separator used in the clinker burning process.The analysis was performed using both CFD and experimental research methods.Very few studies are devoted to such types of cyclone separators,which in addition to their basic functions are also responsible for the technological process.Due to the atypical working conditions of these cyclone separators,they are characterized with a complex geometry,which significantly differs from that of the traditional separators.Furthermore,the evaluation of the accuracy and level of reliability of the two models of turbulence closure—k-e RNG and RSM(RANS),and the LES.The results obtained led to the conclusion that for the lower stage cyclone separators,the LES model proved to be the most accurate(both in the case of forecasting the separation efficiency and pressure drop).The performance parameter(in particular the separation efficiency)values obtained for the RSM model were also characterized by high accuracy.The k-e RNG model was characterized by significantly larger deviations.
基金supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (62311)supported by the National Natural Science Foundation of China(51006040)+1 种基金the Hong Kong Scholar Programthe National Science Fund for Distinguished Young Scholars of China (51125024)
文摘In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enables the study of multicomponent flows with different molecular weights, different viscosities and different Schmidt numbers, but also ensures that the distribution function of each component evolves on the same square lattice without invoking ad- ditional interpolations. Furthermore, the Chapman-Enskog analysis shows that the present model results in the correct hydrodynamic equations, and satisfies the indifferentiability principle. The numerical validation exercises further demonstrate that the favorable performance of the present model.
基金supported by National Natural Science Foundation of China(Grant No.50806031)
文摘Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars,but many phenomena concerned are not well understood,especially in the internal jet flow of the cutting head at the condition of ultra-high pressure.The multiphase flow in the cutting head is numerically simulated to study the abrasive motion mechanism and wear inside the cutting head at the pressure beyond 300 MPa.Visible predictions of the particles trajectories and wear rate in the cutting head are presented.The influences of the abrasive physical properties,size of the jewel orifice and the operating pressure on the trajectories are discussed.Based on the simulation,a wear experiment is carried out under the corresponding pressures.The simulation and experimental results show that the flow in the mixing chamber is composed of the jet core zone and the disturbance zone,both affect the particles trajectories.The mixing efficiency drops with the increase of the abrasive granularity.The abrasive density determines the response of particles to the effects of different flow zones,the abrasive with medium density gives the best general performance.Increasing the operating pressure or using the jewel with a smaller orifice improves the coherency of p articles trajectories but increases the wear rate of the jewel holder at the same time.Walls of the jewel holder,the entrance of the mixing chamber and the convergence part of the mixing tube are subject to wear out.The computational and experimental results give a qualitative consistency which proves that this numerical method can provide a reliable and visible cognition of the flow characteristics of ultra-high pressure abrasive water jet.The investigation is benefit for improving the machining properties of water jet cutting systems and the optimization design of the cutting head.
基金financially supported by the National Natural Science Foundation of China(No.51704062)the Fundamental Research Funds for the Central Universities,China(No.N2025019)。
文摘A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.
基金Supported by the Key Projects of National Natural Science Foundation of China （50736006 9587003-13） the State Key Development Program for Basic Research of China （G1999-0222-08） the National Pandeng Project of China （85-06-1-2）
文摘Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronautical and astronautical,transportation,hydraulic and nuclear engineering. In this paper,advances and re-search needs in fundamental studies of dispersed multiphase flows,including the particle/droplet/bubble dynamics,particle-particle,droplet-droplet and bubble-bubble interactions,gas-particle and bubble-liquid turbulence interac-tions,particle-wall interaction,numerical simulation of dispersed multiphase flows,including Reynolds-averaged modeling(RANS modeling),large-eddy simulation(LES) and direct numerical simulation(DNS) are reviewed. The research results obtained by the present author are also included in this review.
基金The authors would like to acknowledge financial support provided by National Natural Science Foundation of China(No.52074338)We are also grateful to the support of the National Key R&D Program of China(No.2019YFA0708700)+1 种基金National Key Basic Research Program of China(20CX06071A)Bin Yuan would like to thank for the support of Shandong Mountain Tai Scholar Program.Chris Clarkson would like to acknowledge funding support from an NSERC Alliance grant(ALLRP 548652-19)for research related to the topic of this paper.
文摘Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.
文摘Tarbela dam is one of the largest earth filled dam in the world used for power generation and irrigation purposes. Like all reservoirs the sediments inflow in the Tarbela reservoir has resulted in reduction in water storage capacity and is also causing damage to the tunnels, power generating units and ultimately to the plant equipment. This numerical study was performed to predict the flow patterns and characteristics in Tarbela dam. Tunnel 3 and 4 inlets;originally on the bed level were raised in the 3-D model and meshed. Analysis was performed using multiphase flow (water and air) for maximum inflow in the reservoir, i.e., considering summer season and discharging water through different locations, i.e., tunnels and spillways. Pressure, velocities, flow rate and free surface height results obtained were found in good agreement with the analytical and existing results where available. Results show uneven discharge through each gate due to maximum velocity near exits and overall stagnant phenomena of water within the reservoir. Maximum velocity was observed along the spillways outlet. Strong vortex motion was observed near the spillways outlet and tunnel inlets. New design of Tunnel 3 and 4 were suggested to WAPDA in order to decrease the sediment inflow and improvements in design of the spillways were suggested.
基金the National Key R&D Program(No.2023YFB3709900)the National Natural Science Foundation of China(Nos.U22A20171 and 52104343)the High Steel Central(HSC)at North China University of Science and Technology and Yanshan Univ ersity,China。
文摘A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.
