F. Dassi

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, 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, 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.

P. F. Antonietti, F. Dassi, E. Manuzzi

We propose two new strategies based on Machine Learning techniques to handle polyhedral grid refinement, to be possibly employed within an adaptive framework. The first one employs the k-means clustering algorithm to partition the points of the polyhedron to be refined. This strategy is a variation of the well known Centroidal Voronoi Tessellation. The second one employs Convolutional Neural Networks to classify the "shape" of an element so that "ad-hoc" refinement criteria can be defined. This strategy can be used to enhance existing refinement strategies, including the k-means strategy, at a low online computational cost. We test the proposed algorithms considering two families of finite element methods that support arbitrarily shaped polyhedral elements, namely the Virtual Element Method (VEM) and the Polygonal Discontinuous Galerkin (PolyDG) method. We demonstrate that these strategies do preserve the structure and the quality of the underlaying grids, reducing the overall computational cost and mesh complexity.

L. Beirao da Veiga, F. Dassi, G. Vacca

The present paper has two objectives. On one side, we develop and test numerically divergence free Virtual Elements in three dimensions, for variable ``polynomial'' order. These are the natural extension of the two-dimensional divergence free VEM elements, with some modification that allows for a better computational efficiency. We test the element's performance both for the Stokes and (diffusion dominated) Navier-Stokes equation. The second, and perhaps main, motivation is to show that our scheme, also in three dimensions, enjoys an underlying discrete Stokes complex structure. We build a pair of virtual discrete spaces based on general polytopal partitions, the first one being scalar and the second one being vector valued, such that when coupled with our velocity and pressure spaces, yield a discrete Stokes complex.

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.