A hybrid thermal lattice Boltzmann (LB) model is presented to simulate thermal multiphase flows with phase change based on an improved pseudopotential LB approach [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 87, 053301 (2013)]. The present model doe
s not suffer from the spurious term caused by the forcing-term effect, which was encountered in some previous thermal LB models for liquid-vapor phase change. Using the model, the liquid-vapor boiling process is simulated. The boiling curve together with the three boiling stages (nucleate boiling, transition boiling, and film boiling) is numerically reproduced in the LB community for the first time. The numerical results show that the basic features and the fundamental characteristics of boiling heat transfer are well captured, such as the severe fluctuation of transient heat flux in the transition boiling and the feature that the maximum heat transfer coefficient lies at a lower wall superheat than that of the maximum heat flux. Furthermore, the effects of the heating surface wettability on boiling heat transfer are investigated. It is found that an increase in contact angle promotes the onset of boiling but reduces the critical heat flux, and makes the boiling process enter into the film boiling regime at a lower wall superheat, which is consistent with the findings from experimental studies.
In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard D2Q9 lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression wo
rk are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and basic advantages of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.
This paper proposes an improved lattice Boltzmann scheme for incompressible axisymmetric flows. The scheme has the following features. First, it is still within the framework of the standard lattice Boltzmann method using the single-particle density
distribution function and consistent with the philosophy of the lattice Boltzmann method. Second, the source term of the scheme is simple and contains no velocity gradient terms. Owing to this feature, the scheme is easy to implement. In addition, the singularity problem at the axis can be appropriately handled without affecting an important advantage of the lattice Boltzmann method: the easy treatment of boundary conditions. The scheme is tested by simulating Hagen-Poiseuille flow, three-dimensional Womersley flow, Wheeler benchmark problem in crystal growth, and lid-driven rotational flow in cylindrical cavities. It is found that the numerical results agree well with the analytical solutions and/or the results reported in previous studies.
In this brief report, a thermal lattice-Boltzmann (LB) model is presented for axisymmetric thermal flows in the incompressible limit. The model is based on the double-distribution-function LB method, which has attracted much attention since its emerg
ence for its excellent numerical stability. Compared with the existing axisymmetric thermal LB models, the present model is simpler and retains the inherent features of the standard LB method. Numerical simulations are carried out for the thermally developing laminar flows in circular ducts and the natural convection in an annulus between two coaxial vertical cylinders. The Nusselt number obtained from the simulations agrees well with the analytical solutions and/or the results reported in previous studies.
