BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Thomas Takacs

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.DGmath.FAmath.OC

Catalog footprint

Research overview

7works
5topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

7 publications

preprint2022arXiv

Adaptive Refinement for Unstructured T-Splines with Linear Complexity

We present an adaptive refinement algorithm for T-splines on unstructured 2D meshes. While for structured 2D meshes, one can refine elements alternatingly in horizontal and vertical direction, such an approach cannot be generalized directly to unstructured meshes, where no two unique global mesh directions can be assigned. To resolve this issue, we introduce the concept of direction indices, i.e., integers associated to each edge, which are inspired by theory on higher-dimensional structured T-splines. Together with refinement levels of edges, these indices essentially drive the refinement scheme. We combine these ideas with an edge subdivision routine that allows for I-nodes, yielding a very flexible refinement scheme that nicely distributes the T-nodes, preserving global linear independence, analysis-suitability (local linear independence) except in the vicinity of extraordinary nodes, sparsity of the system matrix, and shape regularity of the mesh elements. Further, we show that the refinement procedure has linear complexity in the sense of guaranteed upper bounds on a) the distance between marked and additionally refined elements, and on b) the ratio of the numbers of generated and marked mesh elements.

Open access
preprint2022arXiv

Locally refined quad meshing for linear elasticity problems based on convolutional neural networks

In this paper we propose a method to generate suitably refined finite element meshes using neural networks. As a model problem we consider a linear elasticity problem on a planar domain (possibly with holes) having a polygonal boundary. We impose boundary conditions by fixing the position of a part of the boundary and applying a force on another part of the boundary. The resulting displacement and distribution of stresses depend on the geometry of the domain and on the boundary conditions. When applying a standard Galerkin discretization using quadrilateral finite elements, one usually has to perform adaptive refinement to properly resolve maxima of the stress distribution. Such an adaptive scheme requires a local error estimator and a corresponding local refinement strategy. The overall costs of such a strategy are high. We propose to reduce the costs of obtaining a suitable discretization by training a neural network whose evaluation replaces this adaptive refinement procedure. We set up a single network for a large class of possible domains and boundary conditions and not on a single domain of interest. The computational domain and boundary conditions are interpreted as images, which are suitable inputs for convolution neural networks. We use the U-net architecture and we devise training strategies by dividing the possible inputs into different categories based on their overall geometric complexity. Thus, we compare different training strategies based on varying geometric complexity. One of the advantages of the proposed approach is the interpretation of input and output as images, which do not depend on the underlying discretization scheme. Another is the generalizability and geometric flexibility. The network can be applied to previously unseen geometries, even with different topology and level of detail. Thus, training can easily be extended to other classes of geometries.

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

Convergence of Tikhonov regularization for solving ill-posed operator equations with solutions defined on surfaces

We study Tikhonov regularization for solving ill--posed operator equations where the solutions are functions defined on surfaces. One contribution of this paper is an error analysis of Tikhonov regularization which takes into account perturbations of the surfaces, in particular when the surfaces are approximated by spline surfaces. Another contribution is that we highlight the analysis of regularization for functions with range in vector bundles over surfaces. We also present some practical applications, such as an inverse problem of gravimetry and an imaging problem for denoising vector fields on surfaces, and show the numerical verification.

Open access
preprint2015arXiv

Approximation properties of isogeometric function spaces on singularly parameterized domains

We study approximation error bounds of isogeometric function spaces on a specific type of singularly parameterized domains. In this context an isogeometric function is the composition of a piecewise rational function with the inverse of a piecewise rational geometry parameterization. We consider domains where one edge of the parameter domain is mapped onto one point in physical space. To be more precise, in our configuration the singular patch is derived from a reparameterization of a regular triangular patch. On such a domain one can define an isogeometric function space fulfilling certain regularity criteria that guarantee optimal convergence. The main contribution of this paper is to prove approximation error bounds for the previously defined class of isogeometric discretizations.

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.