Numerical simulation of Buoyancy-Marangoni convection in two superposed immiscible liquid layers with a free surface