The present work presents a numerical study of the evaporation of a sessile liquid droplet deposited on a substrate and subjected to different heating configurations. The physical formulation accounts for evaporation, the Marangoni effect, conductive transfer in the support, radiative heating, and diffusion–convection in the droplet itself. The moving interface is solved using the Arbitrary Lagrangian–Eulerian (ALE) method. Simulations were performed using COMSOL Multiphysics. Different configurations were performed to investigate the effect of the heating conditions on the shape and intensity of the Marangoni circulations. A droplet can be heated by the substrate (different natures and thicknesses were tested) and/or by a heat flux supplied at the top of the droplet. The results show that the Marangoni flow can be controlled by the heating configuration. An upward Marangoni flow was obtained for a heated substrate and a downward Marangoni flow for a flux imposed at the top of the droplet. Using both heat sources generated two vortices with an upward flow from the bottom and a downward flow from the top. The position of the stagnation zone depended on the respective intensities of the heating fluxes. Controlling the circulation in the droplet might have interesting applications, such as the control of the deposition of microparticles in suspension in the liquid, the deposition of the solved constituent, and the enhancement of the evaporation rate.