BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christian Klingenberg

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

ResearcherUnclaimed source record
math.APmath.NANumerical Analysisastro-ph.IMphysics.comp-phastro-ph.COastro-ph.GAastro-ph.SRphysics.flu-dynCell Behaviormath-phmath.MPmath.OC

Catalog footprint

What is connected

24works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

24 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A numerical scheme for the compressible low-Mach number regime of ideal fluid dynamics

Based on the Roe solver a new technique that allows to correctly represent low Mach number flows with a discretization of the compressible Euler equations was proposed in Miczek et al.: New numerical solver for flows at various mach numbers. A&A 576, A50 (2015). We analyze properties of this scheme and demonstrate that its limit yields a discretization of the continuous limit system. Furthermore we perform a linear stability analysis for the case of explicit time integration and study the performance of the scheme under implicit time integration via the evolution of its condition number. A numerical implementation demonstrates the capabilities of the scheme on the example of the Gresho vortex which can be accurately followed down to Mach numbers of ~1e-10 .

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vertical structure of a supernova-driven turbulent magnetized ISM

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 - 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds (< 10 %) embedded within. At |z| ~ 1-2 kpc, transition-temperature (10^5 K) gas accounts for most of the mass and volume, while hot gas dominates at |z| > 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme's robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma beta. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma beta. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512^3 grid cells we identify a narrow inertial range with the expected power-law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate twice as high from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

MHD simulations of a supernova-driven ISM and the warm ionized medium using a positivity preserving ideal MHD scheme

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

0 citations0 reviewsNo code linkedOpen access