NAMay 28, 2019
Finite Element Approximation of Elliptic Homogenization Problems in Nondivergence-FormYves Capdeboscq, Timo Sprekeler, Endre Süli
We use uniform $W^{2,p}$ estimates to obtain corrector results for periodic homogenization problems of the form $A(x/\varepsilon):D^2 u_{\varepsilon} = f$ subject to a homogeneous Dirichlet boundary condition. We propose and rigorously analyze a numerical scheme based on finite element approximations for such nondivergence-form homogenization problems. The second part of the paper focuses on the approximation of the corrector and numerical homogenization for the case of nonuniformly oscillating coefficients. Numerical experiments demonstrate the performance of the scheme.
40.2NAApr 16
A post-processed higher-order multiscale method for nondivergence-form elliptic equationsMoritz Hauck, Roland Maier, Timo Sprekeler
We study the finite element approximation of linear second-order elliptic partial differential equations in nondivergence form with highly heterogeneous diffusion and drift coefficients. A generalized Cordes condition is imposed to guarantee that a suitably renormalized version of the nondivergence-form differential operator is near the Laplacian. Based on a stabilized symmetric formulation for the gradient that enables the use of $H^1$-conforming approximation spaces, we construct a multiscale method following the methodology of the localized orthogonal decomposition with coarse basis functions tailored to the heterogeneous coefficients. We employ a novel post-processing strategy to obtain higher-order convergence rates, overcoming previous limitations imposed by the low regularity of the load functional. Numerical experiments demonstrate the performance of the method.
15.3NAApr 12
Stable localized orthogonal decomposition in Raviart-Thomas spacesPatrick Henning, Hao Li, Timo Sprekeler
This work proposes a computational multiscale method for the mixed formulation of a second-order linear elliptic equation subject to a homogeneous Neumann boundary condition, based on a stable localized orthogonal decomposition (LOD) in Raviart-Thomas finite element spaces. In the spirit of numerical homogenization, the construction provides low-dimensional coarse approximation spaces that incorporate fine-scale information from the heterogeneous coefficients by solving local patch problems on a fine mesh. The resulting numerical scheme is accompanied by a rigorous error analysis, and it is applicable beyond periodicity and scale-separation in spatial dimensions two and three. In particular, this novel realization circumvents the presence of pollution terms observed in a previous LOD construction for elliptic problems in mixed formulation. Finally, various numerical experiments are provided that demonstrate the performance of the method.