Consistent MPFA Discretization for Flow in the Presence of Gravity
A standard practice used in the industry to discretizing the gravity term in the two‐phase Darcy flow equations is to apply an upwind strategy. In this paper, we show that this can give a persistent unphysical flux field and an incorrect pressure distribution. As a solution to this problem, we prese...
Saved in:
Published in | Water resources research Vol. 55; no. 12; pp. 10105 - 10118 |
---|---|
Main Authors | , , , |
Format | Journal Article |
Language | English |
Published |
01.12.2019
|
Online Access | Get full text |
Cover
Loading…
Abstract | A standard practice used in the industry to discretizing the gravity term in the two‐phase Darcy flow equations is to apply an upwind strategy. In this paper, we show that this can give a persistent unphysical flux field and an incorrect pressure distribution. As a solution to this problem, we present a new consistent discretization of flow, termed Gravitationally Consistent Multipoint Flux Approximation (GCMPFA), which is valid for both single‐ and two‐phase flows. The discretization is based on the idea that the gravitational term in the flow equations is treated as part of the discrete flux operator and not as a right‐hand side. Here, the traditional formulation representing pressure as a potential is extended to the case including gravity by introducing an additional set of right‐hand side to the local linear system solved in the MPFA construction, thus obtaining an expression of the fluxes in terms of jumps in cell‐center gravities. Numerical examples showing the convergence of the method are provided for both single‐ and two‐phase flows. For two‐phase flow, we show how our new method is capable of eliminating the unphysical fluxes arising when using a standard upwind scheme, thus converging to the correct pressure distribution.
Key Points
Standard discretization of two‐phase Darcy flow in the presence of gravity is shown to create a persistent unphysical flux field
We present a new consistent discretization of flow, which treats the gravity term as part of the discrete flux operator
Our method is generally second‐order convergent and is capable of eliminating the flux field arising when using a standard discretization |
---|---|
AbstractList | A standard practice used in the industry to discretizing the gravity term in the two‐phase Darcy flow equations is to apply an upwind strategy. In this paper, we show that this can give a persistent unphysical flux field and an incorrect pressure distribution. As a solution to this problem, we present a new consistent discretization of flow, termed Gravitationally Consistent Multipoint Flux Approximation (GCMPFA), which is valid for both single‐ and two‐phase flows. The discretization is based on the idea that the gravitational term in the flow equations is treated as part of the discrete flux operator and not as a right‐hand side. Here, the traditional formulation representing pressure as a potential is extended to the case including gravity by introducing an additional set of right‐hand side to the local linear system solved in the MPFA construction, thus obtaining an expression of the fluxes in terms of jumps in cell‐center gravities. Numerical examples showing the convergence of the method are provided for both single‐ and two‐phase flows. For two‐phase flow, we show how our new method is capable of eliminating the unphysical fluxes arising when using a standard upwind scheme, thus converging to the correct pressure distribution.
Key Points
Standard discretization of two‐phase Darcy flow in the presence of gravity is shown to create a persistent unphysical flux field
We present a new consistent discretization of flow, which treats the gravity term as part of the discrete flux operator
Our method is generally second‐order convergent and is capable of eliminating the flux field arising when using a standard discretization |
Author | Keilegavlen, E. Berre, I. Nordbotten, J. M. Starnoni, M. |
Author_xml | – sequence: 1 givenname: M. orcidid: 0000-0002-8552-6997 surname: Starnoni fullname: Starnoni, M. email: michele.starnoni@uib.no organization: University of Bergen – sequence: 2 givenname: I. orcidid: 0000-0002-0212-7959 surname: Berre fullname: Berre, I. organization: University of Bergen – sequence: 3 givenname: E. orcidid: 0000-0002-0333-9507 surname: Keilegavlen fullname: Keilegavlen, E. organization: University of Bergen – sequence: 4 givenname: J. M. orcidid: 0000-0003-1455-5704 surname: Nordbotten fullname: Nordbotten, J. M. organization: University of Bergen |
BookMark | eNpNj9FKwzAUhoNMcJve-QB5geo5PUnTXM5qpzBxFGWXJW0TjNRUmuKYT-9EL7z64Ofjh2_BZmEIlrFLhCuEVF-ngHpXQSopFydsjlqIRGlFMzYHEJQgaXXGFjG-AaCQmZqzm2II0cfJhok_bssVv_WxHe3kv8zkh8DdMPKyH_bcBz69Wr4dbbShtXxwfD2aTz8dztmpM320F39cspfy7rm4TzZP64ditUkMZTklrhEmaxBllmvTNp2C3AAeV220QIsdyc6BTFuF1jQuPQoyJ42qUSg7dLRk9Pu797091B-jfzfjoUaof-Lr__H1riqqVBAQfQMqDU-s |
CitedBy_id | crossref_primary_10_1016_j_jaerosci_2022_106129 crossref_primary_10_3390_en14102757 crossref_primary_10_1002_nme_6744 |
ContentType | Journal Article |
Copyright | 2019. The Authors. |
Copyright_xml | – notice: 2019. The Authors. |
DBID | 24P WIN |
DOI | 10.1029/2019WR025384 |
DatabaseName | Wiley_OA刊 Wiley Online Library Free Content |
DatabaseTitleList | |
Database_xml | – sequence: 1 dbid: 24P name: Wiley Open Access url: https://authorservices.wiley.com/open-science/open-access/browse-journals.html sourceTypes: Publisher |
DeliveryMethod | fulltext_linktorsrc |
Discipline | Geography Economics |
EISSN | 1944-7973 |
EndPage | 10118 |
ExternalDocumentID | WRCR24303 |
Genre | article |
GroupedDBID | -~X ..I .DC 05W 0R~ 123 1OB 1OC 24P 31~ 33P 3V. 50Y 5VS 6TJ 7WY 7XC 8-1 8CJ 8FE 8FG 8FH 8FL 8G5 8R4 8R5 8WZ A00 A6W AAESR AAHBH AAHHS AAIHA AAIKC AAMNW AANLZ AASGY AAXRX AAYJJ AAYOK AAZKR ABCUV ABJCF ABJNI ABPPZ ABTAH ABUWG ACAHQ ACBWZ ACCFJ ACCZN ACGFO ACGFS ACIWK ACKIV ACNCT ACPOU ACPRK ACXBN ACXQS ADBBV ADEOM ADKYN ADMGS ADOZA ADXAS ADZMN AEEZP AEIGN AENEX AEQDE AETEA AEUYR AFBPY AFGKR AFKRA AFPWT AFRAH AFZJQ AIDBO AIURR AIWBW AJBDE ALMA_UNASSIGNED_HOLDINGS ALUQN ALXUD AMYDB ASPBG ATCPS AVWKF AZFZN AZQEC AZVAB BDRZF BENPR BEZIV BFHJK BGLVJ BHPHI BKSAR BMXJE BPHCQ BRXPI CCPQU CS3 D0L D1J DCZOG DDYGU DPXWK DRFUL DRSTM DU5 DWQXO EBS EJD F5P FEDTE FRNLG G-S GNUQQ GODZA GROUPED_ABI_INFORM_COMPLETE GUQSH HCIFZ HVGLF HZ~ K60 K6~ L6V LATKE LEEKS LITHE LK5 LOXES LUTES LYRES M0C M2O M7R M7S MEWTI MSFUL MSSTM MVM MW2 MXFUL MXSTM MY~ O9- OHT OK1 P-X P2P P2W PALCI PATMY PCBAR PQBIZ PQBZA PQQKQ PROAC PTHSS PYCSY Q2X R.K RIWAO RJQFR ROL SAMSI SUPJJ TAE TN5 TWZ UQL VJK VOH WBKPD WIN WXSBR WYJ XOL XSW YHZ YV5 ZCG ZY4 ZZTAW ~02 ~KM ~OA ~~A |
ID | FETCH-LOGICAL-a3683-fb4a6b115689acbd708a01fb49a941e1d35df052c71eabf2d70583917b715d1f3 |
IEDL.DBID | 24P |
ISSN | 0043-1397 |
IngestDate | Sat Aug 24 01:08:02 EDT 2024 |
IsDoiOpenAccess | true |
IsOpenAccess | true |
IsPeerReviewed | true |
IsScholarly | true |
Issue | 12 |
Language | English |
License | Attribution |
LinkModel | DirectLink |
MergedId | FETCHMERGED-LOGICAL-a3683-fb4a6b115689acbd708a01fb49a941e1d35df052c71eabf2d70583917b715d1f3 |
ORCID | 0000-0003-1455-5704 0000-0002-0333-9507 0000-0002-8552-6997 0000-0002-0212-7959 |
OpenAccessLink | https://onlinelibrary.wiley.com/doi/abs/10.1029%2F2019WR025384 |
PageCount | 14 |
ParticipantIDs | wiley_primary_10_1029_2019WR025384_WRCR24303 |
PublicationCentury | 2000 |
PublicationDate | December 2019 |
PublicationDateYYYYMMDD | 2019-12-01 |
PublicationDate_xml | – month: 12 year: 2019 text: December 2019 |
PublicationDecade | 2010 |
PublicationTitle | Water resources research |
PublicationYear | 2019 |
References | 2007; 106 2011 2019; 77 2002; 6 2006b; 104 2002; 2 2005; 21 1994 2008; 346 2012; 17 2006; 2 2007; 11 1996; 127 1985; 47 1998a; 19 1998b; 19 2010; 116 2006; 44 2000; 160 2019 2007; 6 2008; 24 1998; 2 2011; 47 2007; 44 2006a; 22 |
References_xml | – volume: 22 start-page: 1438 issue: 6 year: 2006a end-page: 1454 article-title: Convergence of multipoint flux approximations on quadrilateral grids publication-title: Numerical Methods for Partial Differential Equations: An International Journal – year: 2011 – volume: 346 start-page: 1007 year: 2008 end-page: 1012 article-title: Convergence of the finite volume MPFA O scheme for heterogeneous anisotropic diffusion problems on general meshes publication-title: Comptes rendus de l'Académie des Sciences – volume: 2 start-page: 325 issue: 4 year: 2002 end-page: 353 article-title: Mathematical and numerical study of an industrial scheme for two‐phase flows in porous media under gravity publication-title: Computational Methods in Applied Mathematics – volume: 116 start-page: 1 issue: 1 year: 2010 end-page: 29 article-title: First‐order convergence of multi‐point flux approximation on triangular grids and comparison with mixed finite element methods publication-title: Numerische Mathematik – volume: 106 start-page: 255 issue: 2 year: 2007 end-page: 288 article-title: Monotonicity of control volume methods publication-title: Numerische Mathematik – volume: 127 start-page: 2 issue: 1 year: 1996 end-page: 14 article-title: Discretization on non‐orthogonal, quadrilateral grids for inhomogeneous, anisotropic media publication-title: Journal of computational physics – volume: 17 start-page: 779 issue: 03 year: 2012 end-page: 793 article-title: Accurate cell‐centered discretizations for modeling multiphase flow in porous media on general hexahedral and simplicial grids publication-title: SPE Journal – volume: 19 start-page: 1717 issue: 5 year: 1998b end-page: 1736 article-title: Discretization on unstructured grids for inhomogeneous, anisotropic media. part II: Discussion and numerical results publication-title: SIAM Journal on Scientific Computing – volume: 44 start-page: 2082 issue: 5 year: 2006 end-page: 2106 article-title: A multipoint flux mixed finite element method publication-title: SIAM Journal on Numerical Analysis – volume: 104 start-page: 317 issue: 3 year: 2006b end-page: 337 article-title: Robust convergence of multi point flux approximation on rough grids publication-title: Numerische Mathematik – volume: 19 start-page: 1700 issue: 5 year: 1998a end-page: 1716 article-title: Discretization on unstructured grids for inhomogeneous, anisotropic media. part i: Derivation of the methods publication-title: SIAM Journal on Scientific Computing – volume: 2 year: 2006 – year: 1994 – volume: 6 start-page: 319 issue: 1 year: 2007 end-page: 346 article-title: A multiscale mortar mixed finite element method publication-title: Multiscale Modeling & Simulation – volume: 47 start-page: 217 issue: 2 year: 1985 end-page: 235 article-title: Two families of mixed finite elements for second order elliptic problems publication-title: Numerische Mathematik – volume: 160 start-page: 1 issue: 1 year: 2000 end-page: 28 article-title: M‐matrix flux splitting for general full tensor discretization operators on structured and unstructured grids publication-title: Journal of Computational Physics – volume: 44 start-page: 1103 issue: 5 year: 2007 end-page: 1119 article-title: A multiscale mortar mixed finite element method for slightly compressible flows in porous media publication-title: Journal of the Korean Mathematical Society – volume: 11 start-page: 333 issue: 4 year: 2007 end-page: 345 article-title: Convergence of a symmetric MPFA method on quadrilateral grids publication-title: Computational Geosciences – volume: 47 issue: 5 year: 2011 article-title: Compositional modeling of three‐phase flow with gravity using higher‐order finite element methods publication-title: Water Resources Research – volume: 24 start-page: 1329 issue: 5 year: 2008 end-page: 1360 article-title: A compact multipoint flux approximation method with improved robustness publication-title: Numerical Methods for Partial Differential Equations: An International Journal – volume: 2 start-page: 259 issue: 4 year: 1998 end-page: 290 article-title: Finite volume discretization with imposed flux continuity for the general tensor pressure equation publication-title: Computational geosciences – volume: 6 start-page: 405 issue: 3‐4 year: 2002 end-page: 432 article-title: An introduction to multipoint flux approximations for quadrilateral grids publication-title: Computational Geosciences – volume: 77 start-page: 1437 issue: 6 year: 2019 end-page: 1452 article-title: Multipoint flux mixed finite element methods for slightly compressible flow in porous media publication-title: Computers & Mathematics with Applications – volume: 6 start-page: 433 issue: 3‐4 year: 2002 end-page: 452 article-title: Unstructured, control‐volume distributed, full‐tensor finite‐volume schemes with flow based grids publication-title: Computational Geosciences – volume: 21 start-page: 1079 issue: 6 year: 2005 end-page: 1098 article-title: On the convergence of the multi‐point flux approximation O‐method: Numerical experiments for discontinuous permeability publication-title: Numerical Methods for Partial Differential Equations: An International Journal – year: 2019 |
SSID | ssj0014567 |
Score | 2.3829334 |
Snippet | A standard practice used in the industry to discretizing the gravity term in the two‐phase Darcy flow equations is to apply an upwind strategy. In this paper,... |
SourceID | wiley |
SourceType | Publisher |
StartPage | 10105 |
Title | Consistent MPFA Discretization for Flow in the Presence of Gravity |
URI | https://onlinelibrary.wiley.com/doi/abs/10.1029%2F2019WR025384 |
Volume | 55 |
hasFullText | 1 |
inHoldings | 1 |
isFullTextHit | |
isPrint | |
link | http://utb.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwlZ3NS8MwFMCDzoNexE_8JgePFpcmbZrjnNYhTMZwbLeSTxxIK1tF_O99Sbsxj95Kmxya15f3a94XQrfCGMmtIBHskGnEEuoipb3P1TkVO-OkCsGYw9d0MGEvs2TWHrj5XJimPsT6wM1rRtivvYJLtWyLDfgamWC5xHQMJptmbBvtANlk_quO2WjtRQA44CsPsyedNvAd5t9vzv7LpcGw5AdovyVC3GtEeIi2bHmEdlcJw0u4bhuVv_8co4fQYBMkU9Z4OMp7-HEOem_rNpsSA4Li_KP6xvMSA9rhUcgu0hZXDj8vpG8UcYIm-dNbfxC1bRAiSdOMRk4xmSp4vzQTUivDu5nsErgrpGDEEkMT47pJrDmxUrkYBiSAPYQrThJDHD1FnbIq7RnC2ukU_vioyzRwENXCUe0MIByz1DLJz9FdWInisyl1UQQXdSyKzeUqpuP-OGZg-y7-N_wS7fkHTTDIFerUiy97DSa9VjdBbr-bnZi_ |
link.rule.ids | 315,786,790,11589,27955,27956,46085,46509 |
linkProvider | Wiley-Blackwell |
linkToHtml | http://utb.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwlV07T8MwELagDGVBPMUbD4xE1LHz8FgKoUBbVVWrskV-ikooRSUI8e85O2lVRrYoOi8-nb_PvrvvELrmWovEcBLACRkHLKI2kMrlXK2VodVWSF-M2R_E3Ql7fo1e6zmnrhem0odYPbi5yPDntQtw9yBdqw04kUyALj4dAWbTlG2iLea04Jy0Mxuu0gjADpJlitlRnbryHdbfrq_-S0w9smS7aKemhLhd-XAPbZhiHzWXHcOf8F1PKn_7OUB3fsImuKYocX-YtfH9DALflHU7JQYOirP3-TeeFRi4HR769iJl8Nzix4VwkyIO0SR7GHe6QT0HIRA0TmlgJROxBOoWp1woqZNWKloE_nLBGTFE00jbVhSqhBghbQgGEfAeksiERJpYeoQaxbwwxwgrq2K48lGbKiBCVHFLldXA4ZihhonkBN34ncg_Kq2L3OeoQ56vb1c-HXVGIQPwO_2f-RVqdsf9Xt57GrycoW1nVFWGnKNGufgyF4Dvpbz0PvwF9VKcJw |
linkToPdf | http://utb.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwlZ1JSwMxFICDVlAv4oq7OXh0sJlklhxr61iXlqFY6m3IigWZljoi_ntf0mmpR2_D8HLJ4-V9ydsQuuZai8RwEsAJGQcsojaQysVcrZWh1VZIn4zZ68fdIXt6i97qBzdXCzPvD7F8cHOW4c9rZ-BTbetmA65HJnguPhqAy6YpW0cbLAZ4cJ2dWb6MIgAcJIsIsyOdOvEd1t-urv7Lpd6xZLtopyZC3JqrcA-tmXIfbS0Khj_hux5U_v5zgO78gE3QTFnhXp61cGcMdm-qupoSA4Li7GPyjcclBrTDua8uUgZPLH6YCTco4hANs_vXdjeoxyAEgsYpDaxkIpZAbnHKhZI6aaaiSeAvF5wRQzSNtG1GoUqIEdKGIBAB9pBEJiTSxNIj1CgnpTlGWFkVw42P2lQBB1HFLVVWA8IxQw0TyQm68TtRTOetLgofog55sbpdxWjQHoQMfN_p_8Sv0GbeyYqXx_7zGdp2MvO8kHPUqGZf5gK8eyUvvQp_AUxVm1A |
openUrl | ctx_ver=Z39.88-2004&ctx_enc=info%3Aofi%2Fenc%3AUTF-8&rfr_id=info%3Asid%2Fsummon.serialssolutions.com&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=Consistent+MPFA+Discretization+for+Flow+in+the+Presence+of+Gravity&rft.jtitle=Water+resources+research&rft.au=Starnoni%2C+M.&rft.au=Berre%2C+I.&rft.au=Keilegavlen%2C+E.&rft.au=Nordbotten%2C+J.+M.&rft.date=2019-12-01&rft.issn=0043-1397&rft.eissn=1944-7973&rft.volume=55&rft.issue=12&rft.spage=10105&rft.epage=10118&rft_id=info:doi/10.1029%2F2019WR025384&rft.externalDBID=10.1029%252F2019WR025384&rft.externalDocID=WRCR24303 |
thumbnail_l | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/lc.gif&issn=0043-1397&client=summon |
thumbnail_m | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/mc.gif&issn=0043-1397&client=summon |
thumbnail_s | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/sc.gif&issn=0043-1397&client=summon |