Heat transfer in the Jeffery-Hamel flow of a yield-stress fluid

Heat transfer in the Jeffery-Hamel flow of a yield-stress fluid that obeys the Bingham-Papanastasiou model is investigated in the present study. Among various rheological models proposed to formulate the flow of yield stress fluids, the Bingham model is usually used to address the effect of yield-stress. A fluid that obeys the Bingham model does not show any shear-thinning or shear-thickening behaviors; this makes the interpretation of the results much easier. Addressing the unyielded regions (regions at which the yield-stress material behaves like a solid) in a yield-stress fluid flow is not a trivial task, and it could be impossible, especially in complex geometries. Unfortunately, the Bingham model itself cannot be used to formulate complex fluid flow problems due to the mathematical difficulties (i.e., addressing the unyielded regions); hence it is usually regularized to avoid these mathematical challenges. One of these regularized models is the Bingham-Papanastasiou model. Using an exponential term enables the Bingham-Papanastasiou model to address the unyielded regions. However, this exponential term results in highly non-linear equations, which are difficult to be handled mathematically. In the present study, the effect of the yield-stress (Bingham number) on the hydrodynamics and heat transfer in the Jeffery-Hamel flow is investigated by using the Bingham-Papanastasiou model. The momentum transfer in the Jeffery-Hamel flow of a Bingham-Papanastasiou fluid is governed by a highly non-linear ordinary differential equation (ODE) which is a boundary value problem. Assuming that the fluid's properties are independent of the temperature and pressure, this highly non-linear ODE is solved numerically. The heat transfer in the flow, with a constant wall temperature boundary condition, is investigated, assuming the flow to be an equilibrium viscous dissipation flow (the fluid's mixed mean temperature variation is negligible along the channel). The resultant linear energy equation is also solved numerically to obtain the temperature distribution across the channel. The results obtained in this work show that the yield-stress increases the wall shear stress, Nusselt number, and Stanton number. In fact, the effect of the yield-stress on the Nusselt number is comparable to the effects of Prandtl and Eckert numbers. It is also found that the yield-stress can prevent the flow separation in a divergent channel. Considering the advantages of using the yield-stress fluids as the working fluids in the Jeffery-Hamel flow, they can be considered as an alternative to nanoparticles which are usually used to enhance the fluid's thermal conductivity and also the Nusselt number.