Researcher profile

Matthias Kirchhart

Matthias Kirchhart contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
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

5 published item(s)

preprint2023arXiv

The Numerical Flow Iteration for the Vlasov-Poisson equation

We present the numerical flow iteration (NuFI) for solving the Vlasov--Poisson equation. In a certain sense specified later herein, NuFI provides infinite resolution of the distribution function. NuFI exactly preserves positivity, all $L^p$-norms, charge, and entropy. Numerical experiments show no energy drift. NuFI is fast, requires several orders of magnitude less memory than conventional approaches, and can very efficiently be parallelised on GPU clusters. Low fidelity simulations provide good qualitative results for extended periods of time and can be computed on low-cost workstations.

preprint2022arXiv

Solidfmm: A highly optimised library of operations on the solid harmonics for use in fast multipole methods

We present solidfmm, a highly optimised C++ library for the solid harmonics as they are needed in fast multipole methods. The library provides efficient, vectorised implementations of the translation operations M2M, M2L, and L2L, and is available as free software. While asymptotically of complexity $O(P^3)$, for all practically relevant expansion orders, the translation operators display an empirical complexity of $O(P^2)$, outperforming the naïve implementation by orders of magnitude.

preprint2021arXiv

Div-Curl Problems and $\mathbf{H}^1$-regular Stream Functions in 3D Lipschitz Domains

We consider the problem of recovering the divergence-free velocity field ${\mathbf U}\in\mathbf{L}^2(Ω)$ of a given vorticity ${\mathbf F}=\mathrm{curl}\,{\mathbf U}$ on a bounded Lipschitz domain $Ω\subset\mathbb{R}^3$. To that end, we solve the "div-curl problem" for a given ${\mathbf F}\in{\mathbf H}^{-1}(Ω)$. The solution is expressed in terms of a vector potential (or stream function) ${\mathbf A}\in{\mathbf H}^1(Ω)$ such that ${\mathbf U}=\mathrm{curl}\,{\mathbf A}$. After discussing existence and uniqueness of solutions and associated vector potentials, we propose a well-posed construction for the stream function. A numerical method based on this construction is presented, and experiments confirm that the resulting approximations display higher regularity than those of another common approach.

preprint2020arXiv

Analytic Integration of the Newton Potential over Cuboids and an Application to Fast Multipole Methods

We present simplified formulae for the analytic integration of the Newton potential of polynomials over boxes in two- and three-dimensional space. These are implemented in an easy-to-use C++ library that allows computations in arbitrary precision arithmetic which is also documented here. We describe how these results can be combined with fast multipole methods for general, non-polynomial data.

preprint2019arXiv

A Particle Method without Remeshing

We propose a simple tweak to a recently developed regularisation scheme for particle methods. This allows us to chose the particle spacing $h$ proportional to the regularisation length $σ$ and achieve optimal error bounds of the form $\mathcal{O}(σ^n)$, $n\in\mathbb{N}$, without any need of remeshing. We prove this result for the linear advection equation but also carry out high-order experiments on the full Navier--Stokes equations. In our experiments the particle methods proved to be highly accurate, long-term stable, and competitive with discontinuous Galerkin methods.