Lourenco Beirao da Veiga

Lourenco Beirao da Veiga, Carlo Lovadina, Alessandro Russo

We analyse the Virtual Element Methods (VEM) on a simple elliptic model problem, allowing for more general meshes than the one typically considered in the VEM literature. For instance, meshes with arbitrarily small edges (with respect to the parent element diameter), can be dealt with. Our general approach applies to different choices of the stability form, including, for example, the "classical" one introduced in [L. Beirao da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, and A. Russo, Basic principles of virtual element methods, Math. Models Methods Appl. Sci. 23 (2013), no. 1, 199-214], and a recent one presented in [Wriggers, P., Rust, W.T., and Reddy, B.D., A virtual element method for contact, submitted for publication]. Finally, we show that the stabilization term can be simplified by dropping the contribution of the internal-to-the-element degrees of freedom. The resulting stabilization form, involving only the boundary degrees of freedom, can be used in the VEM scheme without affecting the stability and convergence properties. The numerical tests are in accordance with the theoretical predictions.

Paola F. Antonietti, Lourenco Beirao da Veiga, Simone Scacchi et al.

In this paper we develop an evolution of the $C^1$ virtual elements of minimal degree for the approximation of the Cahn-Hilliard equation. The proposed method has the advantage of being conforming in $H^2$ and making use of a very simple set of degrees of freedom, namely 3 degrees of freedom per vertex of the mesh. Moreover, although the present method is new also on triangles, it can make use of general polygonal meshes. As a theoretical and practical support, we prove the convergence of the semi-discrete scheme and investigate the performance of the fully discrete scheme through a set of numerical tests.

Edoardo Artioli, Lourenco Beirao da Veiga, Carlo Lovadina et al.

The present work deals with the formulation of a Virtual Element Method (VEM) for two dimensional structural problems. The contribution is split in two parts: in part I, the elastic problem is discussed, while in part II [3] the method is extended to material nonlinearity, considering different inelastic responses of the material. In particular, in part I a standardized procedure for the construction of all the terms required for the implementation of the method in a code is explained. The procedure is initially illustrated for the simplest case of quadrilateral virtual elements with linear approximation of displacement variables on the boundary of the element. Then, the case of polygonal elements with quadratic and, even, higher order interpolation is considered. The construction of the method is detailed, deriving the approximation of the consistent term, the required stabilization term and the loading term for all the considered virtual elements. A wide numerical investigation is performed to assess the performances of the developed virtual elements, considering different number of edges describing the elements and different order of approximations of the unknown field. Numerical results are also compared with the one recovered using the classical finite element method.

Lourenco Beirao da Veiga, Carlo Lovadina, David Mora

A low-order mimetic finite difference (MFD) method for Reissner-Mindlin plate problems is considered. Together with the source problem, the free vibration and the buckling problems are investigated. Full details about the scheme implementation are provided, and the numerical results on several different types of meshes are reported.

Lourenco Beirao da Veiga, Franco Brezzi, Franco Dassi et al.

We numerically validate the Virtual Element Method of order k for general second order elliptic problems with variable coefficients in three dimensions. Moreover, we investigate numerically also the Serendipity version of the VEM (in three dimensions) and the associated computational gain in terms of degrees of freedom.

Lourenco Beirao da Veiga, Luciano Lopez, Giuseppe Vacca

In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in time to integrate the semi-discrete Hamiltonian system. The main characteristic of MFD methods, when applied to stationary problems, is to mimic important properties of the continuous system. This approach, associated with a symplectic method for the time integration yields a full numerical procedure suitable to integrate Hamiltonian problems. A complete theoretical analysis of the method and some numerical simulations are developed in the paper.