Achim Nonnenmacher | EPFL Graph Search

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

A posteriori error estimate in quantities of interest for the finite element heterogeneous multiscale method

Assyr Abdulle, Achim Nonnenmacher

We present an "a posteriori" error analysis in quantities of interest for elliptic homogenization problems discretized by the finite element heterogeneous multiscale method. The multiscale method is based on a macro-to-micro formulation, where the macroscopic physical problem is discretized in a macroscopic finite element space, and the missing macroscopic data are recovered on-the-fly using the solutions of corresponding microscopic problems. We propose a new framework that allows to follow the concept of the (single-scale) dual-weighted residual method at the macroscopic level in order to derive a posteriori error estimates in quantities of interests for multiscale problems. Local error indicators, derived in the macroscopic domain, can be used for adaptive goal-oriented mesh refinement. These error indicators rely only on available macroscopic and microscopic solutions. We further provide a detailed analysis of the data approximation error, including the quadrature errors. Numerical experiments confirm the efficiency of the adaptive method and the effectivity of our error estimates in the quantities of interest. (C) 2013 Wiley Periodicals, Inc.

Wiley-Blackwell2013

Adaptive Finite Element Methods for Multiscale Partial Differential Equations

Achim Nonnenmacher

Engineers rely on efficient simulations that provide them with reliable data in order to make proper engineering design decisions. The purpose of this thesis is to design adaptive numerical methods for multiscale problems in this spirit. We consider elliptic homogenization problems discretized by the finite element heterogeneous multiscale method (FE-HMM). Unlike standard (single-scale) finite element methods, our multiscale discretization scheme relies on coupled macro and micro finite elements. The framework of the HMM allows to design an algorithm that follows the classical finite element structure on the macro level. The fine scales of the multiscale problems are taken into account by replacing the element-wise numerical integration over unknown macroscopic data by a numerical integration over suitably averaged micro solutions. These micro solutions are obtained from micro FE problems on sampling domains within the macro elements. This thesis is divided into two parts. In the first part, we discuss a short and versatile FE implementation of the multiscale algorithm. The implementation is flexible, easy to use and to modify and can handle simplicial or quadrilateral FE and various macro-micro coupling conditions for the constrained micro problems. The implementation of time-dependent problems is also discussed. Numerical examples including three dimensional problems are presented and demonstrate the efficiency and the versatility of the computational strategy. In the second part (the main part of this thesis), we present an a posteriori error analysis for the FE-HMM. The a posteriori analysis enables us to estimate the accuracy of a numerical solution (and therefore its reliability) and further it allows for the design of adaptive numerical methods, which are the most efficient. The crucial component for the design of an adaptive multiscale method is the introduction of appropriate error indicators. As the error indicators depend on macroscopic data (such as the macroscopic diffusion tensor) that are not readily available, we construct error indicators that only depend on the available macro and micro FE solutions, available from previous computations. We provide a posteriori estimates for the upper and lower bound in the energy norm. The corresponding macroscopic mesh refinement strategy is therefore both reliable and efficient. The microscopic mesh is refined simultaneously and – under appropriate assumptions – optimally with the macroscopic mesh. This means that the strategy reduces the macro and micro error at the same rate. In the case of a uniformly oscillating tensor and exact micro computations, the standard a posteriori error estimates for the FEM applied to the homogenized problem are recovered. Numerical experiments confirm the efficiency and reliability of the adaptive multiscale method and demonstrate the optimality of the chosen macro-micro coupling. We extend the adaptive FE-HMM to higher order FE. We further derive a posteriori estimates for the error in quantities of interest that are needed to make certain design decisions; the quantity of interest is represented by a linear functional. We derive and analyze a multiscale counterpart to the classical dual-weighted residual method and design a corresponding goal-oriented adaptive multiscale method. The efficiency of the method is shown in numerical experiments.

EPFL2011

Adaptive Finite Element Methods for Multiscale Partial Differential Equations

Achim Nonnenmacher

EPFL2011