Source author record

Claus-Dieter Munz

Claus-Dieter Munz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA physics.flu-dyn Numerical Analysis physics.comp-ph

Catalog footprint

What is connected

5works

4topics

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.

preprint2022arXiv

Riemann Solvers for Phase Transition in a Compressible Sharp-Interface Method

In this paper, we consider Riemann solvers with phase transition effects based on the Euler-Fourier equation system. One exact and two approximate solutions of the two-phase Riemann problem are obtained by modelling the phase transition process via the theory of classical irreversible thermodynamics. Closure is obtained by appropriate Onsager coefficients for evaporation and condensation. We use the proposed Riemann solvers in a sharp-interface level-set ghost fluid method to couple the individual phases with each other. The Riemann solvers are validated against molecular dynamics data of evaporating Lennard-Jones truncated and shifted fluid. We further study the effects of phase transition on a shock-drop interaction with the novel approximate Riemann solvers.

preprint2020arXiv

A novel regularization strategy for the local discontinuous Galerkin method for level-set reinitialization

In this paper we propose a novel regularization strategy for the local discontinuous Galerkin method to solve the Hamilton-Jacobi equation in the context of level-set reinitialization. The novel regularization idea works in analogy to shock-capturing schemes for discontinuous Galerkin methods, which are based on finite volume sub-cells. In this spirit, the local discontinuous Galerkin method is combined with an upwind/downwind finite volume sub-cell discretization, which is applied in areas of low regularity. To ensure the applicability on unstructured meshes, the finite volume discretization is based on a least squares approach.

preprint2019arXiv

A Parabolic Relaxation Model for the Navier-Stokes-Korteweg Equations

The isothermal Navier-Stokes-Korteweg system is a classical diffuse interface model for compressible two-phase flow. However, the numerical solution faces two major challenges: due to a third-order dispersion contribution in the momentum equations, extended numerical stencils are required for the flux calculation. Furthermore, the equation of state given by a Van-der-Waals law, exhibits non-monotone behaviour in the two-phase state space leading to imaginary eigenvalues of the Jacobian of the first-order fluxes. In this work a lower-order parabolic relaxation model is used to approximate solutions of the classical NSK equations. Whereas an additional parabolic evolution equation for the relaxation variable has to be solved, the system involves no differential operator of higher as second order. The use of a modified pressure function guarantees that the first-order fluxes remain hyperbolic. Altogether, the relaxation system is directly accessible for standard compressible flow solvers. We use the higher-order Discontinuous Galerkin spectral element method as realized in the open source code FLEXI. The relaxation model is validated against solutions of the original NSK model for a variety of 1D and 2D test cases. Three-dimensional simulations of head-on droplet collisions for a range of different collision Weber numbers underline the capability of the approach.

preprint2015arXiv

A sharp interface method for compressible liquid-vapor flow with phase transition and surface tension

The numerical approximation of non-isothermal liquid-vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local multi-phase Riemann problems across the interface. The Riemann solution accounts for the relevant physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the multi-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.

preprint2014arXiv

A Weakly Asymptotic Preserving Low Mach Number Scheme for the Euler Equations of Gas Dynamics

We propose a low Mach number, Godunov-type finite volume scheme for the numerical solution of the compressible Euler equations of gas dynamics. The scheme combines Klein's non-stiff/stiff decomposition of the fluxes (J. Comput. Phys. 121:213-237, 1995) with an explicit/implicit time discretization (Cordier et al., J. Comput. Phys. 231:5685- 5704, 2012) for the split fluxes. This results in a scalar second order partial differential equation (PDE) for the pressure, which we solve by an iterative approximation. Due to our choice of a crucial reference pressure, the stiff subsystem is hyperbolic, and the second order PDE for the pressure is elliptic. The scheme is also uniformly asymptotically consistent. Numerical experiments show that the scheme needs to be stabilized for low Mach numbers. Unfortunately, this affects the asymptotic consistency, which becomes non-uniform in the Mach number, and requires an unduly fine grid in the small Mach number limit. On the other hand, the CFL number is only related to the non-stiff characteristic speeds, independently of the Mach number. Our analytical and numerical results stress the importance of further studies of asymptotic stability in the development of AP (asymptotic preserving) schemes.