“…Conservation of concentration: where div-divergent operator, V =( u , v , w ) – velocity field, ρ – density, p – pressure, μ – dynamic viscosity, β – dimensionless Casson fluid parameter, D / Dt =(∂/∂ t + u ∂/∂ x + v ∂/∂ y + w ∂/∂ z ) – substantial derivative, b =( b x , b y , b z ) – body force, T – temperature of the fluid, C p – effective specific heat, q =− k ∇ T – heat flux, C – concentration of the fluid, j =− D ∇ C –mass flux, k – thermal conductivity, D – diffusion coefficient and t – time. Since the plate is assumed to be infinitely long, all the flow variables are assumed to be function of z and t only (see Mohamad et al , 2018). Hence, the velocity vector is taken as: …”