Source author record

Paola F. Antonietti

Paola F. Antonietti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2026arXiv

A whole-brain model of amyloid beta accumulation and cerebral hypoperfusion in Alzheimer's disease

Accumulation of amyloid beta proteins is a defining feature of Alzheimer's disease, and is usually accompanied by cerebrovascular pathology. Evidence suggests that amyloid beta and cerebrovascular pathology are mutually reinforcing; in particular, amyloid beta suppresses perfusion by constricting capillaries, and hypoperfusion promotes the production of amyloid beta. Here, we propose a whole-brain model coupling amyloid beta and blood vessel through a hybrid model consisting of a reaction-diffusion system for the protein dynamics and porous-medium model of blood flow within and between vascular networks: arterial, capillary and venous. We discretize the resulting parabolic--elliptic system of PDEs by means of a high-order discontinuous Galerkin method in space and an implicit Euler scheme in time. Simulations in realistic brain geometries demonstrate the emergence of multistability, implying that a sufficiently large pathogenic protein seeds is necessary to trigger disease outbreak. Motivated by the "two-hit vascular hypothesis" of Alzheimer's disease that hypoperfusive vascular damage triggers amyloid beta pathology, we also demonstrate that localized hypoperfusion, in response to injury, can destabilize the healthy steady state and trigger brain-wide disease outbreak.

preprint2026arXiv

The Reduced Basis Multigrid scheme for the Virtual Element Method

We present a non-nested W-cycle multigrid scheme for the lowest order Virtual Element Method on polygonal meshes. To avoid the implicit definition of the Virtual Element space, which poses several issues in the computation of intergrid operators that underpin multigrid methods, the proposed scheme uses a fully-conforming auxiliary space constructed by cheaply computing the virtual basis functions via the reduced basis method.

preprint2026arXiv

Virtual Element methods for non-Newtonian shear-thickening fluid flow problems

In this work, we present a comprehensive theoretical analysis for Virtual Element discretizations of incompressible non-Newtonian flows governed by the Carreau-Yasuda constitutive law, in the shear-thickening regime (r > 2) including both degenerate (delta = 0) and non-degenerate (delta > 0) cases. The proposed Virtual Element method features two distinguishing advantages: the construction of an exactly divergence-free discrete velocity field and compatibility with general polygonal meshes. The analysis presented in this work extends a previous work, where only shear-thinning behavior (1 < r < 2) was considered. Indeed, the theoretical analysis of the shear-thickening setting requires several novel analytical tools, including: an inf-sup stability analysis of the discrete velocity-pressure coupling in non-Hilbertian norms, a stabilization term specifically designed to address the nonlinear structure as the exponent r > 2; and the introduction of a suitable discrete norm tailored to the underlying nonlinear constitutive relation. Numerical results demonstrate the practical performance of the proposed formulation.

preprint2025arXiv

Deep Learning Accelerated Algebraic Multigrid Methods for Polytopal Discretizations of Second-Order Differential Problems

Algebraic Multigrid (AMG) methods are state-of-the-art algebraic solvers for partial differential equations. Still, their efficiency depends heavily on the choice of suitable parameters and/or ingredients. Paradigmatic examples include the so-called strong threshold parameter $θ$, which controls the algebraic coarse-grid hierarchy, as well as the smoother, i.e., the relaxation methods used on the fine grid to damp out high-frequency errors. In AMG, since the coarse grids are constructed algebraically (without geometric intuition), the smoother's performance is even more critical. For the linear systems stemming from polytopal discretizations, such as Polytopal Discontinuous Galerkin (PolyDG) and Virtual Element Methods (VEM), AMG sensitivity to such choices is even more critical due to the significant variability of the underlying meshes, which results in algebraic systems with different sparsity patterns. We propose a novel deep learning approach that automatically tunes the strong threshold parameter, as well as the smoother choice in AMG solvers, for linear systems of equations arising from polytopal discretizations, thereby maximizing AMG performance. We interpret the sparse matrix resulting from polytopal discretization as a grayscale image, and by applying pooling, our neural network extracts compact features that preserve the necessary information at a low computational cost. We test various differential problems in both two- and three-dimensional settings, with heterogeneous coefficients and polygonal/polyhedral meshes, and demonstrate that the proposed approach generalizes well. In practice, we demonstrate that we can reduce AMG solver time by up to $27\%$ with minimal changes to existing PolyDG and VEM codes.

preprint2022arXiv

Discontinuous Galerkin approximation of the fully-coupled thermo-poroelastic problem

We present and analyze a discontinuous Galerkin method for the numerical modelling of the non-linear fully-coupled thermo-poroelastic problem. For the spatial discretization, we design a high-order discontinuous Galerkin method on polygonal and polyhedral grids based on a novel four-field formulation of the problem. To handle the non-linear convective transport term in the energy conservation equation we adopt a fixed-point linearization strategy. We perform a robust stability analysis for the linearized semi-discrete problem under mild requirements on the problem data. A priori hp-version error estimates in suitable energy norms are also derived. A complete set of numerical simulations is presented in order to validate the theoretical analysis, to inspect numerically the robustness properties, and to test the capability of the proposed method in a practical scenario inspired by a geothermal problem.

preprint2021arXiv

$C^1$-VEM for some variants of the Cahn-Hilliard equation: a numerical exploration

We consider the $C^1$-Virtual Element Method (VEM) for the conforming numerical approximation of some variants of the Cahn-Hilliard equation on polygonal meshes. In particular, we focus on the discretization of the advective Cahn-Hilliard problem and the Cahn-Hilliard inpainting problem. We present the numerical approximation and several numerical results to assess the efficacy of the proposed methodology.

preprint2016arXiv

Multigrid algorithms for $hp$-version Interior Penalty Discontinuous Galerkin methods on polygonal and polyhedral meshes

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.

preprint2016arXiv

On the Virtual Element Method for Topology Optimization on polygonal meshes: a numerical study

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.

preprint2015arXiv

A $C^1$ virtual element method for the Cahn-Hilliard equation with polygonal meshes

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.

preprint2015arXiv

Discontinuous Galerkin approximation of linear parabolic problems with dynamic boundary conditions

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.

preprint2014arXiv

A Two-Level Method for Mimetic Finite Difference Discretizations of Elliptic Problems

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.

preprint2014arXiv

A uniform additive Schwarz preconditioner for the $hp$-version of Discontinuous Galerkin approximations of elliptic problems

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.

preprint2013arXiv

Stability Analysis for Discontinuous Galerkin approximations of the elastodynamics problem

We consider semi-discrete discontinuous Galerkin approximations of a general elastodynamics problem, in both {\it displacement} and {\it displacement-stress} formulations. We present the stability analysis of all the methods in the natural energy norm and derive optimal a-priori error estimates. For the displacement-stress formulation, schemes preserving the total energy of the system are introduced and discussed. We include some numerical experiments in three dimensions to verify the theory.