Abstract
The Darcy-Forchheimer flow of a Williamson fluid over a Riga plate was analyzed in this paper. Energy and mass equations are modeled with Cattaneo-Christov theory and double stratifications. The governing PDE models are altered into ODE models. These models are numerically solved by MATLAB bvp4c and analytically solved by the homotopy analysis method. The impact of governing flow parameters on fluid velocity, fluid temperature, fluid concentration, skin-friction coefficient, local Nusselt number, and local Sherwood number is scrutinized via graphs and tables. We acknowledged that the speed of the fluid becomes diminishes for more presence of porosity parameter. Also, we noted that the thermal and solutal boundary layer thicknesses are waning due to their corresponding stratification parameters. In addition, the maximum decreasing percentage of skin friction is obtained when the suction/injection parameter varies from 0.0 to 0.4 for Williamson and viscous fluids. The maximum increasing percentage of local Nusselt number occurs when the suction/injection parameter varies from 0.4 to 0.8 for Williamson and viscous fluids.