  In this work, we propose and compare three numerical methods to handle the
one-phase Hele-Shaw problem with surface tension in dimension two by using
three variational approaches in the spirit of the seminal works \cite{Otto,
Gia_Otto}.
