“…A second-order central differencing scheme was used for flow parameters, whereas the turbulent properties were discretized using the van Leer scheme (Van Leer [ 55 ]). The numerical treatment and discretization techniques utilized in this study have already been extensively benchmarked over a wide range of thermal-hydraulic and renewable energy configurations, see Afgan et al [ 56 ], Filippone and Afgan [ 57 ], Guleren et al [ 58 ], Han et al [ 59 ], Wu et al [ 60 ], Wu et al [ 61 ], Wu et al [ 62 ], Kahil et al [ 63 ], Benhamadouche et al [ 64 ], Nguyen et al [ 65 ], Ejeh et al [ 66 ], Ejeh et al [ 67 ], Revell et al [ 68 ], Ahmed et al [ 69 ], Ahmed et al [ 70 ], McNaughton et al [ 71 ] and Abed et al [ 72 ]. The boundary conditions applied in the current study are summarized below:…”