
A MIXED VIRTUAL ELEMENT METHOD FOR THE BOUSSINESQ PROBLEM ON POLYGONAL MESHES
Gabriel N. Gatica, Mauricio Munar, Filánder A. Sequeira
Journal of Computational Mathematics
2021, 39 (3):
392427.
DOI: 10.4208/jcm.2001m20190187
In this work we introduce and analyze a mixed virtual element method (mixedVEM) for the twodimensional stationary Boussinesq problem. The continuous formulation is based on the introduction of a pseudostress tensor depending nonlinearly on the velocity, which allows to obtain an equivalent model in which the main unknowns are given by the aforementioned pseudostress tensor, the velocity and the temperature, whereas the pressure is computed via a postprocessing formula. In addition, an augmented approach together with a fixed point strategy is used to analyze the wellposedness of the resulting continuous formulation. Regarding the discrete problem, we follow the approach employed in a previous work dealing with the NavierStokes equations, and couple it with a VEM for the convectiondiffusion equation modelling the temperature. More precisely, we use a mixedVEM for the scheme associated with the fluid equations in such a way that the pseudostress and the velocity are approximated on virtual element subspaces of ${\Bbb H}$(div) and
H
^{1}, respectively, whereas a VEM is proposed to approximate the temperature on a virtual element subspace of H
^{1}. In this way, we make use of the L
^{2}orthogonal projectors onto suitable polynomial spaces, which allows the explicit integration of the terms that appear in the bilinear and trilinear forms involved in the scheme for the fluid equations. On the other hand, in order to manipulate the bilinear form associated to the heat equations, we define a suitable projector onto a space of polynomials to deal with the fact that the diffusion tensor, which represents the thermal conductivity, is variable. Next, the corresponding solvability analysis is performed using again appropriate fixedpoint arguments. Further, Strangtype estimates are applied to derive the a
priori error estimates for the components of the virtual element solution as well as for the fully computable projections of them and the postprocessed pressure. The corresponding rates of convergence are also established. Finally, several numerical examples illustrating the performance of the mixedVEM scheme and confirming these theoretical rates are presented.
Reference 
Related Articles 
Metrics

