Source author record

Annalisa Buffa

Annalisa Buffa 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
4topics
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 Recovery-Based Error Indicator for Finite Difference Methods

A novel recovery-based error indicator for high-order Finite Difference Methods, based on post-processing of the Finite Difference values is presented. The values obtained on the Finite Difference grid are interpolated into a suitable polynomial Finite Element space. A recovery-based error indicator, with the polynomial-preserving property, is then applied to estimate the gradient error. The performance and accuracy of the proposed error indicator are demonstrated through several numerical experiments, including the two-dimensional Poisson problem solved using second- and fourth-order finite difference schemes. Additional experiments are conducted on elliptic problems with discontinuous coefficients, as well as on the two and three-dimensional wave equation in homogeneous media with second- and fourth-order finite differences, and in heterogeneous media with second-order finite differences.

preprint2022arXiv

An a posteriori error estimator for isogeometric analysis on trimmed geometries

Trimming consists of cutting away parts of a geometric domain, without reconstructing a global parametrization (meshing). It is a widely used operation in computer aided design, which generates meshes that are unfitted with the described physical object. This paper develops an adaptive mesh refinement strategy on trimmed geometries in the context of hierarchical B-spline based isogeometric analysis. A residual a posteriori estimator of the energy norm of the numerical approximation error is derived, in the context of the Poisson equation. The estimator is proven to be reliable, independently of the number of hierarchical levels and of the way the trimmed boundaries cut the underlying mesh. Numerical experiments are performed to validate the presented theory, and to show that the estimator's effectivity index is independent of the size of the active part of the trimmed mesh elements.

preprint2021arXiv

Robust Numerical Integration on Curved Polyhedra Based on Folded Decompositions

We present a novel method to perform numerical integration over curved polyhedra enclosed by high-order parametric surfaces. Such a polyhedron is first decomposed into a set of triangular and/or rectangular pyramids, whose certain faces correspond to the given parametric surfaces. Each pyramid serves as an integration cell with a geometric mapping from a standard parent domain (e.g., a unit cube), where the tensor-product Gauss quadrature is adopted. As no constraint is imposed on the decomposition, certain resulting pyramids may intersect with themselves, and thus their geometric mappings may present negative Jacobian values. We call such cells the folded cells and refer to the corresponding decomposition as a folded decomposition. We show that folded cells do not cause any issues in practice as they are only used to numerically compute certain integrals of interest. The same idea can be applied to planar curved polygons as well. We demonstrate both theoretically and numerically that folded cells can retain the same accuracy as the cells with strictly positive Jacobians. On the other hand, folded cells allow for a much easier and much more flexible decomposition for general curved polyhedra, on which one can robustly compute integrals. In the end, we show that folded cells can flexibly and robustly accommodate real-world complex geometries by presenting several examples in the context of immersed isogeometric analysis, where involved sharp features can be well respected in generating integration cells.

preprint2020arXiv

A minimal stabilization procedure for Isogeometric methods on trimmed geometries

Trimming is a common operation in CAD, and, in its simplest formulation, consists in removing superfluous parts from a geometric entity described via splines (a spline patch). After trimming the geometric description of the patch remains unchanged, but the underlying mesh is unfitted with the physical object. We discuss the main problems arising when solving elliptic PDEs on a trimmed domain. First we prove that, even when Dirichlet boundary conditions are weakly enforced using Nitsche's method, the resulting method suffers lack of stability. Then, we develop novel stabilization techniques based on a modification of the variational formulation, which allow us to recover well-posedness and guarantee accuracy. Optimal a priori error estimates are proven, and numerical examples confirming the theoretical results are provided.

preprint2020arXiv

A projected super-penalty method for the $C^1$-coupling of multi-patch isogeometric Kirchhoff plates

This work focuses on the development of a super-penalty strategy based on the $L^2$-projection of suitable coupling terms to achieve $C^1$-continuity between non-conforming multi-patch isogeometric Kirchhoff plates. In particular, the choice of penalty parameters is driven by the underlying perturbed saddle point problem from which the Lagrange multipliers are eliminated and is performed to guarantee the optimal accuracy of the method. Moreover, by construction, the method does not suffer from locking also on very coarse meshes. We demonstrate the applicability of the proposed coupling algorithm to Kirchhoff plates by studying several benchmark examples discretized by non-conforming meshes. In all cases, we recover the optimal rates of convergence achievable by B-splines where we achieve a substantial gain in accuracy per degree-of-freedom compared to other choices of the penalty parameters.

preprint2019arXiv

A hierarchical approach to the a posteriori error estimation of isogeometric Kirchhoff plates and Kirchhoff-Love shells

This work focuses on the development of a posteriori error estimates for fourth-order, elliptic, partial differential equations. In particular, we propose a novel algorithm to steer an adaptive simulation in the context of Kirchhoff plates and Kirchhoff-Love shells by exploiting the local refinement capabilities of hierarchical B-splines. The method is based on the solution of an auxiliary residual-like variational problem, formulated by means of a space of localized spline functions. This space is characterized by $C^1$ continuous B-splines with compact support on each active element of the hierarchical mesh. We demonstrate the applicability of the proposed estimator to Kirchhoff plates and Kirchhoff-Love shells by studying several benchmark problems which exhibit both smooth and singular solutions. In all cases, we obtain optimal asymptotic rates of convergence for the error measured in the energy norm and an excellent approximation of the true error.

preprint2019arXiv

Adaptive isogeometric analysis on two-dimensional trimmed domains based on a hierarchical approach

The focus of this work is on the development of an error-driven isogeometric framework, capable of automatically performing an adaptive simulation in the context of second- and fourth-order, elliptic partial differential equations defined on two-dimensional trimmed domains. The method is steered by an a posteriori error estimator, which is computed with the aid of an auxiliary residual-like problem formulated onto a space spanned by splines with single element support. The local refinement of the basis is achieved thanks to the use of truncated hierarchical B-splines. We prove numerically the applicability of the proposed estimator to various engineering-relevant problems, namely the Poisson problem, linear elasticity and Kirchhoff-Love shells, formulated on trimmed geometries. In particular, we study several benchmark problems which exhibit both smooth and singular solutions, where we recover optimal asymptotic rates of convergence for the error measured in the energy norm and we observe a substantial increase in accuracy per-degree-of-freedom compared to uniform refinement. Lastly, we show the applicability of our framework to the adaptive shell analysis of an industrial-like trimmed geometry modeled in the commercial software Rhinoceros, which represents the B-pillar of a car.

preprint2016arXiv

A posteriori error estimators for hierarchical B-spline discretizations

In this article we develop function-based a posteriori error estimators for the solution of linear second order elliptic problems considering hierarchical spline spaces for the Galerkin discretization. We prove a global upper bound for the energy error. The theory hinges on some weighted Poincaré type inequalities, where the B-spline basis functions are the weights appearing in the norms. Such inequalities are derived following the lines in [Veeser and Verfürth, 2009], where the case of standard finite elements is considered. Additionally, we present numerical experiments that show the efficiency of the error estimators independently of the degree of the splines used for discretization, together with an adaptive algorithm guided by these local estimators that yields optimal meshes and rates of convergence, exhibiting an excellent performance.

preprint2015arXiv

Adaptive isogeometric methods with hierarchical splines: error estimator and convergence

The problem of developing an adaptive isogeometric method (AIGM) for solving elliptic second-order partial differential equations with truncated hierarchical B-splines of arbitrary degree and different order of continuity is addressed. The adaptivity analysis holds in any space dimensions. We consider a simple residual-type error estimator for which we provide a posteriori upper and lower bound in terms of local error indicators, taking also into account the critical role of oscillations as in a standard adaptive finite element setting. The error estimates are properly combined with a simple marking strategy to define a sequence of admissible locally refined meshes and corresponding approximate solutions. The design of a refine module that preserves the admissibility of the hierarchical mesh configuration between two consectutive steps of the adaptive loop is presented. The contraction property of the quasi-error, given by the sum of the energy error and the scaled error estimator, leads to the convergence proof of the AIGM.

preprint2015arXiv

Complexity of hierarchical refinement for a class of admissible mesh configurations

An adaptive isogeometric method based on $d$-variate hierarchical spline constructions can be derived by considering a refine module that preserves a certain class of admissibility between two consecutive steps of the adaptive loop [6]. In this paper we provide a complexity estimate, i.e., an estimate on how the number of mesh elements grows with respect to the number of elements that are marked for refinement by the adaptive strategy. Our estimate is in the line of the similar ones proved in the finite element context, [3,24].

preprint2015arXiv

New refinable spaces and local approximation estimates for hierarchical splines

We study the local approximation properties in hierarchical spline spaces through multiscale quasi-interpolation operators. This construction suggests the analysis of a subspace of the classical hierarchical spline space (Vuong et al., 2011) which still satisfies the essential properties of the full space. The B-spline basis of such a subspace can be constructed using parent-children relations only, making it well adapted to local refinement algorithms.

preprint2015arXiv

The Influence of Quadrature Errors on Isogeometric Mortar Methods

Mortar methods have recently been shown to be well suited for isogeometric analysis. We review the recent mathematical analysis and then investigate the variational crime introduced by quadrature formulas for the coupling integrals. Motivated by finite element observations, we consider a quadrature rule purely based on the slave mesh as well as a method using quadrature rules based on the slave mesh and on the master mesh, resulting in a non-symmetric saddle point problem. While in the first case reduced convergence rates can be observed, in the second case the influence of the variational crime is less significant.

preprint2014arXiv

Isogeometric mortar methods

The application of mortar methods in the framework of isogeometric analysis is investigated theoretically as well as numerically. For the Lagrange multiplier two choices of uniformly stable spaces are presented, both of them are spline spaces but of a different degree. In one case, we consider an equal order pairing for which a cross point modification based on a local degree reduction is required. In the other case, the degree of the dual space is reduced by two compared to the primal. This pairing is proven to be inf-sup stable without any necessary cross point modification. Several numerical examples confirm the theoretical results and illustrate additional aspects. Keywords: isogeometric analysis, mortar methods, inf-sup stability, cross point modification.