Paola F. Antonietti

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.

P. F. Antonietti, G. Manzini, M. Verani

In this paper we address the numerical approximation of linear fourth-order elliptic problems on polygonal meshes. In particular, we present a novel nonconforming virtual element discretization of arbitrary order of accuracy for biharmonic problems. The approximation space is made of possibly discontinuous functions, thus giving rise to the fully nonconforming virtual element method. We derive optimal error estimates in a suitable (broken) energy norm and present numerical results to assess the validity of the theoretical estimates.

P. F. Antonietti, P. Houston, G. Pennesi et al.

In this article we design and analyze a class of two-level non-overlapping additive Schwarz preconditioners for the solution of the linear system of equations stemming from discontinuous Galerkin discretizations of second-order elliptic partial differential equations on polytopic meshes. The preconditioner is based on a coarse space and a non-overlapping partition of the computational domain where local solvers are applied in parallel. In particular, the coarse space can potentially be chosen to be non-embedded with respect to the finer space; indeed it can be obtained from the fine grid by employing agglomeration and edge coarsening techniques. We investigate the dependence of the condition number of the preconditioned system with respect to the diffusion coefficient and the discretization parameters, i.e., the mesh size and the polynomial degree of the fine and coarse spaces. Numerical examples are presented which confirm the theoretical bounds.

Paola F. Antonietti, Paul Houston, Xiaozhe Hu et al.

In this paper we analyze the convergence properties of two-level and W-cycle multigrid solvers for the numerical solution of the linear system of equations arising from hp-version symmetric interior penalty discontinuous Galerkin discretizations of second-order elliptic partial differential equations on polygonal/polyhedral meshes. We prove that the two-level method converges uniformly with respect to the granularity of the grid and the polynomial approximation degree p, provided that the number of smoothing steps, which depends on p, is chosen sufficiently large. An analogous result is obtained for the W-cycle multigrid algorithm, which is proved to be uniformly convergent with respect to the mesh size, the polynomial approximation degree, and the number of levels, provided the number of smoothing steps is chosen sufficiently large. Numerical experiments are presented which underpin the theoretical predictions; moreover, the proposed theoretical assumptions are not fully satisfied.

Paola F. Antonietti, Francesco Bonaldi, Ilario Mazzieri

We address the spatial discretization of an evolution problem arising from the coupling of viscoelastic and acoustic wave propagation phenomena by employing a discontinuous Galerkin scheme on polygonal and polyhedral meshes. The coupled nature of the problem is ascribed to suitable transmission conditions imposed at the interface between the solid (elastic) and fluid (acoustic) domains. We state and prove a well-posedness result for the strong formulation of the problem, present a stability analysis for the semi-discrete formulation, and finally prove an a priori $hp$-version error estimate for the resulting formulation in a suitable (mesh-dependent) energy norm. We also discuss the time integration scheme employed to obtain the fully discrete system. The convergence results are validated by numerical experiments carried out in a two-dimensional setting.

P. F. Antonietti, G. Pennesi

In this paper we analyse the convergence properties of V-cycle multigrid algorithms for the numerical solution of the linear system of equations arising from discontinuous Galerkin discretization of second-order elliptic partial differential equations on polytopal meshes. Here, the sequence of spaces that stands at the basis of the multigrid scheme is possibly non nested and is obtained based on employing agglomeration with possible edge/face coarsening. We prove that the method converges uniformly with respect to the granularity of the grid and the polynomial approximation degree p, provided that the number of smoothing steps, which depends on p, is chosen sufficiently large.

Matteo Caldana, Paola F. Antonietti, Luca Dede'

Algebraic multigrid (AMG) methods are among the most efficient solvers for linear systems of equations and they are widely used for the solution of problems stemming from the discretization of Partial Differential Equations (PDEs). The most severe limitation of AMG methods is the dependence on parameters that require to be fine-tuned. In particular, the strong threshold parameter is the most relevant since it stands at the basis of the construction of successively coarser grids needed by the AMG methods. We introduce a novel Deep Learning algorithm that minimizes the computational cost of the AMG method when used as a finite element solver. We show that our algorithm requires minimal changes to any existing code. The proposed Artificial Neural Network (ANN) tunes the value of the strong threshold parameter by interpreting the sparse matrix of the linear system as a black-and-white image and exploiting a pooling operator to transform it into a small multi-channel image. We experimentally prove that the pooling successfully reduces the computational cost of processing a large sparse matrix and preserves the features needed for the regression task at hand. We train the proposed algorithm on a large dataset containing problems with a highly heterogeneous diffusion coefficient defined in different three-dimensional geometries and discretized with unstructured grids and linear elasticity problems with a highly heterogeneous Young's modulus. When tested on problems with coefficients or geometries not present in the training dataset, our approach reduces the computational time by up to 30%.

P. F. Antonietti, N. Farenga, E. Manuzzi et al.

Agglomeration-based strategies are important both within adaptive refinement algorithms and to construct scalable multilevel algebraic solvers. In order to automatically perform agglomeration of polygonal grids, we propose the use of Machine Learning (ML) strategies, that can naturally exploit geometrical information about the mesh in order to preserve the grid quality, enhancing performance of numerical methods and reducing the overall computational cost. In particular, we employ the k-means clustering algorithm and Graph Neural Networks (GNNs) to partition the connectivity graph of a computational mesh. Moreover, GNNs have high online inference speed and the advantage to process naturally and simultaneously both the graph structure of mesh and the geometrical information, such as the areas of the elements or their barycentric coordinates. These techniques are compared with METIS, a standard algorithm for graph partitioning, which is meant to process only the graph information of the mesh. We demonstrate that performance in terms of quality metrics is enhanced for ML strategies. Such models also show a good degree of generalization when applied to more complex geometries, such as brain MRI scans, and the capability of preserving the quality of the grid. The effectiveness of these strategies is demonstrated also when applied to MultiGrid (MG) solvers in a Polygonal Discontinuous Galerkin (PolyDG) framework. In the considered experiments, GNNs show overall the best performance in terms of inference speed, accuracy and flexibility of the approach.

Paola F. Antonietti, Matteo Bruggi, Simone Scacchi et al.

It is well known that the solution of topology optimization problems may be affected both by the geometric properties of the computational mesh, which can steer the minimization process towards local (and non-physical) minima, and by the accuracy of the method employed to discretize the underlying differential problem, which may not be able to correctly capture the physics of the problem. In light of the above remarks, in this paper we consider polygonal meshes and employ the virtual element method (VEM) to solve two classes of paradigmatic topology optimization problems, one governed by nearly-incompressible and compressible linear elasticity and the other by Stokes equations. Several numerical results show the virtues of our polygonal VEM based approach with respect to more standard methods.

Paola F. Antonietti, Michele Botti, Ilario Mazzieri

This work is concerned with the analysis of a space-time finite element discontinuous Galerkin method on polytopal meshes (XT-PolydG) for the numerical discretization of wave propagation in coupled poroelastic-elastic media. The mathematical model consists of the low-frequency Biot's equations in the poroelastic medium and the elastodynamics equation for the elastic one. To realize the coupling, suitable transmission conditions on the interface between the two domains are (weakly) embedded in the formulation. The proposed PolydG discretization in space is then coupled with a dG time integration scheme, resulting in a full space-time dG discretization. We present the stability analysis for both the continuous and the semidiscrete formulations, and we derive error estimates for the semidiscrete formulation in a suitable energy norm. The method is applied to a wide set of numerical test cases to verify the theoretical bounds. Examples of physical interest are also presented to investigate the capability of the proposed method in relevant geophysical scenarios.

Paola F. Antonietti, Maurizio Grasselli, Simone Stangalino et al.

In this paper we propose and analyze a Discontinuous Galerkin method for a linear parabolic problem with dynamic boundary conditions. We present the formulation and prove stability and optimal a priori error estimates for the fully discrete scheme. More precisely, using polynomials of degree $p\geq 1$ on meshes with granularity $h$ along with a backward Euler time-stepping scheme with time-step $Δt$, we prove that the fully-discrete solution is bounded by the data and it converges, in a suitable (mesh-dependent) energy norm, to the exact solution with optimal order $h^p + Δt$. The sharpness of the theoretical estimates are verified through several numerical experiments.

P. F. Antonietti, G. Manzini, M. Verani

In this work, we exploit the capability of virtual element methods in accommodating approximation spaces featuring high-order continuity to numerically approximate differential problems of the form $Δ^p u =f$, $p\ge1$. More specifically, we develop and analyze the conforming virtual element method for the numerical approximation of polyharmonic boundary value problems, and prove an abstract result that states the convergence of the method in the energy norm.

P. F. Antonietti, B. Ayuso de Dios, S. Bertoluzza et al.

We propose and study an iterative substructuring method for an h-p Nitsche-type discretization, following the original approach introduced in [Bramble, Pasciack, Schatz (Math Comp. 1986)] for conforming methods. We prove quasi-optimality with respect to the mesh size and the polynomial degree for the proposed preconditioner. Numerical experiments asses the performance of the preconditioner and verify the theory.

Paola F. Antonietti, Mattia Corti, Sergio Gómez et al.

We investigate a two-state conformational conversion system and introduce a novel structure-preserving numerical scheme that couples a local discontinuous Galerkin space discretization with the backward Euler time-integration method. The model is first reformulated in terms of auxiliary variables involving suitable nonlinear transformations, which allow us to enforce positivity and boundedness at the numerical level. Then, we prove a discrete entropy-stability inequality, which we use to show the existence of discrete solutions, as well as to establish the convergence of the scheme by means of some discrete compactness arguments. As a by-product of the theoretical analysis, we also prove the existence of global weak solutions satisfying the system's physical bounds. Numerical results validate the theoretical results and assess the capabilities of the proposed method in practice.

Francesco Daniele, Caterina B. Leimer Saglio, Stefano Pagani et al.

Modeling and simulating how oxygen supply shapes neuronal excitability is crucial for advancing the understanding of brain function in pathological scenarios, such as ischemia. This condition is caused by a reduced blood supply, leading to the deprivation of oxygen and other metabolites; this energy deficit impairs ionic pumps and causes cellular swelling. In this work, this phenomenon is modeled through a volumetric variation law that links cell swelling to local oxygen concentration and the percentage of blood flow reduction. The swelling law links volume changes to local oxygen and the degree of blood-flow depression, providing a simple mechanistic pathway from hypoxia to tortuosity-driven transport impairment. The interplay between oxygen supply and excitability in brain tissue is described by coupling the monodomain model with specific neuronal ionic and metabolic models that characterize ion and metabolite concentration dynamics. The numerical approximation of this coupled multiscale problem is particularly challenging, owing to the presence of sharp and fast-propagating wavefronts and complex geometrical domains. To address these challenges, suitable space- and time-adaptive schemes are employed to capture the action potential dynamics accurately. This multiscale model is discretized in space with the high-order p-adaptive discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) and integrated in time with a Crank-Nicolson scheme. We numerically investigate different pathological scenarios on a two-dimensional idealized square domain and on a realistic geometry, both discretized with a polygonal grid, analyzing how subclinical and severe ischemia can affect brain electrophysiology and metabolic concentrations.

Matteo Caldana, Paola F. Antonietti, Luca Dede'

Finite element-based high-order solvers of conservation laws offer large accuracy but face challenges near discontinuities due to the Gibbs phenomenon. Artificial viscosity is a popular and effective solution to this problem based on physical insight. In this work, we present a physics-informed machine learning algorithm to automate the discovery of artificial viscosity models in a non-supervised paradigm. The algorithm is inspired by reinforcement learning and trains a neural network acting cell-by-cell (the viscosity model) by minimizing a loss defined as the difference with respect to a reference solution thanks to automatic differentiation. This enables a dataset-free training procedure. We prove that the algorithm is effective by integrating it into a state-of-the-art Runge-Kutta discontinuous Galerkin solver. We showcase several numerical tests on scalar and vectorial problems, such as Burgers' and Euler's equations in one and two dimensions. Results demonstrate that the proposed approach trains a model that is able to outperform classical viscosity models. Moreover, we show that the learnt artificial viscosity model is able to generalize across different problems and parameters.

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.

P. F. Antonietti, E. Manuzzi

We propose new strategies to handle polygonal grids refinement based on Convolutional Neural Networks (CNNs). We show that CNNs can be successfully employed to identify correctly the "shape" of a polygonal element so as to design suitable refinement criteria to be possibly employed within adaptive refinement strategies. We propose two refinement strategies that exploit the use of CNNs to classify elements' shape, at a low computational cost. We test the proposed idea considering two families of finite element methods that support arbitrarily shaped polygonal elements, namely Polygonal Discontinuous Galerkin (PolyDG) methods and Virtual Element Methods (VEMs). We demonstrate that the proposed algorithms can greatly improve the performance of the discretization schemes both in terms of accuracy and quality of the underlying grids. Moreover, since the training phase is performed off-line and is independent of the differential model the overall computational costs are kept low.

P. F. Antonietti, L. Mascotto, M. Verani

We present a multigrid algorithm for the solution of the linear systems of equations stemming from the $p-$version of the Virtual Element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed decreasing progressively the polynomial approximation degree of the Virtual Element space, as in standard $p$-multigrid schemes. The construction of the interspace operators relies on auxiliary Virtual Element spaces, where it is possible to compute higher order polynomial projectors. We prove that the multigrid scheme is uniformly convergent, provided the number of smoothing steps is chosen sufficiently large. We also demonstrate that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom that can be employed to accelerate the convergence of classical Krylov-based iterative schemes. Numerical experiments validate the theoretical results.

Paola F. Antonietti, Marco Sarti, Marco Verani et al.

In this paper we design and analyze a uniform preconditioner for a class of high order Discontinuous Galerkin schemes. The preconditioner is based on a space splitting involving the high order conforming subspace and results from the interpretation of the problem as a nearly-singular problem. We show that the proposed preconditioner exhibits spectral bounds that are uniform with respect to the discretization parameters, i.e., the mesh size, the polynomial degree and the penalization coefficient. The theoretical estimates obtained are supported by several numerical simulations.

Paola F. Antonietti, Marco Verani, Ludmil Zikatanov

We propose and analyze a two-level method for mimetic finite difference approximations of second order elliptic boundary value problems. We prove that the two-level algorithm is uniformly convergent, i.e., the number of iterations needed to achieve convergence is uniformly bounded independently of the characteristic size of the underling partition. We also show that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom. Numerical results that validate the theory are also presented.