Computation of Sakiadis flow of an Eyring-Powell rheological fluid from a moving porous surface with a non-Fourier heat flux model

This article examines theoretically and numerically the effect of non-Fourier heat flux on non-Newtonian (Eyring-Powell) Sakiadis convective flow from a moving permeable surface accompanied by a parallel free-stream velocity, as a simulation of polymeric coating processes. The Cattaneo-Christov model is deployed which features thermal relaxation effects as these are important in thermal polymer processing. The physical flow problem is modelled in a Cartesian coordinate system and the governing conservation differential equations and associated boundary conditions are rendered dimensionless by applying suitable transformations. Liquid velocity and thermal distributions are computed considering numerical procedure namely, a shooting method in conjunction with the 5th order Runge-Kutta algorithm (R-K5) executed in a symbolic software. Validation with the three-stage Lobatto IIIA algorithm in MATLAB is included. The impact of key parameters on streamline distributions is also computed. Velocity is increased with increment in Eyring-Powell first parameter for the Sakiadis case whereas it is reduced with Eyring-Powell second parameter for the case where sheet and liquid are inspiring in the similar direction. The special case of Blasius flow is also examined (stationary sheet). For higher injection, there is a solid dampening in the boundary-layer flow for both Sakiadis and Blasius scenarios. With increment in thermal relaxation parameter and Eyring-Powell first and second parameters, temperatures are strongly reduced, and thermal boundary-layer thickness is suppressed. Higher injection at the wall also depletes temperatures. The Cattaneo-Christov heat flux model predicts lower temperature and thermal boundary layer thickness due to the thermal relaxation effect than the classical Fourier model. Flow patterns are displayed through 2D and 3D streamline contour plots and non-Newtonian characteristics are also found to exhibit strong modifications in these plots. The computations are useful in industrial high temperature coating flow designs using non-Newtonian materials on a moving substrate