In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation （FDLBE）. Because source terms that reflect the interaction between phases can be accurately described, the new model is s...In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation （FDLBE）. Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.展开更多
The boiling heat transfer technology with cavity surfaces can provide higher heat flux under lower wall superheat,which is of great significance for the cooling of electronic chips and microelectromechanical devices.I...The boiling heat transfer technology with cavity surfaces can provide higher heat flux under lower wall superheat,which is of great significance for the cooling of electronic chips and microelectromechanical devices.In this paper,the boiling characteristics of the cavity surfaces are investigated based on the lattice Boltzmann(LB)method,focusing on the effects of cavity shapes,sizes,and heater thermal conductivity on the heat transfer performance.The results show that the triangular cavity has the best boiling performance since it has less residual vapor and higher bubble departure frequency than those of the trapezoidal and rectangular cavities.As the cavity size increases,the enhancement of heat transfer by the cavity mouth is suppressed by the heat accumulation effect at the heater bottom.The liquid rewetting process during bubble departure is the reason for the fluctuation of the space-averaged heat flux,and the heater thermal conductivity determines the fluctuation amplitude.The evaporation of liquid in the cavity with high thermal conductivity walls is more intense,resulting in shorter waiting time and higher bubble departure frequency.展开更多
The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invarianc...The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation （LBGKE） induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.展开更多
A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation be...A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation behind Navier–Stokes continuum and momentum equations and constructs an LB equation by rigorously integrating the BGK-Boltzmann equation along characteristics. It is a general theory, supporting most existing LB equations including the standard lattice BGK(LBGK) equation inherited from lattice-gas automata, whose theoretical foundation had been questioned. ACI LB theory also indicates that the characteristic parameter of an LB equation is collision number, depicting the particle-interaction intensity in the time span of the LB equation, instead of the traditionally assumed relaxation time, and the over-relaxation time problem is merely a manifestation of the temporal evolution of equilibrium distribution along characteristics under high collision number, irrelevant to particle kinetics. In ACI LB theory, the temporal evolution of equilibrium distribution along characteristics is the determinant of LB method accuracy and we numerically prove this.展开更多
This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations （NPDEs） in the form ut ＋αuux ＋βu^nuz ＋γuxx ＋δuzxx ＋ζxxxx = 0. This model...This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations （NPDEs） in the form ut ＋αuux ＋βu^nuz ＋γuxx ＋δuzxx ＋ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.展开更多
In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation i...In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation is proposed. With the Taylor and Chapman-Enskog expansion, the nonlinear Klein-Gordon equation is recovered correctly from the lattice Boltzmann equation. The method is applied on some test examples, and the numerical results have been compared with the analytical solutions or the numerical solutions reported in previous studies. The L2, L∞ and Root-Mean-Square (RMS) errors in the solutions show the efficiency of the method computationally.展开更多
We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coeffic...We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.展开更多
A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution functi...A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LB...In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LBM) with the collision and streaming process. The theoretical derivation of lattice Boltzmann model for transient neutron transport problem is proposed for the first time.The fully implicit backward difference scheme is used to ensure the numerical stability, and relaxation time and equilibrium particle distribution function are obtained. To validate the new lattice Boltzmann model, the LBM formulation is tested for a homogenous media with different sources, and both transient and steady-state LBM results get a good agreement with the benchmark solutions.展开更多
Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a...Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.展开更多
A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximat...A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximation for the collision operator, the model becomes fully explicit. In this case, any iterative technique is not needed. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretized accuracy of the LB equation. The numerical results show correct physics of the ocean circulation driven by the double-gyre wind stress with different Reynolds numbers and different spatial resolutions. An intrinsic low-frequency variability of the shallow water model is also found. The wind-driven ocean circulation exhibits subannual and interannual oscillations, which are comparable to those of models in which the conventional numerical methods are used.展开更多
基金The project supported by the National Natural Science Foundation of China(60073044)the State Key Development Programme for Basic Research of China(G1990022207).
文摘In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation （FDLBE）. Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.
基金Project supported by the National Natural Science Foundation of China(Nos.11872083,12172017,12202021)。
文摘The boiling heat transfer technology with cavity surfaces can provide higher heat flux under lower wall superheat,which is of great significance for the cooling of electronic chips and microelectromechanical devices.In this paper,the boiling characteristics of the cavity surfaces are investigated based on the lattice Boltzmann(LB)method,focusing on the effects of cavity shapes,sizes,and heater thermal conductivity on the heat transfer performance.The results show that the triangular cavity has the best boiling performance since it has less residual vapor and higher bubble departure frequency than those of the trapezoidal and rectangular cavities.As the cavity size increases,the enhancement of heat transfer by the cavity mouth is suppressed by the heat accumulation effect at the heater bottom.The liquid rewetting process during bubble departure is the reason for the fluctuation of the space-averaged heat flux,and the heater thermal conductivity determines the fluctuation amplitude.The evaporation of liquid in the cavity with high thermal conductivity walls is more intense,resulting in shorter waiting time and higher bubble departure frequency.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90816013 and 10572083)Shanghai Leading Academic Discipline Project,China (Grant No Y0103)
文摘The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation （LBGKE） induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.
基金Project supported by the National Science and Technology Major Project,China(Grant No.2017ZX06002002)
文摘A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation behind Navier–Stokes continuum and momentum equations and constructs an LB equation by rigorously integrating the BGK-Boltzmann equation along characteristics. It is a general theory, supporting most existing LB equations including the standard lattice BGK(LBGK) equation inherited from lattice-gas automata, whose theoretical foundation had been questioned. ACI LB theory also indicates that the characteristic parameter of an LB equation is collision number, depicting the particle-interaction intensity in the time span of the LB equation, instead of the traditionally assumed relaxation time, and the over-relaxation time problem is merely a manifestation of the temporal evolution of equilibrium distribution along characteristics under high collision number, irrelevant to particle kinetics. In ACI LB theory, the temporal evolution of equilibrium distribution along characteristics is the determinant of LB method accuracy and we numerically prove this.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10661005)Fujian Province Science and Technology Plan Item (Grant No. 2008F5019)
文摘This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations （NPDEs） in the form ut ＋αuux ＋βu^nuz ＋γuxx ＋δuzxx ＋ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.
文摘In this paper, in order to extend the lattice Boltzmann method to deal with more nonlinear equations, a one-dimensional (1D) lattice Boltzmann scheme with an amending function for the nonlinear Klein-Gordon equation is proposed. With the Taylor and Chapman-Enskog expansion, the nonlinear Klein-Gordon equation is recovered correctly from the lattice Boltzmann equation. The method is applied on some test examples, and the numerical results have been compared with the analytical solutions or the numerical solutions reported in previous studies. The L2, L∞ and Root-Mean-Square (RMS) errors in the solutions show the efficiency of the method computationally.
基金The project supported by the Foundation of the Laboratory for Nonlinear Mechanics of Continuous Media,Institute of Mechanics,Chinese Academy of Sciences
文摘We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.
文摘A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金supported by the Foundation of National Key Laboratory of Reactor System Design Technology(No.HT-LW-02-2014003)the State Key Program of National Natural Science of China(No.51436009)
文摘In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LBM) with the collision and streaming process. The theoretical derivation of lattice Boltzmann model for transient neutron transport problem is proposed for the first time.The fully implicit backward difference scheme is used to ensure the numerical stability, and relaxation time and equilibrium particle distribution function are obtained. To validate the new lattice Boltzmann model, the LBM formulation is tested for a homogenous media with different sources, and both transient and steady-state LBM results get a good agreement with the benchmark solutions.
基金Project supported by the National Natural Science Foundation of China(Nos.11572183 and 11272198)
文摘Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.
基金The work was supported by the One Hundred Talents Project of the Chinese Academy of Sciences(Grant No.KCL14014)the Impacts of Ocean-Land-Atmosphere Interactions over the East Asian Mon soon Region on the Climate in China(EAMOLA)(Grant No:ZKCX2-SW-210)the National Outstanding Youth Science Foundation of China(Grant No.40325016).
文摘A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximation for the collision operator, the model becomes fully explicit. In this case, any iterative technique is not needed. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretized accuracy of the LB equation. The numerical results show correct physics of the ocean circulation driven by the double-gyre wind stress with different Reynolds numbers and different spatial resolutions. An intrinsic low-frequency variability of the shallow water model is also found. The wind-driven ocean circulation exhibits subannual and interannual oscillations, which are comparable to those of models in which the conventional numerical methods are used.
