Lattice Boltzmann formulation for conjugate heat transfer in heterogeneous media

Off-lattice Boltzmann simulation of conjugate heat transfer for natural convection in two-dimensional cavities

This study addresses the inadequacy of isothermal wall conditions in predicting accurate flow features and thermal effects in multicomponent systems. A finite-difference characteristic-based off-lattice Boltzmann method (OLBM) with a source term-based conjugate heat transfer (CHT) model is utilized to analyze buoyancy-driven flows in two-dimensional enclosures. The source term-based CHT model [Karani and Huber, Phys. Rev. E 91, 023304 (2015)1539-375510.1103/PhysRevE.91.023304] is extended to handle curved conjugate boundaries. The proposed CHT-OLBM solver is verified using analytical solutions and reference data from the literature. The effects of wall conduction on conjugate natural convection (CNC) problems in square and horizontal annular cavities are systematically examined with a solid wall of a nondimensional thickness of 0.2. For the square cavity problem, the governing parameters considered are 10^{5}≤Gr≤10^{9} and χ=1,5,10, while for the horizontal annulus problem, the governing parameters are taken as 10^{5}≤Ra≤10^{7} and χ=1,5,10,100, where Gr is the Grashof number, Ra is the Rayleigh number, and χ is the thermal conductivity ratio. Qualitative analysis of the simulation results using isotherms and streamlines and quantitative analysis of the local fluid-solid interface temperature, solid wall temperature distribution, Nusselt number profiles, overall Nusselt number, and effective Grashof/Rayleigh number is conducted. The overall heat transfer reduction inside the cavities due to a solid wall is quantified, and correlations are obtained to represent the overall heat transfer inside both enclosures. The findings demonstrate the significance of CHT analysis, and the CHT-OLBM solver developed in this study can successfully investigate steady, unsteady, and chaotic CNC flows.

Three-dimensional lattice Boltzmann flux solver for simulation of fluid-solid conjugate heat transfer problems with curved boundary

A three-dimensional (3D) lattice Boltzmann flux solver is presented in this work for simulation of fluid-solid conjugate heat transfer problems with a curved boundary. In this scheme, the macroscopic governing equations for mass, momentum, and energy conservation are discretized by the finite-volume method, and the numerical fluxes at the cell interface are reconstructed by the local solution of lattice Boltzmann equation. For solving the 3D fluid-solid conjugate heat transfer problems, the density distribution function (D3Q15 model) is utilized to compute the numerical fluxes of continuity and momentum equations, and the total enthalpy distribution function (D3Q7 model) is introduced to calculate the numerical flux of the energy equation. The connections between the macroscopic fluxes and the local solution of the lattice Boltzmann equation are provided by the Chapman-Enskog expansion analysis. As compared with the lattice Boltzmann method, in which the time step and grid spacing are correlated, the local solution of the lattice Boltzmann equation at each cell interface used in the present scheme is independent of each other. As a result, the drawback of the tie-up between the time step and grid spacing can be effectively removed and the developed method applies very well to nonuniform mesh and curved boundaries. To validate the performance of the developed method, the steady and unsteady natural convection in a finned 3D cavity and in a finned 3D annulus are simulated. Numerical results showed that the present scheme can effectively solve the 3D conjugate heat transfer problems with a curved boundary.

Analysis of the typical unified lattice Boltzmann models and a comprehensive multiphase model for convection-diffusion problems in multiphase systems

The present paper analyzes the typical unified lattice Boltzmann (LB) models for different convection-diffusion (CD) problems in multiphase systems. The CD problems in multiphase systems can be roughly classified into three groups: CD problems with a continuous scalar value and a continuous flux, a discontinuous scalar value and a continuous flux, a continuous scalar value and a discontinuous flux. The characteristics of the corresponding unified LB models for the three kinds of CD problems are analyzed and the equivalence between the LB models based on different perspectives or numerical schemes is revealed. Finally, a comprehensive multiphase LB model (CMLBM) capable of solving different isotropic and anisotropic CD problems in multiphase systems is proposed. Four typical CD problems in multiphase systems are calculated to validate the CMLBM; the results show that it performs well against the typical isotropic and anisotropic CD problems in multiphase systems.

Consistent second-order boundary implementations for convection-diffusion lattice Boltzmann method

In this study, an alternative second-order boundary scheme is proposed under the framework of the convection-diffusion lattice Boltzmann (LB) method for both straight and curved geometries. With the proposed scheme, boundary implementations are developed for the Dirichlet, Neumann and linear Robin conditions in a consistent way. The Chapman-Enskog analysis and the Hermite polynomial expansion technique are first applied to derive the explicit expression for the general distribution function with second-order accuracy. Then, the macroscopic variables involved in the expression for the distribution function is determined by the prescribed macroscopic constraints and the known distribution functions after streaming [see the paragraph after Eq. (29) for the discussions of the "streaming step" in LB method]. After that, the unknown distribution functions are obtained from the derived macroscopic information at the boundary nodes. For straight boundaries, boundary nodes are directly placed at the physical boundary surface, and the present scheme is applied directly. When extending the present scheme to curved geometries, a local curvilinear coordinate system and first-order Taylor expansion are introduced to relate the macroscopic variables at the boundary nodes to the physical constraints at the curved boundary surface. In essence, the unknown distribution functions at the boundary node are derived from the known distribution functions at the same node in accordance with the macroscopic boundary conditions at the surface. Therefore, the advantages of the present boundary implementations are (i) the locality, i.e., no information from neighboring fluid nodes is required; (ii) the consistency, i.e., the physical boundary constraints are directly applied when determining the macroscopic variables at the boundary nodes, thus the three kinds of conditions are realized in a consistent way. It should be noted that the present focus is on two-dimensional cases, and theoretical derivations as well as the numerical validations are performed in the framework of the two-dimensional five-velocity lattice model.

Multiple-relaxation-time lattice Boltzmann simulation for flow, mass transfer, and adsorption in porous media

In this paper, to predict the dynamics behaviors of flow and mass transfer with adsorption phenomena in porous media at the representative elementary volume (REV) scale, a multiple-relaxation-time (MRT) lattice Boltzmann (LB) model for the convection-diffusion equation is developed to solve the transfer problem with an unsteady source term in porous media. Utilizing the Chapman-Enskog analysis, the modified MRT-LB model can recover the macroscopic governing equations at the REV scale. The coupled MRT-LB model for momentum and mass transfer is validated by comparing with the finite-difference method and the analytical solution. Moreover, using the MRT-LB method coupled with the linear driving force model, the fluid transfer and adsorption behaviors of the carbon dioxide in a porous fixed bed are explored. The breakthrough curve of adsorption from MRT-LB simulation is compared with the experimental data and the finite-element solution, and the transient concentration distributions of the carbon dioxide along the porous fixed bed are elaborated upon in detail. In addition, the MRT-LB simulation results show that the appearance time of the breakthrough point in the breakthrough curve is advanced as the mass transfer resistance in the linear driving force model increases; however, the saturation point is prolonged inversely.

Role of thermal disequilibrium on natural convection in porous media: Insights from pore-scale study

The present study investigates the role of thermal nonequilibrium on natural convection in a fluid-saturated porous medium heated from below. We conduct high-resolution direct numerical simulation at the pore scale in a two-dimensional regular porous structure by means of the thermal lattice-Boltzmann method (LBM). We perform a combination of linear stability analysis of continuum-scale heat transfer models, and pore-scale and continuum-scale simulations to study the role of thermal conductivity contrasts among phases on natural convection. The comparison of pore-scale lattice-Boltzmann simulations with linear stability analysis reveals that traditional continuum-scale models fail to capture the correct onset of convection, convection mode, and heat transfer when the thermal conductivity of the solid obstacles does not match that of the fluid.

Conjugate heat transfer with the entropic lattice Boltzmann method

A conjugate heat-transfer model is presented based on the two-population entropic lattice Boltzmann method. The present approach relies on the extension of Grad's boundary conditions to the two-population model for thermal flows, as well as on the appropriate exact conjugate heat-transfer condition imposed at the fluid-solid interface. The simplicity and efficiency of the lattice Boltzmann method (LBM), and in particular of the entropic multirelaxation LBM, are retained in the present approach, thus enabling simulations of turbulent high Reynolds number flows and complex wall boundaries. The model is validated by means of two-dimensional parametric studies of various setups, including pure solid conduction, conjugate heat transfer with a backward-facing step flow, and conjugate heat transfer with the flow past a circular heated cylinder. Further validations are performed in three dimensions for the case of a turbulent flow around a heated mounted cube.

Lattice Boltzmann method for convection-diffusion equations with general interfacial conditions

In this work, we propose an interfacial scheme accompanying the lattice Boltzmann method for convection-diffusion equations with general interfacial conditions, including conjugate conditions with or without jumps in heat and mass transfer, continuity of macroscopic variables and normal fluxes in ion diffusion in porous media with different porosity, and the Kapitza resistance in heat transfer. The construction of this scheme is based on our boundary schemes [Huang and Yong, J. Comput. Phys. 300, 70 (2015)JCTPAH0021-999110.1016/j.jcp.2015.07.045] for Robin boundary conditions on straight or curved boundaries. It gives second-order accuracy for straight interfaces and first-order accuracy for curved ones. In addition, the new scheme inherits the advantage of the boundary schemes in which only the current lattice nodes are involved. Such an interfacial scheme is highly desirable for problems with complex geometries or in porous media. The interfacial scheme is numerically validated with several examples. The results show the utility of the constructed scheme and very well support our theoretical predications.

Full Eulerian lattice Boltzmann model for conjugate heat transfer

In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently. The present conjugate interface treatment is validated by several steady and unsteady numerical tests, including pure heat conduction, forced convection, and natural convection problems. Both flat and curved interfaces are also involved. The obtained results show good agreement with the analytical and/or finite volume results.

Counter-extrapolation method for conjugate interfaces in computational heat and mass transfer

In this paper a conjugate interface method is developed by performing extrapolations along the normal direction. Compared to other existing conjugate models, our method has several technical advantages, including the simple and straightforward algorithm, accurate representation of the interface geometry, applicability to any interface-lattice relative orientation, and availability of the normal gradient. The model is validated by simulating the steady and unsteady convection-diffusion system with a flat interface and the steady diffusion system with a circular interface, and good agreement is observed when comparing the lattice Boltzmann results with respective analytical solutions. A more general system with unsteady convection-diffusion process and a curved interface, i.e., the cooling process of a hot cylinder in a cold flow, is also simulated as an example to illustrate the practical usefulness of our model, and the effects of the cylinder heat capacity and thermal diffusivity on the cooling process are examined. Results show that the cylinder with a larger heat capacity can release more heat energy into the fluid and the cylinder temperature cools down slower, while the enhanced heat conduction inside the cylinder can facilitate the cooling process of the system. Although these findings appear obvious from physical principles, the confirming results demonstrates the application potential of our method in more complex systems. In addition, the basic idea and algorithm of the counter-extrapolation procedure presented here can be readily extended to other lattice Boltzmann models and even other computational technologies for heat and mass transfer systems.