Researcher profile

Joachim Schöberl

Joachim Schöberl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 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.