Source author record

Arnold Reusken

Arnold Reusken 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

11works
3topics
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

11 published item(s)

preprint2024arXiv

Analysis of the Taylor-Hood Surface Finite Element Method for the surface Stokes equation

We consider the surface Stokes equation on a smooth closed hypersurface in three-dimensional space. For discretization of this problem a generalization of the surface finite element method (SFEM) of Dziuk-Elliott combined with a Hood-Taylor pair of finite element spaces has been used in the literature. We call this method Hood-Taylor-SFEM. This method uses a penalty technique to weakly satisfy the tangentiality constraint. In this paper we present a discretization error analysis of this method resulting in optimal discretization error bounds in an energy norm. We also address linear algebra aspects related to (pre)conditioning of the system matrix.

preprint2020arXiv

Error analysis of higher order trace finite element methods for the surface Stokes equations

The paper studies a higher order unfitted finite element method for the Stokes system posed on a surface in three-dimensional space. The method employs generalized Taylor-Hood finite element pairs on tetrahedral bulk mesh to discretize the Stokes system on embedded surface. Stability and optimal order convergence results are proved. The proofs include a complete quantification of geometric errors stemming from approximate parametric representation of the surface. Numerical experiments include formal convergence studies and an example of the Kelvin-Helmholtz instability problem on the unit sphere.

preprint2020arXiv

Inf-sup stability of the trace P2-P1 Taylor-Hood elements for surface PDEs

The paper studies a geometrically unfitted finite element method (FEM), known as trace FEM or cut FEM, for the numerical solution of the Stokes system posed on a closed smooth surface. A trace FEM based on standard Taylor-Hood (continuous P2-P1) bulk elements is proposed. A so-called volume normal derivative stabilization, known from the literature on trace FEM, is an essential ingredient of this method. The key result proved in the paper is an inf-sup stability of the trace P2-P1 finite element pair, with the stability constant uniformly bounded with respect to the discretization parameter and the position of the surface in the bulk mesh. Optimal order convergence of a consistent variant of the finite element method follows from this new stability result and interpolation properties of the trace FEM. Properties of the method are illustrated with numerical examples.

preprint2014arXiv

A space-time FEM for PDEs on evolving surfaces

The paper studies a finite element method for computing transport and diffusion along evolving surfaces. The method does not require a parametrization of a surface or an extension of a PDE from a surface into a bulk outer domain. The surface and its evolution may be given implicitly, e.g., as the solution of a level set equation. This approach naturally allows a surface to undergo topological changes and experience local geometric singularities. The numerical method uses space-time finite elements and is provably second order accurate. The paper reviews the method, error estimates and shows results for computing the diffusion of a surfactant on surfaces of two colliding droplets.

preprint2014arXiv

A trace finite element method for a class of coupled bulk-interface transport problems

In this paper we study a system of advection-diffusion equations in a bulk domain coupled to an advection-diffusion equation on an embedded surface. Such systems of coupled partial differential equations arise in, for example, the modeling of transport and diffusion of surfactants in two-phase flows. The model considered here accounts for adsorption-desorption of the surfactants at a sharp interface between two fluids and their transport and diffusion in both fluid phases and along the interface. The paper gives a well-posedness analysis for the system of bulk-surface equations and introduces a finite element method for its numerical solution. The finite element method is unfitted, i.e., the mesh is not aligned to the interface. The method is based on taking traces of a standard finite element space both on the bulk domains and the embedded surface. The numerical approach allows an implicit definition of the surface as the zero level of a level-set function. Optimal order error estimates are proved for the finite element method both in the bulk-surface energy norm and the $L^2$-norm. The analysis is not restricted to linear finite elements and a piecewise planar reconstruction of the surface, but also covers the discretization with higher order elements and a higher order surface reconstruction.

preprint2014arXiv

An Eulerian space-time finite element method for diffusion problems on evolving surfaces

In this paper, we study numerical methods for the solution of partial differential equations on evolving surfaces. The evolving hypersurface in $\Bbb{R}^d$ defines a $d$-dimensional space-time manifold in the space-time continuum $\Bbb{R}^{d+1}$. We derive and analyze a variational formulation for a class of diffusion problems on the space-time manifold. For this variational formulation new well-posedness and stability results are derived. The analysis is based on an inf-sup condition and involves some natural, but non-standard, (anisotropic) function spaces. Based on this formulation a discrete in time variational formulation is introduced that is very suitable as a starting point for a discontinuous Galerkin (DG) space-time finite element discretization. This DG space-time method is explained and results of numerical experiments are presented that illustrate its properties.

preprint2014arXiv

Error analysis of a space-time finite element method for solving PDEs on evolving surfaces

In this paper we present an error analysis of an Eulerian finite element method for solving parabolic partial differential equations posed on evolving hypersurfaces in $\mathbb{R}^d$, $d=2,3$. The method employs discontinuous piecewise linear in time -- continuous piecewise linear in space finite elements and is based on a space-time weak formulation of a surface PDE problem. Trial and test surface finite element spaces consist of traces of standard volumetric elements on a space-time manifold resulting from the evolution of a surface. We prove first order convergence in space and time of the method in an energy norm and second order convergence in a weaker norm. Furthermore, we derive regularity results for solutions of parabolic PDEs on an evolving surface, which we need in a duality argument used in the proof of the second order convergence estimate.

preprint2014arXiv

Optimal preconditioners for Nitsche-XFEM discretizations of interface problems

In the past decade, a combination of unfitted finite elements (or XFEM) with the Nitsche method has become a popular discretization method for elliptic interface problems. This development started with the introduction and analysis of this Nitsche-XFEM technique in the paper [A. Hansbo, P. Hansbo, Comput. Methods Appl. Mech. Engrg. 191 (2002)]. In general, the resulting linear systems have very large condition numbers, which depend not only on the mesh size $h$, but also on how the interface intersects the mesh. This paper is concerned with the design and analysis of optimal preconditioners for such linear systems. We propose an additive subspace preconditioner which is optimal in the sense that the resulting condition number is independent of the mesh size $h$ and the interface position. We further show that already the simple diagonal scaling of the stifness matrix results in a condition number that is bounded by $ch^{-2}$, with a constant $c$ that does not depend on the location of the interface. Both results are proven for the two-dimensional case. Results of numerical experiments in two and three dimensions are presented, which illustrate the quality of the preconditioner.

preprint2013arXiv

A stabilized finite element method for advection-diffusion equations on surfaces

A recently developed Eulerian finite element method is applied to solve advection-diffusion equations posed on hypersurfaces. When transport processes on a surface dominate over diffusion, finite element methods tend to be unstable unless the mesh is sufficiently fine. The paper introduces a stabilized finite element formulation based on the SUPG technique. An error analysis of the method is given. Results of numerical experiments are presented that illustrate the performance of the stabilized method.

preprint2013arXiv

On surface meshes induced by level set functions

The zero level set of a piecewise-affine function with respect to a consistent tetrahedral subdivision of a domain in $\mathbb{R}^3$ is a piecewise-planar hyper-surface. We prove that if a family of consistent tetrahedral subdivions satisfies the minimum angle condition, then after a simple postprocessing this zero level set becomes a consistent surface triangulation which satisfies the maximum angle condition. We treat an application of this result to the numerical solution of PDEs posed on surfaces, using a $P_1$ finite element space on such a surface triangulation. For this finite element space we derive optimal interpolation error bounds. We prove that the diagonally scaled mass matrix is well-conditioned, uniformly with respect to $h$. Furthermore, the issue of conditioning of the stiffness matrix is addressed.

preprint2012arXiv

Reduced Basis A Posteriori Error Bounds for the Instationary Stokes Equations

We present reduced basis approximations and rigorous a posteriori error bounds for the instationary Stokes equations. We shall discuss both a method based on the standard formulation as well as a method based on a penalty approach, which combine techniques developed in our previous work on parametrized saddle point problems with current reduced basis techniques for parabolic problems. The analysis then shows how time integration affects the development of reduced basis a posteriori error bounds as well as the construction of computationally efficient reduced basis approximation spaces. To demonstrate their performance in practice, the methods are applied to a Stokes flow in a two-dimensional microchannel with a parametrized rectangular obstacle; evolution in time is induced by a time-dependent velocity profile on the inflow boundary. Numerical results illustrate (i) the rapid convergence of reduced basis approximations, (ii) the performance of a posteriori error bounds with respect to sharpness, and (iii) computational efficiency.