Finite element framework for describing dynamic wetting phenomena

Finite element framework for describing dynamic wetting phenomena
复制标题

用于描述动态润湿现象的有限元框架

DOI:
10.1002/fld.2603
复制
发表时间:
2011
影响因子:
1.8
通讯作者:
Sprittles J
Sprittles J
中科院分区:
工程技术4区
文献类型:
--
作者:
Sprittles J

文献摘要

相似文献

考虑了动态润湿现象的有限元模拟,需要计算由液体-流体自由表面和液体-固体界面相交所限定的区域中的流动,三相接触线穿过固体。对于这类流动,文献中使用了不同的有限元法(FEM)实现,在某些情况下,这些产生了明显矛盾的结果。在本文中,开发了一个强大的动态润湿流的有限元模拟框架,通过一贯坚持的有限元方法,没有留下任何空间为特设的“可选”的变化,在这些流量的数值处理。所开发的方法使得有可能进行收敛研究,评估所需的空间分辨率,以达到预设的精度,并提供相应的基准计算。这种分析允许识别数值伪影,这在以前被解释为物理效应,并表明,使用边界条件的“强”实现来抑制数值误差在计算域的其他地方产生更大和更少的可检测误差。我们为给定的一组无量纲相似性参数提供了关于数值方案所需的空间分辨率的实用建议,并给出了一个用户友好的分步指南,指定了整个实现,这使得读者可以轻松地重现所有呈现的结果,包括基准计算。它还示出了如何开发的框架容纳概括的数学模型占额外的物理效应,如表面张力梯度。版权所有© 2011约翰威利父子有限公司.
The finite element simulation of dynamic wetting phenomena, requiring the computation of flow in a domain confined by intersecting a liquid–fluid free surface and a liquid–solid interface, with the three‐phase contact line moving across the solid, is considered. For this class of flows, different finite element method (FEM) implementations have been used in the literature, and in some cases, these produced apparently contradictory results. In the present paper, a robust framework for the FEM simulation of dynamic wetting flows is developed, which, by consistently adhering to the FEM methodology, leaves no room forad hoc‘optional’ variations in the numerical handling of these flows. The developed approach makes it possible to conduct a convergence study, assess the spatial resolution required to achieve a preset accuracy and provide the corresponding benchmark calculations. This analysis allows one to identify numerical artefacts, which had previously been interpreted as physical effects, and demonstrates that suppressing numerical errors using a ‘strong’ implementation of a boundary condition creates bigger and less detectable errors elsewhere in the computational domain. We provide practical recommendations on the spatial resolution required by a numerical scheme for a given set of non‐dimensional similarity parameters and give a user‐friendly step‐by‐step guide specifying the entire implementation, which allows the reader to easily reproduce all presented results including the benchmark calculations. It is also shown how the developed framework accommodates generalizations of the mathematical model accounting for additional physical effects, such as gradients in surface tensions. Copyright © 2011 John Wiley & Sons, Ltd.