科研之友
微信
新浪微博
Facebook
分享链接
摘要
We suggest a new positivity-preserving finite volume scheme for anisotropic diffusion problems on arbitrary polygonal grids. The scheme has vertex-centered, edge-midpoint and cell-centered unknowns. The vertex-centered unknowns are primary and have finite volume equations associated with them. The edge-midpoint and cell-centered unknowns are treated as auxiliary ones and are interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. Unlike most existing positivity-preserving schemes, the construction of the scheme is based on a special nonlinear two-point flux approximation that has a fixed stencil and does not require the convex decomposition of the co-normal. In order to solve efficiently the nonlinear systems resulting from the nonlinear scheme, Picard method and its Anderson acceleration are discussed. Numerical experiments demonstrate the second-order accuracy and well positivity of the solution for heterogeneous and anisotropic problems on severely distorted grids. The high efficiency of the Anderson acceleration is also shown on reduction of the number of nonlinear iterations. Moreover, the proposed scheme does not have the so-called numerical heat-barrier issue suffered by most existing cell-centered and hybrid schemes. However, further improvements have to be made if the solution is very close to the machine precision and the mesh distortion is very severe.
- 出版日期2017-9-1
- 单位武汉大学; 中国工程物理研究院; 北京应用物理与计算数学研究所
