BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Marco Verani

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

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

21works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

21 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Agglomeration-based geometric multigrid schemes for the Virtual Element Method

In this paper we analyse the convergence properties of two-level, W-cycle and V-cycle agglomeration-based geometric multigrid schemes for the numerical solution of the linear system of equations stemming from the lowest order $C^0$-conforming Virtual Element discretization of two-dimensional second-order elliptic partial differential equations. The sequence of agglomerated tessellations are nested, but the corresponding multilevel virtual discrete spaces are generally non-nested thus resulting into non-nested multigrid algorithms. We prove the uniform convergence of the two-level method with respect to the mesh size and the uniform convergence of the W-cycle and the V-cycle multigrid algorithms with respect to the mesh size and the number of levels. Numerical experiments confirm the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The conforming virtual element method for polyharmonic and elastodynamics problems: a review

In this paper, we review recent results on the conforming virtual element approximation of polyharmonic and elastodynamics problems. The structure and the content of this review is motivated by three paradigmatic examples of applications: classical and anisotropic Cahn-Hilliard equation and phase field models for brittle fracture, that are briefly discussed in the first part of the paper. We present and discuss the mathematical details of the conforming virtual element approximation of linear polyharmonic problems, the classical Cahn-Hilliard equation and linear elastodynamics problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anisotropic a posteriori error estimate for the Virtual Element Method

We derive an anisotropic a posteriori error estimate for the adaptive conforming Virtual Element approximation of a paradigmatic two-dimensional elliptic problem. In particular, we introduce a quasi-interpolant operator and exploit its approximation results to prove the reliability of the error indicator. We design and implement the corresponding adaptive polygonal anisotropic algorithm. Several numerical tests assess the superiority of the proposed algorithm in comparison with standard polygonal isotropic mesh refinement schemes.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On p-Robust Saturation for hp-AFEM

We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK, REFINE, with ESTIMATE being implemented using the $p$-robust equilibrated flux estimator, and MARK being Dörfler marking. As a refinement strategy we employ $p$-refinement. We investigate the question by which amount the local polynomial degree on any marked patch has to be increase in order to achieve a $p$-independent error reduction. The resulting adaptive method can be turned into an instance optimal $hp$-adaptive method by the addition of a coarsening routine.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive Spectral Galerkin Methods with Dynamic Marking

The convergence and optimality theory of adaptive Galerkin methods is almost exclusively based on the Dörfler marking. This entails a fixed parameter and leads to a contraction constant bounded below away from zero. For spectral Galerkin methods this is a severe limitation which affects performance. We present a dynamic marking strategy that allows for a super-linear relation between consecutive discretization errors, and show exponential convergence with linear computational complexity whenever the solution belongs to a Gevrey approximation class.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence and Optimality of hp-AFEM

We design and analyze an adaptive $hp$-finite element method (hp-AFEM) in dimensions $n=1,2$. The algorithm consists of iterating two routines: hp-NEARBEST finds a near-best $hp$-approximation of the current discrete solution and data to a desired accuracy, and REDUCE improves the discrete solution to a finer but comparable accuracy. The former hinges on a recent algorithm by Binev for adaptive $hp$-approximation, and acts as a coarsening step. We prove convergence and instance optimality.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Contraction and optimality properties of an adaptive Legendre-Galerkin method: the multi-dimensional case

We analyze the theoretical properties of an adaptive Legendre-Galerkin method in the multidimensional case. After the recent investigations for Fourier-Galerkin methods in a periodic box and for Legendre-Galerkin methods in the one dimensional setting, the present study represents a further step towards a mathematically rigorous understanding of adaptive spectral/$hp$ discretizations of elliptic boundary-value problems. The main contribution of the paper is a careful construction of a multidimensional Riesz basis in $H^1$, based on a quasi-orthonormalization procedure. This allows us to design an adaptive algorithm, to prove its convergence by a contraction argument, and to discuss its optimality properties (in the sense of non-linear approximation theory) in certain sparsity classes of Gevrey type.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High order discontinuous Galerkin methods on surfaces

We derive and analyze high order discontinuous Galerkin methods for second-order elliptic problems on implicitely defined surfaces in $\mathbb{R}^{3}$. This is done by carefully adapting the unified discontinuous Galerkin framework of Arnold et al. [2002] on a triangulated surface approximating the smooth surface. We prove optimal error estimates in both a (mesh dependent) energy norm and the $L^2$ norm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shape optimization for Stokes flow: a reference domain approach

In this paper we analyze a shape optimization problem, with Stokes equations as the state problem, defined on a domain with a part of the boundary that is described as the graph of the control function. The state problem formulation is mapped onto a reference domain, which is independent of the control function, and the analysis is mainly led on such domain. The existence of an optimal control function is proved, and optimality conditions are derived. After the analytical inspection of the problem, finite element discretization is considered for both the control function and the state variables, and a priori convergence error estimates are derived. Numerical experiments assess the validity of the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the decay of the inverse of matrices that are sum of Kronecker products

Decay patterns of matrix inverses have recently attracted considerable interest, due to their relevance in numerical analysis, and in applications requiring matrix function approximations. In this paper we analyze the decay pattern of the inverse of banded matrices in the form $S=M \otimes I_n + I_n \otimes M$ where $M$ is tridiagonal, symmetric and positive definite, $I_n$ is the identity matrix, and $\otimes$ stands for the Kronecker product. It is well known that the inverses of banded matrices exhibit an exponential decay pattern away from the main diagonal. However, the entries in $S^{-1}$ show a non-monotonic decay, which is not caught by classical bounds. By using an alternative expression for $S^{-1}$, we derive computable upper bounds that closely capture the actual behavior of its entries. We also show that similar estimates can be obtained when $M$ has a larger bandwidth, or when the sum of Kronecker products involves two different matrices. Numerical experiments illustrating the new bounds are also reported.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive Fourier-Galerkin Methods

We study the performance of adaptive Fourier-Galerkin methods in a periodic box in $\mathbb{R}^d$ with dimension $d\ge 1$. These methods offer unlimited approximation power only restricted by solution and data regularity. They are of intrinsic interest but are also a first step towards understanding adaptivity for the $hp$-FEM. We examine two nonlinear approximation classes, one classical corresponding to algebraic decay of Fourier coefficients and another associated with exponential decay. We study the sparsity classes of the residual and show that they are the same as the solution for the algebraic class but not for the exponential one. This possible sparsity degradation for the exponential class can be compensated with coarsening, which we discuss in detail. We present several adaptive Fourier algorithms, and prove their contraction and optimal cardinality properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Contraction and optimality properties of adaptive Legendre-Galerkin methods: the 1-dimensional case

As a first step towards a mathematically rigorous understanding of adaptive spectral/$hp$ discretizations of elliptic boundary-value problems, we study the performance of adaptive Legendre-Galerkin methods in one space dimension. These methods offer unlimited approximation power only restricted by solution and data regularity. Our investigation is inspired by a similar study that we recently carried out for Fourier-Galerkin methods in a periodic box. We first consider an "ideal" algorithm, which we prove to be convergent at a fixed rate. Next we enhance its performance, consistently with the expected fast error decay of high-order methods, by activating a larger set of degrees of freedom at each iteration. We guarantee optimality (in the non-linear approximation sense) by incorporating a coarsening step. Optimality is measured in terms of certain sparsity classes of the Gevrey type, which describe a (sub-)exponential decay of the best approximation error.

0 citations0 reviewsNo code linkedOpen access