Repository logo
Log in(current)
Repository logoMIT Open ScholarshipDSpace@MIT
  1. Home
  2. MIT Open Access Articles
  3. MIT Open Access Articles
  4. Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form

Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form

Thumbnail Image
Download
Name

10596_2020_9999_ReferencePDF.pdf

Size

1.67 MB

Format

Adobe PDF

Checksum (MD5)

65d51b9aa6e2271739faa04f91cc61d7

Download all files submitted through automated deposit
art_1395197836799702359.zip (1.44 MB)
Author(s)
Jayasinghe, Savithru
•
Darmofal, David L
•
Allmaras, Steven R
•
Dow, Eric
•
Galbraith, Marshall C
Date Issued
October 9, 2020
Publisher
Springer International Publishing
Version
Author's final manuscript
Abstract
Abstract High-order discretizations have become increasingly popular across a wide range of applications, including reservoir simulation. However, the lack of stability and robustness of these discretizations for advection-dominant problems prevent them from being widely adopted. This paper presents work towards improving the stability and robustness of the discontinuous Galerkin (DG) finite element scheme, for advection-dominant two-phase flow problems in particular. A linearized analysis of the two-phase flow equations is used to show that a standard DG discretization of the two-phase flow equations in mass conservation form results in a neutrally stable semi-discrete system in the advection-dominant limit. Furthermore, the analysis is also used to propose additional terms to the DG method which linearly stabilize the discretization. These additional terms are derived by comparing the linearized equations in mass conservation form against an upwinded pressure-saturation form of the equations. Next, a partial differential equation-based artificial viscosity method is proposed for the Buckley-Leverett and two-phase flow equations, as a means of mitigating Gibbs oscillations in high-order discretizations and ensuring convergence to physical solutions. The modified DG method with artificial viscosity is demonstrated on a two-phase flow problem with heterogeneous rock permeabilities, where the high-order discretizations significantly outperform a conventional first-order approach in terms of computational cost required to achieve a given level of error in an output of interest.
MIT Department
Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
Massachusetts Institute of Technology. Gas Turbine Laboratory
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
http://creativecommons.org/licenses/by-nc-sa/4.0/
Persistent DSpace Link
https://hdl.handle.net/1721.1/131988
DOI of Published Version
https://doi.org/10.1007/s10596-020-09999-6
Repository logo
PrivacyPermissionsAccessibilityContact us
Repository logo
Notify us about copyright concerns.