AbstractIn this paper, the unsteady flow of two immiscible fluids with heat transfer is studied numerically with a modified cubic Bâspine Differential Quadrature Method. Generalized Couette flow of two immiscible dusty (fluidâparticle suspension) and pure (Newtonian) fluids are considered through rigid horizontal channels for three separate scenarios: first for nonporous plates with heat transfer, second for porous plates with uniform suction and injection and heat transfer, and third for nonporous plates with interface evolution. The stable liquidâliquid interface is considered for the two immiscible fluids in the first two cases. In the third case, it is assumed that the interface travels from one position to another and may undergo serious deformation; hence the single momentum equation based on the volume of fluid method is combined with the continuum surface approach model, and an interface tracking is proposed. The flow cases are considered to be subjected to three different pressure gradients, of relevance to energy systemsânamely, applied constant, decaying, and periodic pressure gradients. For each case, the coupled partial differential equations are formulated and solved numerically using MCBâDQM to compute the fluids velocities, fluid temperatures, interface evolution. The effects of emerging thermoâfluid parameters, that is, Eckert (dissipation), Reynolds, Prandtl, and Froude numbers, particle concentration parameter, volume fraction parameter, pressure gradient, time, and the ratio of viscosities, densities, thermal conductivities, and specific heats on velocity and temperature characteristics are illustrated through graphs.