We examine graphically the behavior of the dual solution of an unsteady incompressible viscous flow with heat transfer due to a moving porous surface. The mathematical model is developed with the implementation of conservation laws which are then tackled numerically, and the simulations are performed through the Keller–Box technique. Application of the Keller–Box technique allows the non-linear differential equations to transform to a scheme of linear differential equations which are easier to illuminate. Outcomes of velocity and thermal behavior are shown graphically for certain values of the dimensionless parameters. It is obvious from the results that temperature of fluid decreases with the increasing values of dimensionless Prandtl number Pr, whereas the velocity enhances for blowing situations.