Publication Details

Category Text Publication
Reference Category Journals
DOI 10.1007/s12665-017-6970-5
Document Shareable Link
Title (Primary) Comparing global and local implementations of nonlinear complementary problems for the modeling of multi-component two-phase flow with phase change phenomena
Author Huang, Y.; Nagel, T.; Shao, H. ORCID logo
Source Titel Environmental Earth Sciences
Year 2017
Department ENVINF
Volume 76
Issue 18
Page From art. 643
Language englisch
Keywords Phase transition; Multiphase flow; Nonlinear complementary problem; Numerical efficiency and accuracy; OpenGeoSys
UFZ wide themes RU5
Abstract Compositional multiphase flow is considered to be one of the fundamental physical processes in the field of water resources research. The strong nonlinearity and discontinuity emerging from phase transition phenomena pose a serious challenge for numerical modeling. Recently, Lauser et al. (Adv Water Resour 34(8):957–966, 2011) have proposed a numerical scheme, namely the nonlinear complementary problem (NCP), to handle this strong nonlinearity. In this work, the NCP is implemented at both local and global levels of a finite element algorithm. In the former case, the NCP is integrated into the local thermodynamic equilibrium calculation, while in the latter one, it is formulated as one of the governing equations. The two different formulations have been investigated through three well-established benchmarks and analyzed for their efficiency and robustness. It is found that both globally and locally implemented NCP formulations are numerically more efficient and robust in comparison with traditional primary variable switching approach. In homogeneous media, the globally implemented NCP formulation leads to an approximately 20% faster simulation compared to the local NCP. This is because a nested Newton iteration for the local phase state identification can be avoided, and thus, the overall computational resources are saved accordingly. However, for problems involving strongly heterogeneous media, the locally integrated NCP formulation suppresses numerical oscillations and delivers more accurate and robust results, especially at the phase boundary.
Persistent UFZ Identifier https://www.ufz.de/index.php?en=20939&ufzPublicationIdentifier=19394
Huang, Y., Nagel, T., Shao, H. (2017):
Comparing global and local implementations of nonlinear complementary problems for the modeling of multi-component two-phase flow with phase change phenomena
Environ. Earth Sci. 76 (18), art. 643 10.1007/s12665-017-6970-5