“…Therefore, in most cases, analytical solutions are not possible and the equations are often solved numerically. Over the last few decades, various computational approaches have been developed to solve the set of governing PDEs including the boundary value problem for axisymmetric coordinates [32,66,73,120,121], different finite element methods [122,123,124], Monte Carlo methods [125,126,127], finite difference methods [128,129], and the phase field representation of the surface [130,131,132]. Each of these methods has its own advantages and disadvantages and, depending on the complexity of the problem, one or more of them can be implemented.…”