112 research outputs found
A robust a posteriori estimator for the residual-free bubbles method applied to advection-diffusion problems
We develop the a posteriori error analysis for the RFB method, applied to the linear advection-diffusion problem: the numerical error, measured in suitable norms, is estimated in terms of the numerical residual. The robustness is investiged, in the sense that we prove uniform equivalence between a norm of the numerical residual and a particular norm of the error
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
Matrix-free weighted quadrature for a computationally efficient isogeometric -method
The -method is the isogeometric method based on splines (or NURBS, etc.)
with maximum regularity. When implemented following the paradigms of classical
finite element methods, the computational resources required by the method
are prohibitive even for moderate degree. In order to address this issue, we
propose a matrix-free strategy combined with weighted quadrature, which is an
ad-hoc strategy to compute the integrals of the Galerkin system. Matrix-free
weighted quadrature (MF-WQ) speeds up matrix operations, and, perhaps even more
important, greatly reduces memory consumption. Our strategy also requires an
efficient preconditioner for the linear system iterative solver. In this work
we deal with an elliptic model problem, and adopt a preconditioner based on the
Fast Diagonalization method, an old idea to solve Sylvester-like equations. Our
numerical tests show that the isogeometric solver based on MF-WQ is faster than
standard approaches (where the main cost is the matrix formation by standard
Gaussian quadrature) even for low degree. But the main achievement is that,
with MF-WQ, the -method gets orders of magnitude faster by increasing the
degree, given a target accuracy. Therefore, we are able to show the
superiority, in terms of computational efficiency, of the high-degree
-method with respect to low-degree isogeometric discretizations. What we
present here is applicable to more complex and realistic differential problems,
but its effectiveness will depend on the preconditioner stage, which is as
always problem-dependent. This situation is typical of modern high-order
methods: the overall performance is mainly related to the quality of the
preconditioner
Robust isogeometric preconditioners for the Stokes system based on the Fast Diagonalization method
In this paper we propose a new class of preconditioners for the isogeometric
discretization of the Stokes system. Their application involves the solution of
a Sylvester-like equation, which can be done efficiently thanks to the Fast
Diagonalization method. These preconditioners are robust with respect to both
the spline degree and mesh size. By incorporating information on the geometry
parametrization and equation coefficients, we maintain efficiency on
non-trivial computational domains and for variable kinematic viscosity. In our
numerical tests we compare to a standard approach, showing that the overall
iterative solver based on our preconditioners is significantly faster.Comment: 31 pages, 4 figure
Construction of analysis-suitable planar multi-patch parameterizations
Isogeometric analysis allows to define shape functions of global
continuity (or of higher continuity) over multi-patch geometries. The
construction of such -smooth isogeometric functions is a non-trivial
task and requires particular multi-patch parameterizations, so-called
analysis-suitable (in short, AS-) parameterizations, to ensure
that the resulting isogeometric spaces possess optimal approximation
properties, cf. [7]. In this work, we show through examples that it is possible
to construct AS- multi-patch parameterizations of planar domains, given
their boundary. More precisely, given a generic multi-patch geometry, we
generate an AS- multi-patch parameterization possessing the same
boundary, the same vertices and the same first derivatives at the vertices, and
which is as close as possible to this initial geometry. Our algorithm is based
on a quadratic optimization problem with linear side constraints. Numerical
tests also confirm that isogeometric spaces over AS- multi-patch
parameterized domains converge optimally under mesh refinement, while for
generic parameterizations the convergence order is severely reduced
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 -refined meshes
A sparse-grid isogeometric solver
Isogeometric Analysis (IGA) typically adopts tensor-product splines and NURBS
as a basis for the approximation of the solution of PDEs. In this work, we
investigate to which extent IGA solvers can benefit from the so-called
sparse-grids construction in its combination technique form, which was first
introduced in the early 90s in the context of the approximation of
high-dimensional PDEs. The tests that we report show that, in accordance to the
literature, a sparse-grid construction can indeed be useful if the solution of
the PDE at hand is sufficiently smooth. Sparse grids can also be useful in the
case of non-smooth solutions when some a-priori knowledge on the location of
the singularities of the solution can be exploited to devise suitable
non-equispaced meshes. Finally, we remark that sparse grids can be seen as a
simple way to parallelize pre-existing serial IGA solvers in a straightforward
fashion, which can be beneficial in many practical situations.Comment: updated version after revie
The Argyris isogeometric space on unstructured multi-patch planar domains
Multi-patch spline parametrizations are used in geometric design and
isogeometric analysis to represent complex domains. We deal with a particular
class of planar multi-patch spline parametrizations called
analysis-suitable (AS-) multi-patch parametrizations (Collin,
Sangalli, Takacs; CAGD, 2016). This class of parametrizations has to satisfy
specific geometric continuity constraints, and is of importance since it allows
to construct, on the multi-patch domain, isogeometric spaces with optimal
approximation properties. It was demonstrated in (Kapl, Sangalli, Takacs; CAD,
2018) that AS- multi-patch parametrizations are suitable for modeling
complex planar multi-patch domains.
In this work, we construct a basis, and an associated dual basis, for a
specific isogeometric spline space over a given AS-
multi-patch parametrization. We call the space the Argyris
isogeometric space, since it is across interfaces and at all
vertices and generalizes the idea of Argyris finite elements to tensor-product
splines. The considered space is a subspace of the entire
isogeometric space , which maintains the reproduction
properties of traces and normal derivatives along the interfaces. Moreover, it
reproduces all derivatives up to second order at the vertices. In contrast to
, the dimension of does not depend on the domain
parametrization, and admits a basis and dual basis which possess
a simple explicit representation and local support.
We conclude the paper with some numerical experiments, which exhibit the
optimal approximation order of the Argyris isogeometric space and
demonstrate the applicability of our approach for isogeometric analysis
A family of quadrilateral finite elements
We present a novel family of quadrilateral finite elements, which
define global 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 , to all degrees . Thus, we call the family of
finite elements Brenner-Sung quadrilaterals. The proposed 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 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~, and the polynomial degrees
and by employing a splitting into or polynomial
pieces, respectively.
The proposed elements reproduce polynomials of total degree . 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
explicit formulas for the B\'{e}zier or B-spline coefficients of the basis
functions. Numerical experiments by solving the biharmonic equation demonstrate
the potential of the proposed quadrilateral finite element for the
numerical analysis of fourth order problems, also indicating that (for )
the proposed element performs comparable or in general even better than the
Argyris triangle with respect to the number of degrees of freedom
- …