The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated includ...The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to des...Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.展开更多
This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unstead...This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unsteady axisymmetric Navier-Stokes equations for both the drop and the ambient gas flow. The level set method is applied to capturing the liquid/gas interface. Through calculation are obtained four typical breakup modes--oscillation, bag breakup, sheet stripping breakup and shear breakup governed by four non-dimensional numbers which are gas Weber number （Weg）, liquid Reynolds number （Rel）, gas Reynolds number （Reg） and density ratio （γ）. Their effects upon each mode are analyzed. The results indicate that among the four numbers, Weg is of the highest importance with Rel, Reg and γfollowing up. By widening the range of the density ratio up to 1 000, the breakup mode is discovered to be so complicated that a new one called multimode breakup mode turns up. This mode contains the shearing breakup and piercing breakup, which successively happen. The calculation results agree well with what is observed from the experiments.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic s...A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic simulation.The validity of the numerical method is verified by a benchmark problem.The melt interface evolution versus time is captured and the physical quantities such as temperature,velocity and pressure at each time step are obtained with corresponding analysis.A"frozen skin"layer with the thickness increasing versus time during the injection process is found.The fact that the"frozen skin"layer can be reduced by increasing the injection velocity is numerically verified.The fountain flow phenomenon near the melt interface is also captured.Moreover,comparisons with the non-isothermal Newtonian case show that the curvatures of the interface arcs and the pressure contours near the horizontal mid-line of the cavity for the non-Newtonian pseudoplastic case is larger than that for the Newtonian case.The velocity profiles are different at different positions for the non-Newtonian pseudoplastic case,while in the case of Newtonian flow the velocity profiles are parabolic and almost the same at different positions.展开更多
Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous me...Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.展开更多
The mathematical model of mass transfer-induced Marangoni effect is formulated. The drop surface evolution is captured by the level set method, in which the interface is represented by the embedded set of zero level o...The mathematical model of mass transfer-induced Marangoni effect is formulated. The drop surface evolution is captured by the level set method, in which the interface is represented by the embedded set of zero level of a scalar distance function defined in the whole computational domain. Numerical simulation of the Marangoni effect induced by interphase mass transfer to/from deformable single drops in unsteady motion in liquid-liquid extraction systems is performed in a Eulerian axisymmetric reference frame. The occurrence and development of the Marangoni effect are simulated, and the re- sults are in good agreement with the classical theoretical analysis and previous simulation.展开更多
Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eul...Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.展开更多
A series of 2D direct numerical simulations were performed with an accurate level set method for single drop impacts.The adopted ACLS method was validated to be efficient with perfect mass conservation in both normal ...A series of 2D direct numerical simulations were performed with an accurate level set method for single drop impacts.The adopted ACLS method was validated to be efficient with perfect mass conservation in both normal and oblique impacts.A square-root correction for neck bases was modified in accuracy as well as scope of applications.In addition,process of jet formation and evolution was studied to reveal internal dynamics in drop impacts.It's found that pressure gradient and vortex are coexisting and completive reasons for jet topology while the inclined angle has a significant effect on them.Mechanisms of jet formation and evolution are different in the front and back necks.With the help of PDF distribution and correction calculation,a compromise in the competition is observed.This work lays a solid foundation for further studies of dynamics in gas-liquid flows.展开更多
Cold atmospheric plasmas(CAPs)have attracted considerable interest in the field of plasma medicine.Generated reactive species such as hydroxyl(OH)species play an important role in applications of CAPs.Transportation o...Cold atmospheric plasmas(CAPs)have attracted considerable interest in the field of plasma medicine.Generated reactive species such as hydroxyl(OH)species play an important role in applications of CAPs.Transportation of OH species towards the target and distribution of these OH species in the plasma plume play an important role in the applications of plasma medicine.In the present work,a computational model was built to simulate the transportation and distribution of OH species in CAP discharges,which was based on the level set method to dynamically track the propagation of plasma carrier gas in air.A reaction term was incorporated for the OH species.The OH species tended to diffuse around the main stream of the carrier gas,and thus covered larger radial and axial distances.A CAP discharge onto a skin layer led to the largest accumulation of OH species at the central part of the exposed area.The distribution of OH species on the skin was asymmetric,which agreed with experiments.The computational model itself and the obtained results would be useful for future development of plasma medicine.展开更多
Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was ...Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was first validated and improved by several test cases of a still droplet.The five distinct flow regimes of the kerosene-water system,previously identified in the experiments from Zhao et al.,were reproduced.The quantitative and qualitative agreement between the simulations and the experimental data show the effectiveness of the numerical method.The roles of the interfacial tension and contact angle on the flow patterns and shapes of droplets were discussed and highlighted according to the numerical results based on the improved two-phase LB model.This work demonstrated that the developed LBM simulator is a viable tool to study immiscible two-phase flows in microchannels,and such a tool could provide tangible guidance for the design of various microfluidic devices that involve immiscible multi-phase flows.展开更多
A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Sto...A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.展开更多
We propose an efficient numerical method for the simulation of the twophase flows with moving contact lines in three dimensions.The mathematical model consists of the incompressible Navier-Stokes equations for the two...We propose an efficient numerical method for the simulation of the twophase flows with moving contact lines in three dimensions.The mathematical model consists of the incompressible Navier-Stokes equations for the two immiscible fluids with the standard interface conditions,the Navier slip condition along the solid wall,and a contact angle condition(Ren et al.(2010)[28]).In the numerical method,the governing equations for the fluid dynamics are coupledwith an advection equation for a level-set function.The latter models the dynamics of the fluid interface.Following the standard practice,the interface conditions are taken into account by introducing a singular force on the interface in themomentum equation.This results in a single set of governing equations in the whole fluid domain.Similarly,the contact angle condition is imposed by introducing a singular force,which acts in the normal direction of the contact line,into theNavier slip condition.The newboundary condition,which unifies the Navier slip condition and the contact angle condition,is imposed along the solid wall.The model is solved using the finite difference method.Numerical results are presented for the spreading of a droplet on both homogeneous and inhomogeneous solid walls,as well as the dynamics of a droplet on an inclined plate under gravity.展开更多
The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses...The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses.At first,the strong governing formulation of these microstructures is briefly constructed,and the relevant weak formulation is then deduced based on the virtual work theorem.Afterwards,a subdividing approach,which cuts the elements intercepted by the interfaces between distinct phases automatically,is established with the aid of the level set method(LSM),and the discrete counterpart of the governing formula is obtained by combining the weak formulation derived and a discretized model.To be noted,two silent merits are found when the elaborated strategy is applied:(1) the continents constituting the microstructures of different scales can be arbitrarily-shaped and conveniently reproduced;(2) the periodic boundary condition commonly employed can be enforced on the external surfaces of representative unit cells(RUCs) with no difficulty.Besides,a boundary value problem(BVP) involving a simplified HTE nanostructure is designed,analytically solved,and hereafter applied to verify the correctness of the proposed strategy.It is observed that both the displacement and stress predictions by the computational approach are in good agreement with the relevant analytical results irrespective of the material combinations applied.Eventually,discussions are made on the influence of material organizations of both the 2D and 3D HTE microstructures at the ultrastructural and repeated rod levels,and some concluding remarks are drawn.展开更多
Electrowetting has been proposed as a technique for manipulating dropletssurrounded by air or oil. In this paper, we discuss the modeling and simulation of thedroplet fission process between two parallel plates inside...Electrowetting has been proposed as a technique for manipulating dropletssurrounded by air or oil. In this paper, we discuss the modeling and simulation of thedroplet fission process between two parallel plates inside an electrowetting on dielectric (EWOD) device. Since the gap between the plates is small, we use the two-phaseHele-Shaw flow as a model. While there are several high order methods around, suchas the immersed interface methods [1, 2], we decide to use two first-order methods forsimplicity. A ghost-fluid (GF) method is employed to solve the governing equationsand a local level set method is used to track the drop interface. For comparison purposes, the same set of two-phase Hele-Shaw equations are also solved directly usingthe immersed boundary (IB) method. Numerical results are consistent with experimental observations reported in the literature.展开更多
We propose a new semi-implicit level set approach to a class of curvature dependent flows.The method generalizes a recent algorithm proposed for the motion by mean curvature where the interface is updated by solving t...We propose a new semi-implicit level set approach to a class of curvature dependent flows.The method generalizes a recent algorithm proposed for the motion by mean curvature where the interface is updated by solving the Rudin-Osher-Fatemi(ROF)model for image regularization.Our proposal is general enough so that one can easily extend and apply the method to other curvature dependent motions.Since the derivation is based on a semi-implicit time discretization,this suggests that the numerical scheme is stable even using a time-step significantly larger than that of the corresponding explicit method.As an interesting application of the numerical approach,we propose a new variational approach for extracting limit cycles in dynamical systems.The resulting algorithm can automatically detect multiple limit cycles staying inside the initial guess with no condition imposed on the number nor the location of the limit cycles.Further,we also propose in this work an Eulerian approach based on the level set method to test if the limit cycles are stable or unstable.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
基金Project(21406141)supported by the National Natural Science Foundation of ChinaProject(20141078)supported by the Scientific Research Starting Foundation for Doctors of Liaoning Province,China+1 种基金Project(L2014060)supported by the Foundation of Department of Education of Liaoning Province,ChinaProject(157B21)supported by the Scientific Research Starting Foundation for Doctors of Shenyang Aerospace University,China
文摘The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
基金National High Technology Research and Development Program of China (863program) (2006AA04Z140)National Natural Science Foundation of China (NSFC) (50605024)
文摘Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.
文摘This article studies numerically a familiar important phenomenon in spray combustion which is deformation and breakup of liquid drops in gas flow. The SIMPLER method is used to solve the two-dimensional （2D） unsteady axisymmetric Navier-Stokes equations for both the drop and the ambient gas flow. The level set method is applied to capturing the liquid/gas interface. Through calculation are obtained four typical breakup modes--oscillation, bag breakup, sheet stripping breakup and shear breakup governed by four non-dimensional numbers which are gas Weber number （Weg）, liquid Reynolds number （Rel）, gas Reynolds number （Reg） and density ratio （γ）. Their effects upon each mode are analyzed. The results indicate that among the four numbers, Weg is of the highest importance with Rel, Reg and γfollowing up. By widening the range of the density ratio up to 1 000, the breakup mode is discovered to be so complicated that a new one called multimode breakup mode turns up. This mode contains the shearing breakup and piercing breakup, which successively happen. The calculation results agree well with what is observed from the experiments.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
基金Supported by the National Natural Science Foundation of China（10871159） the National Basic Research Program of China（2005CB321704）
文摘A non-isothermal injection molding process for a non-Newtonian viscous pseudoplastic fluid is simulated.A conservative interface capturing technique and the flow field solving method are coupled to perform a dynamic simulation.The validity of the numerical method is verified by a benchmark problem.The melt interface evolution versus time is captured and the physical quantities such as temperature,velocity and pressure at each time step are obtained with corresponding analysis.A"frozen skin"layer with the thickness increasing versus time during the injection process is found.The fact that the"frozen skin"layer can be reduced by increasing the injection velocity is numerically verified.The fountain flow phenomenon near the melt interface is also captured.Moreover,comparisons with the non-isothermal Newtonian case show that the curvatures of the interface arcs and the pressure contours near the horizontal mid-line of the cavity for the non-Newtonian pseudoplastic case is larger than that for the Newtonian case.The velocity profiles are different at different positions for the non-Newtonian pseudoplastic case,while in the case of Newtonian flow the velocity profiles are parabolic and almost the same at different positions.
基金funded by Key Laboratory of Coal-based CO_2 Capture and Geological Storage,Jiangsu Province,ChinaUS Advanced Coal Technology Consortium(No.2013 DFB60140-08)
文摘Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.
基金Supported by the National Natural Science Foundation of China (Grant Nos. 20490206 & 20576133, 50404009)PetroChina and the National Basic Research Program of China (Grant Nos. 2004CB619203, 2004CB619203)
文摘The mathematical model of mass transfer-induced Marangoni effect is formulated. The drop surface evolution is captured by the level set method, in which the interface is represented by the embedded set of zero level of a scalar distance function defined in the whole computational domain. Numerical simulation of the Marangoni effect induced by interphase mass transfer to/from deformable single drops in unsteady motion in liquid-liquid extraction systems is performed in a Eulerian axisymmetric reference frame. The occurrence and development of the Marangoni effect are simulated, and the re- sults are in good agreement with the classical theoretical analysis and previous simulation.
文摘Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.
基金Supported by the National Natural Science Foundation of China(91541202,51276163)
文摘A series of 2D direct numerical simulations were performed with an accurate level set method for single drop impacts.The adopted ACLS method was validated to be efficient with perfect mass conservation in both normal and oblique impacts.A square-root correction for neck bases was modified in accuracy as well as scope of applications.In addition,process of jet formation and evolution was studied to reveal internal dynamics in drop impacts.It's found that pressure gradient and vortex are coexisting and completive reasons for jet topology while the inclined angle has a significant effect on them.Mechanisms of jet formation and evolution are different in the front and back necks.With the help of PDF distribution and correction calculation,a compromise in the competition is observed.This work lays a solid foundation for further studies of dynamics in gas-liquid flows.
基金funded by National Natural Science Foundation of China (Nos. U1632145, 81573093 and 81227902)funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) and Jiangsu Provincial Key Laboratory of Radiation Medicine and Protection, China Postdoctoral Science Foundation (No. 2016M592584)Strategic Research Grant 7004641 from City University of Hong Kong
文摘Cold atmospheric plasmas(CAPs)have attracted considerable interest in the field of plasma medicine.Generated reactive species such as hydroxyl(OH)species play an important role in applications of CAPs.Transportation of OH species towards the target and distribution of these OH species in the plasma plume play an important role in the applications of plasma medicine.In the present work,a computational model was built to simulate the transportation and distribution of OH species in CAP discharges,which was based on the level set method to dynamically track the propagation of plasma carrier gas in air.A reaction term was incorporated for the OH species.The OH species tended to diffuse around the main stream of the carrier gas,and thus covered larger radial and axial distances.A CAP discharge onto a skin layer led to the largest accumulation of OH species at the central part of the exposed area.The distribution of OH species on the skin was asymmetric,which agreed with experiments.The computational model itself and the obtained results would be useful for future development of plasma medicine.
基金supported by Corning Incorporated, the National Natural Science Foundation of China (20990224, 20976177)National Science Fund for Distinguished Young Scholars (21025627)the National Basic Research Program of China (2009CB623406)
文摘Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was first validated and improved by several test cases of a still droplet.The five distinct flow regimes of the kerosene-water system,previously identified in the experiments from Zhao et al.,were reproduced.The quantitative and qualitative agreement between the simulations and the experimental data show the effectiveness of the numerical method.The roles of the interfacial tension and contact angle on the flow patterns and shapes of droplets were discussed and highlighted according to the numerical results based on the improved two-phase LB model.This work demonstrated that the developed LBM simulator is a viable tool to study immiscible two-phase flows in microchannels,and such a tool could provide tangible guidance for the design of various microfluidic devices that involve immiscible multi-phase flows.
基金Projects(51304145,51301118,51304152)supported by the National Natural Science Foundation of ChinaProject(2013JQ7016)supported by the Natural Science Foundation of Shannxi Province,China+1 种基金Project(2013T002)supported by the Science Foundation of Taiyuan University of Technology,ChinaProject(2013JK0904)supported by Shannxi Provincial Education Department,China
文摘A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.
基金partially supported by Singapore MOE AcRF grants(R-146-000-285-114,R-146-000-327-112)NSFC(NO.11871365)supported by the National Natural Science Foundation of China(NO.12071190).
文摘We propose an efficient numerical method for the simulation of the twophase flows with moving contact lines in three dimensions.The mathematical model consists of the incompressible Navier-Stokes equations for the two immiscible fluids with the standard interface conditions,the Navier slip condition along the solid wall,and a contact angle condition(Ren et al.(2010)[28]).In the numerical method,the governing equations for the fluid dynamics are coupledwith an advection equation for a level-set function.The latter models the dynamics of the fluid interface.Following the standard practice,the interface conditions are taken into account by introducing a singular force on the interface in themomentum equation.This results in a single set of governing equations in the whole fluid domain.Similarly,the contact angle condition is imposed by introducing a singular force,which acts in the normal direction of the contact line,into theNavier slip condition.The newboundary condition,which unifies the Navier slip condition and the contact angle condition,is imposed along the solid wall.The model is solved using the finite difference method.Numerical results are presented for the spreading of a droplet on both homogeneous and inhomogeneous solid walls,as well as the dynamics of a droplet on an inclined plate under gravity.
基金supported by the National Natural Science Foundation of China(Grant Nos.51535010,51305362,11372260&11572266)the Fundamental Research Funds for the Central Universities(Grant Nos.2682014BR016&2682016CX024)the China Scholar Council
文摘The present paper aims to develop an automatical strategy for generating accurate different-scale microstructures of human tooth enamels(HTEs),and to elaborate a numerical scheme for simulating their elastic responses.At first,the strong governing formulation of these microstructures is briefly constructed,and the relevant weak formulation is then deduced based on the virtual work theorem.Afterwards,a subdividing approach,which cuts the elements intercepted by the interfaces between distinct phases automatically,is established with the aid of the level set method(LSM),and the discrete counterpart of the governing formula is obtained by combining the weak formulation derived and a discretized model.To be noted,two silent merits are found when the elaborated strategy is applied:(1) the continents constituting the microstructures of different scales can be arbitrarily-shaped and conveniently reproduced;(2) the periodic boundary condition commonly employed can be enforced on the external surfaces of representative unit cells(RUCs) with no difficulty.Besides,a boundary value problem(BVP) involving a simplified HTE nanostructure is designed,analytically solved,and hereafter applied to verify the correctness of the proposed strategy.It is observed that both the displacement and stress predictions by the computational approach are in good agreement with the relevant analytical results irrespective of the material combinations applied.Eventually,discussions are made on the influence of material organizations of both the 2D and 3D HTE microstructures at the ultrastructural and repeated rod levels,and some concluding remarks are drawn.
基金The work was supported in part by Chinese NSF Project 10431030,NSERC and MITACS(Canada).
文摘Electrowetting has been proposed as a technique for manipulating dropletssurrounded by air or oil. In this paper, we discuss the modeling and simulation of thedroplet fission process between two parallel plates inside an electrowetting on dielectric (EWOD) device. Since the gap between the plates is small, we use the two-phaseHele-Shaw flow as a model. While there are several high order methods around, suchas the immersed interface methods [1, 2], we decide to use two first-order methods forsimplicity. A ghost-fluid (GF) method is employed to solve the governing equationsand a local level set method is used to track the drop interface. For comparison purposes, the same set of two-phase Hele-Shaw equations are also solved directly usingthe immersed boundary (IB) method. Numerical results are consistent with experimental observations reported in the literature.
基金The work of Leung was supported in part by the RGC under Grant 605612。
文摘We propose a new semi-implicit level set approach to a class of curvature dependent flows.The method generalizes a recent algorithm proposed for the motion by mean curvature where the interface is updated by solving the Rudin-Osher-Fatemi(ROF)model for image regularization.Our proposal is general enough so that one can easily extend and apply the method to other curvature dependent motions.Since the derivation is based on a semi-implicit time discretization,this suggests that the numerical scheme is stable even using a time-step significantly larger than that of the corresponding explicit method.As an interesting application of the numerical approach,we propose a new variational approach for extracting limit cycles in dynamical systems.The resulting algorithm can automatically detect multiple limit cycles staying inside the initial guess with no condition imposed on the number nor the location of the limit cycles.Further,we also propose in this work an Eulerian approach based on the level set method to test if the limit cycles are stable or unstable.
基金Supported by the National Natural Science Foundation of China (No.20490206) and the Special Funds for Major State Basic Research Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
基金the National Natural Science Foundation of China (No.20490206) the Special Funds for Major State BasicResearch Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
