AbstractIn this paper, the heat transfer properties in the three-dimensional (3D) magnetized with the Darcy-Forchheimer flow over a shrinking surface of the $$Cu + Al_{2} O_{3} /$$
C
u
+
A
l
2
O
3
/
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 $$F_{s} \le F_{sc}$$
F
s
≤
F
sc
, $$M \ge M_{C}$$
M
≥
M
C
,$$S \ge S_{C} ,$$
S
≥
S
C
,
and $$K_{C} \le K$$
K
C
≤
K
whereas no solution exists in the ranges of $$F_{s} > F_{sc}$$
F
s
>
F
sc
, $$M < M_{c}$$
M
<
M
c
, $$S < S_{c}$$
S
<
S
c
, and $$K_{C} > K$$
K
C
>
K
. Further, a reduction in the rate of heat transfer is noticed with a rise in the parameter of the copper solid volume fraction.