A numerical study is presented of the free surface deformation and Marangoni convection in immiscible droplets positioned by an electrostatic field and heated by laser beams under microgravity. The boundary element and the weighted residuals methods are applied to iteratively solve for the electric field distribution and for the unknown free surface shapes, while the Galerkin finite element method for the thermal and fluid flow field in both the transient and steady states. Results show that the inner interface demarking the two immiscible fluids in an electrically conducting droplet maintains its sphericity in microgravity. The free surface of the droplet, however, deforms into an oval shape in an electric field, owing to the pulling action of the normal component of the Maxwell stress. The thermal and fluid flow distributions are rather complex in an immiscible droplet, with conduction being the main mechanism for the thermal transport. The non-uniform temperature along the free surface induces the flow in the outer layer, whereas the competition between the interfacial surface tension gradient and the inertia force in the outer layer is responsible for the flows in the inner core and near the immiscible interface. As the droplet cools into an undercooled state, surface radiation causes a reversal of the surface temperature gradients along the free surface, which in turn reverses the surface tension driven flow in the outer layer. The flow near the interfacial region, on the other hand, is driven by a complimentary mechanism between the interfacial and the inertia forces during the time when the thermal gradient on the free surface has been reversed while that on the interface has not yet. After the completion of the interfacial thermal gradient reversal, however, the interfacial flows are largely driven by the inertia forces of the outer layer fluid.