Convergence of the CEM-GMsFEM for Stokes flows in heterogeneous perforated domains

From MaRDI portal
Publication:2226262



Abstract: In this paper, we consider the incompressible Stokes flow problem in a perforated domain and employ the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) to solve this problem. The proposed method provides a flexible and systematical approach to construct crucial divergence-free multiscale basis functions for approximating the displacement field. These basis functions are constructed by solving a class of local energy minimization problems over the eigenspaces that contain local information on the heterogeneities. These multiscale basis functions are shown to have the property of exponential decay outside the corresponding local oversampling regions. By adapting the technique of oversampling, the spectral convergence of the method with error bounds related to the coarse mesh size is proved.


Let \(\Omega\subset\mathbb R^d\) be a bounded domain and \(B_{\epsilon}\) be a set of perforations within this domain. The set of perforations \(B_{\epsilon}\) is assumed to be a union of connected circular disks. Each of these disks is of diameter of order \(0<\epsilon\ll \operatorname{diam}(\Omega)\) and might be separated by a distance also of order \(\epsilon\). The work develops the so-called CEM-GMs FEM-method (constraint energy minimizing generalized multiscale finite element method) for studying the following Stokes problem in the theory of incompressible fluid mechanics, namely \begin{align*} -\Delta\mathbf{u}+\nabla p&=\mathbf{f} \text{ in } \Omega^{\epsilon},\\ \nabla\cdot\mathbf{u}&=0 \text{ in } \Omega^{\epsilon},\\ \mathbf{u}&=\mathbf{g} \text{ on } \partial\Omega^{\epsilon}\cap\partial\Omega, \\ \mathbf{u}&=\mathbf{0} \text{ on } \partial\Omega^{\epsilon}\cap B_{\epsilon}. \end{align*} The authors reformulate the problem as a variational one using classical Lebesgue spaces \(L^2_0(\Omega^{\epsilon})\) with zero mean, Sobolev spaces \(H^1_0(\Omega^{\epsilon})\) with vanishing trace and their vector variant \(\mathbf{V}_0=(H^1_0(\Omega^{\epsilon}))^d\). Find \((\mathbf{u},p)\in\mathbf{V}_0\times L^2_0(\Omega^{\epsilon})\) such that \begin{align*} a(\mathbf{u},\mathbf{v})-b(\mathbf{v},p)&=(\mathbf{f},\mathbf{v}) \text{ for all } \mathbf{v}\in\mathbf{V}_0,\\ b(\mathbf{u},q)&=0 \text{ for all } q\in L^2_0(\Omega^{\epsilon}), \end{align*} where \[ a(\mathbf{u},\mathbf{v})=\int\limits_{\Omega^{\epsilon}}\nabla\mathbf{u}\cdot\nabla\mathbf{v}dx,\quad b(\mathbf{u},q)=\int\limits_{\Omega^{\epsilon}}q\nabla\cdot\mathbf{u}dx. \] After this reformulation, they introduce auxiliary spaces and the construction of multiscale basis functions for pressure using relaxed constraint energy minimization. The multiscale basis functions are constructed by solving a class of local spectral problems and constrained minimization problems. Convergence results are given in a separate section.



Cites work









This page was built for publication: Convergence of the CEM-GMsFEM for Stokes flows in heterogeneous perforated domains

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2226262)