Investigation of boundary stresses on MHD flow in a convergent/divergent channel: An analytical and numerical study

This study investigates the boundary stresses prescribed by the traction boundary condition on an incompressible magnetohydrodynamic Jeffery-Hamel flow of viscous fluid. The governing partial differential equations are transformed to a nonlinear ordinary differential equation (ODE) by the transformation obtained from continuity equation. We have modeled traction boundary conditions to solve the problem. We have not imposed flow symmetry condition and applied prescribed flow rate condition. We have computed analytical and numerical solutions of the problem. We have developed modified Adomian decomposition method based on Daun-Rach Approach (DRA) for general nonlinear third order ODE subject to Robin’s and integral boundary conditions. The governing equation is solved using the developed scheme. The Mathematica routine, NDSolve, is used to obtain numerical solution by first approximating the constant flow rate condition by Gauss-Lobatto integration formula of four and five points as well as by trapezoidal rule. Analytical and numerical results are compared and found in good agreement. The slip and no-slip scenarios are observed both for inertial and non-inertial flows. The back flow occurs both for convergent and divergent channel geometries. Flow bifurcations and boundary layer get strengthened due to presence of boundary stresses. The symmetric and asymmetric natures of flow are also captured due to presence of boundary stresses.