“…The method of manufactured solutions is adopted [20,21,24,25]; hence, we consider the model equations with the functions f w (y, t), f ϕ (y, t) and f θ (y, t) added to the right hand side of (1), ( 2) and (3), respectively, where Further, we impose the initial and boundary conditions w(y, 0) = cos(πy), ϕ(y, 0) = e −y 2 , θ(y, 0) = 1 (14) and…”