Speed limit measures are ubiquitous due to the complexity of the road environment,which can be supplied with the help of vehicle to everything(V2X)communication technology.Therefore,the influence of speed limit on tra...Speed limit measures are ubiquitous due to the complexity of the road environment,which can be supplied with the help of vehicle to everything(V2X)communication technology.Therefore,the influence of speed limit on traffic system will be investigated to construct a two-lane lattice model accounting for the speed limit effect during the lane change process under V2X environment.Accordingly,the stability condition and the mKdV equation are closely associated with the speed limit effect through theory analysis.Moreover,the evolution of density and hysteresis loop is simulated to demonstrate the positive role of the speed limit effect on traffic stability in the cases of strong reaction intensity and high limited speed.展开更多
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.展开更多
An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D...An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.展开更多
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.展开更多
A novel lattice hydrodynamic model is proposed by integrating the cooperative deviation of density and optimal flux under vehicle to X(V2X) environment. According to the theoretical analysis, the stability conditions ...A novel lattice hydrodynamic model is proposed by integrating the cooperative deviation of density and optimal flux under vehicle to X(V2X) environment. According to the theoretical analysis, the stability conditions and the mKdV equations affected by the cooperative deviation of traffic information are explored. And the density wave, hysteresis loop and energy consumption of the traffic flow have been investigated via numerical simulation. The results indicate that the cooperative deviation of density and optimal flux can effectively alleviate the traffic congestion. More importantly, our new consideration can reduce fuel consumption and exhaust emission under the V2X environment.展开更多
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.展开更多
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.展开更多
In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtai...In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtained by using the linear stability theory. The result shows that considering the flow difference effect leads to stabilization of the system compared with the original lattice hydrodynamic model. The jamming transitions among the freely moving phase, the coexisting phase, and the uniform congested phase are studied by nonlinear analysis. The modified KdV equation near the critical point is derived to describe the traffic jam, and kink -antikink soliton solutions related to the traffic density waves are obtained. The simulation results are consistent with the theoretical analysis for the new model.展开更多
A new lattice Boltzmann model for compressible perfect gas is proposed. The numerical example shows that it can be used to simulate shock wave and contact discontinuity. The results are comparable with those obtained ...A new lattice Boltzmann model for compressible perfect gas is proposed. The numerical example shows that it can be used to simulate shock wave and contact discontinuity. The results are comparable with those obtained by traditional methods. The ratio of specific heats gamma may be chosen according to the requirement of problems.展开更多
This paper proposes two lattice traffic models by taking into account the drivers' delay in response. The lattice versions of the hydrodynamic model are described by the differential-difference equation and differenc...This paper proposes two lattice traffic models by taking into account the drivers' delay in response. The lattice versions of the hydrodynamic model are described by the differential-difference equation and difference-difference equation, respectively. The stability conditions for the two models are obtained by using the linear stability theory. The modified KdV equation near the critical point is derived to describe the traffic jam by using the reductive perturbation method, and the kink-antikink soliton solutions related to the traffic density waves are obtained. The results show that the drivers' delay in sensing headway plays an important role in jamming transition.展开更多
We present a new multi-anticipation lattice hydrodynamic model based on the traffic anticipation effect in the real world. Applying the linear stability theory, we obtain the linear stability condition of the model. T...We present a new multi-anticipation lattice hydrodynamic model based on the traffic anticipation effect in the real world. Applying the linear stability theory, we obtain the linear stability condition of the model. Through nonlinear analysis, we derive the modified Korteweg-de Vries equation to describe the propagating behaviour of a traffic density wave near the critical point. The good agreement between the simulation results and the analytical results shows that the stability of traffic flow can be enhanced when the multi^anticipation effect is considered.展开更多
By introducing the traffic anticipation effect in the real world into the original lattice hydrodynamic model, we present a new anticipation effect lattice hydrodynamic (AELH) model, and obtain the linear stability ...By introducing the traffic anticipation effect in the real world into the original lattice hydrodynamic model, we present a new anticipation effect lattice hydrodynamic (AELH) model, and obtain the linear stability condition of the model by applying the linear stability theory. Through nonlinear analysis, we derive the Burgers equation and Korteweg-de Vries (KdV) equation, to describe the propagating behaviour of traffic density waves in the stable and the metastable regions, respectively. The good agreement between simulation results and analytical results shows that the stability of traffic flow can be enhanced when the anticipation effect is considered.展开更多
Lattice structures with excellent physical properties have attracted great research interest.In this paper,a novel volume parametric modeling method based on the skeleton model is proposed for the construction of thre...Lattice structures with excellent physical properties have attracted great research interest.In this paper,a novel volume parametric modeling method based on the skeleton model is proposed for the construction of threedimensional lattice structures.The skeleton model is divided into three types of nodes.And the corresponding algorithms are utilized to construct diverse types of volume parametric nodes.The unit-cell is assembled with distinct nodes according to the geometric features.The final lattice structure is created by the periodic arrangement of unit-cells.Several different types of volume parametric lattice structures are constructed to prove the stability and applicability of the proposed method.The quality is assessed in terms of the value of the Jacobian matrix.Moreover,the volume parametric lattice structures are tested with the isogeometric analysis to verify the feasibility of integration of modeling and simulation.展开更多
A three-dimensional(3D) lattice model for predicting the rheological behavior of asphalt mixtures was presented.In this model asphalt mixtures were described as a two-phase composite material consisting of asphalt san...A three-dimensional(3D) lattice model for predicting the rheological behavior of asphalt mixtures was presented.In this model asphalt mixtures were described as a two-phase composite material consisting of asphalt sand and coarse aggregates distributed randomly.Asphalt sand was regarded as a viscoelastic material and aggregates as an elastic material.The rheological response of asphalt mixture subjected to different constant stresses was simulated.The calibrated overall creep strain shows a good approximation to experimental results.展开更多
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 Xue's lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is o...Based on Xue's lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries (mKdV) equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.展开更多
Based on the stability analysis with no linearization and expansion, it is argued that instability in the lattice BGK model is originated from the linear relaxation hypothesis of collision in the model. The hypothesis...Based on the stability analysis with no linearization and expansion, it is argued that instability in the lattice BGK model is originated from the linear relaxation hypothesis of collision in the model. The hypothesis stands up only when the deviation from the local equilibrium is weak. In this case the computation is absolutely stable for real fluids. But for flows of high Reynolds number, this hypothesis is violated and then instability takes place physically. By performing a transformation a quantified stability criteria is put forward without those approximation. From the criteria a sufficient condition for stability can be obtained and serve as an estimation of the limited Reynolds number as high as possible.展开更多
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 single-phase lattice Boltzmann model with modified surface tension is developed in this paper to solve the problem of high-density-ratio free surface flow.The computational efficiency and accuracy are both enhanced....A single-phase lattice Boltzmann model with modified surface tension is developed in this paper to solve the problem of high-density-ratio free surface flow.The computational efficiency and accuracy are both enhanced.The restriction to the relaxation factor (which needs to be smaller than 1) is circumvented by the new surface tension algebra,due to its rational physical nature compared with the treatment of Xing,Buther and Yang in their paper (Comp.Mater.Sci.,2007,39(2):282-290).The proposed stable surface tension scheme is applied to simulate the free deformation of a square droplet with surface tension effect and the process of a droplet impinging on a liquid film.The numerical solution for free deformation of a droplet agrees well with thermodynamic principles,and also achieves high accuracy in comparison with Xing,et al.'s model.Three typical impinging modes are successfully obtained with the new scheme,and another particular mode found by Wang and Chen is also successfully simulated.The evolutions of liquid crown agree well with the power law related to time.展开更多
We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined ...We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined by generating accurate data for lattices with L=8, 10, 12, 15, 20, 25, 30, 40 and 50. The updated spin algorithm we employ has the advantages of both a Metropolis algorithm and a single-update method. Our study indicates that the transition is continuous at Тc=3.6403(2). A convincing finite-size scaling analysis of the model yields ν=0.9995(21), β/ν=0.12400(17), γ/v=1.75223(22), γ^1/ν=1.7555(22), α/ν=0.00077(420) (scaling) and α/ν=0.0010(42) (hyperscaling). The present scheme yields more accurate estimates for all the critical exponents than the Monte Carlo method, and our estimates are shown to be in excellent agreement with their predicted values.展开更多
基金Project supported by the Guangxi Natural Science Foundation,China(Grant No.2022GXNSFDA035080)the Central Government Guidance Funds for Local Scientific and Technological Development,China(Grant No.Guike ZY22096024)the National Natural Science Foundation,China(Grant No.61963008).
文摘Speed limit measures are ubiquitous due to the complexity of the road environment,which can be supplied with the help of vehicle to everything(V2X)communication technology.Therefore,the influence of speed limit on traffic system will be investigated to construct a two-lane lattice model accounting for the speed limit effect during the lane change process under V2X environment.Accordingly,the stability condition and the mKdV equation are closely associated with the speed limit effect through theory analysis.Moreover,the evolution of density and hysteresis loop is simulated to demonstrate the positive role of the speed limit effect on traffic stability in the cases of strong reaction intensity and high limited speed.
基金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.
基金the Youth Program of the National Natural Science Foundation of China(Grant Nos.11972272,12072246,and 12202331)the National Key Project,China(Grant No.GJXM92579)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JQ-028)。
文摘An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.
基金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.
基金Project supported by the National Natural Science Foundation of China (Grant No.61963008)Guangxi Natural Science Foundation (Grant No.2022GXNSFDA035080)+1 种基金Guangxi Innovation-Driven Development Special Fund Project (Grant No.GUIKEAA19254034-3)Doctor Scientific Research Startup Project Foundation of Guangxi Normal University,China (Grant No.2018BQ007)。
文摘A novel lattice hydrodynamic model is proposed by integrating the cooperative deviation of density and optimal flux under vehicle to X(V2X) environment. According to the theoretical analysis, the stability conditions and the mKdV equations affected by the cooperative deviation of traffic information are explored. And the density wave, hysteresis loop and energy consumption of the traffic flow have been investigated via numerical simulation. The results indicate that the cooperative deviation of density and optimal flux can effectively alleviate the traffic congestion. More importantly, our new consideration can reduce fuel consumption and exhaust emission under the V2X environment.
基金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(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.
基金Project supported by the National Basic Research Program of China (Grant No. G2006CB705500)the National Natural Science Foundation of China (Grant Nos. 70501004,70701004 and 70631001)Program for New Century Excellent Talents in University(Grant No. NCET-07-0057)
文摘In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtained by using the linear stability theory. The result shows that considering the flow difference effect leads to stabilization of the system compared with the original lattice hydrodynamic model. The jamming transitions among the freely moving phase, the coexisting phase, and the uniform congested phase are studied by nonlinear analysis. The modified KdV equation near the critical point is derived to describe the traffic jam, and kink -antikink soliton solutions related to the traffic density waves are obtained. The simulation results are consistent with the theoretical analysis for the new model.
基金The project supported by the National Natural Science Foundation of China
文摘A new lattice Boltzmann model for compressible perfect gas is proposed. The numerical example shows that it can be used to simulate shock wave and contact discontinuity. The results are comparable with those obtained by traditional methods. The ratio of specific heats gamma may be chosen according to the requirement of problems.
基金Project supported by the National Basic Research Program of China (Grant No 2006CB705500)the National Natural Science Foundation of China (Grant No 10532060)+1 种基金the Natural Science Foundation of Ningbo (Grant Nos 2008A610022 and 2007A610050)K. C. Wang Magna Fund in Ningbo University, China
文摘This paper proposes two lattice traffic models by taking into account the drivers' delay in response. The lattice versions of the hydrodynamic model are described by the differential-difference equation and difference-difference equation, respectively. The stability conditions for the two models are obtained by using the linear stability theory. The modified KdV equation near the critical point is derived to describe the traffic jam by using the reductive perturbation method, and the kink-antikink soliton solutions related to the traffic density waves are obtained. The results show that the drivers' delay in sensing headway plays an important role in jamming transition.
基金Project supported by the National High Technology and Development Program of China (Grant No. 511-0910-1031)Doctoral Fund of Ministry of Education of China (Grant No. 20090191110022)
文摘We present a new multi-anticipation lattice hydrodynamic model based on the traffic anticipation effect in the real world. Applying the linear stability theory, we obtain the linear stability condition of the model. Through nonlinear analysis, we derive the modified Korteweg-de Vries equation to describe the propagating behaviour of a traffic density wave near the critical point. The good agreement between the simulation results and the analytical results shows that the stability of traffic flow can be enhanced when the multi^anticipation effect is considered.
基金supported by the Fundamental Research Funds for the Central Universities (Grant No. CDJZR11170002)the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20090191110022)
文摘By introducing the traffic anticipation effect in the real world into the original lattice hydrodynamic model, we present a new anticipation effect lattice hydrodynamic (AELH) model, and obtain the linear stability condition of the model by applying the linear stability theory. Through nonlinear analysis, we derive the Burgers equation and Korteweg-de Vries (KdV) equation, to describe the propagating behaviour of traffic density waves in the stable and the metastable regions, respectively. The good agreement between simulation results and analytical results shows that the stability of traffic flow can be enhanced when the anticipation effect is considered.
基金supported by the National Nature Science Foundation of China under Grant No.52075340.
文摘Lattice structures with excellent physical properties have attracted great research interest.In this paper,a novel volume parametric modeling method based on the skeleton model is proposed for the construction of threedimensional lattice structures.The skeleton model is divided into three types of nodes.And the corresponding algorithms are utilized to construct diverse types of volume parametric nodes.The unit-cell is assembled with distinct nodes according to the geometric features.The final lattice structure is created by the periodic arrangement of unit-cells.Several different types of volume parametric lattice structures are constructed to prove the stability and applicability of the proposed method.The quality is assessed in terms of the value of the Jacobian matrix.Moreover,the volume parametric lattice structures are tested with the isogeometric analysis to verify the feasibility of integration of modeling and simulation.
基金Project(10672063) supported by the National Natural Science Foundation of China
文摘A three-dimensional(3D) lattice model for predicting the rheological behavior of asphalt mixtures was presented.In this model asphalt mixtures were described as a two-phase composite material consisting of asphalt sand and coarse aggregates distributed randomly.Asphalt sand was regarded as a viscoelastic material and aggregates as an elastic material.The rheological response of asphalt mixture subjected to different constant stresses was simulated.The calibrated overall creep strain shows a good approximation to experimental results.
基金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 High Technology Research and Development Program of China(Grant No.511-0910-1031)
文摘Based on Xue's lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries (mKdV) equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.
文摘Based on the stability analysis with no linearization and expansion, it is argued that instability in the lattice BGK model is originated from the linear relaxation hypothesis of collision in the model. The hypothesis stands up only when the deviation from the local equilibrium is weak. In this case the computation is absolutely stable for real fluids. But for flows of high Reynolds number, this hypothesis is violated and then instability takes place physically. By performing a transformation a quantified stability criteria is put forward without those approximation. From the criteria a sufficient condition for stability can be obtained and serve as an estimation of the limited Reynolds number as high as possible.
文摘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.
基金supported by the National Natural Science Foundation of China (Grant Nos.10625210,50609020 and 10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission (Grant No.J50501)
文摘A single-phase lattice Boltzmann model with modified surface tension is developed in this paper to solve the problem of high-density-ratio free surface flow.The computational efficiency and accuracy are both enhanced.The restriction to the relaxation factor (which needs to be smaller than 1) is circumvented by the new surface tension algebra,due to its rational physical nature compared with the treatment of Xing,Buther and Yang in their paper (Comp.Mater.Sci.,2007,39(2):282-290).The proposed stable surface tension scheme is applied to simulate the free deformation of a square droplet with surface tension effect and the process of a droplet impinging on a liquid film.The numerical solution for free deformation of a droplet agrees well with thermodynamic principles,and also achieves high accuracy in comparison with Xing,et al.'s model.Three typical impinging modes are successfully obtained with the new scheme,and another particular mode found by Wang and Chen is also successfully simulated.The evolutions of liquid crown agree well with the power law related to time.
基金Project supported partially by Guangdong Natural Science Foundation (GDNSF) of China (Grant No 07300793)One of authors(Loan Mushtaq) was partially supported by the Guangdong Ministry of Education,China
文摘We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined by generating accurate data for lattices with L=8, 10, 12, 15, 20, 25, 30, 40 and 50. The updated spin algorithm we employ has the advantages of both a Metropolis algorithm and a single-update method. Our study indicates that the transition is continuous at Тc=3.6403(2). A convincing finite-size scaling analysis of the model yields ν=0.9995(21), β/ν=0.12400(17), γ/v=1.75223(22), γ^1/ν=1.7555(22), α/ν=0.00077(420) (scaling) and α/ν=0.0010(42) (hyperscaling). The present scheme yields more accurate estimates for all the critical exponents than the Monte Carlo method, and our estimates are shown to be in excellent agreement with their predicted values.