Hostname: page-component-76fb5796d-22dnz Total loading time: 0 Render date: 2024-04-26T14:38:12.934Z Has data issue: false hasContentIssue false
  • English
  • Français

Monotone Finite Difference Schemes for Anisotropic Diffusion Problems via Nonnegative Directional Splittings

Part of: Partial differential equations, initial value and time-dependent initial-boundary value problems Partial differential equations, boundary value problems

Published online by Cambridge University Press:  01 February 2016

Cuong Ngo  and
Weizhang Huang
Cuong Ngo
Affiliation:
Department of Mathematics, University of Kansas, Lawrence, KS 66045, USA.
Weizhang Huang*
Affiliation:
Department of Mathematics, University of Kansas, Lawrence, KS 66045, USA.
*
*Corresponding author. Email addresses:cngo@ku.edu (C. Ngo), whuang@ku.edu (W. Huang)
Get access

Abstract

Nonnegative directional splittings of anisotropic diffusion operators in the divergence form are investigated. Conditions are established for nonnegative directional splittings to hold in a neighborhood of an arbitrary interior point. The result is used to construct monotone finite difference schemes for the boundary value problem of anisotropic diffusion operators. It is shown that such a monotone scheme can be constructed if the underlying diffusion matrix is continuous on the closure of the physical domain and symmetric and uniformly positive definite on the domain, the mesh spacing is sufficiently small, and the size of finite difference stencil is sufficiently large. An upper bound for the stencil size is obtained, which is determined completely by the diffusion matrix. Loosely speaking, the more anisotropic the diffusion matrix is, the larger stencil is required. An exception is the situation with a strictly diagonally dominant diffusion matrix where a three-by-three stencil is sufficient for the construction of a monotone finite difference scheme. Numerical examples are presented to illustrate the theoretical findings.

Keywords

Anisotropic diffusionanisotropic coefficientdiscrete maximum principlemonotone schemefinite differencenonnegative directional splitting

MSC classification

Secondary: 65N06: Finite difference methods 65M06: Finite difference methods
Type
Research Article
Information
Communications in Computational Physics , Volume 19 , Issue 2 , February 2016 , pp. 473 - 495
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Barles, G. and Souganidis, P.. Convergence of approximation schemes for fully nonlinear second order equations. Asymptotic Analysis, 4(3):271–283,1991.Google Scholar
[2]Bramble, J. H., Hubbard, B., and Thomée, V.. Convergence estimates for essentially positive type discrete dirichlet problems. Math. Comp., 23:695–709,1969.CrossRefGoogle Scholar
[3]Ciarlet, P.. Discrete maximum principle for finite difference operators. Aequationes Math., 4:338–352,1970.CrossRefGoogle Scholar
[4]Ciarlet, P. and Raviart, P.-A.. Maximum principle and uniform convergence for the finite element method. Comput. Meth. Appl. Mech. Eng., 2:17–31,1973.Google Scholar
[5]Crumpton, P. I., Shaw, G.J., and Ware, A. R. Discretisation and multigrid solution of elliptic equations with mixed derivative terms and strongly discontinuous coefficients. J. Comput.Phys., 116:343–358,1995.Google Scholar
[6]Drǎgǎnescu, A., Dupont, T. R, and Scott, L. R.. Failure of the discrete maximum principle for an elliptic finite element problem. Math. Comp., 74:1–23,2004.Google Scholar
[7]Ertekin, T., Abou-Kassem, J. H., and King, G. R.. Basic Applied Reservoir Simulation. SPE textbook series, Vol. 7, Richardson, Texas, 2001.CrossRefGoogle Scholar
[8]Greenspan, D. and Jain, P. C.. On non-negative difference analogues of elliptic differential equations. J. Franklin Inst., 279:360–365,1965.Google Scholar
[9]Gűnter, S. and Lackner, K.. A mixed implicit-explicit finite difference scheme for heat transport in magnetised plasmas. J. Comput. Phys., 228:282–293,2009.Google Scholar
[10]Günter, S., Lackner, K., and Tichmann, C.. Finite element and higher order difference formulations for modelling heat transport in magnetised plasmas. J. Comput. Phys., 226:2306–2316, 2007.CrossRefGoogle Scholar
[11]Günter, S., Yu, Q., Kruger, J., and Lackner, K.. Modelling of heat transport in magnetised plasmas using non-aligned coodinates. J. Comput. Phys., 209:354–370,2005.Google Scholar
[12]Karras, D. A. and Mertzios, G. B.. New PDE-based methods for image enhancement using SOM and Bayesian inference in various discretization schemes. Meas. Sci. Technol., 20:104012, 2009.CrossRefGoogle Scholar
[13]Kuzmin, D.,Shashkov, J. J., and Svyatskiy, D.. A constrained finite element method satisfying the discrete maximum principle for anisotropic diffusion problems. J. Comput. Phys., 228:3448–3463,2009.Google Scholar
[14]Li, X. P. and Huang, W.. An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems. J. Comput. Phys., 229:8072–8094,2010.(arXiv:1003.4530).Google Scholar
[15]Li, X. P. and Huang, W.. Maximum principle for the finite element solution of time dependent anisotropic diffusion problems. Numer. Meth. P. D. E., 29:1963–1985,2013. (arXiv:1209.5657).Google Scholar
[16]Li, X. P., Svyatskiy, D., and Shashkov, M.. Mesh adaptation and discrete maximum principle for 2d anisotropic diffusion problems. Technical report, LANL, 2007. Final Report of the Summer Student Program.Google Scholar
[17]Lipnikov, K., Shashkov, M., Svyatskiy, D., and Vassilevski, Y.. Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes. J. Comput. Phys., 227:492–512,2007.CrossRefGoogle Scholar
[18]Liska, R. and Shashkov, M.. Enforcing the discrete maximum principle for linear finite element solutions of second-order elliptic problems. Comm. Comput. Phys., 3:852–877,2008.Google Scholar
[19]Mlacnik, M.J. and Durlofsky, L.J.. Unstructured grid optimizatioin for improved monotonicity of discrete solutions of elliptic equations with highly anisotropic cofficients. J. Comput. Phys., 216:337–361,2006.CrossRefGoogle Scholar
[20]Mlacnik, M.J. and Durlofsky, L.J.. Unstructured grid optimization for improved monotonicity of discrete solutions of elliptic equations with highly anisotropic coefficients. J. Comput. Phys., 216:337–361,2006.Google Scholar
[21]Motzkin, T. S. and Wasow, W.. On the approximation of linear elliptic differential equations by difference equations with positive coefficients. J. Math. Physics, 31:253–259,1953.CrossRefGoogle Scholar
[22]Mrázek, P. and Navara, M.. Consistent positive directional splitting of anisotropic diffusion. In Likar, B., editor, Proc. of Computer Vision Winter Workshop, pages 37–48, Bled, Slovenia, February 2001.Google Scholar
[23]Oberman, A.. Convergent difference schemes for degenerate elliptic and parabolic equations: Hamilton-jacobi equations and free boundary problems. SIAM J. Numer. Anal, 44(2):879–895, 2006.Google Scholar
[24]Potier, C. L.. Schéma volumes finis monotone pour des opérateurs de diffusion fortement anisotropes sur des maillages de triangles non structurés. C. R. Math. Acad. Sci. Paris, 341:787–792,2005.Google Scholar
[25]Potier, C. L.. A nonlinear finite volume scheme satisfying maximum and minimum principles for diffusion operators. Int. J. Finite Vol., 6:20,2009.Google Scholar
[26]Potier, C. L.. Un schéma linéaire vérifiant le principe du maximum pour des opérateurs de diffusion très anisotropes sur des maillages déformés. C. R. Math. Acad. Sci. Paris, 347:105–110,2009.Google Scholar
[27]Sharma, P. and Hammett, G. W.. Preserving monotonicity in anisotropic diffusion. J. Comput. Phys., 227:123–142,2007.Google Scholar
[28]Stoyan, G.. On a maximum principle for matrices, and on conservation of monotonicity, with applications to discretization methods. ZAMM, 62:375–381,1982.Google Scholar
[29]Stoyan, G.. On maximum principles for monotone matrices. Lin. Alg. Appl., 78:147–161,1986.Google Scholar
[30]Varga, R.. Matrix Iterative Analysis. Prentice-Hall, New Jersey, 1962.Google Scholar
[31]Weickert, J.. Anisotropic Diffusion in Image Processing. B.G Teubner Stuttgart, 1998.Google Scholar
[32]Yuan, G. and Sheng, Z.. Monotone finite volume schemes for diffusion equations on polygonal meshes. J. Comput. Phys., 227:6288–6312,2008.CrossRefGoogle Scholar
[33]Zhang, Y., Zhang, X., and Shu, C.-W.. Maximum-principle-satisfying second order discontinuous Galerkin schemes for convection-diffusion equations on triangular meshes. J. Comput. Phys., 234:295–316,2013.Google Scholar