Hydrodynamic and thermosolutal analysis of MHD ternary hybrid nanofluids in BFS configuration: high-order FEM for drag reduction and thermal enhancement

Purpose This study aims to numerically investigate the combined effects of thermosolutal convection, magnetohydrodynamics (MHD) and radiative heat transfer in a backward-facing step (BFS) channel filled with a ternary hybrid ferrofluid suspension (Cu–Fe 3 O 4 –CoFe 2 O 4 /water) modeled as a Casson fluid. The primary objective is to analyze how key parameters, such as the Reynolds number ( Re ), Hartmann number ( Ha ), Lewis number ( Le ) and obstacle positioning, influence hydrodynamic forces (drag and lift coefficients), heat and mass transfer and flow stability. The study aims to provide actionable insights for optimizing thermal management systems, enhancing microfluidic device performance and advancing biomedical applications involving hybrid nanofluids and non-Newtonian fluids. Design/methodology/approach The governing equations for mass, momentum, energy and solute transport are solved using a high-order finite element method (FEM), with nonlinearities addressed via Newton’s method. Time integration is carried out using a nonstationary scheme based on the backward differentiation formula (BDF). The model accounts for magnetohydrodynamic (MHD) effects, thermal radiation and the rheological behavior of Casson fluid. The numerical implementation is validated against experimental data and benchmark solutions prior to performing the simulations. Findings Key results show that the ternary hybrid nanofluid enhances heat transfer, with a 1.03% increase in the Nusselt number, while the Casson fluid reduces drag and stabilizes flow reattachment. Increasing Re enlarges recirculation zones but decreases the drag coefficient ( C D ) by 95%. In contrast, higher Ha increases C D by 92% due to Lorentz forces. Obstacle positioning significantly alters hydrodynamic forces, with minimal C D at y 0 = 0.7 H and maximum shear-induced drag at y 0 = 1.3 H . The lift coefficient ( C L ) transitions nonmonotonically with x 0 , and magnetic fields redistribute pressure, amplifying C L . Originality/value This