Researcher profile

Claudio Canuto

Claudio Canuto contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 published item(s)

preprint2022arXiv

A posteriori error analysis and adaptivity for a VEM discretization of the Navier-Stokes equations

We consider the Virtual Element method (VEM) introduced by Beirão da Veiga, Lovadina and Vacca in 2016 for the numerical solution of the steady, incompressible Navier-Stokes equations; the method has arbitrary order $k \geq 2$ and guarantees divergence-free velocities. For such discretization, we develop a residual-based a posteriori error estimator, which is a combination of standard terms in VEM analysis (residual terms, data oscillation, and VEM stabilization), plus some other terms originated by the VEM discretization of the nonlinear convective term. We show that a linear combination of the velocity and pressure errors is upper-bounded by a multiple of the estimator (reliability). We also establish some efficiency results, involving lower bounds of the error. Some numerical tests illustrate the performance of the estimator and of its components while refining the mesh uniformly, yielding the expected decay rate. At last, we apply an adaptive mesh refinement strategy to the computation of the low-Reynolds flow around a square cylinder inside a channel.

preprint2022arXiv

Variational Physics Informed Neural Networks: the role of quadratures and test functions

In this work we analyze how quadrature rules of different precisions and piecewise polynomial test functions of different degrees affect the convergence rate of Variational Physics Informed Neural Networks (VPINN) with respect to mesh refinement, while solving elliptic boundary-value problems. Using a Petrov-Galerkin framework relying on an inf-sup condition, we derive an a priori error estimate in the energy norm between the exact solution and a suitable high-order piecewise interpolant of a computed neural network. Numerical experiments confirm the theoretical predictions and highlight the importance of the inf-sup condition. Our results suggest, somehow counterintuitively, that for smooth solutions the best strategy to achieve a high decay rate of the error consists in choosing test functions of the lowest polynomial degree, while using quadrature formulas of suitably high precision.

preprint2020arXiv

Analysis of the Shifted Boundary Method for the Poisson Problem in General Domains

The shifted boundary method (SBM) is an approximate domain method for boundary value problems, in the broader class of unfitted/embedded/immersed methods. It has proven to be quite efficient in handling problems with complex geometries, ranging from Poisson to Darcy, from Navier-Stokes to elasticity and beyond. The key feature of the SBM is a {\it shift} in the location where Dirichlet boundary conditions are applied - from the true to a surrogate boundary - and an appropriate modification (again, a {\it shift}) of the value of the boundary conditions, in order to reduce the consistency error. In this paper we provide a sound analysis of the method in smooth and non-smooth domains, highlighting the influence of geometry and distance between exact and surrogate boundaries upon the convergence rate. Without loss of generality, we consider the Poisson problem with Dirichlet boundary conditions as a model and we first detail a procedure to obtain the crucial shifting between the surrogate and the true boundaries. Next, we give a sufficient condition for the well-posedness and stability of the discrete problem. The behavior of the consistency error arising from shifting the boundary conditions is thoroughly analyzed, for smooth boundaries and for boundaries with corners and edges. The convergence rate is proven to be optimal in the energy norm, and is further enhanced in the $L^2$-norm.

preprint2020arXiv

Efficient computation of bifurcation diagrams with a deflated approach to reduced basis spectral element method

The majority of the most common physical phenomena can be described using partial differential equations (PDEs). However, they are very often characterized by strong nonlinearities. Such features lead to the coexistence of multiple solutions studied by the bifurcation theory. Unfortunately, in practical scenarios, one has to exploit numerical methods to compute the solutions of systems of PDEs, even if the classical techniques are usually able to compute only a single solution for any value of a parameter when more branches exist. In this work we implemented an elaborated deflated continuation method, that relies on the spectral element method (SEM) and on the reduced basis (RB) one, to efficiently compute bifurcation diagrams with more parameters and more bifurcation points. The deflated continuation method can be obtained combining the classical continuation method and the deflation one: the former is used to entirely track each known branch of the diagram, while the latter is exploited to discover the new ones. Finally, when more than one parameter is considered, the efficiency of the computation is ensured by the fact that the diagrams can be computed during the online phase while, during the offline one, one only has to compute one-dimensional diagrams. In this work, after a more detailed description of the method, we will show the results that can be obtained using it to compute a bifurcation diagram associated with a problem governed by the Navier-Stokes equations.