A mixed discontinuous Galerkin method with symmetric stress for Brinkman problem based on the velocity–pseudostress formulation
Yanxia Qian,Shuonan Wu,Fei Wang
DOI: https://doi.org/10.1016/j.cma.2020.113177
IF: 6.588
2020-08-01
Computer Methods in Applied Mechanics and Engineering
Abstract:<p>The Brinkman equations can be regarded as a combination of the Stokes and Darcy equations which model transitions between the fast flow in channels (governed by Stokes equations) and the slow flow in porous media (governed by Darcy's law). The numerical challenge for this model is the designing of a numerical scheme which is stable for both the Stokes-dominated (high permeability) and the Darcy-dominated (low permeability) equations. In this paper, we solve the Brinkman model in <span class="math"><math>n</math></span> dimensions (<span class="math"><math>n=2,3</math></span>) by using the mixed discontinuous Galerkin (MDG) method, which meets this challenge. This MDG method is based on the pseudostress–velocity formulation and uses a discontinuous piecewise polynomial pair <span class="math"><math>P̲k+1S</math></span>-<span class="math"><math>Pk</math></span><span class="math"><math>(k≥0)</math></span>, where the stress field is symmetric. The main unknowns are the pseudostress and the velocity, whereas the pressure is easily recovered through a simple postprocessing. A key step in the analysis is to establish the parameter-robust inf–sup stability through specific parameter-dependent norms at both continuous and discrete levels. Therefore, the stability results presented here are uniform with respect to the permeability. Thanks to the parameter-robust stability analysis, we obtain optimal error estimates for the stress in broken <span class="math"><math>H̲(div)</math></span>-norm and velocity in <span class="math"><math>L2</math></span>-norm. Furthermore, the <span class="math"><math>L̲2</math></span> error estimate for pseudostress is derived under certain conditions. Finally, numerical experiments are provided to support the theoretical results and to show the robustness, accuracy, and flexibility of the MDG method.</p>
mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications