A computational method for sharp interface advection

Royal Society Open Science - Tập 3 Số 11 - Trang 160405 - 2016
Johan Roenby1, Henrik Bredmose2, Hrvoje Jasak3,4
1DHI, Department of Ports and Offshore Technology, Agern Allé 5, 2970 Hørsholm, Denmark
2Department of Wind Energy, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
3University of Zagreb, Faculty of Mechanical Engineering and Naval Architecture, Ivana Lučića 5, Zagreb, Croatia
4Wikki Ltd, 459 Southbank House, SE1 7SJ, London, UK

Tóm tắt

We devise a numerical method for passive advection of a surface, such as the interface between two incompressible fluids, across a computational mesh. The method is called isoAdvector, and is developed for general meshes consisting of arbitrary polyhedral cells. The algorithm is based on the volume of fluid (VOF) idea of calculating the volume of one of the fluids transported across the mesh faces during a time step. The novelty of the isoAdvector concept consists of two parts. First, we exploit an isosurface concept for modelling the interface inside cells in a geometric surface reconstruction step. Second, from the reconstructed surface, we model the motion of the face–interface intersection line for a general polygonal face to obtain the time evolution within a time step of the submerged face area. Integrating this submerged area over the time step leads to an accurate estimate for the total volume of fluid transported across the face. The method was tested on simple two-dimensional and three-dimensional interface advection problems on both structured and unstructured meshes. The results are very satisfactory in terms of volume conservation, boundedness, surface sharpness and efficiency. The isoAdvector method was implemented as an OpenFOAM ® extension and is published as open source.

Từ khóa


Tài liệu tham khảo

10.1017/CBO9780511975264

10.1016/S0021-9991(03)00298-5

The OpenFOAM Foundation. OpenFOAM. See www.openfoam.org.

10.1063/1.168744

10.1016/0021-9991(81)90145-5

10.1088/1749-4699/5/1/014016

10.1016/j.jcp.2007.06.033

10.1002/fld.1776

10.1002/fld.1775

Ivey CB Moin P. 2012 Conservative volume of fluid advection method on unstructured grids in three dimensions. Center for Turbulence Research Annual Research Briefs.

Maric T Marschall H Bothe D. 2013 voFoam – A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM. (http://arxiv.org/abs/1305.3417) [physics].

10.1002/fld.3968

10.1016/j.compfluid.2014.02.001

Roenby J. isoAdvector. See www.github.com/isoadvector.

Kitware Sandia National Labs and CSimSoft ‘ParaView.’ See www.paraview.org.

Muzaferija S Peric M Sames P Schellin T. 1998 A two-fluid Navier-Stokes solver to simulate water entry. In Proceedings of the 22nd symposium on naval hydrodynamics Washington DC pp. 277–289. Washington DC: National Acadamies Press.

10.1006/jcph.1999.6276

Nielsen KB. 2003 Numerical prediction of green water loads on ships. Ph.D. thesis Technical University of Denmark Department of Mechanical Engineering.

10.1007/s10915-012-9614-7

10.1016/j.jcp.2008.12.031

10.1016/j.jcp.2012.07.019

10.1006/jcph.2000.6510

10.1006/jcph.1998.5906

10.1002/(SICI)1097-0363(19980815)28:2<357::AID-FLD750>3.0.CO;2-D

10.1002/(SICI)1097-0363(19970415)24:7<671::AID-FLD508>3.0.CO;2-9

10.1002/fld.228

10.1137/0733033

10.1016/j.jcp.2011.04.040

10.1016/j.compfluid.2005.09.003

10.1016/j.compstruc.2004.04.024

Roenby J Bredmose H Jasak H Uroic T Vallier A. 2016 Data from: A computational method for sharp interface advection. Dryad Digital Repository. (http://dx.doi.org/10.5061/dryad.66840)

Thông tin
Thông tin xuất bản

Royal Society Open Science

Tập 3 Số 11

160405

Thông tin tác giả