Researcher profile

Christian Klingenberg

Christian Klingenberg contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
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

17 published item(s)

preprint2026arXiv

BGK model for rarefied gas in a bounded domain

We study the Bathnagar-Gross-Krook (BGK) equation in a smooth bounded domain featuring a diffusive reflection boundary condition with general collision frequency. We prove that the BGK equation admits a unique global solution with an exponential convergence rate if the initial condition is a small perturbation around the global Maxwellian in the $L^\infty$ space. For the proof, we utilize the dissipative nature from the linearized BGK operator and establish an $L^2$ coercive estimate. Next, we derive the a priori estimate by obtaining an $L^\infty$ bound on the nonlinear operator; this requires a delicate analysis to manage its intrinsic nonlinear structure. Finally, we establish the $L^\infty$ stability estimate and introduce sequential arguments for the nonlinear BGK operator, thereby concluding both well-posedness and positivity.

preprint2024arXiv

Reconstructing the kinetic chemotaxis kernel using macroscopic data: well-posedness and ill-posedness

Bacterial motion is steered by external stimuli (chemotaxis), and the motion described on the mesoscopic scale is uniquely determined by a parameter $K$ that models velocity change response from the bacteria. This parameter is called chemotaxis kernel. In a practical setting, it is inferred by experimental data. We deploy a PDE-constrained optimization framework to perform this reconstruction using velocity-averaged, localized data taken in the interior of the domain. The problem can be well-posed or ill-posed depending on the data preparation and the experimental setup. In particular, we propose one specific design that guarantees numerical reconstructability and local convergence. This design is adapted to the discretization of $K$ in space and decouples the reconstruction of local values of $K$ into smaller cell problems, opening up parallelization opportunities. Numerical evidences support the theoretical findings.

preprint2022arXiv

BGK model for multi-component gases near a global Maxwellian

In this paper, we establish the existence of the unique global-in-time classical solutions to the multi-component BGK model suggested in \cite{mixmodel} when the initial data is a small perturbation of global equilibrium. For this, we carefully analyze the dissipative nature of the linearized multi-component relaxation operator, and observe that the partial dissipation from the intra-species and the inter-species linearized relaxation operators are combined in a complementary manner to give rise to the desired dissipation estimate of the model We also observe that the convergence rate of the distribution function increases as the momentum-energy interchange rate between the different components of the gas increases.

preprint2022arXiv

Exact solution and the multidimensional Godunov scheme for the acoustic equations

The acoustic equations derived as a linearization of the Euler equations are a valuable system for studies of multi-dimensional solutions. Additionally they possess a low Mach number limit analogous to that of the Euler equations. Aiming at understanding the behaviour of the multi-dimensional Godunov scheme in this limit, first the exact solution of the corresponding Cauchy problem in three spatial dimensions is derived. The appearance of logarithmic singularities in the exact solution of the 4-quadrant Riemann Problem in two dimensions is discussed. The solution formulae are then used to obtain the multidimensional Godunov finite volume scheme in two dimensions. It is shown to be superior to the dimensionally split upwind/Roe scheme concerning its domain of stability and ability to resolve multi-dimensional Riemann problems. It is shown experimentally and theoretically that despite taking into account multi-dimensional information it is, however, not able to resolve the low Mach number limit.

preprint2022arXiv

Well-balanced Central Scheme for the System of MHD Equations with Gravitational source term

A well-balanced second order finite volume central scheme for the magnetohydrodynamic (MHD) equations with gravitational source term is developed in this paper. The scheme is an unstaggered central scheme that evolves the numerical solution on a single grid and avoids solving Riemann problems at the cell interfaces using ghost staggered cells. A subtraction technique is used on the conservative variables with the support of a known steady state in order to manifest the well-balanced property of the scheme. The divergence-free constraint of the magnetic field is satisfied after applying the constrained transport method (CTM) for unstaggered central schemes at the end of each time-step by correcting the components of the magnetic field. The robustness of the proposed scheme is verified on a list of numerical test cases from the literature.

preprint2021arXiv

Euler system with a polytropic equation of state as a vanishing viscosity limit

We consider the Euler system of gas dynamics endowed with the incomplete equation of state relating the internal energy to the mass density and the pressure. We show that any sufficiently smooth solution can be recovered as a vanishing viscosity - heat conductivity limit of the Navier--Stokes--Fourier system with a properly defined temperature. The result is unconditional in the case of the Navier type (slip) boundary conditions and extends to the no-slip condition for the velocity under some extra hypotheses of Kato's type concerning the behavior of the fluid in the boundary layer.

preprint2020arXiv

BGK model of the multi-species Uehling-Uhlenbeck equation

We propose a BGK model of the quantum Boltzmann equation for gas mixtures. We also provide a sufficient condition that guarantees the existence of equilibrium coefficients so that the model shares the same conservation laws and $H$-theorem with the quantum Boltzmann equation. Unlike the classical BGK for gas mixtures, the equilibrium coefficinets of the local equilibiriums for quantum multi-species gases are defined through highly nonlinear relations that are not explicitly solvable. We verify in a unified way that such nonlinear relations uniquely determine the equilibrium coefficients under the condition, leading to the well-definedness of our model.

preprint2020arXiv

Entropy Stable Numerical Fluxes for Compressible Euler Equations which are Suitable for All Mach Numbers

We propose two novel two-state approximate Riemann solvers for the compressible Euler equations which are provably entropy dissipative and suitable for the simulation of low Mach numbers. What is new, is that one of our two methods in addition is provably kinetic energy stable. Both methods are based on the entropy satisfying and kinetic energy consistent methods of Chandrashekar (2013). The low Mach number compliance is achieved by rescaling some speed of sound terms in the diffusion matrix in the spirit of Li & Gu (2008). In numerical tests we demonstrate the low Mach number compliance and the entropy stability of the proposed fluxes.

preprint2020arXiv

High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws

We introduce a general framework for the construction of well-balanced finite volume methods for hyperbolic balance laws. We use the phrase well-balancing in a broader sense, since our proposed method can be applied to exactly follow any solution of any system of hyperbolic balance laws in multiple spatial dimensions and not only time independent solutions. The solution has to be known a priori, either as an analytical expression or as discrete data. The proposed framework modifies the standard finite volume approach such that the well-balancing property is obtained and in case the method is high order accurate, this is maintained under our modification. We present numerical tests for the compressible Euler equations with and without gravity source term and with different equations of state, and for the equations of compressible ideal magnetohydrodynamics.

preprint2020arXiv

Multilevel Monte Carlo Finite Volume Methods for Random Conservation Laws with Discontinuous Flux

We consider conservation laws with discontinuous flux where the initial datum, the flux function, and the discontinuous spatial dependency coefficient are subject to randomness. We establish a notion of random adapted entropy solutions to these equations and prove well-posedness provided that the spatial dependency coefficient is piecewise constant with finitely many discontinuities. In particular, the setting under consideration allows the flux to change across finitely many points in space whose positions are uncertain. We propose a single- and multilevel Monte Carlo method based on a finite volume approximation for each sample. Our analysis includes convergence rate estimates of the resulting Monte Carlo and multilevel Monte Carlo finite volume methods as well as error versus work rates showing that the multilevel variant outperforms the single-level method in terms of efficiency. We present numerical experiments motivated by two-phase reservoir simulations for reservoirs with varying geological properties.

preprint2020arXiv

On oscillatory solutions to the complete Euler system

The Euler system in fluid dynamics is a model of a compressible inviscid fluid incorporating the three basic physical principles: Conservation of mass, momentum, and energy. We show that the Cauchy problem is basically ill-posed for the $L^\infty$-initial data in the class of weak entropy solutions. As a consequence, there are infinitely many measure-valued solutions for a vast set of initial data. Finally, using the concept of relative energy, we discuss a singular limit problem for the measure-valued solutions, where the Mach and Froude number are proportional to a small parameter.

preprint2019arXiv

An all speed second order IMEX relaxation scheme for the Euler equations

We present an implicit-explicit finite volume scheme for the Euler equations. We start from the non-dimensionalised Euler equations where we split the pressure in a slow and a fast acoustic part. We use a Suliciu type relaxation model which we split in an explicit part, solved using a Godunov-type scheme based on an approximate Riemann solver, and an implicit part where we solve an elliptic equation for the fast pressure. The relaxation source terms are treated projecting the solution on the equilibrium manifold. The proposed scheme is positivity preserving with respect to the density and internal energy and asymptotic preserving towards the incompressible Euler equations. For this first order scheme we give a second order extension which maintains the positivity property. We perform numerical experiments in 1D and 2D to show the applicability of the proposed splitting and give convergence results for the second order extension.

preprint2019arXiv

An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity

We present an implicit-explicit well-balanced finite volume scheme for the Euler equations with a gravitational source term which is able to deal also with low Mach flows. To visualize the different scales we use the non-dimensionalized equations on which we apply a pressure splitting and a Suliciu relaxation. On the resulting model, we apply a splitting of the flux into a linear implicit and an non-linear explicit part that leads to a scale independent time-step. The explicit step consists of a Godunov type method based on an approximative Riemann solver where the source term is included in the flux formulation. We develop the method for a first order scheme and give an extension to second order. Both schemes are designed to be well-balanced, preserve the positivity of density and internal energy and have a scale independent diffusion. We give the low Mach limit equations for well-prepared data and show that the scheme is asymptotic preserving. These properties are numerically validated by various test cases.

preprint2019arXiv

High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes

We present a new family of very high order accurate direct Arbitrary-Lagrangian-Eulerian (ALE) Finite Volume (FV) and Discontinuous Galerkin (DG) schemes for the solution of nonlinear hyperbolic PDE systems on moving 2D Voronoi meshes that are regenerated at each time step and which explicitly allow topology changes in time. The Voronoi tessellations are obtained from a set of generator points that move with the local fluid velocity. We employ an AREPO-type approach, which rapidly rebuilds a new high quality mesh rearranging the element shapes and neighbors in order to guarantee a robust mesh evolution even for vortex flows and very long simulation times. The old and new Voronoi elements associated to the same generator are connected to construct closed space--time control volumes, whose bottom and top faces may be polygons with a different number of sides. We also incorporate degenerate space--time sliver elements, needed to fill the space--time holes that arise because of topology changes. The final ALE FV-DG scheme is obtained by a redesign of the fully discrete direct ALE schemes of Boscheri and Dumbser, extended here to moving Voronoi meshes and space--time sliver elements. Our new numerical scheme is based on the integration over arbitrary shaped closed space--time control volumes combined with a fully-discrete space--time conservation formulation of the governing PDE system. In this way the discrete solution is conservative and satisfies the GCL by construction. Numerical convergence studies as well as a large set of benchmarks for hydrodynamics and magnetohydrodynamics (MHD) demonstrate the accuracy and robustness of the proposed method. Our numerical results clearly show that the new combination of very high order schemes with regenerated meshes with topology changes lead to substantial improvements compared to direct ALE methods on conforming meshes.

preprint2019arXiv

Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in 2D

The question of well- and ill-posedness of entropy admissible solutions to the multi-dimensional systems of conservation laws has been studied recently in the case of isentropic Euler equations. In this context special initial data were considered, namely the 1D Riemann problem which is extended trivially to a second space dimension. It was shown that there exist infinitely many bounded entropy admissible weak solutions to such a 2D Riemann problem for isentropic Euler equations, if the initial data give rise to a 1D self-similar solution containing a shock. In this work we study such a 2D Riemann problem for the full Euler system in two space dimensions and prove the existence of infinitely many bounded entropy admissible weak solutions in the case that the Riemann initial data give rise to the 1D self-similar solution consisting of two shocks and possibly a contact discontinuity.