Coupling lattice Boltzmann model for simulation of thermal flows on standard lattices

Color-gradient-based phase-field equation for multiphase flow

In this paper, the underlying problem with the color-gradient (CG) method in handling density-contrast fluids is explored. It is shown that the CG method is not fluid invariant. Based on nondimensionalizing the CG method, a phase-field interface-capturing model is proposed which tackles the difficulty of handling density-contrast fluids. The proposed formulation is developed for incompressible, immiscible two-fluid flows without phase-change phenomena, and a solver based on the lattice Boltzmann method is proposed. Coupled with an available robust hydrodynamic solver, a binary fluid flow package that handles fluid flows with high density and viscosity contrasts is presented. The macroscopic and lattice Boltzmann equivalents of the formulation, which make the physical interpretation of it easier, are presented. In contrast to existing color-gradient models where the interface-capturing equations are coupled with the hydrodynamic ones and include the surface tension forces, the proposed formulation is in the same spirit as the other phase-field models such as the Cahn-Hilliard and the Allen-Cahn equations and is solely employed to capture the interface advected due to a flow velocity. As such, similarly to other phase-field models, a so-called mobility parameter comes into play. In contrast, the mobility is not related to the density field but a constant coefficient. This leads to a formulation that avoids individual speed of sound for the different fluids. On the lattice Boltzmann solver side, two separate distribution functions are adopted to solve the formulation, and another one is employed to solve the Navier-Stokes equations, yielding a total of three equations. Two series of numerical tests are conducted to validate the accuracy and stability of the model, where we compare simulated results with available analytical and numerical solutions, and good agreement is observed. In the first set the interfacial evolution equations are assessed, while in the second set the hydrodynamic effects are taken into account.

Generalized equilibria for color-gradient lattice Boltzmann model based on higher-order Hermite polynomials: A simplified implementation with central moments

We propose generalized equilibria of a three-dimensional color-gradient lattice Boltzmann model for two-component two-phase flows using higher-order Hermite polynomials. Although the resulting equilibrium distribution function, which includes a sixth-order term on the velocity, is computationally cumbersome, its equilibrium central moments (CMs) are velocity-independent and have a simplified form. Numerical experiments show that our approach, as in Wen et al. [Phys. Rev. E 100, 023301 (2019)2470-004510.1103/PhysRevE.100.023301] who consider terms up to third order, improves the Galilean invariance compared to that of the conventional approach. Dynamic problems can be solved with high accuracy at a density ratio of 10; however, the accuracy is still limited to a density ratio of 1000. For lower density ratios, the generalized equilibria benefit from the CM-based multiple-relaxation-time model, especially at very high Reynolds numbers, significantly improving the numerical stability.

Discrete unified gas kinetic scheme for continuum compressible flows

In this paper, a discrete unified gas kinetic scheme (DUGKS) is proposed for continuum compressible gas flows based on the total energy kinetic model [Guo et al., Phys. Rev. E 75, 036704 (2007)1539-375510.1103/PhysRevE.75.036704]. The proposed DUGKS can be viewed as a special finite-volume lattice Boltzmann method for the compressible Navier-Stokes equations in the double distribution function formulation, in which the mass and momentum transport are described by the kinetic equation for a density distribution function (g), and the energy transport is described by the other one for an energy distribution function (h). To recover the full compressible Navier-Stokes equations exactly, the corresponding equilibrium distribution functions g^{eq} and h^{eq} are expanded as Hermite polynomials up to third and second orders, respectively. The velocity spaces for the kinetic equations are discretized according to the seventh and fifth Gauss-Hermite quadratures. Consequently, the computational efficiency of the present DUGKS can be much improved in comparison with previous versions using more discrete velocities required by the ninth Gauss-Hermite quadrature.

Mesoscopic kinetic approach for studying nonequilibrium hydrodynamic and thermodynamic effects of shock wave, contact discontinuity, and rarefaction wave in the unsteady shock tube

This paper presents a detailed description of a molecular velocity distribution-based mesoscopic kinetic approach that enables a better understanding of various nonequilibrium hydrodynamic and thermodynamic effects in shock waves, contact discontinuities, and rarefaction waves. This builds on the mesoscopic kinetic approach in a previous investigation into regular reflection shocks by further addressing the mesoscopic physical meaning of kinetic moments from the view of kinetics and the implications of the magnitude and sign of nonequilibrium kinetic moments. To deepen understanding of nonequilibrium effects, this work focuses on the one-dimensional unsteady shock tube problem, which contains the typical and essential features of the discontinuous flows, and has no interference of two-dimensional flow direction. The approach uses a lattice Boltzmann method to solve the flow field, and describes nonequilibrium effects through the nonequilibrium kinetic moments of molecular velocity distribution functions. The mechanism of nonequilibrium effect in discontinuous flows is further probed. This work develops the mesoscopic kinetic approach and clarifies the mesoscopic physics of shock waves, contact discontinuities, and rarefaction waves.

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.

A Multiple-Grid Lattice Boltzmann Method for Natural Convection under Low and High Prandtl Numbers

A multi-distribution lattice Boltzmann Bhatnagar–Gross–Krook (BGK) model with a multiple-grid lattice Boltzmann (MGLB) model is proposed to efficiently simulate natural convection over a wide range of Prandtl numbers. In this method, different grid sizes and time steps for heat transfer and fluid flow equations are chosen. The model is validated against natural convection in a square cavity, since extensive benchmark solutions are available for that problem. The proposed method can resolve the computational difficulty in simulating problems with very different time scales, in particular, when using extremely low or high Prandtl numbers. The technique can also enhance computational speed and stability while keeping the simplicity of the BGK method. Compared with the conventional lattice Boltzmann method, the simulation time can be reduced up to one-tenth of the time while maintaining the accuracy in an acceptable range. The proposed model can be extended to other lattice Boltzmann collision models and three-dimensional cases, making it a great candidate for large-scale simulations.

Achieving thermodynamic consistency in a class of free-energy multiphase lattice Boltzmann models

The free-energy lattice Boltzmann (LB) model is one of the major multiphase models in the LB community. The present study is focused on a class of free-energy LB models in which the divergence of thermodynamic pressure tensor or its equivalent form expressed by the chemical potential is incorporated into the LB equation via a forcing term. Although this class of free-energy LB models may be thermodynamically consistent at the continuum level, it suffers from thermodynamic inconsistency at the discrete lattice level owing to numerical errors [Guo et al., Phys. Rev. E 83, 036707 (2010)10.1103/PhysRevE.83.036707]. The numerical error term mainly includes two parts: one comes from the discrete gradient operator and the other can be identified in a high-order Chapman-Enskog analysis. In this paper, we propose an improved scheme to eliminate the thermodynamic inconsistency of the aforementioned class of free-energy LB models. The improved scheme is constructed by modifying the equation of state of the standard LB equation, through which the discretization of ∇(ρc_{s}^{2}) is no longer involved in the force calculation and then the numerical errors can be significantly reduced. Numerical simulations are subsequently performed to validate the proposed scheme. The numerical results show that the improved scheme is capable of eliminating the thermodynamic inconsistency and can significantly reduce the spurious currents in comparison with the standard forcing-based free-energy LB model.

Lattice Boltzmann simulations of thermal flows beyond the Boussinesq and ideal-gas approximations

In this work, the recent lattice Boltzmann model with self-tuning equation of state (EOS) [R. Huang et al., J. Comput. Phys. 392, 227 (2019)]JCTPAH0021-999110.1016/j.jcp.2019.04.044 is improved in three aspects to simulate the thermal flows beyond the Boussinesq and ideal-gas approximations. First, an improved scheme is proposed to eliminate the additional cubic terms of velocity, which can significantly improve the numerical accuracy. Second, a local scheme is proposed to calculate the density gradient instead of the conventional finite-difference scheme. Third, a scaling factor is introduced into the lattice sound speed, which can be adjusted to effectively enhance numerical stability. The thermal Couette flow of a nonattracting rigid-sphere fluid, which is described by the Carnahan-Starling EOS, is first simulated, and the better performance of the present improvements on the numerical accuracy and stability is demonstrated. As a further application, the turbulent Rayleigh-Bénard convection in a supercritical fluid slightly above its critical point, which is described by the van der Waals EOS, is successfully simulated by the present lattice Boltzmann model. The piston effect of the supercritical fluid is successfully captured, which induces a fast and homogeneous increase of the temperature in the bulk region, and the time evolution from the initiation of heating to the final turbulent state is analyzed in detail and divided into five stages.

Simulation of Boiling Heat Transfer at Different Reduced Temperatures with an Improved Pseudopotential Lattice Boltzmann Method

The pseudopotential Lattice Boltzmann Method has attracted much attention in the recent years for the simulation of boiling heat transfer. Many studies have been published recently for the simulation of the bubble cycle (nucleation, growth and departure from a heated surface). This paper puts forward two-dimensional simulations of bubble nucleation, growth and departure using an improved pseudopotential Lattice Boltzmann Model from the literature at different reduced temperatures, Tr=0.76 and Tr=0.86. Two different models using the Bhatnagar–Gross–Krook (BGK) and the Multiple-Relaxation-Time (MRT) collision operators with appropriate forcing schemes are used. The results for pool boiling show that the bubbles exhibit axial symmetry during growth and departure. Numerical results of departure diameter and release period for pool boiling are compared against empirical correlations from the literature by varying the gravitational acceleration. Reasonable agreement is observed. Nucleate boiling trends with heat flux are also captured by the simulations. Numerical results of flow boiling simulations are compared by varying the Reynolds number for both reduced temperatures with the MRT model. It was found that the departure diamenter and release period decreases with the increase of the Reynolds number. These results are a direct effect of the drag force. Proper conclusions are commented at the end of the paper.

Lattice Boltzmann Modeling of Condensation Heat Transfer on Downward-Facing Surfaces with Different Wettabilities

The model of vapor condensation heat transfer on downward-facing surfaces with different wettabilities is built by a two-dimensional (2D) lattice Boltzmann method. Dynamic evolution of condensate microdroplets on different wettability surfaces is simulated and the influence on heat transfer performance is analyzed. Moreover, the mechanism of a heterogeneous wettability surface enhancing condensation heat transfer is explored by investigating the condensate behaviors in the process of condensation. The numerical results indicate that as the contact angle of the homogeneous wettability surface increases, the initial nucleation time of the condensate is prolonged, while the departure time of the condensate is reduced significantly. The temperature adjacent to the gas-liquid interface, especially in the three-phase contact line region, is much higher than elsewhere due to the release of latent heat during condensation. Coalescence and detachment behaviors of condensate droplets cause the average heat flux to fluctuate locally with time. For the hybrid wettability surface, if the proportion of hydrophobic regions is small, the condensation heat transfer performance will be deteriorated. However, increasing the hydrophobic-hydrophilic ratio has a positive effect on enhancing heat transfer. It is found that a critical hydrophobic-hydrophilic ratio exists to optimize the heat transfer performance. For the gradient wettability surface, directional migration induced by capillary force facilitates the removal of condensate droplets, thereby enhancing the condensation heat transfer. Furthermore, a larger wetting gradient benefits to further improve the heat transfer performance. The results are valuable for optimally designing the heat transfer enhancement of vapor condensation on functionalized surfaces with heterogeneous wettability.

Improved three-dimensional color-gradient lattice Boltzmann model for immiscible two-phase flows

In this paper, an improved three-dimensional color-gradient lattice Boltzmann (LB) model is proposed for simulating immiscible two-phase flows. Compared with the previous three-dimensional color-gradient LB models, which suffer from the lack of Galilean invariance and considerable numerical errors in many cases owing to the error terms in the recovered macroscopic equations, the present model eliminates the error terms and therefore improves the numerical accuracy and enhances the Galilean invariance. To validate the proposed model, numerical simulations are performed. First, the test of a moving droplet in a uniform flow field is employed to verify the Galilean invariance of the improved model. Subsequently, numerical simulations are carried out for the layered two-phase flow and three-dimensional Rayleigh-Taylor instability. It is shown that, using the improved model, the numerical accuracy can be significantly improved in comparison with the color-gradient LB model without the improvements. Finally, the capability of the improved color-gradient LB model for simulating dynamic two-phase flows at a relatively large density ratio is demonstrated via the simulation of droplet impact on a solid surface.

Lattice Boltzmann model for compressible flows on standard lattices: Variable Prandtl number and adiabatic exponent

A lattice Boltzmann model for compressible flows on standard lattices is developed and analyzed. A consistent two-population thermal lattice Boltzmann is used which allows a variable Prandtl number and a variable adiabatic exponent, and appropriate correction terms are introduced into the kinetic equations to compensate for deviations in the hydrodynamic limit. Using the concept of a shifted lattice, the model is extended to supersonic flows involving shock waves, and the shock-vortex interaction problem is simulated to show the accuracy of the proposed model. Numerical results demonstrate that the proposed model is a viable candidate for compressible flow simulations.

Improved thermal lattice Boltzmann model for simulation of liquid-vapor phase change

In this paper, an improved thermal lattice Boltzmann (LB) model is proposed for simulating liquid-vapor phase change, which is aimed at improving an existing thermal LB model for liquid-vapor phase change [S. Gong and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012)10.1016/j.ijheatmasstransfer.2012.04.037]. First, we emphasize that the replacement of ∇·(λ∇T)/∇·(λ∇T)ρc_{V}ρc_{V} with ∇·(χ∇T) is an inappropriate treatment for diffuse interface modeling of liquid-vapor phase change. Furthermore, the error terms ∂_{t_{0}}(Tv)+∇·(Tvv), which exist in the macroscopic temperature equation recovered from the previous model, are eliminated in the present model through a way that is consistent with the philosophy of the LB method. Moreover, the discrete effect of the source term is also eliminated in the present model. Numerical simulations are performed for droplet evaporation and bubble nucleation to validate the capability of the model for simulating liquid-vapor phase change. It is shown that the numerical results of the improved model agree well with those of a finite-difference scheme. Meanwhile, it is found that the replacement of ∇·(λ∇T)/∇·(λ∇T)ρc_{V}ρc_{V} with ∇·(χ∇T) leads to significant numerical errors and the error terms in the recovered macroscopic temperature equation also result in considerable errors.

Multiple-relaxation-time color-gradient lattice Boltzmann model for simulating two-phase flows with high density ratio

In this paper we propose a color-gradient lattice Boltzmann (LB) model for simulating two-phase flows with high density ratio and high Reynolds number. The model applies a multirelaxation-time (MRT) collision operator to enhance the stability of the simulation. A source term, which is derived by the Chapman-Enskog analysis, is added into the MRT LB equation so that the Navier-Stokes equations can be exactly recovered. Also, a form of the equilibrium density distribution function is used to simplify the source term. To validate the proposed model, steady flows of a static droplet and the layered channel flow are first simulated with density ratios up to 1000. Small values of spurious velocities and interfacial tension errors are found in the static droplet test, and improved profiles of velocity are obtained by the present model in simulating channel flows. Then, two cases of unsteady flows, Rayleigh-Taylor instability and droplet splashing on a thin film, are simulated. In the former case, the density ratio of 3 and Reynolds numbers of 256 and 2048 are considered. The interface shapes and spike and bubble positions are in good agreement with the results of previous studies. In the latter case, the droplet spreading radius is found to obey the power law proposed in previous studies for the density ratio of 100 and Reynolds number up to 500.

Multiblock approach for the passive scalar thermal lattice Boltzmann method

A multiblock approach for the passive scalar thermal lattice Boltzmann method (TLBM) with multiple-relaxation-time collision scheme is proposed based on the Chapman-Enskog analysis. The interaction between blocks is executed in the moment space directly and an external force term is considered. Theoretical analysis shows that all the nonequilibrium parts of the nonconserved moments should be rescaled, while the nonequilibrium parts of the conserved moments can be calculated directly. Moreover, a local scheme based on the pseudoparticles for computing heat flux is proposed with no need to calculate temperature gradient based on the finite-difference scheme. In order to validate the multiblock approach and local scheme for computing heat flux, thermal Couette flow with wall injection is simulated and good results are obtained, which show that the adoption of the multiblock approach does not deteriorate the convergence rate of TLBM and the local scheme for computing heat flux has second-order convergence rate. Further application of the present approach is the simulation of natural convection in a square cavity with the Rayleigh number up to 109.

Consistent two-population lattice Boltzmann model for thermal flows

Theory of two-population lattice Boltzmann equations for thermal flow simulations is revisited. The present approach makes use of a consistent division of the conservation laws between the two lattices, where mass and the momentum are conserved quantities on the first lattice, and the energy is conserved quantity of the second lattice. The theory of such a division is developed, and the advantage of energy conservation in the model construction is demonstrated in detail. The present fully local lattice Boltzmann theory is specified on the standard lattices for the simulation of thermal flows. Extension to the subgrid entropic lattice Boltzmann formulation is also given. The theory is validated with a set of standard two-dimensional simulations including planar Couette flow and natural convection in two dimensions.