This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.

Search for Publication

Year(s) from:  to 
Keywords (separated by spaces):

Fully implicit Finite Element method for the modeling of free surface flows with surface tension effect

Aymen Laadhari, Gabor Székely
International Journal for Numerical Methods in Engineering
February 2017, in press


A new approach based on the use of the Newton and level set methods allows to follow the motion of interfaces with surface tension immersed in an incompressible Newtonian fluid. Our method features the use of a high order fully implicit time integration scheme that circumvents the stability issues related to the explicit discretization of the capillary force when capillary effects dominate. A strategy based on a consistent Newton-Raphson linearization is introduced, and performances are enhanced by using an exact Newton variant which guarantees a third-order convergence behavior without requiring second order derivatives. The problem is approximated by mixed finite elements, while the anisotropic adaptive mesh refinements enable us to increase the computational accuracy. Numerical investigations of the convergence properties and comparisons with benchmark results provide evidence regarding the efficacy of the methodology. The robustness of the method is tested with respect to the standard explicit method, and stability is maintained for significantly larger time steps compared to those allowed by the stability condition.

Link to publisher's page
Download in pdf format
  author = {Aymen Laadhari and Gabor Székely},
  title = {Fully implicit Finite Element method for the modeling of free surface flows with surface tension effect},
  journal = {International Journal for Numerical Methods in Engineering},
  year = {2017},
  month = {February},
  pages = {},
  volume = {},
  number = {},
  keywords = {Free surface flow; Surface tension; Cubically-convergent Newton method; Finite element method; Anisotropic mesh adaptation; Adaptive time stepping},
  note = {in press}