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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectan...The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.展开更多
A phase-field-based lattice Boltzmann model is proposed for the interface capturing of multi-phase flows based on the conservative Allen–Cahn equation(ACE).By adopting the improved form of a relaxation matrix and an ...A phase-field-based lattice Boltzmann model is proposed for the interface capturing of multi-phase flows based on the conservative Allen–Cahn equation(ACE).By adopting the improved form of a relaxation matrix and an equilibrium distribution function,the time derivative?t(φ)induced by recovering the diffusion term in ACE is eliminated.The conducted Chapman–Enskog analysis demonstrates that the correct conservative ACE is recovered.Four benchmark cases including Zalesak’s disk rotation,vortex droplet,droplet impact on thin film,and Rayleigh–Taylor instability are investigated to validate the proposed model.The numerical results indicate that the proposed model can accurately describe the complex interface deformation.展开更多
Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different des...Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different descriptions of the flow depending on the scale that we are examining. Here we discuss multiphase flows from a multiscale perspective. Those include both how DNS are providing insight and understanding for modeling of scales much larger than the ＂dominant scale＂ （defined where surface tension, viscous forces or inertia are important）, as well as how DNS are often limited by the need to resolve processes taking place on much smaller scales. Both problems can be cast into a language introduced for general classes of multiscale problems and reveal that while the classification may be new, the issues are not.展开更多
The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carre...The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.展开更多
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.展开更多
In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high ...In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.展开更多
A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We ...A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell.The modification procedure can be used on the GFM with any base schemes.In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization.The level set method is used to capture the interface.Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.展开更多
This study investigates the effects of multiphase internal flows that consider hydrate phase transitions on the parametric stability of marine risers.A numerical model of the multiphase internal flow that considers a ...This study investigates the effects of multiphase internal flows that consider hydrate phase transitions on the parametric stability of marine risers.A numerical model of the multiphase internal flow that considers a hydrate phase transition is established.The model first solves the flow parameters and subsequently obtains the natural frequencies of risers with different gas intake ratios.The stability charts of marine risers with different gas intake ratios are plotted by applying Floquet theory,and the effects of the gas intake ratio on the instability and vibration response of the risers are identified.The natural frequency increases with an increase in the gas intake ratio;thus,instability zones move to higher frequency ranges in the stability charts.As the increasing gas intake ratio reduces the damping effect of the Coriolis force,the critical amplitude of the heave in the unstable region decreases,especially when hydrodynamic damping is not considered.As a result,higher-order unstable regions are excited.When in an unstable region,the vibration response curve of a riser with a high gas intake ratio excited by parametric resonance diverges quickly due to parametric resonance.展开更多
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.展开更多
This paper aims to tackle the calculation efficiency problem raised in the cavitation-flow simulation of the aviation centrifugal pump due to the fading-away interface resulting from the dissipation of numerics used i...This paper aims to tackle the calculation efficiency problem raised in the cavitation-flow simulation of the aviation centrifugal pump due to the fading-away interface resulting from the dissipation of numerics used in the phase-change control equation for unstructured-grid multiphase flow,and due to the limitation of flow time-step in whole flow regimes,the control equation of vapor–liquid two-phase flow considering cavitation mass transport is established firstly,modifying the momentum equation by introducing the surface tension,and adding the artificial convective flow to the phase equation to solve the numerical dissipation problem.Secondly,in consideration of the local time step principle and based on the multi-dimensional general limiter algorithm with explicit solutions under the OpenFOAM platform,a solution method of steady-state VOF (Volume of Fluid) model considering cavitation two-phase change is constructed,and the feasibility of this method is verified by NACA hydrofoil and NASA flat plate inducer.Finally,based on the platform developed,the cavitation performance of an aviation centrifugal pump inducer is analyzed.The research results show that the error of the calculated cavitation pressure distribution for NACA hydrofoil between the simulation test and the experimental-test is less than 5%,and the maximum error of calculated cavitation number at pump head dropping for NASA high-speed flat plate inducer between the simulation test and the experimental-test is 2.1%.The cavitation area observed in the simulation test is the same as that obtained in the high-speed photography test.Based on the OpenFOAM simulation method,the position of pump head dropping of the fuel centrifugal pump can be accurately captured.The error of the calculated cavitation number at pump head dropping between the simulation test and the experimental test is about 3.7%,showing high calculation accuracy.展开更多
A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubb...A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubbles,and the effects of different temperatures,anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.%butanedinitrile alloy during the solidification process.The model adopts a multiphase field model to simulate the growth of dendrites,calculates the growth motions of dendrites based on the interfacial solute equilibrium;and adopts a lattice Boltzmann model(LBM)based on the Shan-Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase,which includes the interaction between solid-liquid-gas phases.The simulation results show that during the directional growth of columnar dendrites,bubbles first precipitate out slowly at the very bottom of the dendrites,and then rise up due to the different solid-liquid densities and pressure differences.The bubbles will interact with the dendrite in the process of flow migration,such as extrusion,overflow,fusion and disappearance.In the case of wide gaps in the dendrite channels,bubbles will fuse to form larger irregular bubbles,and in the case of dense channels,bubbles will deform due to the extrusion of dendrites.In the simulated region,as the dendrites converge and diverge,the bubbles precipitate out of the dendrites by compression and diffusion,which also causes physical phenomena such as fusion and spillage of the bubbles.These results reveal the physical mechanisms of bubble nucleation,growth and kinematic evolution during solidification and interaction with dendrite growth.展开更多
The present simulation investigates the multiphase cavitating flow around an underwater projectile.Based on the Homogeneous Equilibrium Flow assumption,a mixture model is applied to simulate the multiphase cavitating ...The present simulation investigates the multiphase cavitating flow around an underwater projectile.Based on the Homogeneous Equilibrium Flow assumption,a mixture model is applied to simulate the multiphase cavitating flow including ventilated cavitation caused by air injection as well as natural cavitation that forms in a region where the pressure of liquid falls below its vapor pressure. The transport equation cavitating model is applied.The calculations are executed based on a suite of CFD code.The hydrodynamics characteristics of flow field under the interaction of natural cavitation and ventilated cavitation is analyzed.The results indicate that the ventilated cavitation number is under a combined effect of the natural cavitation number and gas flow rate in the multiphase cavitating flows.展开更多
During perfusion culture,the growth of bone tissues in the scaffold was closely related to the locations of initial adhered cells and their density.In this study,the fluid mechanical responses of Voronoi-lattice scaff...During perfusion culture,the growth of bone tissues in the scaffold was closely related to the locations of initial adhered cells and their density.In this study,the fluid mechanical responses of Voronoi-lattice scaffolds and initial adhered cells on scaffolds were quantitatively investigated.Multiphase fluid-structure interaction(FSI)model was verified by comparing with the results of Diamond scaffolds culture in the literature.Fluid mechanical responses of Voronoi-lattice scaffolds and cells were analyzed by multiphase FSI model.Regression equations were established by response surface method(RSM)to determine relationships between structural design factors of Voronoi-lattice scaffolds and fluid mechanical response parameters of scaffolds and cells.The results showed that the percentage of adhered cells and the locations of initial adhered cells obtained by multiphase FSI model of Diamond scaffolds had the same trend with that obtained by perfusion culture.Regression equations established based on RSM could well predict the fluid mechanical response parameters of Voronoi-scaffolds and cells.The multiphase FSI model closely related the densities of cells and the locations of adhered cells to bone tissue growth.The model could provide a certain theoretical basis for constructing and culturing engineered bone tissues in vitro perfusion.展开更多
Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettabi...Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettability alteration by surfactants. Although the interacting capillary bundle(ICB) model shows potential in characterizing imbibition rates in different pores during wettability alteration, the existing ICB models neglect the influence of wettability and viscosity ratio on the imbibition behavior, making it difficult to accurately describe the oil-water imbibition behavior within the porous media. In this work,a new ICB mathematical model is established by introducing pressure balance without assuming the position of the leading front to comprehensively describe the imbibition behavior in a porous medium under different conditions, including gas-liquid spontaneous imbibition and oil-water imbibition.When the pore size distribution of a tight rock is known, this new model can predict the changes of water saturation during the displacement process in the tight rock, and also determine the imbibition rate in pores of different sizes. The water saturation profiles obtained from the new model are validated against the waterflooding simulation results from the CMG, while the imbibition rates calculated by the model are validated against the experimental observations of gas-liquid spontaneous imbibition. The good match above indicates the newly proposed model can show the water saturation profile at a macroscopic scale while capture the underlying physics of the multiphase flow in a porous medium at a microscopic scale. Simulation results obtained from this model indicate that both wettability and viscosity ratio can affect the sequence of fluid imbibition into pores of different sizes during the multiphase flow, where less-viscous wetting fluid is preferentially imbibed into larger pores while more-viscous wetting fluid tends to be imbibed into smaller pores. Furthermore, this model provides an avenue to calculate the imbibition rate in pores of different sizes during wettability alteration and capture the non-Darcy effect in micro-and nano-scale pores.展开更多
基金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.
基金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.
文摘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.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(No.11425210)the Fundamental Research Funds for the Central Universities(No.WK2090050025)
文摘The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFB1504301)the Australian Research Council(Grant No.DE160101098)。
文摘A phase-field-based lattice Boltzmann model is proposed for the interface capturing of multi-phase flows based on the conservative Allen–Cahn equation(ACE).By adopting the improved form of a relaxation matrix and an equilibrium distribution function,the time derivative?t(φ)induced by recovering the diffusion term in ACE is eliminated.The conducted Chapman–Enskog analysis demonstrates that the correct conservative ACE is recovered.Four benchmark cases including Zalesak’s disk rotation,vortex droplet,droplet impact on thin film,and Rayleigh–Taylor instability are investigated to validate the proposed model.The numerical results indicate that the proposed model can accurately describe the complex interface deformation.
文摘Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different descriptions of the flow depending on the scale that we are examining. Here we discuss multiphase flows from a multiscale perspective. Those include both how DNS are providing insight and understanding for modeling of scales much larger than the ＂dominant scale＂ （defined where surface tension, viscous forces or inertia are important）, as well as how DNS are often limited by the need to resolve processes taking place on much smaller scales. Both problems can be cast into a language introduced for general classes of multiscale problems and reveal that while the classification may be new, the issues are not.
文摘The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.
基金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.
文摘In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.
基金supported by NSFC 10531080,10972230,and 973 project 2005CB321703 and 2010CB731505supported by ARO grant W911NF-08-1-0520 and NSF grant DMS-0809086.
文摘A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell.The modification procedure can be used on the GFM with any base schemes.In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization.The level set method is used to capture the interface.Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.
基金funded by the National Natural Sci-ence Foundation of China(No.U2006226)the National Key Research and Development Program of China(No.2016YFC0303800)the National Natural Science Foundation of China(No.51579245)。
文摘This study investigates the effects of multiphase internal flows that consider hydrate phase transitions on the parametric stability of marine risers.A numerical model of the multiphase internal flow that considers a hydrate phase transition is established.The model first solves the flow parameters and subsequently obtains the natural frequencies of risers with different gas intake ratios.The stability charts of marine risers with different gas intake ratios are plotted by applying Floquet theory,and the effects of the gas intake ratio on the instability and vibration response of the risers are identified.The natural frequency increases with an increase in the gas intake ratio;thus,instability zones move to higher frequency ranges in the stability charts.As the increasing gas intake ratio reduces the damping effect of the Coriolis force,the critical amplitude of the heave in the unstable region decreases,especially when hydrodynamic damping is not considered.As a result,higher-order unstable regions are excited.When in an unstable region,the vibration response curve of a riser with a high gas intake ratio excited by parametric resonance diverges quickly due to parametric resonance.
基金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.
基金supported by the National Science and Technology Major Project, China (No. J2019-V-0016-0111)the Aviation Engine and Gas Turbine Basic Science Center Project, China (No. P2022-B-V-003-001)+3 种基金the Defense Industrial Technology Development Program, China (No. JCKY2022607C002)the AECC Industry University Cooperation Project, China (No. HFZL2022CXY013)the National Natural Science Foundation of China (No. 52372396)the Key R&D Project in Shaanxi Province, China (No. 2021GXLH-01-16)。
文摘This paper aims to tackle the calculation efficiency problem raised in the cavitation-flow simulation of the aviation centrifugal pump due to the fading-away interface resulting from the dissipation of numerics used in the phase-change control equation for unstructured-grid multiphase flow,and due to the limitation of flow time-step in whole flow regimes,the control equation of vapor–liquid two-phase flow considering cavitation mass transport is established firstly,modifying the momentum equation by introducing the surface tension,and adding the artificial convective flow to the phase equation to solve the numerical dissipation problem.Secondly,in consideration of the local time step principle and based on the multi-dimensional general limiter algorithm with explicit solutions under the OpenFOAM platform,a solution method of steady-state VOF (Volume of Fluid) model considering cavitation two-phase change is constructed,and the feasibility of this method is verified by NACA hydrofoil and NASA flat plate inducer.Finally,based on the platform developed,the cavitation performance of an aviation centrifugal pump inducer is analyzed.The research results show that the error of the calculated cavitation pressure distribution for NACA hydrofoil between the simulation test and the experimental-test is less than 5%,and the maximum error of calculated cavitation number at pump head dropping for NASA high-speed flat plate inducer between the simulation test and the experimental-test is 2.1%.The cavitation area observed in the simulation test is the same as that obtained in the high-speed photography test.Based on the OpenFOAM simulation method,the position of pump head dropping of the fuel centrifugal pump can be accurately captured.The error of the calculated cavitation number at pump head dropping between the simulation test and the experimental test is about 3.7%,showing high calculation accuracy.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52161002,51661020,and 11364024)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology of China(Grant No.J201304).
文摘A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubbles,and the effects of different temperatures,anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.%butanedinitrile alloy during the solidification process.The model adopts a multiphase field model to simulate the growth of dendrites,calculates the growth motions of dendrites based on the interfacial solute equilibrium;and adopts a lattice Boltzmann model(LBM)based on the Shan-Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase,which includes the interaction between solid-liquid-gas phases.The simulation results show that during the directional growth of columnar dendrites,bubbles first precipitate out slowly at the very bottom of the dendrites,and then rise up due to the different solid-liquid densities and pressure differences.The bubbles will interact with the dendrite in the process of flow migration,such as extrusion,overflow,fusion and disappearance.In the case of wide gaps in the dendrite channels,bubbles will fuse to form larger irregular bubbles,and in the case of dense channels,bubbles will deform due to the extrusion of dendrites.In the simulated region,as the dendrites converge and diverge,the bubbles precipitate out of the dendrites by compression and diffusion,which also causes physical phenomena such as fusion and spillage of the bubbles.These results reveal the physical mechanisms of bubble nucleation,growth and kinematic evolution during solidification and interaction with dendrite growth.
基金supported by the National Natural Science Foundation of China(Grant No:10832007)Shanghai Leading Academic Discipline Project(Project No. B206)
文摘The present simulation investigates the multiphase cavitating flow around an underwater projectile.Based on the Homogeneous Equilibrium Flow assumption,a mixture model is applied to simulate the multiphase cavitating flow including ventilated cavitation caused by air injection as well as natural cavitation that forms in a region where the pressure of liquid falls below its vapor pressure. The transport equation cavitating model is applied.The calculations are executed based on a suite of CFD code.The hydrodynamics characteristics of flow field under the interaction of natural cavitation and ventilated cavitation is analyzed.The results indicate that the ventilated cavitation number is under a combined effect of the natural cavitation number and gas flow rate in the multiphase cavitating flows.
基金supported by the National Natural Science Foundation of China(Grant No.12272029).
文摘During perfusion culture,the growth of bone tissues in the scaffold was closely related to the locations of initial adhered cells and their density.In this study,the fluid mechanical responses of Voronoi-lattice scaffolds and initial adhered cells on scaffolds were quantitatively investigated.Multiphase fluid-structure interaction(FSI)model was verified by comparing with the results of Diamond scaffolds culture in the literature.Fluid mechanical responses of Voronoi-lattice scaffolds and cells were analyzed by multiphase FSI model.Regression equations were established by response surface method(RSM)to determine relationships between structural design factors of Voronoi-lattice scaffolds and fluid mechanical response parameters of scaffolds and cells.The results showed that the percentage of adhered cells and the locations of initial adhered cells obtained by multiphase FSI model of Diamond scaffolds had the same trend with that obtained by perfusion culture.Regression equations established based on RSM could well predict the fluid mechanical response parameters of Voronoi-scaffolds and cells.The multiphase FSI model closely related the densities of cells and the locations of adhered cells to bone tissue growth.The model could provide a certain theoretical basis for constructing and culturing engineered bone tissues in vitro perfusion.
基金financially supported by the General Program Grant from the National Natural Science Foundation of China(52274051 and 52174045)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-01)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(51521063)。
文摘Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettability alteration by surfactants. Although the interacting capillary bundle(ICB) model shows potential in characterizing imbibition rates in different pores during wettability alteration, the existing ICB models neglect the influence of wettability and viscosity ratio on the imbibition behavior, making it difficult to accurately describe the oil-water imbibition behavior within the porous media. In this work,a new ICB mathematical model is established by introducing pressure balance without assuming the position of the leading front to comprehensively describe the imbibition behavior in a porous medium under different conditions, including gas-liquid spontaneous imbibition and oil-water imbibition.When the pore size distribution of a tight rock is known, this new model can predict the changes of water saturation during the displacement process in the tight rock, and also determine the imbibition rate in pores of different sizes. The water saturation profiles obtained from the new model are validated against the waterflooding simulation results from the CMG, while the imbibition rates calculated by the model are validated against the experimental observations of gas-liquid spontaneous imbibition. The good match above indicates the newly proposed model can show the water saturation profile at a macroscopic scale while capture the underlying physics of the multiphase flow in a porous medium at a microscopic scale. Simulation results obtained from this model indicate that both wettability and viscosity ratio can affect the sequence of fluid imbibition into pores of different sizes during the multiphase flow, where less-viscous wetting fluid is preferentially imbibed into larger pores while more-viscous wetting fluid tends to be imbibed into smaller pores. Furthermore, this model provides an avenue to calculate the imbibition rate in pores of different sizes during wettability alteration and capture the non-Darcy effect in micro-and nano-scale pores.
