BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christian Engwer

Christian Engwer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and ScienceMathematical SoftwareCell BehaviorDistributed, Parallel, and Cluster ComputingNeurons and Cognition

Catalog footprint

What is connected

11works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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)

preprint2026arXiv

Domain-of-dependence-stabilized cut-cell discretizations of linear kinetic models with summation-by-parts properties

We employ the summation-by-parts (SBP) framework to extend the recent domain-of-dependence (DoD) stabilization for cut cells to linear kinetic models in diffusion scaling. Numerical methods for these models are challenged by increased stiffness for small scaling parameters and the necessity of asymptotics preservation regarding a parabolic limit equation. As a prototype model, we consider the telegraph equation in one spatial dimension subject to periodic boundary conditions with an asymptotic limit given by the linear heat equation. We provide a general semidiscrete stability result for this model when spatially discretized by arbitrary periodic (upwind) SBP operators and formally prove that the fully discrete scheme is asymptotic preserving. Moreover, we prove that DoD with central numerical fluxes leads to periodic SBP operators. Furthermore, we show that adapting the upwind DoD scheme yields periodic upwind SBP operators. Consequently, the DoD stabilization possesses the desired properties considered in the first part of this work and thus leads to a stable and asymptotic preserving scheme for the telegraph equation. We back our theoretical results with numerical simulations and demonstrate the applicability of this cut-cell stabilization for implicit time integration in the heat equation limit.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Vectorized Parallel in Time methods for low-order discretizations with application to Porous Media problems

High order methods have shown great potential to overcome performance issues of simulations of partial differential equations (PDEs) on modern hardware, still many users stick to low-order, matrix-based simulations, in particular in porous media applications. Heterogeneous coefficients and low regularity of the solution are reasons not to employ high order discretizations. We present a new approach for the simulation of instationary PDEs that allows to partially mitigate the performance problems. By reformulating the original problem we derive a parallel in time time integrator that increases the arithmetic intensity and introduces additional structure into the problem. By this it helps accelerate matrix-based simulations on modern hardware architectures. Based on a system for multiple time steps we will formulate a matrix equation that can be solved using vectorized solvers like Block Krylov methods. The structure of this approach makes it applicable for a wide range of linear and nonlinear problems. In our numerical experiments we present some first results for three different PDEs, a linear convection-diffusion equation, a nonlinear diffusion-reaction equation and a realistic example based on the Richards' equation.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

DoD stabilization for higher-order advection in two dimensions

When solving time-dependent hyperbolic conservation laws on cut cell meshes one has to overcome the small cell problem: standard explicit time stepping is not stable on small cut cells if the time step is chosen with respect to larger background cells. The domain of dependence (DoD) stabilization is designed to solve this problem in a discontinuous Galerkin framework. It adds a penalty term to the space discretization that restores proper domains of dependency. In this contribution we introduce the DoD stabilization for solving the advection equation in 2d with higher order. We show an $L^2$ stability result for the stabilized semi-discrete scheme for arbitrary polynomial degrees $p$ and provide numerical results for convergence tests indicating orders of $p+1$ in the $L^1$ norm and between $p+\frac 1 2$ and $p+1$ in the $L^{\infty}$ norm.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

DoD Stabilization of linear hyperbolic PDEs on general cut-cell meshes

Standard numerical methods for hyperbolic PDEs require for stability a CFL-condition which implies that the time step size depends on the size of the elements of the mesh. On cut-cell meshes, elements can become arbitrarily small and thus the time step size cannot take the size of small cut-cells into account but has to be chosen based on the background mesh elements. A remedy for this is the so called DoD (domain of dependence) stabilization for which several favorable theoretical and numerical properties have been shown in one and two space dimensions. Up to now the method is restricted to stabilization of cut-cells with exactly one inflow and one outflow face, i.e. triangular cut-cells with a no-flow face. We extend the DoD stabilization to cut-cells with multiple in- and out-flow faces by properly considering the flow distribution inside the cut-cell. We further prove L2-stability for the semi-discrete formulation in space and present numerical results to validate the proposed extension.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Hardware-aware and Stable Orthogonalization Framework

The orthogonalization process is an essential building block in Krylov space methods, which takes up a large portion of the computational time. Commonly used methods, like the Gram-Schmidt method, consider the projection and normalization separately and store the orthogonal base explicitly. We consider the problem of orthogonalization and normalization as a QR decomposition problem on which we apply known algorithms, namely CholeskyQR and TSQR. This leads to methods that solve the orthogonlization problem with reduced communication costs, while maintaining stability and stores the orthogonal base in a locally orthogonal representation. Furthermore, we discuss the novel method as a framework which allows us to combine different orthogonalization algorithms and use the best algorithm for each part of the hardware. After the formulation of the methods, we show their advantageous performance properties based on a performance model that takes data transfers within compute nodes as well as message passing between compute nodes into account. The theoretic results are validated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Estimating the extent of glioblastoma invasion

Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection-diffusion-reaction equations. The patient-wise parametrisation of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are generally missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalised approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalisation is introduced via an analytical solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The DUNE Framework: Basic Concepts and Recent Developments

This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A discontinuous Galerkin Method for the EEG Forward Problem using the Subtraction Approach

In order to perform electroencephalography (EEG) source reconstruction, i.e., to localize the sources underlying a measured EEG, the electric potential distribution at the electrodes generated by a dipolar current source in the brain has to be simulated, which is the so-called EEG forward problem. To solve it accurately, it is necessary to apply numerical methods that are able to take the individual geometry and conductivity distribution of the subject's head into account. In this context, the finite element method (FEM) has shown high numerical accuracy with the possibility to model complex geometries and conductive features, e.g., white matter conductivity anisotropy. In this article, we introduce and analyze the application of a discontinuous Galerkin (DG) method, a finite element method that includes features of the finite volume framework, to the EEG forward problem. The DG-FEM approach fulfills the conservation property of electric charge also in the discrete case, making it attractive for a variety of applications. Furthermore, as we show, this approach can alleviate modeling inaccuracies that might occur in head geometries when using classical FE methods, e.g., so-called "skull leakage effects", which may occur in areas where the thickness of the skull is in the range of the mesh resolution. Therefore, we derive a DG formulation of the FEM subtraction approach for the EEG forward problem and present numerical results that highlight the advantageous features and the potential benefits of the proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces

Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Geometric Integration Over Irregular Domains with topologic Guarantees

Implicitly described domains are a well established tool in the simulation of time dependent problems, e.g. using level-set methods. In order to solve partial differential equations on such domains, a range of numerical methods was developed, e.g. the Immersed Boundary method, Unfitted Finite Element or Unfitted discontinuous Galerkin methods, eXtended or Generalised Finite Element methods, just to name a few. Many of these methods involve integration over cut-cells or their boundaries, as they are described by sub-domains of the original level-set mesh. We present a new algorithm to geometrically evaluate the integrals over domains described by a first-order, conforming level-set function. The integration is based on a polyhedral reconstruction of the implicit geometry, following the concepts of the Marching Cubes algorithm. The algorithm preserves various topological properties of the implicit geometry in its polyhedral reconstruction, making it suitable for Finite Element computations. Numerical experiments show second order accuracy of the integration. An implementation of the algorithm is available as free software, which allows for an easy incorporation into other projects. The software is in productive use within the DUNE framework.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A numerically stable a posteriori error estimator for reduced basis approximations of elliptic equations

The Reduced Basis (RB) method is a well established method for the model order reduction of problems formulated as parametrized partial differential equations. One crucial requirement for the application of RB schemes is the availability of an a posteriori error estimator to reliably estimate the error introduced by the reduction process. However, straightforward implementations of standard residual based estimators show poor numerical stability, rendering them unusable if high accuracy is required. In this work we propose a new algorithm based on representing the residual with respect to a dedicated orthonormal basis, which is both easy to implement and requires little additional computational overhead. A numerical example is given to demonstrate the performance of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access