This paper aims at showing that to prescribe a flow rate at the inlet section of a vertical channel with heated walls leads to surprising and counterintuitive physical solutions, especially when the problem is modeled as elliptical. Such an approach can give rise to the onset of recirculation cells in the entry region while the heat transfer is slightly increased under the influence of the buoyancy force. We suggest an alternative model based on more realistic boundary conditions based on a prescribed total pressure at the inlet and a fixed pressure at the outlet sections. In this case, the pressure and buoyancy forces act effectively in the same direction and, the concept of buoyancy aiding convection makes sense. The numerical results emphasize the large differences between solutions based on prescribed inlet velocity and those obtained with the present pressure-based boundary conditions.
List of symbols aThermal diffusivity (m 2 s -1 )c z Stretching parameters, Eq. 13 D Plate spacing (m) D h Hydraulic diameter, D h = 2D (m) g Gravitational acceleration (m s -2 ) Gr H Grashof number based on H, Gr H = gb 0 DTH 3 /m 0 2 h Heat transfer coefficient (W m -2 K -1 ) H Channel height (m) k Thermal conductivity (W m -1 K -1 ) L Channel length in the spanwise direction (m) _ m Mass flow rate (kg s -1 ) n x , n z Numbers of grid points in x-and z-directions p Pressure (Pa) p s Pressure at the outlet section (Pa) Pr Prandtl number, Pr = m 0 /a 0 Q Heat flux (W) Q en Enthalpy heat flux (W) Q 2w Convective heat flux along the two channel walls (W) Re Reynolds number based on D h , Re = w 0 D h /m 0 Ri Richardson number, Ri = Gr/Re 2 S c Area of the channel cross section, S c = DL (m 2 ) t Time (s) T Temperature (K) u, w Velocity components (m s -1 ) x, z Coordinates (m) Greeks b Coefficient of thermal expansion, b = 1/T 0 (K -1 ) DT Temperature difference, DT = (T h -T 0 ) (K) l Dynamic viscosity (Pa s) m Kinematic viscosity (m 2 s -1 ) q Density (kg m -3 ) h Dimensionless temperature ratio, h = (T -T 0 )/DT s Dimensionless time Subscripts a, b Analytical solutions h Hot wall H Quantity based on channel height nc Natural convection 0 Inlet section Superscripts -Averaged quantity * Dimensionless quantity