Upscaling for the Laplace problem using a discontinuous Galerkin method

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:240

Article

Barucq, H.²  Chaumont-Frelet, T.²  Diaz, T. J.²  Peron, V.^1,2 

[1] Univ Pau & Pays Adour, F-64010 Pau, France

[2] INRIA Bordeaux Sud Ouest, MAGIQUE 3D, Bordeaux, France

关键词: Upscaling;    Laplace problem;    Multiscale methods;    Discontinuous Galerkin;    Interior penalty;   

DOI  :  10.1016/j.cam.2012.05.025

来源: Elsevier

【 摘 要 】

Scientists and engineers generally tackle problems that include multiscale effects and that are thus difficult to solve numerically. The main difficulty is to capture both the fine and the coarse scales to get an accurate numerical solution. Indeed, the computations are generally performed by using numerical schemes based on grids. But the stability and thus the accuracy of the numerical method depends on the size of the grid which must be refined drastically in the case of very fine scales. That implies huge computational costs and in particular the limitations of the memory capacity are often reached. It is thus necessary to use numerical methods that are able to capture the fine scale effects with computations on coarse meshes. Operator-based upscaling is one of them and we present a first attempt to adapt that technique to a Discontinuous Galerkin Method (DGM). We consider the Laplace problem as a benchmark and we compare the performance of the resulting numerical scheme with the classical one using Lagrange finite elements. The comparison involves both an accuracy analysis and a complexity calculus. This work shows that there is an interest of combining DGM with upscaling. (c) 2012 Elsevier B.V. All rights reserved.

【 授权许可】

Free