Heat Transfer in Porous Microchannels with Second-Order Slipping Boundary Conditions

A. A. Avramenko,Y. Y. Kovetska,I. V. Shevchuk,A. I. Tyrinov,V. I. Shevchuk
DOI: https://doi.org/10.1007/s11242-019-01300-3
IF: 3.6097
2019-05-28
Transport in Porous Media
Abstract:The paper presents results of a study of forced convection in a vertical flat and circular microchannels incorporating porous medium under slip boundary conditions of the first and second orders. The problem was solved analytically and compared with numerical simulations based on the lattice Boltzmann method. Effects of porosity and slip velocity on velocity and temperature profiles were investigated. Behavior of the normalized Nusselt number as a function of the Knudsen and Prandtl numbers, as well as the parameter <em class="EmphasisTypeItalic">M</em> characterizing porosity of the medium in the microchannel, was also elucidated. Computations indicated that a decrease in porosity (an increase in the parameter <em class="EmphasisTypeItalic">M</em>) causes a decrease in the velocity and temperature jumps on the wall, which contributes to the increase in the Nusselt number. Considering the effect of the second-order slip boundary conditions leads to the decrease in the velocity jump on the wall, when the coefficient <em class="EmphasisTypeItalic">A</em><sub>2</sub> (this coefficient considers the second-order slip boundary conditions) changes from negative to positive values. The heat transfer rate at high Prandtl numbers increases with the increasing Knudsen number, because of the improved thermal interaction of the flow with the channel wall. Given the second-order boundary conditions, the effect of parameter <em class="EmphasisTypeItalic">A</em><sub>2</sub> was not observed at small Prandtl numbers (<em class="EmphasisTypeItalic">Pr</em> ≤ 1). For <em class="EmphasisTypeItalic">A</em><sub>2</sub> &lt; 0, the second-order boundary conditions cause an increase in the normalized Nusselt number, whereas for <em class="EmphasisTypeItalic">A</em><sub>2</sub> &gt; 0 the normalized Nusselt number decreases in comparison with the case of <em class="EmphasisTypeItalic">A</em><sub>2</sub> = 0 (the first-order boundary conditions). The results by the analytical solution and the lattice Boltzmann method are in a good agreement with each other for <em class="EmphasisTypeItalic">A</em><sub>2</sub> ≥ 0, with the differences ≤ 1%. For <em class="EmphasisTypeItalic">A</em><sub>2</sub> &lt; 0, the differences between two models are more noticeable. For <em class="EmphasisTypeItalic">A</em><sub>2</sub> ≥ 0, predictions by the analytical model lie higher than the numerical results, whereas for <em class="EmphasisTypeItalic">A</em><sub>2</sub> &lt; 0, numerical simulations exceed the analytical solution.
engineering, chemical
What problem does this paper attempt to address?