Toggle Main Menu Toggle Search

Open Access padlockePrints

A novel Godunov-type scheme for free-surface flows with artificial compressibility

Lookup NU author(s): Dr Shannon Leakey, Dr Vassilis Glenis, Dr Caspar HewettORCiD



This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND).


Variable density incompressible flows are governed by parabolic equations. The artificial compressibility method makes these equations hyperbolic-type, which means that they can be solved using techniques developed for compressible flows, such as Godunov-type schemes. While the artificial compressibility method is well-established, its application to variable density flows has been largely neglected in the literature. This paper harnesses recent advances in the wider field by applying a more robust Riemann solver and a more easily parallelisable time discretisation to the variable density equations than previously. We also develop a new method for calculating the pressure gradient as part of the second-order reconstruction step. Based on a rearrangement of the momentum equation and an exploitation of the other gradients and source terms, the new pressure gradient calculation automatically captures the pressure gradient discontinuity at the free surface. Benchmark tests demonstrate the improvements gained by this robust Riemann solver and new pressure gradient calculation.

Publication metadata

Author(s): Leakey S, Glenis V, Hewett CJM

Publication type: Article

Publication status: Published

Journal: Computer Methods in Applied Mechanics and Engineering

Year: 2022

Volume: 393

Print publication date: 01/04/2022

Online publication date: 07/03/2022

Acceptance date: 14/02/2022

Date deposited: 03/05/2022

ISSN (electronic): 0045-7825

Publisher: Elsevier


DOI: 10.1016/j.cma.2022.114763

ePrints DOI: 10.57711/8k54-g218


Altmetrics provided by Altmetric