Source author record

Matthias Kirchhart

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

ResearcherUnclaimed source record

math.NA Numerical Analysis math.AP Mathematical Software physics.comp-ph

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.