This study aims to perform a numerical simulation of the boundary flow with
the characteristic Sakiadis flow of the magnetohydrodynamic (MHD) Jeffrey
fluid under the Cattaneo-Christov heat flux model over the horizontal plate.
The similarity transformation for the local similarity solution was used to
reduce the set of governing equations to non-linear ordinary differential
equations. The equations were solved by using ?dsolve? command with the
numeric option for the boundary value problem in Maple. Simulations have
been carried out for different values of the relaxation to retardation
times, the Deborah number, the magnetic field parameter, the heat flux
relaxation time, the Prandtl number, and the Schmidt parameter. A
comparative study of the numerical results from the previously published
paper with the present result for the dimensionless velocity gradient over
the horizontal plate shows excellent agreement. It has been found that the
growth of the Deborah number leads to the dimensionless velocity gradient
enhancement, while the increment of the relaxation to retardation times
parameter and the magnetic field parameter indicates the opposite trend. The
heat transfer rate noticeably decreased with an increment in the Prandtl
number and thermal relaxation time at the fluid regime. Also, fluid
concentration decreases with larger values of the Schmidt parameter.