The MHD oscillatory flow of two immiscible, viscous liquids in a porous channel with heat transfer is the subject of this investigation. The two liquid layers with different viscosities flow in both regions. The analytical expressions for velocity and temperature distribution have been derived by solving the governing flow equations using the regular perturbation method. The effects of various parameters on the velocity, temperature, and Nusselt number have been shown graphically, and numerical values of skin friction and flow rate are presented in tabular form and discussed. According to our analysis, the mass flux reduces as the magnetic field strength rises. While the temperature of the liquid enhances with an increase in the Eckert number and the Prandtl number, the temperature distribution rises with a decrease in the thermal conductivity ratio. To validate the results, the analytical solutions are compared with the fourth-order numerical Runge–Kutta method coupled with the shooting approach, and the results are found to be in excellent agreement.