Artificial viscosity model to mitigate numerical artefacts at fluid interfaces with surface tension
File(s)
Published version
Author(s)
Denner, F
Evrard, F
Serfaty, R
van Wachem, B
Type
Journal Article
Abstract
The numerical onset of parasitic and spurious artefacts in the vicinity of uid interfaces with surface tension is
an important and well-recognised problem with respect to the accuracy and numerical stability of interfacial
ow simulations. Issues of particular interest are spurious capillary waves, which are spatially underresolved
by the computational mesh yet impose very restrictive time-step requirements, as well as parasitic currents,
typically the result of a numerically unbalanced curvature evaluation. We present an arti cial viscosity
model to mitigate numerical artefacts at surface-tension-dominated interfaces without adversely a ecting
the accuracy of the physical solution. The proposed methodology computes an additional interfacial shear
stress term, including an interface viscosity, based on the local ow data and uid properties that reduces the
impact of numerical artefacts and dissipates underresolved small scale interface movements. Furthermore,
the presented methodology can be readily applied to model surface shear viscosity, for instance to simulate
the dissipative e ect of surface-active substances adsorbed at the interface. The presented analysis of
numerical test cases demonstrates the e cacy of the proposed methodology in diminishing the adverse
impact of parasitic and spurious interfacial artefacts on the convergence and stability of the numerical
solution algorithm as well as on the overall accuracy of the simulation results.
an important and well-recognised problem with respect to the accuracy and numerical stability of interfacial
ow simulations. Issues of particular interest are spurious capillary waves, which are spatially underresolved
by the computational mesh yet impose very restrictive time-step requirements, as well as parasitic currents,
typically the result of a numerically unbalanced curvature evaluation. We present an arti cial viscosity
model to mitigate numerical artefacts at surface-tension-dominated interfaces without adversely a ecting
the accuracy of the physical solution. The proposed methodology computes an additional interfacial shear
stress term, including an interface viscosity, based on the local ow data and uid properties that reduces the
impact of numerical artefacts and dissipates underresolved small scale interface movements. Furthermore,
the presented methodology can be readily applied to model surface shear viscosity, for instance to simulate
the dissipative e ect of surface-active substances adsorbed at the interface. The presented analysis of
numerical test cases demonstrates the e cacy of the proposed methodology in diminishing the adverse
impact of parasitic and spurious interfacial artefacts on the convergence and stability of the numerical
solution algorithm as well as on the overall accuracy of the simulation results.
Date Issued
2016-11-15
Date Acceptance
2016-11-14
Citation
Computers & Fluids, 2016, 143, pp.59-72
ISSN
0045-7930
Publisher
Elsevier
Start Page
59
End Page
72
Journal / Book Title
Computers & Fluids
Volume
143
Copyright Statement
© 2016 The Authors. Published by Elsevier Ltd.
This is an open access article under the CC BY license. (http://creativecommons.org/licenses/by/4.0/)
This is an open access article under the CC BY license. (http://creativecommons.org/licenses/by/4.0/)
Sponsor
PETROLEO BRASILEIRO S. A. PETROBRAS
Engineering & Physical Science Research Council (EPSRC)
Grant Number
0050.0077680.12.2
EP/M021556/1
Subjects
Applied Mathematics
0102 Applied Mathematics
0915 Interdisciplinary Engineering
0913 Mechanical Engineering
Publication Status
Published