Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form
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...
| Main Authors: | , , , , |
|---|---|
| Other Authors: | |
| Format: | Article |
| Language: | English |
| Published: |
Springer International Publishing
2021
|
| Online Access: | https://hdl.handle.net/1721.1/131988 |
| _version_ | 1811096952067063808 |
|---|---|
| author | Jayasinghe, Savithru Darmofal, David L Allmaras, Steven R Dow, Eric Galbraith, Marshall C |
| author2 | Massachusetts Institute of Technology. Department of Aeronautics and Astronautics |
| author_facet | Massachusetts Institute of Technology. Department of Aeronautics and Astronautics Jayasinghe, Savithru Darmofal, David L Allmaras, Steven R Dow, Eric Galbraith, Marshall C |
| author_sort | Jayasinghe, Savithru |
| collection | MIT |
| description | 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. |
| first_indexed | 2024-09-23T16:51:59Z |
| format | Article |
| id | mit-1721.1/131988 |
| institution | Massachusetts Institute of Technology |
| language | English |
| last_indexed | 2024-09-23T16:51:59Z |
| publishDate | 2021 |
| publisher | Springer International Publishing |
| record_format | dspace |
| spelling | mit-1721.1/1319882023-01-20T20:29:40Z Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form Jayasinghe, Savithru Darmofal, David L Allmaras, Steven R Dow, Eric Galbraith, Marshall C Massachusetts Institute of Technology. Department of Aeronautics and Astronautics Massachusetts Institute of Technology. Gas Turbine Laboratory 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. 2021-09-20T17:41:17Z 2021-09-20T17:41:17Z 2020-10-09 2021-01-19T04:24:34Z Article http://purl.org/eprint/type/JournalArticle https://hdl.handle.net/1721.1/131988 en https://doi.org/10.1007/s10596-020-09999-6 Creative Commons Attribution-Noncommercial-Share Alike http://creativecommons.org/licenses/by-nc-sa/4.0/ Springer Nature Switzerland AG application/pdf Springer International Publishing Springer International Publishing |
| spellingShingle | Jayasinghe, Savithru Darmofal, David L Allmaras, Steven R Dow, Eric Galbraith, Marshall C Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title | Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title_full | Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title_fullStr | Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title_full_unstemmed | Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title_short | Upwinding and artificial viscosity for robust discontinuous Galerkin schemes of two-phase flow in mass conservation form |
| title_sort | upwinding and artificial viscosity for robust discontinuous galerkin schemes of two phase flow in mass conservation form |
| url | https://hdl.handle.net/1721.1/131988 |
| work_keys_str_mv | AT jayasinghesavithru upwindingandartificialviscosityforrobustdiscontinuousgalerkinschemesoftwophaseflowinmassconservationform AT darmofaldavidl upwindingandartificialviscosityforrobustdiscontinuousgalerkinschemesoftwophaseflowinmassconservationform AT allmarasstevenr upwindingandartificialviscosityforrobustdiscontinuousgalerkinschemesoftwophaseflowinmassconservationform AT doweric upwindingandartificialviscosityforrobustdiscontinuousgalerkinschemesoftwophaseflowinmassconservationform AT galbraithmarshallc upwindingandartificialviscosityforrobustdiscontinuousgalerkinschemesoftwophaseflowinmassconservationform |