Source author record

Thomas Takacs

Thomas Takacs 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

7works
5topics
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

7 published item(s)

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.

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.

preprint2020arXiv

A family of $C^1$ quadrilateral finite elements

We present a novel family of $C^1$ quadrilateral finite elements, which define global $C^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\geq 6$, to all degrees $p \geq 3$. Thus, we call the family of $C^1$ finite elements Brenner-Sung quadrilaterals. The proposed $C^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 $C^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~$5$, and the polynomial degrees $p=3$ and $p=4$ by employing a splitting into $3\times3$ or $2\times2$ polynomial pieces, respectively. The proposed elements reproduce polynomials of total degree $p$. 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\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 $C^1$ quadrilateral finite element for the numerical analysis of fourth order problems, also indicating that (for $p=5$) the proposed element performs comparable or in general even better than the Argyris triangle with respect to the number of degrees of freedom.

preprint2016arXiv

Analysis-suitable $G^1$ multi-patch parametrizations for $C^1$ isogeometric spaces

One key feature of isogeometric analysis is that it allows smooth shape functions. Indeed, when isogeometric spaces are constructed from $p$-degree splines (and extensions, such as NURBS), they enjoy up to $C^{p-1}$ continuity within each patch. However, global continuity beyond $C^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 $C^0$ at the patch interface. On such domains we study the $h$-refinement of $C^1$-continuous isogeometric spaces. These spaces in general do not have optimal approximation properties. The reason is that the $C^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) $C^1$ splines. This is the starting point of our study. We introduce the class of analysis-suitable $G^1$ geometry parametrizations, which includes piecewise bilinear parametrizations. We then analyze the structure of $C^1$ isogeometric spaces over analysis-suitable $G^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 $C^1$ isogeometric spaces is prevented.

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.

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.

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 $h$-refined meshes.