In this paper, the effects of chemical reaction on free convective flow of electrically conducting and viscous incompressible immiscible fluids are analyzed. The coupled nonlinear equations governing the heat and mass transfer are solved analytically and numerically with appropriate boundary conditions for each fluid and the solutions have been matched at the interface. The analytical solutions are solved by using regular perturbation method valid for small values of perturbation parameter and numerically by using finite difference method. The numerical results for various values of thermal Grashof number, mass Grashof number, Hartman number, viscosity ratio, width ratio, conductivity ratio, and chemical reaction parameter have been presented graphically in the presence and in the absence of electric field load parameter. In addition, the closed form expression for volumetric flow rate, Nusselt number, species concentration, and total heat rate added to the flow is also analyzed. The solutions obtained by finite difference method and perturbation method agree very well to the order of 10 −4 for small values of perturbation parameter. C⃝ 2014 Wiley Periodicals, Inc. Heat Trans Asian Res, 44(7): 657-680, 2015; Published online in Wiley Online Library (wileyonlinelibrary.com/journal/htj).