A monotonicity-preserving higher-order accurate finite-volume method for Kapila's two-fluid flow model

R. de Böck (Corresponding author), A.S. Tijsseling, Barry Koren

Research output: Contribution to journalArticleAcademicpeer-review

2 Downloads (Pure)

Abstract

In preparation of the study of liquefied natural gas (LNG) sloshing in ships and vehicles, we model and numerically analyze compressible two-fluid flow. We consider a five-equation two-fluid flow model, assuming velocity and pressure continuity across two-fluid interfaces, with a separate equation to track the interfaces. The system of partial differential equations is hyperbolic and quasi-conservative. It is discretized in space with a tailor-made third-order accurate finite-volume method, employing an HLLC approximate Riemann solver. The third-order accuracy is obtained through spatial reconstruction with a limiter function, for which a novel formulation is presented. The non-homogeneous term is handled in a way consistent with the HLLC treatment of the convection operator. We study the one-dimensional case of a liquid column impacting onto a gas pocket entrapped at a solid wall. It mimics the impact of a breaking wave in an LNG containment system, where a gas pocket is entrapped at the tank wall below the wave crest. Furthermore, the impact of a shock wave on a gas bubble containing the heavy gas R22, immersed in air, is simulated in two dimensions and compared with experimental results. The numerical scheme is shown to be higher-order accurate in space and capable of capturing the important characteristics of compressible two-fluid flow.

Original languageEnglish
Article number104272
Number of pages16
JournalComputers and Fluids
Volume193
DOIs
Publication statusPublished - 30 Oct 2019

Fingerprint

Finite volume method
Flow of fluids
Liquefied natural gas
Gases
Liquid sloshing
Limiters
Shock waves
Partial differential equations
Ships
Fluids
Liquids
Air

Keywords

  • Finite-volume method
  • Hyperbolic system
  • Limiter functions
  • MUSCL
  • Source-term treatment
  • Two-fluid model
  • Wave impacts

Cite this

@article{998e4ca4383b4902b2291b179da9b809,
title = "A monotonicity-preserving higher-order accurate finite-volume method for Kapila's two-fluid flow model",
abstract = "In preparation of the study of liquefied natural gas (LNG) sloshing in ships and vehicles, we model and numerically analyze compressible two-fluid flow. We consider a five-equation two-fluid flow model, assuming velocity and pressure continuity across two-fluid interfaces, with a separate equation to track the interfaces. The system of partial differential equations is hyperbolic and quasi-conservative. It is discretized in space with a tailor-made third-order accurate finite-volume method, employing an HLLC approximate Riemann solver. The third-order accuracy is obtained through spatial reconstruction with a limiter function, for which a novel formulation is presented. The non-homogeneous term is handled in a way consistent with the HLLC treatment of the convection operator. We study the one-dimensional case of a liquid column impacting onto a gas pocket entrapped at a solid wall. It mimics the impact of a breaking wave in an LNG containment system, where a gas pocket is entrapped at the tank wall below the wave crest. Furthermore, the impact of a shock wave on a gas bubble containing the heavy gas R22, immersed in air, is simulated in two dimensions and compared with experimental results. The numerical scheme is shown to be higher-order accurate in space and capable of capturing the important characteristics of compressible two-fluid flow.",
keywords = "Finite-volume method, Hyperbolic system, Limiter functions, MUSCL, Source-term treatment, Two-fluid model, Wave impacts",
author = "{de B{\"o}ck}, R. and A.S. Tijsseling and Barry Koren",
year = "2019",
month = "10",
day = "30",
doi = "10.1016/j.compfluid.2019.104272",
language = "English",
volume = "193",
journal = "Computers & Fluids",
issn = "0045-7930",
publisher = "Elsevier",

}

TY - JOUR

T1 - A monotonicity-preserving higher-order accurate finite-volume method for Kapila's two-fluid flow model

AU - de Böck, R.

AU - Tijsseling, A.S.

AU - Koren, Barry

PY - 2019/10/30

Y1 - 2019/10/30

N2 - In preparation of the study of liquefied natural gas (LNG) sloshing in ships and vehicles, we model and numerically analyze compressible two-fluid flow. We consider a five-equation two-fluid flow model, assuming velocity and pressure continuity across two-fluid interfaces, with a separate equation to track the interfaces. The system of partial differential equations is hyperbolic and quasi-conservative. It is discretized in space with a tailor-made third-order accurate finite-volume method, employing an HLLC approximate Riemann solver. The third-order accuracy is obtained through spatial reconstruction with a limiter function, for which a novel formulation is presented. The non-homogeneous term is handled in a way consistent with the HLLC treatment of the convection operator. We study the one-dimensional case of a liquid column impacting onto a gas pocket entrapped at a solid wall. It mimics the impact of a breaking wave in an LNG containment system, where a gas pocket is entrapped at the tank wall below the wave crest. Furthermore, the impact of a shock wave on a gas bubble containing the heavy gas R22, immersed in air, is simulated in two dimensions and compared with experimental results. The numerical scheme is shown to be higher-order accurate in space and capable of capturing the important characteristics of compressible two-fluid flow.

AB - In preparation of the study of liquefied natural gas (LNG) sloshing in ships and vehicles, we model and numerically analyze compressible two-fluid flow. We consider a five-equation two-fluid flow model, assuming velocity and pressure continuity across two-fluid interfaces, with a separate equation to track the interfaces. The system of partial differential equations is hyperbolic and quasi-conservative. It is discretized in space with a tailor-made third-order accurate finite-volume method, employing an HLLC approximate Riemann solver. The third-order accuracy is obtained through spatial reconstruction with a limiter function, for which a novel formulation is presented. The non-homogeneous term is handled in a way consistent with the HLLC treatment of the convection operator. We study the one-dimensional case of a liquid column impacting onto a gas pocket entrapped at a solid wall. It mimics the impact of a breaking wave in an LNG containment system, where a gas pocket is entrapped at the tank wall below the wave crest. Furthermore, the impact of a shock wave on a gas bubble containing the heavy gas R22, immersed in air, is simulated in two dimensions and compared with experimental results. The numerical scheme is shown to be higher-order accurate in space and capable of capturing the important characteristics of compressible two-fluid flow.

KW - Finite-volume method

KW - Hyperbolic system

KW - Limiter functions

KW - MUSCL

KW - Source-term treatment

KW - Two-fluid model

KW - Wave impacts

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

U2 - 10.1016/j.compfluid.2019.104272

DO - 10.1016/j.compfluid.2019.104272

M3 - Article

AN - SCOPUS:85072562798

VL - 193

JO - Computers & Fluids

JF - Computers & Fluids

SN - 0045-7930

M1 - 104272

ER -