Premium
Extension theorems for Stokes and Lamé equations for nearly incompressible media and their applications to numerical solution of problems with highly discontinuous coefficients
Author(s) -
Bakhvalov N. S.,
Knyazev A. V.,
Parashkevov R. R.
Publication year - 2002
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.259
Subject(s) - preconditioner , mathematics , mathematical analysis , compressibility , finite element method , domain (mathematical analysis) , linear system , physics , aerospace engineering , thermodynamics , engineering
We prove extension theorems in the norms described by Stokes and Lamé operators for the three‐dimensional case with periodic boundary conditions. For the Lamé equations, we show that the extension theorem holds for nearly incompressible media, but may fail in the opposite limit, i.e. for case of absolutely compressible media. We study carefully the latter case and associate it with the Cosserat problem. Extension theorems serve as an important tool in many applications, e.g. in domain decomposition and fictitious domain methods, and in analysis of finite element methods. We consider an application of established extension theorems to an efficient iterative solution technique for the isotropic linear elasticity equations for nearly incompressible media and for the Stokes equations with highly discontinuous coefficients. The iterative method involves a special choice for an initial guess and a preconditioner based on solving a constant coefficient problem. Such preconditioner allows the use of well‐known fast algorithms for preconditioning. Under some natural assumptions on smoothness and topological properties of subdomains with small coefficients, we prove convergence of the simplest Richardson method uniform in the jump of coefficients. For the Lamé equations, the convergence is also uniform in the incompressible limit. Our preliminary numerical results for two‐dimensional diffusion problems show fast convergence uniform in the jump and in the mesh size parameter. Copyright © 2002 John Wiley & Sons, Ltd.