Schemes of high-order accuracy for the multi-dimensional heat conduction equation

A.A. Samarskii

DOI: https://doi.org/10.1016/0041-5553(63)90104-6

1963-01-01

USSR Computational Mathematics and Mathematical Physics

Abstract:In [1] an economical scheme is put forward for the heat conduction equation ∂u∂t = ∑α=1p ∂2u∂xα2, p⩽4. with accuracy O(h4 + τ2), where h is the step of the space net, and τ the time step. It is a three-layer scheme.In this paper (Section 2–4) we examine two-layer schemes of O(h4 + τ2) (i.e. of order (4,2)) which are suitable for p ⩽ 3.∗∗ They are put into effect with the aid of a number of splitting algorithms or alternating direction algorithms which involve practically the same number of operations as the corresponding algorithms of O(f2 + τ2) (see [2]–[6]). It is shown that these schemes are absolutely stable and that for any values of γ = τh2 they converge in the mean at a rate O(h4 + τ2). In Section 5 we consider a three-layer scheme of order (4,2) which is more economical than the scheme of [1] and is suitable for p ⩽ 4. This scheme is absolutely stable and converges for any γ (the scheme of [1] converges for γ ⩾ γ0 = const. > 0).A scheme of order (4, 2) for the equation ∂u∂t = ∑α,beta;=12 aαβ∂2u∂xα∂xβ, ∑α,beta;=12 aαβξαξβ ⩾ C1 ∑α=12 ξα2, aαβ = const. is given in Section 6.In Section 7 we examine two-layer schemes of the same order of accuracy for the equation with variable coefficients, c(x,t) ∂u∂t = ∑α=1p ∂∂xαkα(x,t) ∂∂xα + ƒ(x,t), x = (x1…xp), p⩾2.

English Else