A preconditioned implicit free-surface capture scheme for large density ratio on tetrahedral grids

Xin Lv, Qingping Zou, D. E. Reeve, Yong Zhao

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

We present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems.

Original languageEnglish
Pages (from-to)215-248
Number of pages34
JournalCommunications in Computational Physics
Volume11
Issue number1
DOIs
Publication statusPublished - Jan 2012
Externally publishedYes

Fingerprint

grids
preconditioning
Runge-Kutta method
eccentrics
linear systems
fluids
matrices

Keywords

  • Free surface
  • Implicit method
  • Level set
  • Restarted GMRES
  • Tetrahedral grid
  • Unstructured finite volume method
  • VOF

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)

Cite this

A preconditioned implicit free-surface capture scheme for large density ratio on tetrahedral grids. / Lv, Xin; Zou, Qingping; Reeve, D. E.; Zhao, Yong.

In: Communications in Computational Physics, Vol. 11, No. 1, 01.2012, p. 215-248.

Research output: Contribution to journalArticle

@article{39bcb155d6234e41919f3a928b240d6c,
title = "A preconditioned implicit free-surface capture scheme for large density ratio on tetrahedral grids",
abstract = "We present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems.",
keywords = "Free surface, Implicit method, Level set, Restarted GMRES, Tetrahedral grid, Unstructured finite volume method, VOF",
author = "Xin Lv and Qingping Zou and Reeve, {D. E.} and Yong Zhao",
year = "2012",
month = "1",
doi = "10.4208/cicp.170510.290311a",
language = "English",
volume = "11",
pages = "215--248",
journal = "Communications in Computational Physics",
issn = "1815-2406",
publisher = "Global Science Press",
number = "1",

}

TY - JOUR

T1 - A preconditioned implicit free-surface capture scheme for large density ratio on tetrahedral grids

AU - Lv, Xin

AU - Zou, Qingping

AU - Reeve, D. E.

AU - Zhao, Yong

PY - 2012/1

Y1 - 2012/1

N2 - We present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems.

AB - We present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems.

KW - Free surface

KW - Implicit method

KW - Level set

KW - Restarted GMRES

KW - Tetrahedral grid

KW - Unstructured finite volume method

KW - VOF

UR - http://www.scopus.com/inward/record.url?scp=84055190657&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84055190657&partnerID=8YFLogxK

U2 - 10.4208/cicp.170510.290311a

DO - 10.4208/cicp.170510.290311a

M3 - Article

VL - 11

SP - 215

EP - 248

JO - Communications in Computational Physics

JF - Communications in Computational Physics

SN - 1815-2406

IS - 1

ER -