The interior penalty virtual element method for the biharmonic problem
Jikun Zhao, Shipeng Mao, Bei Zhang, Fei Wang
Abstract
In this paper, an interior penalty virtual element method (IPVEM) is developed for solving the biharmonic problem on polygonal meshes. By modifying the existing <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper H squared"> <mml:semantics> <mml:msup> <mml:mi>H</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">H^2</mml:annotation> </mml:semantics> </mml:math> </inline-formula> -conforming virtual element, an <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper H Superscript 1"> <mml:semantics> <mml:msup> <mml:mi>H</mml:mi> <mml:mn>1</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">H^1</mml:annotation> </mml:semantics> </mml:math> </inline-formula> -nonconforming virtual element is obtained with the same degrees of freedom as the usual <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper H Superscript 1"> <mml:semantics> <mml:msup> <mml:mi>H</mml:mi> <mml:mn>1</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">H^1</mml:annotation> </mml:semantics> </mml:math> </inline-formula> -conforming virtual element, such that it locally has <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper H squared"> <mml:semantics> <mml:msup> <mml:mi>H</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">H^2</mml:annotation> </mml:semantics> </mml:math> </inline-formula> -regularity on each polygon in meshes. To enforce the <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper C Superscript 1"> <mml:semantics> <mml:msup> <mml:mi>C</mml:mi> <mml:mn>1</mml:mn> </mml:msup> <mml:annotation encoding="application/x-tex">C^1</mml:annotation> </mml:semantics> </mml:math> </inline-formula> continuity of the solution, an interior penalty formulation is adopted. Hence, this new numerical scheme can be regarded as a combination of the virtual element space and discontinuous Galerkin scheme. Compared with the existing methods, this approach has some advantages in reducing the degree of freedom and capability of handling hanging nodes. The well-posedness and optimal convergence of the IPVEM are proven in a mesh-dependent norm. We also derive a sharp estimate of the condition number of the linear system associated with IPVEM. Some numerical results are presented to verify the optimal convergence of the IPVEM and the sharp estimate of the condition number of the discrete problem. Besides, in the numerical test, the IPVEM has a good performance in computational accuracy by contrast with the other VEMs solving the biharmonic problem.