Numerical simulation for MHD Williamson fluid utilizing modified Darcy’s law