L. Beirão da Veiga

L. Beirão da Veiga, C. Lovadina, D. Mora

We present a Virtual Element Method (VEM) for possibly nonlinear elastic and inelastic problems, mainly focusing on a small deformation regime. The numerical scheme is based on a low-order approximation of the displacement field, as well as a suitable treatment of the displacement gradient. The proposed method allows for general polygonal and polyhedral meshes, it is efficient in terms of number of applications of the constitutive law, and it can make use of any standard black-box constitutive law algorithm. Some theoretical results have been developed for the elastic case. Several numerical results within the 2D setting are presented, and a brief discussion on the extension to large deformation problems is included.

L. Beirão da Veiga, F. Dassi, A. Russo et al.

The Virtual Element Method (VEM) is a well-established framework for solving partial differential equations on polygonal and polyhedral meshes. In this paper, we introduce a novel hybrid VEM that integrates both conforming and nonconforming virtual spaces. We apply this formulation to a three-dimensional linear elasticity problem, providing rigorous theoretical analysis to demonstrate optimal convergence rates. Furthermore, we explore the extension of this approach to domains with curved boundaries.

L. Beirão da Veiga, C. Lovadina, G. Vacca

A family of Virtual Element Methods for the 2D Navier-Stokes equations is proposed and analysed. The schemes provide a discrete velocity field which is point-wise divergence-free. A rigorous error analysis is developed, showing that the methods are stable and optimally convergent. Several numerical tests are presented, confirming the theoretical predictions. A comparison with some mixed finite elements is also performed.

L. Beirão da Veiga, A. Buffa, C. Lovadina et al.

We present a new isogeometric method for the discretization of the Reissner-Mindlin plate bending problem. The proposed scheme follows a recent theoretical framework that makes possible to construct a space of smooth discrete deflections $W_h$ and a space of smooth discrete rotations $\Rots_h$ such that the Kirchhoff contstraint is exactly satisfied at the limit. Therefore we obtain a formulation which is natural from the theoretical/mechanical viewpoint and locking free by construction.

L. Beirão da Veiga, F. Brezzi, F. Dassi et al.

We consider, as a simple model problem, the application of Virtual Element Methods (VEM) to the linear Magnetostatic three-dimensional problem in the formulation of F. Kikuchi. In doing so, we also introduce new serendipity VEM spaces, where the serendipity reduction is made only on the faces of a general polyhedral decomposition (assuming that internal degrees of freedom could be more easily eliminated by static condensation). These new spaces are meant, more generally, for the combined approximation of $H^1$-conforming ($0$-forms), $H({\rm {\bf curl}})$-conforming ($1$-forms), and $H({\rm div})$-conforming ($2$-forms) functional spaces in three dimensions, and they would surely be useful for other problems and in more general contexts.

L. Beirão da Veiga, D. Mora, G. Vacca

In the present paper, we investigate the underlying Stokes complex structure of the Virtual Element Method for Stokes and Navier--Stokes introduced in previous papers by the same authors, restricting our attention to the two dimensional case. We introduce a Virtual Element space $Φ_h \subset H^2(Ω)$ and prove that the triad $\{Φ_h, V_h, Q_h\}$ (with $V_h$ and $Q_h$ denoting the discrete velocity and pressure spaces) is an exact Stokes complex. Furthermore, we show the computability of the associated differential operators in terms of the adopted degrees of freedom and explore also a different discretization of the convective trilinear form. The theoretical findings are supported by numerical tests.

L. Beirão da Veiga, F. Brezzi, F. Dassi et al.

We give here a simplified presentation of the lowest order Serendipity Virtual Element method, and show its use for the numerical solution of linear magneto-static problems in three dimensions. The method can be applied to very general decompositions of the computational domain (as is natural for Virtual Element Methods) and uses as unknowns the (constant) tangential component of the magnetic field $\mathbf{H}$ on each edge, and the vertex values of the Lagrange multiplier $p$ (used to enforce the solenoidality of the magnetic induction $\mathbf{B}=μ\mathbf{H}$). In this respect the method can be seen as the natural generalization of the lowest order Edge Finite Element Method (the so-called "first kind Nédélec" elements) to polyhedra of almost arbitrary shape, and as we show on some numerical examples it exhibits very good accuracy (for being a lowest order element) and excellent robustness with respect to distortions.

L. Beirão da Veiga, A. Russo, G. Vacca

In this paper we initiate the investigation of Virtual Elements with curved faces. We consider the case of a fixed curved boundary in two dimensions, as it happens in the approximation of problems posed on a curved domain or with a curved interface. While an approximation of the domain with polygons leads, for degree of accuracy $k \geq 2$, to a sub-optimal rate of convergence, we show (both theoretically and numerically) that the proposed curved VEM lead to an optimal rate of convergence.

L. Beirão da Veiga, F. Dassi, A. Russo

We develop a numerical assessment of the Virtual Element Method for the discretization of a diffusion-reaction model problem, for higher "polynomial" order k and three space dimensions. Although the main focus of the present study is to illustrate some h-convergence tests for different orders k, we also hint on other interesting aspects such as structured polyhedral Voronoi meshing, robustness in the presence of irregular grids, sensibility to the stabilization parameter and convergence with respect to the order k.

L. Beirão da Veiga, A. Chernov, L. Mascotto et al.

In the present paper we initiate the study of $hp$ Virtual Elements. We focus on the case with uniform polynomial degree across the mesh and derive theoretical convergence estimates that are explicit both in the mesh size $h$ and in the polynomial degree $p$ in the case of finite Sobolev regularity. Exponential convergence is proved in the case of analytic solutions. The theoretical convergence results are validated in numerical experiments. Finally, an initial study on the possible choice of local basis functions is included.

L. Beirão da Veiga, F. Brezzi, L. D. Marini et al.

We consider the discretization of a boundary value problem for a general linear second-order elliptic operator with smooth coefficients using the Virtual Element approach. As in [59] the problem is supposed to have a unique solution, but the associated bilinear form is not supposed to be coercive. Contrary to what was previously done for Virtual Element Methods (as for instance in [9]), we use here, in a systematic way, the L^2-projection operators as designed in [1]. In particular, the present method does not reduce to the original Virtual Element Method of [9] for simpler problems as the classical Laplace operator (apart from the lowest order cases). Numerical experiments show the accuracy and the robustness of the method, and they show as well that a simple-minded extension of the method in [9] to the case of variable coefficients produces, in general, sub-optimal results.