首页|Positive Scharfetter-Gummel finite volume method for convection-diffusion equations on polygonal meshes
Positive Scharfetter-Gummel finite volume method for convection-diffusion equations on polygonal meshes
扫码查看
点击上方二维码区域，可以放大扫码查看
原文链接
NSTL
Elsevier
In this paper, we develop and study a fully implicit positive finite volume scheme that allows an accurate approximation of the nonlinear highly anisotropic convection-diffusion equations on almost arbitrary girds. The key idea is to relate the oscillatory fluxes, including the convective ones, to the normal monotonic diffusive flux thanks to a technique used in the Scharfetter-Gummel discretizations. Then, we obtain a nonlinear two-point-like scheme with positive coefficients on primal and dual meshes. We check that the structure of the scheme naturally ensures the nonnegativity of the approximate solutions. We also establish energy estimates, which leads to a proof of existence of the numerical solutions. This analytical study is accompanied with a series of numerical results and simulations. They highlight the fulfillment of the discrete maximum principle, the optimal accuracy of our scheme, and its robustness with respect to the mesh and to high ratios of anisotropy. (C) 2022 Elsevier Inc. All rights reserved.
Degenerate parabolic equationsFinite volume methodScharfetter-Gummel schemeGeneric meshesPositivityCoercivitySecond order accuracyNUMERICAL-ANALYSISSCHEMEDISCRETIZATIONAPPROXIMATIONCONVERGENCE2ND-ORDEROPERATORS
NSTL
Elsevier
