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

Hendrik Ranocha

Hendrik Ranocha contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phphysics.flu-dynDistributed, Parallel, and Cluster ComputingMathematical SoftwarePerformance

Trust snapshot

Quick read

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

13 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

GPU-Accelerated Energy-Conserving Methods for the Hyperbolized Serre-Green-Naghdi Equations in 2D

We develop energy-conserving numerical methods for a two-dimensional hyperbolic approximation of the Serre-Green-Naghdi equations with variable bathymetry for both periodic and reflecting boundary conditions. The hyperbolic formulation avoids the costly inversion of an elliptic operator present in the classical model. Our schemes combine split forms with summation-by-parts (SBP) operators to construct semidiscretizations that conserve the total water mass and the total energy. We provide analytical proofs of these conservation properties and also verify them numerically. While the framework is general, our implementation focuses on second-order finite-difference SBP operators. The methods are implemented in Julia for CPU and GPU architectures (AMD and NVIDIA) and achieve substantial speedups on modern accelerators. We validate the approach through convergence studies based on solitary-wave and manufactured-solution tests, and by comparisons to analytical, experimental, and existing numerical results. All source code to reproduce our results is available online.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Note on Numerical Fluxes Conserving a Member of Harten's One-Parameter Family of Entropies for the Compressible Euler Equations

Entropy-conserving numerical fluxes are a cornerstone of modern high-order entropy-dissipative discretizations of conservation laws. In addition to entropy conservation, other structural properties mimicking the continuous level such as pressure equilibrium and kinetic energy preservation are important. This note proves that there are no numerical fluxes conserving (one of) Harten's entropies for the compressible Euler equations that also preserve pressure equilibria and have a density flux independent of the pressure. This is in contrast to fluxes based on the physical entropy, where even kinetic energy preservation can be achieved in addition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing

We present Trixi.jl, a Julia package for adaptive high-order numerical simulations of hyperbolic partial differential equations. Utilizing Julia's strengths, Trixi.jl is extensible, easy to use, and fast. We describe the main design choices that enable these features and compare Trixi.jl with a mature open source Fortran code that uses the same numerical methods. We conclude with an assessment of Julia for simulation-focused scientific computing, an area that is still dominated by traditional high-performance computing languages such as C, C++, and Fortran.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the entropy projection and the robustness of high order entropy stable discontinuous Galerkin schemes for under-resolved flows

High order entropy stable schemes provide improved robustness for computational simulations of fluid flows. However, additional stabilization and positivity preserving limiting can still be required for variable-density flows with under-resolved features. We demonstrate numerically that entropy stable DG methods which incorporate an "entropy projection" are less likely to require additional limiting to retain positivity for certain types of flows. We conclude by investigating potential explanations for this observed improvement in robustness.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reinterpretation and Extension of Entropy Correction Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization

For the general class of residual distribution (RD) schemes, including many finite element (such as continuous/discontinuous Galerkin) and flux reconstruction methods, an approach to construct entropy conservative/ dissipative semidiscretizations by adding suitable correction terms has been proposed by Abgrall (J.~Comp.~Phys. 372: pp. 640--666, 2018). In this work, the correction terms are characterized as solutions of certain optimization problems and are adapted to the SBP-SAT framework, focusing on discontinuous Galerkin methods. Novel generalizations to entropy inequalities, multiple constraints, and kinetic energy preservation for the Euler equations are developed and tested in numerical experiments. For all of these optimization problems, explicit solutions are provided. Additionally, the correction approach is applied for the first time to obtain a fully discrete entropy conservative/dissipative RD scheme. Here, the application of the deferred correction (DeC) method for the time integration is essential. This paper can be seen as describing a systematic method to construct structure preserving discretization, at least for the considered example.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Broad Class of Conservative Numerical Methods for Dispersive Wave Equations

We develop a general framework for designing conservative numerical methods based on summation by parts operators and split forms in space, combined with relaxation Runge-Kutta methods in time. We apply this framework to create new classes of fully-discrete conservative methods for several nonlinear dispersive wave equations: Benjamin-Bona-Mahony (BBM), Fornberg-Whitham, Camassa-Holm, Degasperis-Procesi, Holm-Hone, and the BBM-BBM system. These full discretizations conserve all linear invariants and one nonlinear invariant for each system. The spatial semidiscretizations include finite difference, spectral collocation, and both discontinuous and continuous finite element methods. The time discretization is essentially explicit, using relaxation Runge-Kutta methods. We implement some specific schemes from among the derived classes, and demonstrate their favorable properties through numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fully-Discrete Explicit Locally Entropy-Stable Schemes for the Compressible Euler and Navier-Stokes Equations

Recently, relaxation methods have been developed to guarantee the preservation of a single global functional of the solution of an ordinary differential equation. Here, we generalize this approach to guarantee local entropy inequalities for finitely many convex functionals (entropies) and apply the resulting methods to the compressible Euler and Navier-Stokes equations. Based on the unstructured $hp$-adaptive SSDC framework of entropy conservative or dissipative semidiscretizations using summation-by-parts and simultaneous-approximation-term operators, we develop the first discretizations for compressible computational fluid dynamics that are primary conservative, locally entropy stable in the fully discrete sense under a usual CFL condition, explicit except for the parallelizable solution of a single scalar equation per element, and arbitrarily high-order accurate in space and time. We demonstrate the accuracy and the robustness of the fully-discrete explicit locally entropy-stable solver for a set of test cases of increasing complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Relaxation Runge-Kutta Methods for Hamiltonian Problems

The recently-introduced relaxation approach for Runge-Kutta methods can be used to enforce conservation of energy in the integration of Hamiltonian systems. We study the behavior of implicit and explicit relaxation Runge-Kutta methods in this context. We find that, in addition to their useful conservation property, the relaxation methods yield other improvements. Experiments show that their solutions bear stronger qualitative similarity to the true solution and that the error grows more slowly in time. We also prove that these methods are superconvergent for a certain class of Hamiltonian systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards Green Computing: A Survey of Performance and Energy Efficiency of Different Platforms using OpenCL

When considering different hardware platforms, not just the time-to-solution can be of importance but also the energy necessary to reach it. This is not only the case with battery powered and mobile devices but also with high-performance parallel cluster systems due to financial and practical limits on power consumption and cooling. Recent developments in hard- and software have given programmers the ability to run the same code on a range of different devices giving rise to the concept of heterogeneous computing. Many of these devices are optimized for certain types of applications. To showcase the differences and give a basic outlook on the applicability of different architectures for specific problems, the cross-platform OpenCL framework was used to compare both time- and energy-to-solution. A large set of devices ranging from ARM processors to server CPUs and consumer and enterprise level GPUs has been used with different benchmarking testcases taken from applied research applications. While the results show the overall advantages of GPUs in terms of both runtime and energy efficiency compared to CPUs, ARM devices show potential for certain applications in massively parallel systems. This study also highlights how OpenCL enables the use of the same codebase on many different systems and hardware platforms without specific code adaptations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Analysis of Artificial Dissipation of Explicit and Implicit Time-Integration Methods

Stability is an important aspect of numerical methods for hyperbolic conservation laws and has received much interest. However, continuity in time is often assumed and only semidiscrete stability is studied. Thus, it is interesting to investigate the influence of explicit and implicit time integration methods on the stability of numerical schemes. If an explicit time integration method is applied, spacially stable numerical schemes for hyperbolic conservation laws can result in unstable fully discrete schemes. Focusing on the explicit Euler method (and convex combinations thereof), undesired terms in the energy balance trigger this phenomenon and introduce an erroneous growth of the energy over time. In this work, we study the influence of artificial dissipation and modal filtering in the context of discontinuous spectral element methods to remedy these issues. In particular, lower bounds on the strength of both artificial dissipation and modal filtering operators are given and an adaptive procedure to conserve the (discrete) $\mathbf{L}_2$ norm of the numerical solution in time is derived. This might be beneficial in regions where the solution is smooth and for long time simulations. Moreover, this approach is used to study the connections between explicit and implicit time integration methods and the associated energy production. By adjusting the adaptive procedure, we demonstrate that filtering in explicit time integration methods is able to mimic the dissipative behavior inherent in implicit time integration methods. This contribution leads to a better understanding of existing algorithms and numerical techniques, in particular the application of artificial dissipation as well as modal filtering in the context of numerical methods for hyperbolic conservation laws together with the selection of explicit or implicit time integration methods.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Discrete Vector Calculus and Helmholtz Hodge Decomposition for Classical Finite Difference Summation by Parts Operators

In this article, discrete variants of several results from vector calculus are studied for classical finite difference summation by parts operators in two and three space dimensions. It is shown that existence theorems for scalar/vector potentials of irrotational/solenoidal vector fields cannot hold discretely because of grid oscillations, which are characterised explicitly. This results in a non-vanishing remainder associated to grid oscillations in the discrete Helmholtz Hodge decomposition. Nevertheless, iterative numerical methods based on an interpretation of the Helmholtz Hodge decomposition via orthogonal projections are proposed and applied successfully. In numerical experiments, the discrete remainder vanishes and the potentials converge with the same order of accuracy as usual in other first order partial differential equations. Motivated by the successful application of the Helmholtz Hodge decomposition in theoretical plasma physics, applications to the discrete analysis of magnetohydrodynamic (MHD) wave modes are presented and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On Strong Stability of Explicit Runge-Kutta Methods for Nonlinear Semibounded Operators

Explicit Runge-Kutta methods are classical and widespread techniques in the numerical solution of ordinary differential equations (ODEs). Considering partial differential equations, spatial semidiscretisations can be used to obtain systems of ODEs that are solved subsequently, resulting in fully discrete schemes. However, certain stability investigations of high-order methods for hyperbolic conservation laws are often conducted only for the semidiscrete versions. Here, strong stability (also known as monotonicity) of explicit Runge-Kutta methods for ODEs with nonlinear and semibounded (also known as dissipative) operators is investigated. Contrary to the linear case, it is proven that many strong stability preserving (SSP) schemes of order two or greater are not strongly stable for general smooth and semibounded nonlinear operators. Additionally, it is shown that there are first order accurate explicit SSP Runge-Kutta methods that are strongly stable (monotone) for semibounded (dissipative) and Lipschitz continuous operators.

0 citations0 reviewsNo code linkedOpen access