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
    0 references
    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
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    0 references
    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
    0 references