BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christian Engwer

Christian Engwer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisComputational Engineering, Finance, and ScienceMathematical SoftwareCell BehaviorDistributed, Parallel, and Cluster Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
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

7 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