Influence of a Darcy-Forchheimer porous medium on the flow of a radiative magnetized rotating hybrid nanofluid over a shrinking surface

Sumera Dero; Hisamuddin Shaikh; Ghulam Hyder Talpur; Ilyas Khan; Sayer O Alharbim; Mulugeta Andualem

doi:10.1038/s41598-021-03470-x

Influence of a Darcy-Forchheimer porous medium on the flow of a radiative magnetized rotating hybrid nanofluid over a shrinking surface

Sci Rep. 2021 Dec 20;11(1):24257. doi: 10.1038/s41598-021-03470-x.

Authors

Sumera Dero¹, Hisamuddin Shaikh², Ghulam Hyder Talpur³, Ilyas Khan⁴, Sayer O Alharbim⁵, Mulugeta Andualem⁶

Affiliations

¹ Institute of Mathematics and Computer Science, University of Sindh, Jamshoro, 76080, Pakistan.
² Department of Mathematics, Shah Abdul Latif University Khairpur Mirs, Sindh, Pakistan.
³ Department of Statistics, University of Sindh, Jamshoro, Pakistan.
⁴ Department of Mathematics, College of Science Al-Zulfi, Majmah University, Al-Majmah, 11952, Saudi Arabia. i.said@mu.edu.sa.
⁵ Department of Mathematics, College of Science Al-Zulfi, Majmah University, Al-Majmah, 11952, Saudi Arabia. so.alharbi@mu.edu.sa.
⁶ Department of Mathematics, Bonga University, Bonga, Ethiopia. mulugetaandualem4@gmail.com.

Abstract

In this paper, the heat transfer properties in the three-dimensional (3D) magnetized with the Darcy-Forchheimer flow over a shrinking surface of the [Formula: see text] water hybrid nanofluid with radiation effect were studied. Valid linear similarity variables convert the partial differential equations (PDEs) into the ordinary differential equations (ODEs). With the help of the shootlib function in the Maple software, the generalized model in the form of ODEs is numerically solved by the shooting method. Shooting method can produce non-unique solutions when correct initial assumptions are suggested. The findings are found to have two solutions, thereby contributing to the introduction of a stability analysis that validates the attainability of first solution. Stability analysis is performed by employing if bvp4c method in MATLAB software. The results show limitless values of dual solutions at many calculated parameters allowing the turning points and essential values to not exist. Results reveal that the presence of dual solutions relies on the values of the porosity, coefficient of inertia, magnetic, and suction parameters for the specific values of the other applied parameters. Moreover, it has been noted that dual solutions exist in the ranges of [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text] whereas no solution exists in the ranges of [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]. Further, a reduction in the rate of heat transfer is noticed with a rise in the parameter of the copper solid volume fraction.