Source author record

Joachim Schöberl

Joachim Schöberl 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

16works
6topics
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

16 published item(s)

preprint2026arXiv

Generalizing Riemann curvature to Regge metrics

In this paper, we propose a generalization of the Riemann curvature tensor on manifolds (of dimension two or higher) endowed with a Regge metric. Specifically, while all components of the metric tensor are assumed to be smooth within elements of a triangulation of the manifold, they need not be smooth across element interfaces, where only continuity of the tangential components is assumed. While linear derivatives of the metric can be generalized to Schwartz distributions, similarly generalizing the classical Riemann curvature tensor, a nonlinear second-order derivative of the metric, requires more care. We propose a generalization that combines the classical angle defect and jumps in the second fundamental form across element interfaces, and argue its correctness. Specifically, if a piecewise smooth metric approximates a globally smooth metric, then our generalized Riemann curvature tensor approximates the classical Riemann curvature tensor associated with the latter. Moreover, we show that if the metric approximation converges at some rate in a mesh-dependent norm equivalent to the $L^2$ norm, then the curvature approximation converges in the negative Sobolev space $H^{-2}$, the dual space of $H^2_0$, at the same rate, under additional assumptions. By appropriate contractions of the generalized Riemann curvature tensor, this work also provides generalizations of scalar curvature, the Ricci curvature tensor, and the Einstein tensor in any dimension.

preprint2023arXiv

High-order projection-based upwind method for implicit large eddy simulation

We assess the ability of three different approaches based on high-order discontinuous Galerkin methods to simulate under-resolved turbulent flows. The capabilities of the mass conserving mixed stress method as structure resolving large eddy simulation solver are examined. A comparison of a variational multiscale model to no-model or an implicit model approach is presented via numerical results. In addition, we present a novel approach for turbulent modeling in wall-bounded flows. This new technique provides a more accurate representation of the actual subgrid scales in the near wall region and gives promising results for highly under-resolved flow problems. In this paper, the turbulent channel flow and periodic hill flow problem are considered as benchmarks for our simulations.

preprint2022arXiv

A conforming auxiliary space preconditioner for the mass conserving mixed stress method

We are studying the efficient solution of the system of linear equation stemming from the mass conserving mixed stress (MCS) method discretization of the Stokes equations. To that end we perform static condensation to arrive at a system for the pressure and velocity unknowns. An auxiliary space preconditioner for the positive definite velocity block makes use of efficient and scalable solvers for conforming Finite Element spaces of low order and is analyzed with emphasis placed on the polynomial degree of the discretization. Numerical experiments demonstrate the potential of this approach and the efficiency of the implementation.

preprint2022arXiv

Hybridized Discontinuous Galerkin Methods for a Multiple Network Poroelasticity Model with Medical Applications

The quasi-static multiple network poroelastic theory (MPET) model, first introduced in the context of geomechanics, has recently found new applications in medicine. In practice, the parameters in the MPET equations can vary over several orders of magnitude which makes their stable discretization and fast solution a challenging task. Here, a new efficient parameter-robust hybridized discontinuous Galerkin method, which also features fluid mass conservation, is proposed for the MPET model. Its stability analysis which is crucial for the well-posedness of the discrete problem is performed and cost-efficient fast parameter-robust preconditioners are derived. We present a series of numerical computations for a 4-network MPET model of a human brain which support the performance of the new algorithms.

preprint2022arXiv

Primal and mixed finite element formulations for the relaxed micromorphic model

The classical Cauchy continuum theory is suitable to model highly homogeneous materials. However, many materials, such as porous media or metamaterials, exhibit a pronounced microstructure. As a result, the classical continuum theory cannot capture their mechanical behaviour without fully resolving the underlying microstructure. In terms of finite element computations, this can be done by modelling the entire body, including every interior cell. The relaxed micromorphic continuum offers an alternative method by instead enriching the kinematics of the mathematical model. The theory introduces a microdistortion field, encompassing nine extra degrees of freedom for each material point. The corresponding elastic energy functional contains the gradient of the displacement field, the microdistortion field and its Curl (the micro-dislocation). Therefore, the natural spaces of the fields are $[\mathit{H}^1]^3$ for the displacement and $[\mathit{H}(\mathrm{curl})]^3$ for the microdistortion, leading to unusual finite element formulations. In this work we describe the construction of appropriate finite elements using Nédélec and Raviart-Thomas subspaces, encompassing solutions to the orientation problem and the discrete consistent coupling condition. Further, we explore the numerical behaviour of the relaxed micromorphic model for both a primal and a mixed formulation. The focus of our benchmarks lies in the influence of the characteristic length $L_\mathrm{c}$ and the correlation to the classical Cauchy continuum theory.

preprint2020arXiv

A matrix-free Discontinuous Galerkin method for the time dependent Maxwell equations in unbounded domains

A Discontinuous Galerkin (DG) Finite Element Method (FEM) approach for the 3D time dependent Maxwell equations in unbounded domains is presented. The method (implemented in the FEM library NGsolve) is based on the covariant transformation of a modal orthogonal polynomial basis, originally defined on a reference simplex. The approach leads to an explicit time stepping scheme for which the mass matrix to be inverted is at most $d\!\times\!d$ block-diagonal (in $d\!=\!2,3$ spatial dimensions) while the matrix which discretizes the curl operators on the right-hand side of the system is a small reference matrix, independent from geometric properties of mesh elements. Furthermore, we show that the introduced optimizations are preserved when unbounded domains are also included in the formulation through a complex-stretching based approach.

preprint2020arXiv

A Reduced Basis Method For Fractional Diffusion Operators II

We present a novel numerical scheme to approximate the solution map $s\mapsto u(s) := \mathcal{L}^{-s}f$ to partial differential equations involving fractional elliptic operators. Reinterpreting $\mathcal{L}^{-s}$ as interpolation operator allows us to derive an integral representation of $u(s)$ which includes solutions to parametrized reaction-diffusion problems. We propose a reduced basis strategy on top of a finite element method to approximate its integrand. Unlike prior works, we deduce the choice of snapshots for the reduced basis procedure analytically. Avoiding further discretization, the integral is interpreted in a spectral setting to evaluate the surrogate directly. Its computation boils down to a matrix approximation $L$ of the operator whose inverse is projected to a low-dimensional space, where explicit diagonalization is feasible. The universal character of the underlying $s$-independent reduced space allows the approximation of $(u(s))_{s\in(0,1)}$ in its entirety. We prove exponential convergence rates and confirm the analysis with a variety of numerical examples. Further improvements are proposed in the second part of this investigation to avoid inversion of $L$. Instead, we directly project the matrix to the reduced space, where its negative fractional power is evaluated. A numerical comparison with the predecessor highlights its competitive performance.

preprint2020arXiv

Structure aware Runge-Kutta time stepping for spacetime tents

We introduce a new class of Runge-Kutta type methods suitable for time stepping to propagate hyperbolic solutions within tent-shaped spacetime regions. Unlike standard Runge-Kutta methods, the new methods yield expected convergence properties when standard high order spatial (discontinuous Galerkin) discretizations are used. After presenting a derivation of nonstandard order conditions for these methods, we show numerical examples of nonlinear hyperbolic systems to demonstrate the optimal convergence rates. We also report on the discrete stability properties of these methods applied to linear hyperbolic equations.

preprint2019arXiv

Divergence-free tangential finite element methods for incompressible flows on surfaces

In this work we consider the numerical solution of incompressible flows on two-dimensional manifolds. Whereas the compatibility demands of the velocity and the pressure spaces are known from the flat case one further has to deal with the approximation of a velocity field that lies only in the tangential space of the given geometry. Abandoning $H^1$-conformity allows us to construct finite elements which are -- due to an application of the Piola transformation -- exactly tangential. To reintroduce continuity (in a weak sense) we make use of (hybrid) discontinuous Galerkin techniques. To further improve this approach, $H(\operatorname{div}_Γ)$-conforming finite elements can be used to obtain exactly divergence-free velocity solutions. We present several new finite element discretizations. On a number of numerical examples we examine and compare their qualitative properties and accuracy.

preprint2019arXiv

The Hellan-Herrmann-Johnson Method for Nonlinear Shells

In this paper we derive a new finite element method for nonlinear shells. The Hellan-Herrmann-Johnson (HHJ) method is a mixed finite element method for fourth order Kirchhoff plates. It uses convenient Lagrangian finite elements for the vertical deflection, and introduces sophisticated finite elements for the moment tensor. In this work we present a generalization of this method to nonlinear shells, where we allow finite strains and large rotations. The geometric interpretation of degrees of freedom allows a straight forward discretization of structures with kinks. The performance of the proposed elements is demonstrated by means of several established benchmark examples.

preprint2018arXiv

Computational micromagnetics with Commics

We present our open-source Python module Commics for the study of the magnetization dynamics in ferromagnetic materials via micromagnetic simulations. It implements state-of-the-art unconditionally convergent finite element methods for the numerical integration of the Landau-Lifshitz-Gilbert equation. The implementation is based on the multiphysics finite element software Netgen/NGSolve. The simulation scripts are written in Python, which leads to very readable code and direct access to extensive post-processing. Together with documentation and example scripts, the code is freely available on GitLab.

preprint2016arXiv

Divergence-free Reconstruction Operators for Pressure-Robust Stokes Discretizations With Continuous Pressure Finite Elements

Classical inf-sup stable mixed finite elements for the incompressible (Navier-)Stokes equations are not pressure-robust, i.e., their velocity errors depend on the continuous pressure. However, a modification only in the right hand side of a Stokes discretization is able to reestablish pressure-robustness, as shown recently for several inf-sup stable Stokes elements with discontinuous discrete pressures. In this contribution, this idea is extended to low and high order Taylor-Hood and mini elements, which have continuous discrete pressures. For the modification of the right hand side a velocity reconstruction operator is constructed that maps discretely divergence-free test functions to exactly divergence-free ones. The reconstruction is based on local $H(\mathrm{div})$-conforming flux equilibration on vertex patches, and fulfills certain orthogonality properties to provide consistency and optimal a-priori error estimates. Numerical examples for the incompressible Stokes and Navier-Stokes equations confirm that the new pressure-robust Taylor-Hood and mini elements converge with optimal order and outperform significantly the classical versions of those elements when the continuous pressure is comparably large.

preprint2016arXiv

High order exactly divergence-free Hybrid Discontinuous Galerkin Methods for unsteady incompressible flows

In this paper we present an efficient discretization method for the solution of the unsteady incompressible Navier-Stokes equations based on a high order (Hybrid) Discontinuous Galerkin formulation. The crucial component for the efficiency of the discretization method is the disctinction between stiff linear parts and less stiff non-linear parts with respect to their temporal and spatial treatment. Exploiting the flexibility of operator-splitting time integration schemes we combine two spatial discretizations which are tailored for two simpler sub-problems: a corresponding hyperbolic transport problem and an unsteady Stokes problem. For the hyperbolic transport problem a spatial discretization with an Upwind Discontinuous Galerkin method and an explicit treatment in the time integration scheme is rather natural and allows for an efficient implementation. The treatment of the Stokes part involves the solution of linear systems. In this case a discretization with Hybrid Discontinuous Galerkin methods is better suited. We consider such a discretization for the Stokes part with two important features: H(div)-conforming finite elements to garantuee exactly divergence-free velocity solutions and a projection operator which reduces the number of globally coupled unknowns. We present the method, discuss implementational aspects and demonstrate the performance on two and three dimensional benchmark problems.

preprint2016arXiv

Polynomial robust stability analysis for $H(\textrm{div})$-conforming finite elements for the Stokes equations

In this work we consider a discontinuous Galerkin method for the discretization of the Stokes problem. We use $H(\textrm{div})$-conforming finite elements as they provide major benefits such as exact mass conservation and pressure-independent error estimates. The main aspect of this work lies in the analysis of high order approximations. We show that the considered method is uniformly stable with respect to the polynomial order $k$ and provides optimal error estimates $ \| \boldsymbol{u} - \boldsymbol{u}_h \|_{1_h} + \| Π^{Q_h}p-p_h \| \le c \left( h/k \right)^s \| \boldsymbol{u} \|_{s+1} $. To derive those estimates, we prove a $k$-robust LBB condition. This proof is based on a polynomial $H^2$-stable extension operator. This extension operator itself is of interest for the numerical analysis of $C^0$-continuous discontinuous Galerkin methods for $4^{th}$ order problems.

preprint2012arXiv

A uniform inf--sup condition with applications to preconditioning

A uniform inf-sup condition related to a parameter dependent Stokes problem is established. Such conditions are intimately connected to the construction of uniform preconditioners for the problem, i.e., preconditioners which behave uniformly well with respect to variations in the model parameter as well as the discretization parameter. For the present model, similar results have been derived before, but only by utilizing extra regularity ensured by convexity of the domain. The purpose of this paper is to remove this artificial assumption. As a byproduct of our analysis, in the two dimensional case we also construct a new projection operator for the Taylor-Hood element which is uniformly bounded in $L^2$ and commutes with the divergence operator. This construction is based on a tight connection between a subspace of the Taylor-Hood velocity space and the lowest order Nedelec edge element.

preprint2011arXiv

High order Curl-conforming Hardy space infinite elements for exterior Maxwell problems

A construction of prismatic Hardy space infinite elements to discretize wave equations on unbounded domains $Ω$ in $H^1_{loc}(Ω)$, $H_{loc}(curl;Ω)$ and $H_{loc}(div;Ω)$ is presented. As our motivation is to solve Maxwell's equations we take care that these infinite elements fit into the discrete de Rham diagram, i.e. they span discrete spaces, which together with the exterior derivative form an exact sequence. Resonance as well as scattering problems are considered in the examples. Numerical tests indicate super-algebraic convergence in the number of additional unknowns per degree of freedom on the coupling boundary that are required to realize the Dirichlet to Neumann map.