The investigation of Darcy-Rayleigh convection with uniform heat flux Model A

The paper investigates the effect of uniform heat flux Model A on Darcy-Rayleigh convection with a free surface on the top. The assumption of Model A as thermal boundary conditions is used at the lower surface, while the upper one is only fixed to isothermal conditions. Both approaches of energy equations have been taken into consideration to account for the local thermal non-equilibrium. Linear stability analysis has been performed in the condition where the perturbation terms are evaluated in the form of plane waves. Otherwise, the eigenvalue problem is either dealt with the analytical solution by a dispersion relation or numerically by employing the Range–Kutta solver with the shooting method. According to the results, the transition of the critical values from isothermal to isoflux boundary conditions may emerge in one of the local thermal non-equilibrium assumptions. On the other hand, infinite values may appear in the critical Rayleigh number at the local thermal equilibrium regime. Overall, the thermal instability becomes extremely dominant when the configuration approximates the model of gas saturated with metal foam which is one of the cases of the local thermal non-equilibrium.