“…Among the most popular are the surface mesh method [3] that discretizes the free boundary with a finite-element mesh known as the interface segregating different phases, and the diffuse-interface method [61] that assumes a continuous variation of the phase within the interfacial regions described by a smooth field (spacedependent function) known as the phase field. In recent decades, there has been a multitude of research in free-boundary NLCs with fairly diverse modelling approaches, both in free boundary and in NLC [28,66,36,16,63,30,41,11,20,2,37]. The common practice is to design energy functionals with respect to the shape and the NLC order parameter, and find the stable configuration with energy minimization methods.…”