Stability, convergence and order of the extrapolations of the residual smoothing scheme in energy norm (Q643860)
From MaRDI portal
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Stability, convergence and order of the extrapolations of the residual smoothing scheme in energy norm |
scientific article |
Statements
Stability, convergence and order of the extrapolations of the residual smoothing scheme in energy norm (English)
0 references
2 November 2011
0 references
The authors consider an abstract differential equation in a Hilbert space having in view application of the method of lines to parabolic equations discretized by spectral or finite element methods. The residual smoothing scheme consists in preconditioning the discrete time derivative of the explicit Euler method by a suitable operator. This, like the elliptic operator, is assumed to be self-adjoint and non-negative, further, it contains a free (stability) parameter. After introducing order relations for such operators and proving that the one-step operator has energy norm not exceeding 1, they investigate properties of Richardson extrapolation. The resulting operator is shown to be stable in the energy norm for sufficiently large values of the stability parameter and to assure convergence with any fixed order \(k\) in time. Comparing with results of the alternating directions implicit (ADI) time method we can say that about fifty years ago cheaper ways of preconditioning parabolic equations were found and shown to converge, see Yanenko (1959), Bramble, D'yakonov and Samarskij -- where the Russian papers all were translated to English [e.g., \textit{N. N. Yanenko}, The method of fractional steps. The solution of problems of mathematical physics in several variables. Berlin-Heidelberg-New York: Springer Verlag (1971; Zbl 0209.47103)]. True, the operators were imagined to arise from finite difference methods, and higher order was reached for \(d\)-dimensional parallelepipedons only. However, the stability theory of Samarskij was assuming essentially just the same properties like the present paper: self-adjointness, spectral equivalence of elliptic operator and preconditioner, and resulted convergence in energy norm (a German translation appeared in 1984 of his work from the seventies, see \textit{A. A. Samarskij} [Theorie der Differenzenverfahren. Mathematik und ihre Anwendungen in Physik und Technik, Bd. 40. Leipzig: Akademische Verlagsgesellschaft Geest \& Portig K.-G. (1984; Zbl 0543.65067)]).
0 references
residual smoothing scheme
0 references
parabolic equations
0 references
convergence
0 references
energy norm
0 references
spectral method
0 references
abstract differential equation
0 references
Hilbert space
0 references
method of lines
0 references
finite element methods
0 references
preconditioning
0 references
Euler method
0 references
Richardson extrapolation
0 references
stability
0 references
finite difference methods
0 references