Search

Mixed convection in a inclined cavity has not been investigated by LBM in case of an imposed non zero heat flux. This type of boundary condition, representing very usual situations in physical world, is not simple to model in lattice Boltzmann schemes. In effect, the only boundary condition able to simulate an imposed temperature and an imposed heat flux at a boundary has been presented by D’Orazio et al. in previous works, where the boundary was at rest. In this work, laminar mixed convective heat transfer in two-dimensional rectangular inclined driven cavity is studied numerically by means of a double population thermal Lattice Boltzmann method. The counter-slip internal energy density boundary condition, able to simulate an imposed heat flux at the wall, is applied. Through the top moving lid the heat flux enters the cavity and it leaves the system through the bottom wall; side walls are adiabatic. Results are analyzed over a range of the Richardson numbers and tilting angles of the enclosure, encompassing the dominating forced convection, mixed convection, and dominating natural convection flow regimes. The results show that, as expected, heat transfer rate increases as increases the inclination angle, but this effect is significant for higher Richardson numbers, when buoyancy forces dominate the problem; for horizontal cavity, average Nusselt number decreases with the increase of Richardson number because of the stratified field configuration. This study shows that the counter-slip internal energy density boundary condition can be effectively used to simulate heat transfer phenomena also in case of moving walls and it makes the Lattice Boltzmann Method able to simulate a wide class of cooling process where a given thermal power must be removed.

A thermal lattice BGK model with doubled populations, together with a new boundary condition for temperature and heat flux, is proposed to simulate the two-dimensional natural convection flow in a cavity. Numerical results for the problem of insulated horizontal walls and vertical walls at different temperatures are presented; they are found to be in good agreement with those of previous works. In addition, the problem of insulated horizontal walls and vertical walls heated and cooled by means of an imposed heat flux is discussed. Numerical results are presented as compared with theoretical values and they are found to be in satisfactory agreement. Therefore, the method can capture fundamental behaviors in thermal flows of engineering interest.

Laminar forced convection heat transfer of water-Cu nanofluids in a microchannel is studied using the double population Thermal Lattice Boltzmann method (TLBM). The entering flow is at a lower temperature compared to the microchannel walls. The middle section of the microchannel is heated with a constant and uniform heat flux, simulated by means of the counter slip thermal energy boundary condition. Simulations are performed for nanoparticle volume fractions equal to 0.00%, 0.02% and 0.04% and slip coefficient equal to 0.001, 0.01 and 0.1. Reynolds number is equal to 1, 10 and 50.The model predictions are found to be in good agreement with earlier studies. Streamlines, isotherms, longitudinal variations of Nusselt number and slip velocity as well as velocity and temperature profiles for different cross sections are presented. The results indicate that LBM can be used to simulate forced convection for the nanofluid micro flows. They show that the microchannel performs better heat transfers at higher values of the Reynolds number. For all values of the Reynolds considered in this study, the average Nusselt number increases slightly as the solid volume fraction increases and the slip coefficient increases. The rate of this increase is more significant at higher values of the Reynolds number.

Thermal boundary conditions for a doubled-populations BGK model are introduced and numerically demonstrated. The unknown thermal distribution functions at the boundary are assumed to be equilibrium distribution functions, with a counter-slip internal energy density which is determined consistently with Dirichlet and/or Neumann boundary constraints. The hydrodynamic boundary conditions are adapted to situations of engineering interest, and viscous heating effects are taken in account. The method is used to simulate channel flows; numerical results and theoretical solutions are found in satisfactory agreement for both hydrodynamic and thermal fields.