Numerical investigation of viscoelastic boundary layer in forced convection flow on surface under prescribed heat flux

The paper intends to analyze numerically the effect of forced convection with viscoelastic boundary on porous media subjected to constant heat flux on the surface. The modified Navier–Stokes equations in nondimensional forms are formulated and modeled with Darcy–Forchheimer–Brinkman to solve the equations using bvp4c with MATLAB package. The numerical results are found for local Nusselt number, shear stress, temperature distribution, and finally the velocity of the boundary layer over a horizontal plate. The analysis includes the influence of several parameters, such as Forchheimer, porous media, Nusselt number, and viscoelastic on distributions of velocity and temperature. Also, the effects comprise the friction and heat transfer coefficients. The numerical results under constant heat flux indicate an increase in shear stress and boundary layer velocity as the Darcy parameter increases and at the same time, the values of temperature and Nusselt number are decreased. A significant increase was observed for the values of Nusselt number, temperature, and shear stress under large values of the Forchheimer parameter. The numerical results show a good agreement with available previous works.