On a Finite Difference Analogue of an Elliptic Boundary Problem which is Neither Diagonally Dominant Nor of Non‐negative Type
On a Finite Difference Analogue of an Elliptic Boundary Problem which is Neither Diagonally Dominant Nor of Non‐negative Type
复制标题
非对角占优非非负型椭圆边界问题的有限差分模拟
DOI:
10.1002/sapm1964431117
复制
发表时间:
1964
期刊:
影响因子:
--
通讯作者:
B. Hubbard
中科院分区:
文献类型:
--
作者:
J. Bramble;B. Hubbard
1. Introduction. In the usual study of the discretIzation error resulting from a, pproximatmg boundary problems for ellIptic equations by finite difference methods the maximum prmCIple plays a central role. In 1930 S. Gerschgorin (14) gave a method for estimatmg the order of convergence of the solutIOn to a certam class of fimte difference analogues to the solutIOn of the DirIchlet problem for ellIptic equatIOns. The matrIX of the resultmg system of SImultaneous lmear equatIOns possesses the property of diagonal dominance, ie the sum of the absolute values of the off-dIagonal elements in each row does not exceed the magmtude of the dIagonal element. Furthermore It satIsfies the condItIon that the dIagonal elements are all posItive and the off-dIagonal elements are non-POSItIve (2.10). Following this, others (3),(8),(9),(12),(20),(21) have extended the results of Gerschgorm within the framework of these condItIOns. Recently the authors (5),(6) gave a theorem on the formulatIon of fimte dIfference analogues of the Dmchlet problem for elliptIC equatIOns III whICh the propertIes (2.10) were relaxed near the boundary. As examples of the theorem certain hIgher order fimte dIfference analogues were dIscussed and shown to have convergence properties prevIously unexpected. The question of a maXImum prmciple for the entIre problem was CIrcumvented however and has been dealt WIth only recently by M Rockoff (18). Hence It IS clear that the suffiCIent condI-tIons (2 10) are not necessaryOne purpose of this paper is to give a specific example in which conditions (2.10) are violated at every interIOr point but for which a maximum principle still holds. This is done in section 2 where an O (h') approximation to the two point boundary problem (2.1) is studIed (h is the mesh constant). The interesting fact that the approximation to the operator near the boundary need be only O (h2) without destroymg the overall accuracy of the problem is shown to be true.