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.展开更多
A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving speci...A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.展开更多
A multicomponent multiphase（MCMP） pseudopotential lattice Boltzmann（LB） model with large liquid–gas density ratios is proposed for simulating the wetting phenomena. In the proposed model, two layers of neighborin...A multicomponent multiphase（MCMP） pseudopotential lattice Boltzmann（LB） model with large liquid–gas density ratios is proposed for simulating the wetting phenomena. In the proposed model, two layers of neighboring nodes are adopted to calculate the fluid–fluid cohesion force with higher isotropy order. In addition, the different-time-step method is employed to calculate the processes of particle propagation and collision for the two fluid components with a large pseudoparticle mass contrast. It is found that the spurious current is remarkably reduced by employing the higher isotropy order calculation of the fluid–fluid cohesion force. The maximum spurious current appearing at the phase interfaces is evidently influenced by the magnitudes of fluid–fluid and fluid–solid interaction strengths, but weakly affected by the time step ratio.The density ratio analyses show that the liquid–gas density ratio is dependent on both the fluid–fluid interaction strength and the time step ratio. For the liquid–gas flow simulations without solid phase, the maximum liquid–gas density ratio achieved by the present model is higher than 1000：1. However, the obtainable maximum liquid–gas density ratio in the solid–liquid–gas system is lower. Wetting phenomena of droplets contacting smooth/rough solid surfaces and the dynamic process of liquid movement in a capillary tube are simulated to validate the proposed model in different solid–liquid–gas coexisting systems. It is shown that the simulated intrinsic contact angles of droplets on smooth surfaces are in good agreement with those predicted by the constructed LB formula that is related to Young＇s equation. The apparent contact angles of droplets on rough surfaces compare reasonably well with the predictions of Cassie＇s law. For the simulation of liquid movement in a capillary tube, the linear relation between the liquid–gas interface position and simulation time is observed, which is identical to the analytical prediction. The simulation results regarding the wetting phenomena of droplets on smooth/rough surfaces and the dynamic process of liquid movement in the capillary tube demonstrate the quantitative capability of the proposed model.展开更多
The lattice Boltzmann model with coupled chemical reaction was proposed to simulate the ion exchange process of rare earth leaching and verified by comparison with both empirical correlation of mass transfer coefficie...The lattice Boltzmann model with coupled chemical reaction was proposed to simulate the ion exchange process of rare earth leaching and verified by comparison with both empirical correlation of mass transfer coefficient and unreacted-core shrinking model. By simulation, the zonation phenomenon of leaching reagent in the leaching column was presented, and the breakthrough curve of leaching reagent was obtained. When t=50 s, there existed the saturated and exchange zones, and the leaching reagent concentration decreased gradually from 20 to 9.3 g/L. In accordance with the breakthrough curve, the breakthrough capacity of ion-type rare earth ore and the adsorbed ion concentration of leaching reagent were derived, the time of t=25 s was the breakthrough point of ammonium ion in leaching reagent and the breakthrough capacity of the rare earth ore was 125 g/L. Besides, the chemical kinetics parameters used for the solute transfer process of rare earth leaching were obtained by the simulation and then were used to determine the rate-limiting steps of rare earth leaching process.展开更多
The search for the development of a reliable mathematical model for understanding bubble dynamics behavior is an ongoing endeavor.A long list of complex phenomena underlies the physics of this problem.In the past deca...The search for the development of a reliable mathematical model for understanding bubble dynamics behavior is an ongoing endeavor.A long list of complex phenomena underlies the physics of this problem.In the past decades,the lattice Boltzmann method has emerged as a promising tool to address such complexities.In this regard,we have applied a 121-velocity multiphase lattice Boltzmann model to an asymmetric cluster of bubbles in an acoustic field.A problem as a benchmark is studied to check the consistency and applicability of the model.The problem of interest is to study the deformation and coalescence phenomena in bubble cluster dynamics,as well as the screening effect on an acoustic multibubble medium.It has been observed that the LB model is able to simulate the combination of the three aforementioned phenomena for a bubble cluster as a whole and for every individual bubble in the cluster.展开更多
A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection condit...A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.展开更多
The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received increasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time ...The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received increasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time Shan–Chen model is built to study the cavitation bubble collapse. Using this model, the cavitation phenomena induced by density perturbation are simulated to obtain the coexistence densities at certain temperature and to demonstrate the Young–Laplace equation. Then, the cavitation bubble collapse near a curved rigid wall and the consequent high-speed jet towards the wall are simulated. Moreover, the influences of initial pressure difference and bubble-wall distance on the cavitation bubble collapse are investigated.展开更多
A coupled lattice Boltzmann （LB） model with second-order accuracy is applied to the reduced-gravity, shallow water, 2.5-layer model for wind-driven double-gyre ocean circulation. By introducing the secondorder integ...A coupled lattice Boltzmann （LB） model with second-order accuracy is applied to the reduced-gravity, shallow water, 2.5-layer model for wind-driven double-gyre ocean circulation. By introducing the secondorder integral approximation for the collision operator, the model becomes fully explicit. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretization accuracy of the LB equation. The feature of the multiple equilibria solutions is found in the numerical experiments under different Reynolds numbers based on this LB scheme. With the Reynolds number increasing from 3000 to 4000, the solution of this model is destabilized from the anti-symmetric double-gyre solution to the subtropic gyre solution and then to the subpolar gyre solution. The transitions between these equilibria states are also found in some parameter ranges. The time-dependent variability of the circulation based on this LB simulation is also discussed for varying viscosity regimes. The flow of this model exhibits oscillations with different timescales varying from subannual to interannual. The corresponding statistical oscillation modes are obtained by spectral analysis. By analyzing the spatiotemporal structures of these modes, it is found that the subannual oscillation with a 9-month period originates from the barotropic Rossby basin mode, and the interarmual oscillations with periods ranging from 1.5 years to 4.6 years originate from the recirculation gyre modes, which include the barotropic and the baroclinic recirculation gyre modes.展开更多
The two-dimensional Rayleigh-Taylor Instability(RTI)under multi-mode perturbation in compressible flow is probed via the Discrete Boltzmann Modeling(DBM)with tracers.The distribution of tracers provides clear boundari...The two-dimensional Rayleigh-Taylor Instability(RTI)under multi-mode perturbation in compressible flow is probed via the Discrete Boltzmann Modeling(DBM)with tracers.The distribution of tracers provides clear boundaries between light and heavy fluids in the position space.Besides,the position-velocity phase space offers a new perspective for understanding the flow behavior of RTI with intuitive geometrical correspondence.The effects of viscosity,acceleration,compressibility,and Atwood number on the mixing of material and momentum and the mean nonequilibrium strength at the interfaces are investigated separately based on both the mixedness defined by the tracers and the non-equilibrium strength defined by the DBM.The mixedness increases with viscosity during early stage but decreases with viscosity at the later stage.Acceleration,compressibility,and Atwood number show enhancement effects on mixing based on different mechanisms.After the system relaxes from the initial state,the mean non-equilibrium strength at the interfaces presents an initially increasing and then declining trend,which is jointly determined by the interface length and the macroscopic physical quantity gradient.We conclude that the four factors investigated all significantly affect early evolution behavior of an RTI system,such as the competition between interface length and macroscopic physical quantity gradient.The results contribute to the understanding of the multi-mode RTI evolutionary mechanism and the accompanied kinetic effects.展开更多
There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity sch...There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.展开更多
In this paper,we present a Cole-Hopf transformation based lattice Boltzmann(LB)model for solving one-dimensional Burgers'equation,and compared to available LB models,the effect of nonlinear convection term can be ...In this paper,we present a Cole-Hopf transformation based lattice Boltzmann(LB)model for solving one-dimensional Burgers'equation,and compared to available LB models,the effect of nonlinear convection term can be eliminated.Through Chapman-Enskog analysis,it can be found that the converted diffusion equation based on the Cole-Hopf transformation can be recovered correctly from present LB model.Some numerical tests are also performed to validate the present LB model,and the numerical results show that,similar to previous LB models,the present model also has a second-order convergence rate in space,but it is more accurate than the previous ones.展开更多
The dynamics of the cavitation bubble collapse is a fundamental issue for the bubble collapse application and prevention. In the present work, the modified forcing scheme for the pseudopotential multi-relaxation-time ...The dynamics of the cavitation bubble collapse is a fundamental issue for the bubble collapse application and prevention. In the present work, the modified forcing scheme for the pseudopotential multi-relaxation-time lattice Boltzmann model developed by Li Q et al. [ Li Q, Luo K H and Li X J 2013 Phys. Rev. E 87 053301] is adopted to develop a cavitation bubble collapse model. In the respects of coexistence curves and Laplace law verification, the improved pseudopotential multi-relaxation-time lattice Boltzmann model is investigated. It is found that the thermodynamic consistency and surface tension are independent of kinematic viscosity. By homogeneous and heterogeneous cavitation simulation, the ability of the present model to describe the cavitation bubble development as well as the cavitation inception is verified. The bubble collapse between two parallel walls is simulated. The dynamic process of a collapsing bubble is consistent with the results from experiments and simulations by other numerical methods. It is demonstrated that the present pseudopotential multirelaxation-time lattice Boltzmann model is applicable and efficient, and the lattice Boltzmann method is an alternative tool for collapsing bubble modeling.展开更多
Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic con...Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic consistency by tuning a parameter related with the particle interaction range, is adopted to achieve desired stability and density ratio. The bubble collapse near rough solid wall was simulated by the improved MRT pseudopotential LB model. The mechanism of bubble collapse is studied by investigating the bubble profiles, pressure field and velocity field evolution. The eroding effects of collapsing bubble are analyzed in details. It is found that the process and the effect of the interaction between bubble collapse and rough solid wall are affected seriously by the geometry of solid boundary. At the same time, it demonstrates that the MRT pseudopotential LB model is a potential tool for the investigation of the interaction mechanism between the collapsing bubble and complex geometry boundary.展开更多
A relaxation-rate formula is presented for the entropic lattice Boltzmann model(ELBM) - a discrete kinetic theory for hydrodynamics. The simple formula not only guarantees the discrete time H-theorem but also gives fu...A relaxation-rate formula is presented for the entropic lattice Boltzmann model(ELBM) - a discrete kinetic theory for hydrodynamics. The simple formula not only guarantees the discrete time H-theorem but also gives full consideration to the consistency with hydrodynamics. The relaxation rate calculated with the formula effectively characterizes the drastic changes of the flow fields. By using this formula, the computational cost of the ELBM is significantly reduced and the model now can be efficiently used for a broad range of applications including high Reynolds number flows.展开更多
Lattice Boltzmann method is one of the widely used in multiphase fluid flow.However,the two main disadvantages of this method are the instability of numerical calculations due to the large density ratio of two phases ...Lattice Boltzmann method is one of the widely used in multiphase fluid flow.However,the two main disadvantages of this method are the instability of numerical calculations due to the large density ratio of two phases and impossibility of the temperature distribution to be fed back into the velocity distribution function when the temperature is simulated.Based on the combination prescribed by Inamuro,the large density ratio two-phase flow model and thermal model makes the density ratio of the model simulation to be increased to 2778:1 by optimizing the interface distribution function of two-phase which improves the accuracy of differential format.The phase transition term is added as source term into the distribution function controlling two phase order parameters to describe the temperature effect on the gas-liquid phase transition.The latent heat generated from the phase change is also added as a source term into the temperature distribution function which simulates the movement of the flow under the common coupling of density,velocity,pressure and temperature.The density and the temperature distribution of single bubble are simulated.Comparison of the simulation results with experimental results indicates a good agreement pointing out the effectiveness of the improved model.展开更多
The lattice Boltzmann method (LBM) is a numerical simplification of the Boltzmann equation of the kinetic theory of gases that describes fluid motions by tracking the evolution of the particle velocity distribution fu...The lattice Boltzmann method (LBM) is a numerical simplification of the Boltzmann equation of the kinetic theory of gases that describes fluid motions by tracking the evolution of the particle velocity distribution function based on linear streaming with nonlinear collision. To verify the reliability and accuracy of the model simulation of incompressible fluid flow, we program to simulate two-dimensional Poiseuille flow which has the analytical solution by using the three mode: D2Q9 model, He-Luo model, Guo model (D2G9). The tests show that the Guo model gives better results. So, in the article, the Guo model is used to stimulate the jet flow field, then the result of which is compared to the result from the experience. The research of this article is the first step to make use of the LBM on the aeroacoutics of jet flow and to provide a theoretical basis on jet aeroacoustics for further study.展开更多
We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-ric...We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.展开更多
Based on the lattice Boltzmann method（LBM）,an improved pseudo-potential model,combined with a method of adding force term,is used to simulate the two-phase flows caused by a liquid droplet impacting on a liquid film...Based on the lattice Boltzmann method（LBM）,an improved pseudo-potential model,combined with a method of adding force term,is used to simulate the two-phase flows caused by a liquid droplet impacting on a liquid film.In this model,the different phases are treated as one fluid,and the interfaces between the vapor and liquid phases can be obtained by density value of the fluid.This variant of the LBM allows one to obtain the densities of vapor and liquid with high accuracy.The model is validated by an example of phase separation.The early stage of the impact of droplet on liquid film is simulated,and the results are qualitatively consistent with physical phenomena.展开更多
基金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 National Natural Science Foundation of China under Grant Nos.11772064,and 11502117CAEP Foundation under Grant No.CX2019033+1 种基金the opening project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)Science Challenge Project under Grant No.JCKY2016212A501
文摘A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.
基金supported by the National Natural Science Foundation of China(Grant Nos.51371051 and 51306037)the Scientific Research Foundation of Graduate School of Southeast University,China(Grant No.YBJJ1627)
文摘A multicomponent multiphase（MCMP） pseudopotential lattice Boltzmann（LB） model with large liquid–gas density ratios is proposed for simulating the wetting phenomena. In the proposed model, two layers of neighboring nodes are adopted to calculate the fluid–fluid cohesion force with higher isotropy order. In addition, the different-time-step method is employed to calculate the processes of particle propagation and collision for the two fluid components with a large pseudoparticle mass contrast. It is found that the spurious current is remarkably reduced by employing the higher isotropy order calculation of the fluid–fluid cohesion force. The maximum spurious current appearing at the phase interfaces is evidently influenced by the magnitudes of fluid–fluid and fluid–solid interaction strengths, but weakly affected by the time step ratio.The density ratio analyses show that the liquid–gas density ratio is dependent on both the fluid–fluid interaction strength and the time step ratio. For the liquid–gas flow simulations without solid phase, the maximum liquid–gas density ratio achieved by the present model is higher than 1000：1. However, the obtainable maximum liquid–gas density ratio in the solid–liquid–gas system is lower. Wetting phenomena of droplets contacting smooth/rough solid surfaces and the dynamic process of liquid movement in a capillary tube are simulated to validate the proposed model in different solid–liquid–gas coexisting systems. It is shown that the simulated intrinsic contact angles of droplets on smooth surfaces are in good agreement with those predicted by the constructed LB formula that is related to Young＇s equation. The apparent contact angles of droplets on rough surfaces compare reasonably well with the predictions of Cassie＇s law. For the simulation of liquid movement in a capillary tube, the linear relation between the liquid–gas interface position and simulation time is observed, which is identical to the analytical prediction. The simulation results regarding the wetting phenomena of droplets on smooth/rough surfaces and the dynamic process of liquid movement in the capillary tube demonstrate the quantitative capability of the proposed model.
基金supported by the National Natural Science Foundation of China(51674125,51776212,91434113)National Key Basic Research Program of China(2015CB251402)Chinese Academy of Sciences(QYZDB-SSW-SYS029)and Outstanding Doctoral Dissertation Project Fund of JXUST(YB2016001)
文摘The lattice Boltzmann model with coupled chemical reaction was proposed to simulate the ion exchange process of rare earth leaching and verified by comparison with both empirical correlation of mass transfer coefficient and unreacted-core shrinking model. By simulation, the zonation phenomenon of leaching reagent in the leaching column was presented, and the breakthrough curve of leaching reagent was obtained. When t=50 s, there existed the saturated and exchange zones, and the leaching reagent concentration decreased gradually from 20 to 9.3 g/L. In accordance with the breakthrough curve, the breakthrough capacity of ion-type rare earth ore and the adsorbed ion concentration of leaching reagent were derived, the time of t=25 s was the breakthrough point of ammonium ion in leaching reagent and the breakthrough capacity of the rare earth ore was 125 g/L. Besides, the chemical kinetics parameters used for the solute transfer process of rare earth leaching were obtained by the simulation and then were used to determine the rate-limiting steps of rare earth leaching process.
文摘The search for the development of a reliable mathematical model for understanding bubble dynamics behavior is an ongoing endeavor.A long list of complex phenomena underlies the physics of this problem.In the past decades,the lattice Boltzmann method has emerged as a promising tool to address such complexities.In this regard,we have applied a 121-velocity multiphase lattice Boltzmann model to an asymmetric cluster of bubbles in an acoustic field.A problem as a benchmark is studied to check the consistency and applicability of the model.The problem of interest is to study the deformation and coalescence phenomena in bubble cluster dynamics,as well as the screening effect on an acoustic multibubble medium.It has been observed that the LB model is able to simulate the combination of the three aforementioned phenomena for a bubble cluster as a whole and for every individual bubble in the cluster.
基金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,China(Grant No.J201304).
文摘A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11674173,81627802,11474161,11374155 and 11474001the Qing Lan Project
文摘The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received increasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time Shan–Chen model is built to study the cavitation bubble collapse. Using this model, the cavitation phenomena induced by density perturbation are simulated to obtain the coexistence densities at certain temperature and to demonstrate the Young–Laplace equation. Then, the cavitation bubble collapse near a curved rigid wall and the consequent high-speed jet towards the wall are simulated. Moreover, the influences of initial pressure difference and bubble-wall distance on the cavitation bubble collapse are investigated.
基金The work is supported by the ＂100 Talent project＂ of Chinese Academy of Sciences （Grant No. KCL14014） the National 0utstanding Youth Science Foundation of China （Grant No. 40325016）.
文摘A coupled lattice Boltzmann （LB） model with second-order accuracy is applied to the reduced-gravity, shallow water, 2.5-layer model for wind-driven double-gyre ocean circulation. By introducing the secondorder integral approximation for the collision operator, the model becomes fully explicit. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretization accuracy of the LB equation. The feature of the multiple equilibria solutions is found in the numerical experiments under different Reynolds numbers based on this LB scheme. With the Reynolds number increasing from 3000 to 4000, the solution of this model is destabilized from the anti-symmetric double-gyre solution to the subtropic gyre solution and then to the subpolar gyre solution. The transitions between these equilibria states are also found in some parameter ranges. The time-dependent variability of the circulation based on this LB simulation is also discussed for varying viscosity regimes. The flow of this model exhibits oscillations with different timescales varying from subannual to interannual. The corresponding statistical oscillation modes are obtained by spectral analysis. By analyzing the spatiotemporal structures of these modes, it is found that the subannual oscillation with a 9-month period originates from the barotropic Rossby basin mode, and the interarmual oscillations with periods ranging from 1.5 years to 4.6 years originate from the recirculation gyre modes, which include the barotropic and the baroclinic recirculation gyre modes.
基金supported by the National Natural Science Foundation of China(under Grant No.12172061)the Opening Project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)under Grant No.KFJJ21-16 MFoundation of Laboratory of Computational Physics。
文摘The two-dimensional Rayleigh-Taylor Instability(RTI)under multi-mode perturbation in compressible flow is probed via the Discrete Boltzmann Modeling(DBM)with tracers.The distribution of tracers provides clear boundaries between light and heavy fluids in the position space.Besides,the position-velocity phase space offers a new perspective for understanding the flow behavior of RTI with intuitive geometrical correspondence.The effects of viscosity,acceleration,compressibility,and Atwood number on the mixing of material and momentum and the mean nonequilibrium strength at the interfaces are investigated separately based on both the mixedness defined by the tracers and the non-equilibrium strength defined by the DBM.The mixedness increases with viscosity during early stage but decreases with viscosity at the later stage.Acceleration,compressibility,and Atwood number show enhancement effects on mixing based on different mechanisms.After the system relaxes from the initial state,the mean non-equilibrium strength at the interfaces presents an initially increasing and then declining trend,which is jointly determined by the interface length and the macroscopic physical quantity gradient.We conclude that the four factors investigated all significantly affect early evolution behavior of an RTI system,such as the competition between interface length and macroscopic physical quantity gradient.The results contribute to the understanding of the multi-mode RTI evolutionary mechanism and the accompanied kinetic effects.
基金sponsored by the National Natural Science Foundation of China under Grant No.52206101Shanghai Sailing Program under Grant No.20YF1431200the Experiments for Space Exploration Program and the Qian Xuesen Laboratory,China Academy of Space Technology under Grant No.TKTSPY-2020-01-01.
文摘There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.
基金Supported by the National Natural Science Foundation of China under Grant No.51576079
文摘In this paper,we present a Cole-Hopf transformation based lattice Boltzmann(LB)model for solving one-dimensional Burgers'equation,and compared to available LB models,the effect of nonlinear convection term can be eliminated.Through Chapman-Enskog analysis,it can be found that the converted diffusion equation based on the Cole-Hopf transformation can be recovered correctly from present LB model.Some numerical tests are also performed to validate the present LB model,and the numerical results show that,similar to previous LB models,the present model also has a second-order convergence rate in space,but it is more accurate than the previous ones.
基金supported by the National Natural Science Foundation of China(Grant Nos.11274092 and 1140040119)the Natural Science Foundation of Jiangsu Province,China(Grant No.SBK2014043338)
文摘The dynamics of the cavitation bubble collapse is a fundamental issue for the bubble collapse application and prevention. In the present work, the modified forcing scheme for the pseudopotential multi-relaxation-time lattice Boltzmann model developed by Li Q et al. [ Li Q, Luo K H and Li X J 2013 Phys. Rev. E 87 053301] is adopted to develop a cavitation bubble collapse model. In the respects of coexistence curves and Laplace law verification, the improved pseudopotential multi-relaxation-time lattice Boltzmann model is investigated. It is found that the thermodynamic consistency and surface tension are independent of kinematic viscosity. By homogeneous and heterogeneous cavitation simulation, the ability of the present model to describe the cavitation bubble development as well as the cavitation inception is verified. The bubble collapse between two parallel walls is simulated. The dynamic process of a collapsing bubble is consistent with the results from experiments and simulations by other numerical methods. It is demonstrated that the present pseudopotential multirelaxation-time lattice Boltzmann model is applicable and efficient, and the lattice Boltzmann method is an alternative tool for collapsing bubble modeling.
文摘Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic consistency by tuning a parameter related with the particle interaction range, is adopted to achieve desired stability and density ratio. The bubble collapse near rough solid wall was simulated by the improved MRT pseudopotential LB model. The mechanism of bubble collapse is studied by investigating the bubble profiles, pressure field and velocity field evolution. The eroding effects of collapsing bubble are analyzed in details. It is found that the process and the effect of the interaction between bubble collapse and rough solid wall are affected seriously by the geometry of solid boundary. At the same time, it demonstrates that the MRT pseudopotential LB model is a potential tool for the investigation of the interaction mechanism between the collapsing bubble and complex geometry boundary.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11471185,11801030,and 11861131004)
文摘A relaxation-rate formula is presented for the entropic lattice Boltzmann model(ELBM) - a discrete kinetic theory for hydrodynamics. The simple formula not only guarantees the discrete time H-theorem but also gives full consideration to the consistency with hydrodynamics. The relaxation rate calculated with the formula effectively characterizes the drastic changes of the flow fields. By using this formula, the computational cost of the ELBM is significantly reduced and the model now can be efficiently used for a broad range of applications including high Reynolds number flows.
基金supported by the National Natural Science Foundation of China (51609131)Shandong Provincial Natural Science Foundation of China (ZR2017MEE031)+1 种基金 Weihai Science and Technology Development Plan (2017GNS18)Shandong Provincial Higher Educational Science and Technology Foundation of China (J16LA61)
文摘Lattice Boltzmann method is one of the widely used in multiphase fluid flow.However,the two main disadvantages of this method are the instability of numerical calculations due to the large density ratio of two phases and impossibility of the temperature distribution to be fed back into the velocity distribution function when the temperature is simulated.Based on the combination prescribed by Inamuro,the large density ratio two-phase flow model and thermal model makes the density ratio of the model simulation to be increased to 2778:1 by optimizing the interface distribution function of two-phase which improves the accuracy of differential format.The phase transition term is added as source term into the distribution function controlling two phase order parameters to describe the temperature effect on the gas-liquid phase transition.The latent heat generated from the phase change is also added as a source term into the temperature distribution function which simulates the movement of the flow under the common coupling of density,velocity,pressure and temperature.The density and the temperature distribution of single bubble are simulated.Comparison of the simulation results with experimental results indicates a good agreement pointing out the effectiveness of the improved model.
文摘The lattice Boltzmann method (LBM) is a numerical simplification of the Boltzmann equation of the kinetic theory of gases that describes fluid motions by tracking the evolution of the particle velocity distribution function based on linear streaming with nonlinear collision. To verify the reliability and accuracy of the model simulation of incompressible fluid flow, we program to simulate two-dimensional Poiseuille flow which has the analytical solution by using the three mode: D2Q9 model, He-Luo model, Guo model (D2G9). The tests show that the Guo model gives better results. So, in the article, the Guo model is used to stimulate the jet flow field, then the result of which is compared to the result from the experience. The research of this article is the first step to make use of the LBM on the aeroacoutics of jet flow and to provide a theoretical basis on jet aeroacoustics for further study.
基金supported by the Strategic Program of Chinese Academy of Sciences(Grant No.XDB10030400)the Hundred Talent Program of Chinese Academy of Sciences(Grant No.Y323081C01)the National Natural Science Foundation of China(Grant No.51439008)
文摘We propose a boundary scheme for addressing multi-mechanism flow in a porous medium in slip and early transition flow regimes, which is frequently encountered in shale gas reservoirs. Micro-gaseous flow in organic-rich shale involves a complex flow mechanism. A self-developed boundary scheme that combines the non-equilibrium extrapolation scheme and the combined diffusive reflection and bounce-back scheme(half-way DBB) to embed the Langmuir slip boundary into the single-relaxation-time lattice Boltzmann method(SRT-LBM) enables us to describe this process, namely, the coupling effect of micro-gaseous flow and surface diffusion in organic-rich nanoscale pores. The present LBM model comes with the careful consideration of the local Knudsen number, local pressure gradient, viscosity correction model, and regularization procedure to account for the rarefied gas flows in irregular pores. Its validity and accuracy are verified by several benchmarking cases, and the calculated results by this boundary scheme accord well with our analytical solutions.This boundary scheme shows a higher accuracy than the existing studies. Additionally, a subiteration strategy is presented to tackle the coupled micro-gaseous flow and surface diffusion, which necessitates the iteration process matching of these two mechanisms. The multi-mechanism flow in the self-developed irregular pores is also numerically investigated and analyzed over a wide range of parameters. The results indicate that the present model can effectively capture the coupling effect of micro-gaseous flow and surface diffusion in a tree-like porous medium.
基金Project supported by the National Natrual Science Foundation of China (Grant No.10872123)
文摘Based on the lattice Boltzmann method（LBM）,an improved pseudo-potential model,combined with a method of adding force term,is used to simulate the two-phase flows caused by a liquid droplet impacting on a liquid film.In this model,the different phases are treated as one fluid,and the interfaces between the vapor and liquid phases can be obtained by density value of the fluid.This variant of the LBM allows one to obtain the densities of vapor and liquid with high accuracy.The model is validated by an example of phase separation.The early stage of the impact of droplet on liquid film is simulated,and the results are qualitatively consistent with physical phenomena.
