High-order hybridizable discontinuous Galerkin formulation with fully implicit temporal schemes for the simulation of two-phase flow through porous media
Visualitza/Obre
Cita com:
hdl:2117/342787
Tipus de documentArticle
Data publicació2021-07-30
EditorJohn Wiley & Sons
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
ProjecteINTEGRACION DE NUEVAS METODOLOGIAS PARA LA PLANIFICACION Y ANALISIS MEDIOAMBIENTAL DE YACIMIENTOS PETROLIFEROS (MINECO-CTM2014-55014-C3-3-R)
SIMULACION IN VIVO DEL EFECTO DE LA HIPOXIA Y LA DOSIS DEL FARMACO EN EL CRECIMIENTO DEL GLIOBLASTOMA (AEI-PGC2018-097257-B-C33)
SIMULACION IN VIVO DEL EFECTO DE LA HIPOXIA Y LA DOSIS DEL FARMACO EN EL CRECIMIENTO DEL GLIOBLASTOMA (AEI-PGC2018-097257-B-C33)
Abstract
We present a memory-efficient high-order hybridizable discontinuous Galerkin (HDG) formulation coupled with high-order fully implicit Runge-Kutta schemes for immiscible and incompressible two-phase flow through porous media. To obtain the same high-order accuracy in space and time, we propose using high-order temporal schemes that allow using large time steps. Therefore, we require unconditionally stable temporal schemes for any combination of element size, polynomial degree and time step. Specifically, we use the Radau IIA and Gauss-Legendre schemes, which are unconditionally stable, achieve high-order accuracy with few stages, and do not suffer order reduction in this problem. To reduce the memory footprint of coupling these spatial and temporal high-order schemes, we rewrite the non-linear system. In this way, we achieve a better sparsity pattern of the Jacobian matrix and less coupling between stages. Furthermore, we propose a fix-point iterative method to further reduce the memory consumption. The saturation solution may present sharp fronts. Thus, the high-order approximation may contain spurious oscillations. To reduce them, we introduce artificial viscosity. We detect the elements with high-oscillations using a computationally efficient shock sensor obtained from the saturation solution and the post-processed saturation of HDG. Finally, we present several examples to assess the capabilities of our formulation.
CitacióCosta, A.; Ruiz-Gironès, E.; Sarrate, J. High-order hybridizable discontinuous Galerkin formulation with fully implicit temporal schemes for the simulation of two-phase flow through porous media. "International journal for numerical methods in engineering", 30 Juliol 2021, vol. 122, núm. 14, p. 3583-3612.
ISSN0029-5981
Versió de l'editorhttps://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.6674
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
CASE_ELOI.pdf | 2,906Mb | Visualitza/Obre |