Loading…

A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes

We suggest a new positivity-preserving discrete duality finite volume (DDFV) scheme for anisotropic diffusion problems on general polygonal meshes. Like existing DDFV schemes, this scheme is built on the primary and dual meshes and has finite volume (FV) equations for both vertex and cell-centered u...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2018-11, Vol.372, p.773-798
Main Authors: Su, Shuai, Dong, Qiannan, Wu, Jiming
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
cited_by cdi_FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53
cites cdi_FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53
container_end_page 798
container_issue
container_start_page 773
container_title Journal of computational physics
container_volume 372
creator Su, Shuai
Dong, Qiannan
Wu, Jiming
description We suggest a new positivity-preserving discrete duality finite volume (DDFV) scheme for anisotropic diffusion problems on general polygonal meshes. Like existing DDFV schemes, this scheme is built on the primary and dual meshes and has finite volume (FV) equations for both vertex and cell-centered unknowns. The transpose of the coefficient matrix for the FV equations of the cell-centered unknowns is an M-matrix while that for the vertex unknowns is not an M-matrix but a symmetric and positive definite matrix. By employing a certain truncation technique for the vertex unknowns, the positivity-preserving property for both categories of unknowns is guaranteed. Local conservation is strictly maintained for the cell-centered unknowns and conditionally maintained for the vertex unknowns. Since the FV equations of the vertex unknowns can be solved independently, the two sets of FV equations are decoupled. In contrast to existing nonlinear positivity-preserving schemes, the new scheme requires no nonlinear iterations for linear problems. For nonlinear problems, the positivity-preserving mechanism of the new scheme is decoupled from its nonlinear iteration so that any nonlinear solver can be adopted. Moreover, the positivity of the discrete solution is proved and the well-posedness is analyzed rigorously for linear problems. In numerical experiments, the new scheme is examined extensively and compared with two positivity-preserving cell-centered schemes and a nonlinear DDFV scheme. Numerical results show that the present DDFV scheme achieves second-order accuracy and preserves the positivity of the solution for heterogeneous and anisotropic problems on severely distorted grids. The high efficiency of the scheme is also demonstrated by the comparison of computation time and number of nonlinear iterations.
doi_str_mv 10.1016/j.jcp.2018.06.052
format article
fullrecord <record><control><sourceid>proquest_cross</sourceid><recordid>TN_cdi_proquest_journals_2116622890</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><els_id>S0021999118304303</els_id><sourcerecordid>2116622890</sourcerecordid><originalsourceid>FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53</originalsourceid><addsrcrecordid>eNp9UE1LxDAQDaLguvoDvAU8t05ik7Z4EvELBC96Dm0yXVO6Tc20C_sX_NVG1rOnecO89-bxGLsUkAsQ-rrPezvlEkSVg85BySO2ElBDJkuhj9kKQIqsrmtxys6IegCoVFGt2Pcdd2jDMg3oeDM6PgXys9_5eZ9NEQnjzo8b7jzZiDNytzRDuvHOjz6tuzAsW-RkPzGNLsTk4SnMMUzeJlXXLeTDyKcY2gG3xBPe4IixGdKnYb8JY0JbpE-kc3bSNQPhxd9cs4_Hh_f75-z17enl_u41szdSzZmUTotWtYBWaYVNJWUNhS1tqaoSoHWtqIpSdlCWFptELXQLXaFqbYWqnbpZs6uDbwr1tSDNpg9LTDnISCG0lrKqIbHEgWVjIIrYmSn6bRP3RoD5rdz0JlVufis3oE2qPGluDxpM8XceoyHrcbTofEQ7Gxf8P-ofntqMkQ</addsrcrecordid><sourcetype>Aggregation Database</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>2116622890</pqid></control><display><type>article</type><title>A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes</title><source>ScienceDirect Freedom Collection</source><creator>Su, Shuai ; Dong, Qiannan ; Wu, Jiming</creator><creatorcontrib>Su, Shuai ; Dong, Qiannan ; Wu, Jiming</creatorcontrib><description>We suggest a new positivity-preserving discrete duality finite volume (DDFV) scheme for anisotropic diffusion problems on general polygonal meshes. Like existing DDFV schemes, this scheme is built on the primary and dual meshes and has finite volume (FV) equations for both vertex and cell-centered unknowns. The transpose of the coefficient matrix for the FV equations of the cell-centered unknowns is an M-matrix while that for the vertex unknowns is not an M-matrix but a symmetric and positive definite matrix. By employing a certain truncation technique for the vertex unknowns, the positivity-preserving property for both categories of unknowns is guaranteed. Local conservation is strictly maintained for the cell-centered unknowns and conditionally maintained for the vertex unknowns. Since the FV equations of the vertex unknowns can be solved independently, the two sets of FV equations are decoupled. In contrast to existing nonlinear positivity-preserving schemes, the new scheme requires no nonlinear iterations for linear problems. For nonlinear problems, the positivity-preserving mechanism of the new scheme is decoupled from its nonlinear iteration so that any nonlinear solver can be adopted. Moreover, the positivity of the discrete solution is proved and the well-posedness is analyzed rigorously for linear problems. In numerical experiments, the new scheme is examined extensively and compared with two positivity-preserving cell-centered schemes and a nonlinear DDFV scheme. Numerical results show that the present DDFV scheme achieves second-order accuracy and preserves the positivity of the solution for heterogeneous and anisotropic problems on severely distorted grids. The high efficiency of the scheme is also demonstrated by the comparison of computation time and number of nonlinear iterations.</description><identifier>ISSN: 0021-9991</identifier><identifier>EISSN: 1090-2716</identifier><identifier>DOI: 10.1016/j.jcp.2018.06.052</identifier><language>eng</language><publisher>Cambridge: Elsevier Inc</publisher><subject>Algorithms ; Anisotropy ; Computational physics ; DDFV ; Decoupled algorithm ; Diffusion ; Diffusion problems ; Discrete element method ; Finite volume method ; Iterative methods ; Mathematical analysis ; Mathematical problems ; Matrix methods ; Positivity-preserving scheme ; Well posed problems ; Well-posedness</subject><ispartof>Journal of computational physics, 2018-11, Vol.372, p.773-798</ispartof><rights>2018 Elsevier Inc.</rights><rights>Copyright Elsevier Science Ltd. Nov 1, 2018</rights><lds50>peer_reviewed</lds50><woscitedreferencessubscribed>false</woscitedreferencessubscribed><citedby>FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53</citedby><cites>FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53</cites></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><link.rule.ids>314,780,784,27924,27925</link.rule.ids></links><search><creatorcontrib>Su, Shuai</creatorcontrib><creatorcontrib>Dong, Qiannan</creatorcontrib><creatorcontrib>Wu, Jiming</creatorcontrib><title>A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes</title><title>Journal of computational physics</title><description>We suggest a new positivity-preserving discrete duality finite volume (DDFV) scheme for anisotropic diffusion problems on general polygonal meshes. Like existing DDFV schemes, this scheme is built on the primary and dual meshes and has finite volume (FV) equations for both vertex and cell-centered unknowns. The transpose of the coefficient matrix for the FV equations of the cell-centered unknowns is an M-matrix while that for the vertex unknowns is not an M-matrix but a symmetric and positive definite matrix. By employing a certain truncation technique for the vertex unknowns, the positivity-preserving property for both categories of unknowns is guaranteed. Local conservation is strictly maintained for the cell-centered unknowns and conditionally maintained for the vertex unknowns. Since the FV equations of the vertex unknowns can be solved independently, the two sets of FV equations are decoupled. In contrast to existing nonlinear positivity-preserving schemes, the new scheme requires no nonlinear iterations for linear problems. For nonlinear problems, the positivity-preserving mechanism of the new scheme is decoupled from its nonlinear iteration so that any nonlinear solver can be adopted. Moreover, the positivity of the discrete solution is proved and the well-posedness is analyzed rigorously for linear problems. In numerical experiments, the new scheme is examined extensively and compared with two positivity-preserving cell-centered schemes and a nonlinear DDFV scheme. Numerical results show that the present DDFV scheme achieves second-order accuracy and preserves the positivity of the solution for heterogeneous and anisotropic problems on severely distorted grids. The high efficiency of the scheme is also demonstrated by the comparison of computation time and number of nonlinear iterations.</description><subject>Algorithms</subject><subject>Anisotropy</subject><subject>Computational physics</subject><subject>DDFV</subject><subject>Decoupled algorithm</subject><subject>Diffusion</subject><subject>Diffusion problems</subject><subject>Discrete element method</subject><subject>Finite volume method</subject><subject>Iterative methods</subject><subject>Mathematical analysis</subject><subject>Mathematical problems</subject><subject>Matrix methods</subject><subject>Positivity-preserving scheme</subject><subject>Well posed problems</subject><subject>Well-posedness</subject><issn>0021-9991</issn><issn>1090-2716</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2018</creationdate><recordtype>article</recordtype><recordid>eNp9UE1LxDAQDaLguvoDvAU8t05ik7Z4EvELBC96Dm0yXVO6Tc20C_sX_NVG1rOnecO89-bxGLsUkAsQ-rrPezvlEkSVg85BySO2ElBDJkuhj9kKQIqsrmtxys6IegCoVFGt2Pcdd2jDMg3oeDM6PgXys9_5eZ9NEQnjzo8b7jzZiDNytzRDuvHOjz6tuzAsW-RkPzGNLsTk4SnMMUzeJlXXLeTDyKcY2gG3xBPe4IixGdKnYb8JY0JbpE-kc3bSNQPhxd9cs4_Hh_f75-z17enl_u41szdSzZmUTotWtYBWaYVNJWUNhS1tqaoSoHWtqIpSdlCWFptELXQLXaFqbYWqnbpZs6uDbwr1tSDNpg9LTDnISCG0lrKqIbHEgWVjIIrYmSn6bRP3RoD5rdz0JlVufis3oE2qPGluDxpM8XceoyHrcbTofEQ7Gxf8P-ofntqMkQ</recordid><startdate>20181101</startdate><enddate>20181101</enddate><creator>Su, Shuai</creator><creator>Dong, Qiannan</creator><creator>Wu, Jiming</creator><general>Elsevier Inc</general><general>Elsevier Science Ltd</general><scope>AAYXX</scope><scope>CITATION</scope><scope>7SC</scope><scope>7SP</scope><scope>7U5</scope><scope>8FD</scope><scope>JQ2</scope><scope>L7M</scope><scope>L~C</scope><scope>L~D</scope></search><sort><creationdate>20181101</creationdate><title>A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes</title><author>Su, Shuai ; Dong, Qiannan ; Wu, Jiming</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2018</creationdate><topic>Algorithms</topic><topic>Anisotropy</topic><topic>Computational physics</topic><topic>DDFV</topic><topic>Decoupled algorithm</topic><topic>Diffusion</topic><topic>Diffusion problems</topic><topic>Discrete element method</topic><topic>Finite volume method</topic><topic>Iterative methods</topic><topic>Mathematical analysis</topic><topic>Mathematical problems</topic><topic>Matrix methods</topic><topic>Positivity-preserving scheme</topic><topic>Well posed problems</topic><topic>Well-posedness</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>Su, Shuai</creatorcontrib><creatorcontrib>Dong, Qiannan</creatorcontrib><creatorcontrib>Wu, Jiming</creatorcontrib><collection>CrossRef</collection><collection>Computer and Information Systems Abstracts</collection><collection>Electronics &amp; Communications Abstracts</collection><collection>Solid State and Superconductivity Abstracts</collection><collection>Technology Research Database</collection><collection>ProQuest Computer Science Collection</collection><collection>Advanced Technologies Database with Aerospace</collection><collection>Computer and Information Systems Abstracts – Academic</collection><collection>Computer and Information Systems Abstracts Professional</collection><jtitle>Journal of computational physics</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>Su, Shuai</au><au>Dong, Qiannan</au><au>Wu, Jiming</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes</atitle><jtitle>Journal of computational physics</jtitle><date>2018-11-01</date><risdate>2018</risdate><volume>372</volume><spage>773</spage><epage>798</epage><pages>773-798</pages><issn>0021-9991</issn><eissn>1090-2716</eissn><abstract>We suggest a new positivity-preserving discrete duality finite volume (DDFV) scheme for anisotropic diffusion problems on general polygonal meshes. Like existing DDFV schemes, this scheme is built on the primary and dual meshes and has finite volume (FV) equations for both vertex and cell-centered unknowns. The transpose of the coefficient matrix for the FV equations of the cell-centered unknowns is an M-matrix while that for the vertex unknowns is not an M-matrix but a symmetric and positive definite matrix. By employing a certain truncation technique for the vertex unknowns, the positivity-preserving property for both categories of unknowns is guaranteed. Local conservation is strictly maintained for the cell-centered unknowns and conditionally maintained for the vertex unknowns. Since the FV equations of the vertex unknowns can be solved independently, the two sets of FV equations are decoupled. In contrast to existing nonlinear positivity-preserving schemes, the new scheme requires no nonlinear iterations for linear problems. For nonlinear problems, the positivity-preserving mechanism of the new scheme is decoupled from its nonlinear iteration so that any nonlinear solver can be adopted. Moreover, the positivity of the discrete solution is proved and the well-posedness is analyzed rigorously for linear problems. In numerical experiments, the new scheme is examined extensively and compared with two positivity-preserving cell-centered schemes and a nonlinear DDFV scheme. Numerical results show that the present DDFV scheme achieves second-order accuracy and preserves the positivity of the solution for heterogeneous and anisotropic problems on severely distorted grids. The high efficiency of the scheme is also demonstrated by the comparison of computation time and number of nonlinear iterations.</abstract><cop>Cambridge</cop><pub>Elsevier Inc</pub><doi>10.1016/j.jcp.2018.06.052</doi><tpages>26</tpages></addata></record>
fulltext fulltext
identifier ISSN: 0021-9991
ispartof Journal of computational physics, 2018-11, Vol.372, p.773-798
issn 0021-9991
1090-2716
language eng
recordid cdi_proquest_journals_2116622890
source ScienceDirect Freedom Collection
subjects Algorithms
Anisotropy
Computational physics
DDFV
Decoupled algorithm
Diffusion
Diffusion problems
Discrete element method
Finite volume method
Iterative methods
Mathematical analysis
Mathematical problems
Matrix methods
Positivity-preserving scheme
Well posed problems
Well-posedness
title A decoupled and positivity-preserving discrete duality finite volume scheme for anisotropic diffusion problems on general polygonal meshes
url http://sfxeu10.hosted.exlibrisgroup.com/loughborough?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2024-12-24T17%3A56%3A29IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_cross&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=A%20decoupled%20and%20positivity-preserving%20discrete%20duality%20finite%20volume%20scheme%20for%20anisotropic%20diffusion%20problems%20on%20general%20polygonal%20meshes&rft.jtitle=Journal%20of%20computational%20physics&rft.au=Su,%20Shuai&rft.date=2018-11-01&rft.volume=372&rft.spage=773&rft.epage=798&rft.pages=773-798&rft.issn=0021-9991&rft.eissn=1090-2716&rft_id=info:doi/10.1016/j.jcp.2018.06.052&rft_dat=%3Cproquest_cross%3E2116622890%3C/proquest_cross%3E%3Cgrp_id%3Ecdi_FETCH-LOGICAL-c325t-22d61b5b0ec565ea822904c7c758700bdb18472f077cea61b46b0f4596c159d53%3C/grp_id%3E%3Coa%3E%3C/oa%3E%3Curl%3E%3C/url%3E&rft_id=info:oai/&rft_pqid=2116622890&rft_id=info:pmid/&rfr_iscdi=true