BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Giancarlo Sangalli

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Science

Catalog footprint

Research overview

8works
3topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

8 publications

preprint2022arXiv

A IETI-DP method for discontinuous Galerkin discretizations in Isogeometric Analysis with inexact local solvers

We construct solvers for an isogeometric multi-patch discretization, where the patches are coupled via a discontinuous Galerkin approach, which allows the consideration of discretizations that do not match on the interfaces. We solve the resulting linear system using a Dual-Primal IsogEometric Tearing and Interconnecting (IETI-DP) method. We are interested in solving the arising patch-local problems using iterative solvers since this allows the reduction of the memory footprint. We solve the patch-local problems approximately using the Fast Diagonalization method, which is known to be robust in the grid size and the spline degree. To obtain the tensor structure needed for the application of the Fast Diagonalization method, we introduce an orthogonal splitting of the local function spaces. We present a convergence theory that confirms that the condition number of the preconditioned system only grows poly-logarithmically with the grid size. The numerical experiments confirm this finding. Moreover, they show that the convergence of the overall solver only mildly depends on the spline degree. We observe a mild reduction of the computational times and a significant reduction of the memory requirements in comparison to standard IETI-DP solvers using sparse direct solvers for the local subproblems. Furthermore, the experiments indicate good scaling behavior on distributed memory machines.

Open access
preprint2022arXiv

Combining the Morris Method and Multiple Error Metrics to Assess Aquifer Characteristics and Recharge in the Lower Ticino Basin, in Italy

Groundwater flow model accuracy is often limited by the uncertainty in model parameters that characterize aquifer properties and aquifer recharge. Aquifer properties such as hydraulic conductivity can have an uncertainty spanning orders of magnitude. Meanwhile, parameters used to configure model boundary conditions can introduce additional uncertainty. In this study, the Morris Method sensitivity analysis is performed on multiple quantities of interest to assess the sensitivity of a steady-state groundwater flow model to uncertain input parameters. The Morris Method determines which of these parameters are less influential on model outputs. Uninfluential parameters can be set constant during subsequent parameter optimization to reduce computational expense. Combining multiple quantities of interest (e.g., RMSE, groundwater fluxes) when performing both the Morris Method and parameter optimization offers a more complete assessment of groundwater models, providing a more reliable and physically consistent estimate of uncertain parameters. The parameter optimization procedure also provides us an estimate of the residual uncertainty in the parameter values, resulting in a more complete estimate of the remaining uncertainty. By employing such techniques, the current study was able to estimate the aquifer hydraulic conductivity and recharge rate due to rice field irrigation in a groundwater basin in Northern Italy, revealing that a significant proportion of surficial aquifer recharge (approximately 81-94%) during the later summer is due to the flood irrigation practices applied to these fields.

Open access
preprint2021arXiv

Weighted quadrature for hierarchical B-splines

We present weighted quadrature for hierarchical B-splines to address the fast formation of system matrices arising from adaptive isogeometric Galerkin methods with suitably graded hierarchical meshes. By exploiting a local tensor-product structure, we extend the construction of weighted rules from the tensor-product to the hierarchical spline setting. The proposed algorithm has a computational cost proportional to the number of degrees of freedom and advantageous properties with increasing spline degree. To illustrate the performance of the method and confirm the theoretical estimates, a selection of 2D and 3D numerical tests is provided.

Open access
preprint2020arXiv

A domain decomposition method for Isogeometric multi-patch problems with inexact local solvers

In Isogeometric Analysis, the computational domain is often described as multi-patch, where each patch is given by a tensor product spline/NURBS parametrization. In this work we propose a FETI-like solver where local inexact solvers exploit the tensor product structure at the patch level. To this purpose, we extend to the isogeometric framework the so-called All-Floating variant of FETI, that allows us to use the Fast Diagonalization method at the patch level. We construct then a preconditioner for the whole system and prove its robustness with respect to the local mesh-size hh and patch-size HH (i.e., we have scalability). Our numerical tests confirm the theory and also show a favourable dependence of the computational cost of the method from the spline degree pp.

Open access
preprint2020arXiv

A family of C1C^1 quadrilateral finite elements

We present a novel family of C1C^1 quadrilateral finite elements, which define global C1C^1 spaces over a general quadrilateral mesh with vertices of arbitrary valency. The elements extend the construction by (Brenner and Sung, J. Sci. Comput., 2005), which is based on polynomial elements of tensor-product degree p≥6p\geq 6, to all degrees p≥3p \geq 3. Thus, we call the family of C1C^1 finite elements Brenner-Sung quadrilaterals. The proposed C1C^1 quadrilateral can be seen as a special case of the Argyris isogeometric element of (Kapl, Sangalli and Takacs, CAGD, 2019). The quadrilateral elements possess similar degrees of freedom as the classical Argyris triangles. Just as for the Argyris triangle, we additionally impose C2C^2 continuity at the vertices. In this paper we focus on the lower degree cases, that may be desirable for their lower computational cost and better conditioning of the basis: We consider indeed the polynomial quadrilateral of (bi-)degree~55, and the polynomial degrees p=3p=3 and p=4p=4 by employing a splitting into 3×33\times3 or 2×22\times2 polynomial pieces, respectively. The proposed elements reproduce polynomials of total degree pp. We show that the space provides optimal approximation order. Due to the interpolation properties, the error bounds are local on each element. In addition, we describe the construction of a simple, local basis and give for p∈{3,4,5}p\in\{3,4,5\} explicit formulas for the Bézier or B-spline coefficients of the basis functions. Numerical experiments by solving the biharmonic equation demonstrate the potential of the proposed C1C^1 quadrilateral finite element for the numerical analysis of fourth order problems, also indicating that (for p=5p=5) the proposed element performs comparable or in general even better than the Argyris triangle with respect to the number of degrees of freedom.

Open access
preprint2016arXiv

Analysis-suitable G1G^1 multi-patch parametrizations for C1C^1 isogeometric spaces

One key feature of isogeometric analysis is that it allows smooth shape functions. Indeed, when isogeometric spaces are constructed from pp-degree splines (and extensions, such as NURBS), they enjoy up to Cp−1C^{p-1} continuity within each patch. However, global continuity beyond C0C^0 on so-called multi-patch geometries poses some significant difficulties. In this work, we consider planar multi-patch domains that have a parametrization which is only C0C^0 at the patch interface. On such domains we study the hh-refinement of C1C^1-continuous isogeometric spaces. These spaces in general do not have optimal approximation properties. The reason is that the C1C^1-continuity condition easily over-constrains the solution which is, in the worst cases, fully locked to linears at the patch interface. However, recent studies by Kapl et al. have given numerical evidence that optimal convergence occurs for bilinear two-patch geometries and cubic (or higher degree) C1C^1 splines. This is the starting point of our study. We introduce the class of analysis-suitable G1G^1 geometry parametrizations, which includes piecewise bilinear parametrizations. We then analyze the structure of C1C^1 isogeometric spaces over analysis-suitable G1G^1 parametrizations and, by theoretical results and numerical testing, discuss their approximation properties. We also consider examples of geometry parametrizations that are not analysis-suitable, showing that in this case optimal convergence of C1C^1 isogeometric spaces is prevented.

Open access
preprint2016arXiv

Isogeometric preconditioners based on fast solvers for the Sylvester equation

We consider large linear systems arising from the isogeometric discretization of the Poisson problem on a single-patch domain. The numerical solution of such systems is considered a challenging task, particularly when the degree of the splines employed as basis functions is high. We consider a preconditioning strategy which is based on the solution of a Sylvester-like equation at each step of an iterative solver. We show that this strategy, which fully exploits the tensor structure that underlies isogeometric problems, is robust with respect to both mesh size and spline degree, although it may suffer from the presence of complicated geometry or coefficients. We consider two popular solvers for the Sylvester equation, a direct one and an iterative one, and we discuss in detail their implementation and efficiency for 2D and 3D problems on single-patch or conforming multi-patch NURBS geometries. Numerical experiments for problems with different domain geometries are presented, which demonstrate the potential of this approach.

Open access
preprint2015arXiv

Unstructured spline spaces for isogeometric analysis based on spline manifolds

Based on spline manifolds we introduce and study a mathematical framework for analysis-suitable unstructured B-spline spaces. In this setting the parameter domain has a manifold structure, which allows for the definition of function spaces that have a tensor-product structure locally, but not globally. This includes configurations such as B-splines over multi-patch domains with extraordinary points, analysis-suitable unstructured T-splines, or more general constructions. Within this framework, we generalize the concept of dual-compatible B-splines, which was originally developed for structured T-splines. This allows us to prove the key properties that are needed for isogeometric analysis, such as linear independence and optimal approximation properties for hh-refined meshes.

Open access

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.