VOF evaluation of the surface tension by using variational representation and Galerkin interpolation projection

被引:6
作者
Chirco, L. [1 ]
Da Via, R. [1 ]
Manservisi, S. [1 ]
机构
[1] Univ Bologna, Lab Montecuccolino, DIN, Via Colli 16, I-40136 Bologna, Italy
关键词
Navier-Stokes equations; Finite element methods; VOF; Two-phase flows; Capillary force; FINITE-ELEMENT-METHOD; LEVEL SET APPROACH; 2-PHASE; VOLUME; TRACKING; FORCE; FORMULATION; SIMULATION; INTERFACES; ALGORITHM;
D O I
10.1016/j.jcp.2019.06.036
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this work we propose a variational approach with cell-to-point Galerkin projections for studying two-phase interface advection problems dominated by surface tension. A Volume Of Fluid (VOF) algorithm is used for tracking and locating the evolution of the two-phase interface on a Cartesian grid and a finite element numerical scheme for solving the velocity-pressure state. The velocity field that drives the evolution of this interface is computed from the weak form of the Navier-Stokes equation where the surface tension force is represented in variational form by the continuous surface force (CSF) and continuous surface stress (CSS) methods. Standard numerical approaches solve the strong form of the Navier-Stokes equations and define the CSS term by taking the divergence of the surface tension tensor. This computation of the divergence term results in a singular force which is difficult to compute when the grid is refined since the tensor is computed in a discontinuous cell-by-cell way. In this work we use the variational formulation of the Navier-Stokes equation and avoid differentiation. The tensor, which is a function of the unit normal, is evaluated over regular Sobolev spaces by using a cell-to-point Galerkin projection. This allows a regular piece-wise continuous representation of the surface tensor and the unit normal based on the VOF reconstruction. In standard approaches the CSF surface force is computed by using the curvature, which is the divergence of the unit normal. In this paper we recover the curvature with point-wise Galerkin projection avoiding direct differentiation. Tests on convergence for two and three-dimension in the static and dynamical cases are reported to show the correct representation in the desired spaces. This method is also natural for coupling non uniform grid computation of the fluid with Cartesian grid of the VOF algorithm. (C) 2019 Elsevier Inc. All rights reserved.
引用
收藏
页码:537 / 562
页数:26
相关论文
共 38 条
[1]   Interface tracking with dynamically-redistributed surface markers in unstructured quadrangular grids [J].
Aubert, F. ;
Aulisa, E. ;
Manservisi, S. ;
Scardovelli, R. .
COMPUTERS & FLUIDS, 2006, 35 (10) :1332-1343
[2]   A novel representation of the surface tension force for two-phase flow with reduced spurious currents [J].
Aulisa, E. ;
Manservisi, S. ;
Scardovelli, R. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2006, 195 (44-47) :6239-6257
[3]   A mixed markers and volume-of-fluid method for the reconstruction and advection of interfaces in two-phase and free-boundary flows [J].
Aulisa, E ;
Manservisi, S ;
Scardovelli, R .
JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 188 (02) :611-639
[4]   A Multilevel Domain Decomposition Approach for Studying Coupled Flow Applications [J].
Aulisa, E. ;
Cervone, A. ;
Manservisi, S. ;
Seshaiyer, P. .
COMMUNICATIONS IN COMPUTATIONAL PHYSICS, 2009, 6 (02) :319-341
[5]   VOFI - A library to initialize the volume fraction scalar field [J].
Bna, S. ;
Manservisi, S. ;
Scardovelli, R. ;
Yecko, P. ;
Zaleski, S. .
COMPUTER PHYSICS COMMUNICATIONS, 2016, 200 :291-299
[6]   A CONTINUUM METHOD FOR MODELING SURFACE-TENSION [J].
BRACKBILL, JU ;
KOTHE, DB ;
ZEMACH, C .
JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 100 (02) :335-354
[7]   Variational formulations for surface tension, capillarity and wetting [J].
Buscaglia, Gustavo C. ;
Ausas, Roberto F. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2011, 200 (45-46) :3011-3025
[8]   A projection method for coupling two-phase VOF and fluid structure interaction simulations [J].
Cerroni, Daniele ;
DaVia, Roberto ;
Manservisi, Sandro .
JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 354 :646-671
[9]   An optimal constrained approach for divergence-free velocity interpolation and multilevel VOF method [J].
Cervone, Antonio ;
Manservisi, Sandro ;
Scardovelli, Ruben .
COMPUTERS & FLUIDS, 2011, 47 (01) :101-114
[10]   An extended finite element method for two-phase fluids [J].
Chessa, J ;
Belytschko, T .
JOURNAL OF APPLIED MECHANICS-TRANSACTIONS OF THE ASME, 2003, 70 (01) :10-17