By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous me...By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore, the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.展开更多
A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear progr...A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear programming. Compared with the traditional two-phase method, it has advantages of saving the memories and reducing the computational efforts.展开更多
The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper,...The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on t...Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.展开更多
In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori...In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori discrete L2-error estimates.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm...In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm, and then use the final point in the first stage as a new initial point to turn to a projected semismooth asymptotically newton method for fast convergence.展开更多
The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase...The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.展开更多
Most numerical transient flow models that consider dynamic friction employ a finite differences approach or the method of characteristics. These models assume a single fluid (water only) with constant density and pres...Most numerical transient flow models that consider dynamic friction employ a finite differences approach or the method of characteristics. These models assume a single fluid (water only) with constant density and pressure wave velocity. But when transient flow modeling attempts to integrate the presence of air, which produces a variable density and pressure-wave velocity, the resolution scheme becomes increasingly complex. Techniques such as finite volumes are often used to improve the quality of results because of their conservative form. This paper focuses on a resolution technique for unsteady friction using the Godunov approach in a finite volume method employing single-equivalent twophase flow equations. The unsteady friction component is determined by taking into account local and convective instantaneous accelerations and the sign of both convective acceleration and velocity values. The approach is used to reproduce a set of transient flow experiments reported in the literature, and good agreement between simulated and experimental results is found.展开更多
In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial ...In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.展开更多
Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleto...Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives（normals and curvatures） and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.展开更多
We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eu...We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.展开更多
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.展开更多
In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The la...In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.展开更多
An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solutio...An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.展开更多
At the late stage of continuous casting（CC）ladle teeming,sink vortex can suck the liquid slag into tundish,and cause negative influences on the cleanliness of molten steel.To address this issue,a twophase fluid mech...At the late stage of continuous casting（CC）ladle teeming,sink vortex can suck the liquid slag into tundish,and cause negative influences on the cleanliness of molten steel.To address this issue,a twophase fluid mechanical modeling method for ladle teeming was proposed.Firstly,a dynamic model for vortex suction process was built,and the profiles of vortex flow field were acquired.Then,based on the level set method（LSM）,a two-phase 3Dinterface coupling model for slag entrapment was built.Finally,in combination with high-order essentially non-oscillatory（ENO）and total variation diminishing（TVD）methods,a LSM-based numerical solution method was proposed to obtain the 3Dcoupling evolution regularities in vortex suction process.Numerical results show that the vortex with higher kinetic energy can form an expanded sandglass-shape region with larger slag fraction and lower rotating velocity;there is a pressure oscillation phenomenon at the vortex penetration state,which is caused by the energy shock of two-phase vortex penetration coupling.展开更多
The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influence...The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.展开更多
The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） p...The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.展开更多
From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow t...From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow transient is calculated and it agrees well with theexperimental result. It is shown that the method is satisfactory for engineering demand.展开更多
文摘By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore, the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.
文摘A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear programming. Compared with the traditional two-phase method, it has advantages of saving the memories and reducing the computational efforts.
基金Sponsored by the National Natural Science Foundation of China (10672080)
文摘The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
文摘Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.
基金Supported by National Natural Science Foundation of ChinaChina Key Project for Basic Researches
文摘In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori discrete L2-error estimates.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
文摘In this paper, we proposed a spectral gradient-Newton two phase method for constrained semismooth equations. In the first stage, we use the spectral projected gradient to obtain the global convergence of the algorithm, and then use the final point in the first stage as a new initial point to turn to a projected semismooth asymptotically newton method for fast convergence.
文摘The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.
文摘Most numerical transient flow models that consider dynamic friction employ a finite differences approach or the method of characteristics. These models assume a single fluid (water only) with constant density and pressure wave velocity. But when transient flow modeling attempts to integrate the presence of air, which produces a variable density and pressure-wave velocity, the resolution scheme becomes increasingly complex. Techniques such as finite volumes are often used to improve the quality of results because of their conservative form. This paper focuses on a resolution technique for unsteady friction using the Godunov approach in a finite volume method employing single-equivalent twophase flow equations. The unsteady friction component is determined by taking into account local and convective instantaneous accelerations and the sign of both convective acceleration and velocity values. The approach is used to reproduce a set of transient flow experiments reported in the literature, and good agreement between simulated and experimental results is found.
文摘In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10972010,11028206,11371069,11372052,11402029,and 11472060)the Science and Technology Development Foundation of China Academy of Engineering Physics(CAEP),China(Grant No.2014B0201030)the Defense Industrial Technology Development Program of China(Grant No.B1520132012)
文摘Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives（normals and curvatures） and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.
基金supported by the National Natural Science Foundation of China(Grant No.11572062)the Fundamental Research Funds for the Central Universities,China(Grant No.CDJZR13248801)+2 种基金the Program for Changjiang Scholars and Innovative Research Team in University,China(Grant No.IRT13043)Key Laboratory of Functional Crystals and Laser Technology,TIPCChinese Academy of Sciences
文摘We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10872222 and 50921063)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20110191110037)the Fundamental Research Funds for the Central Universities,China(Grant Nos.CDJXS11240011 and CDJXS10241103)
文摘In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.
文摘An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.
基金supported by NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization(U1509212)National Natural Science Foundation of China(51375446)Zhejiang Provincial Natural Science Foundation for Distinguished Young Scientists(LR16E050001)
文摘At the late stage of continuous casting（CC）ladle teeming,sink vortex can suck the liquid slag into tundish,and cause negative influences on the cleanliness of molten steel.To address this issue,a twophase fluid mechanical modeling method for ladle teeming was proposed.Firstly,a dynamic model for vortex suction process was built,and the profiles of vortex flow field were acquired.Then,based on the level set method（LSM）,a two-phase 3Dinterface coupling model for slag entrapment was built.Finally,in combination with high-order essentially non-oscillatory（ENO）and total variation diminishing（TVD）methods,a LSM-based numerical solution method was proposed to obtain the 3Dcoupling evolution regularities in vortex suction process.Numerical results show that the vortex with higher kinetic energy can form an expanded sandglass-shape region with larger slag fraction and lower rotating velocity;there is a pressure oscillation phenomenon at the vortex penetration state,which is caused by the energy shock of two-phase vortex penetration coupling.
基金Financial support from NSFC(Grant No.41572303,4151001059,41101008)Key Projects in the National Science & Technology Pillar Program(2014BAL05B01)CAS "Light of West China" Program
文摘The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.
文摘The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.
文摘From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow transient is calculated and it agrees well with theexperimental result. It is shown that the method is satisfactory for engineering demand.
