Premium
Numerical treatment of degenerate diffusion equations via Feller's boundary classification, and applications
Author(s) -
Cacio Emanuela,
Cohn Stephen E.,
Spigler Renato
Publication year - 2012
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.20657
Subject(s) - mathematics , degenerate energy levels , diffusion , boundary (topology) , mathematical analysis , physics , quantum mechanics , thermodynamics
A numerical method is devised to solve a class of linear boundary‐value problems for one‐dimensional parabolic equations degenerate at the boundaries. Feller theory, which classifies the nature of the boundary points, is used to decide whether boundary conditions are needed to ensure uniqueness, and, if so, which ones they are. The algorithm is based on a suitable preconditioned implicit finite‐difference scheme, grid, and treatment of the boundary data. Second‐order accuracy, unconditional stability, and unconditional convergence of solutions of the finite‐difference scheme to a constant as the time‐step index tends to infinity are further properties of the method. Several examples, pertaining to financial mathematics, physics, and genetics, are presented for the purpose of illustration. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011