**Analysis and Numerics of sharp and diffuse interface models for droplet dynamics**

## Description

For the numerical simulation of droplet dynamic processes the liquid phase and the surrounding gaseous fluid have been previously modeled incompressible. But if we consider the interaction of droplets with high-speed flows, droplets near the critical point or even evaporation processes of highly volatile substances, compressible effects can no longer be neglected. In this subproject fully compressible models for the dynamics of free single droplets and droplet ensembles are developed and studied analytically and numerically. The sharp interface approach (SI) and the diffuse interface approach (DI) are considered comparatively.

The SI-approach describes the dynamics in time dependent regions by the hydrodynamic equations for the two fluids. These equations are coupled with jump conditions, so that an overall free boundary value problem is set. For the numerical approximation algorithms are developed, which, in particular, take into account the jumps in the equations and the curvature of the surface. These methods are incorporated in a ghost fluid approach in the code FS3D of the SFB-TRR 75 for solving the SI model numerically.

The SI-approach is invalid when the topology changes – e.g. in case of drop division – and questionable as modelling approach in the transition to supercritical states. Alternatively in phase transition problems phase field models (DI-models) are used very succesfully. Within the DI-models the location of the interface is not exactly given by a jump in the state variables, but smeared out over a small area, such that a system of evolution equations is valid in the whole area. The thickness of the boundary layer is controlled by a regularization parameter. In this subproject a fully compressible DI-model for any state equations will be developed and validated with tools of asymptotic analysis. The focus is on the numerical solution of the DI-model in several dimensions. Therefore an hp-adaptive discontinuous Galerkin method is developed on unstructured grids and implemented. The numerical method is used for the simulation of droplet experiments.

## Team

## Prof. Dr. rer. nat. habil. Christian Rohde |
Director subarea A3 | This email address is being protected from spambots. You need JavaScript enabled to view it. | +49 711 685-65524 | ||

## Jim Magiera, M.Sc. |
A3 | This email address is being protected from spambots. You need JavaScript enabled to view it. | +49 711 685-67646 |

## Publications

#### 2017

Fechter, S.; Zeiler, C.; Munz, C.-D.; Rohde, C.: *A sharp interface method for compressible liquid-vapor flow with phase transition and surface tension*,

J. Comput. Phys., 336, 347-374, 2017.

Chertock, A.; Degond, P. , Neusser, J. *An Asymptotic-Preserving **Method for a Relaxation of the Navier-Stokes-Korteweg Equations*,

Journal of Computational Physics*, *2017*,** *335, 387-403

Magiera, J.; Rohde, C.: *A particle-based multiscale solver for compressible liquid-vapor flow*,

erscheint in: Klingenberg, C.; Westdickenberg, M. (eds.), Hyperbolic Problems. Springer Proc. Math. Stat., Springer, Cham, 2017.

Fechter, S.; Zeiler, C.; Munz, C.-D.; Rohde, C.: *Approximate Riemann solver for compressible liquid vapor flow with phase transition and surface tension*,

erscheint bei Comput. Fluids, http://doi.org/10.1016/j.compfluid.2017.03.026.

#### 2016

Rohde, C., Zeiler, C.: *On Riemann Solvers and Kinetic Relations for Isothermal Two-Phase Flows with Surface Tension*,

eingereicht bei ZAMP, 2016.

B. Kabil, C. Rohde.*Persistence of undercompressive phase boundaries for isothermal Euler equations including configurational forces and surface tension,*

Mathematical Methods in the Applied Sciences, 2016, 39, 5409-5426.

Schleper, V.: *A HLL-type Riemann solver for two-phase flow with surface forces and phase transitions,*

Appl. Numer. Math., 2016, 108, 256-270

Sharanya, V.; Raja Sekhar, G. P. & Rohde, C.: *Bed of polydisperse viscous spherical drops under thermocapillary effects*,

Zeitschrift für angewandte Mathematik und Physik, 2016, 67, 1-17

Dragomirescu, I.; Eisenschmidt, K.; Rohde, C. , Weigand, B.*Perturbation solutions for the finite radially symmetric Stefan problem,*

Inter. J. Thermal Sci., 104, 386-395, 2016.

#### 2015

Rohde, C.; Zeiler, C.: *A relaxation Riemann solver for compressible two-phase flow with phase transition and surface tension*,

Appl. Numer. Math., 95, 267-279, 2015

J. Neusser, C. Rohde, V. Schleper. *Relaxation of the Navier–Stokes–Korteweg equations for compressible two-phase flow with phase transition,*

Int. J. Numer. Meth. Fluids (2015), 79, 615-639.

#### 2014

Fechter, S.; Zeiler, C.; Munz, C.-D.; Rohde, C.: *Simulation of compressible multi-phase flows at extreme ambient conditions using a Discontinuous-Galerkin method*,

ILASS-Europe, 26th European Conference on Liquid Atomization and Spray Systems, 2014.

Engel, P.; Viorel, A.; Rohde, C.*A Low-Order Approximation for Viscous-Capillary Phase Transition Dynamics*,

Port. Math*., * 2014*, 70*, 319-344.

Corli, A.; Rohde, C. , Schleper, V. *Parabolic approximations of **diusive-dispersive equations,*

J. Math. Anal. Appl., 2014, 414,773-798.

Eymard, R.,Schleper, V. *Study of a numerical scheme for miscible**two-phase ow in porous media*,

Numer. Meth. Part. D. E., 2014, 30, 723-748.

Maria Wiebe.

*Ein Sharp-Interface-Ansatz fr Phasenbergangsprobleme **(A Sharp-Interface Approach for Phase Transition Problems)*,

Master Thesis, Institute of Applied Analysis and Numerical Simulation, University of Stuttgart, (2014).

Kabil, B.; Rohde, C.: *The influence of surface tension and configurational forces on the stability of liquid-vapor interfaces*,

Nonlinear Analysis-Theor., 107, 63-75, 2014.

**2013**

K. Eisenschmidt, P. Rauschenberger, C. Rohde, B. Weigand. *Modelling of freezing processes in super-cooled droplets on sub-grid scale*,

ILASS-Europe 2013, European Conference on Liquid Atomization and Spray Systems, Chania, Greece, 2013

Fechter, S., Jägle, F., Schleper, V.* Exact and approximate Riemann solvers at phase boundaries*

Computers & Fluids, 2013, 75, 112-126.

**2012**

Chalons, C., Coquel, F., Engel, P., Rohde, C.*Fast relaxation solvers for hyperbolic-elliptic transition problems*

SIAM J. Sci. Comput. 34 (2012), A1753-A1776

Colombo, R. M., Schleper, V.*Two-phase flow: Non-smooth well-posedness and the compressible to incompressible limit*

Nonlinear Anal. Real World Appl. 13 (2012), 2195-2213.

Corli, A., Rohde, C.*Singular limits for a parabolic regularization of scalar conservation laws*

J. Differential Equations 253 (2012), 1399-1421.

Dreyer, W., Giesselmann, J., Kraus, C., Rohde, C.*Asymptotic analysis for Korteweg models*

Interfaces Free Bound. 14 (2012), 105-143.

Schleper, V.*On the coupling of compressible and incompressible fluids*

Numerical Methods for Hyperbolic Equations (eds. E. Vazquez-Cendon, A. Hidalgo, P. Garcia-Navarro, L. Cea), (2012), 8 Seiten, Taylor & Francis.

Fechter, S., Jagle, F.; Schleper, V.,*A multiscale algorithm for compressible liquid-vapor flow with surface tension*

preprint 2012

Jagle, F.; Zeiler, C.; Rohde, C.*Exact and approximate Riemannsolvers at phase boundaries*,

Preprint 2012.

Eymard, R., Schleper, V*Study of a numerical scheme for miscible two-phase flow*

Numerical Meth- ods for Partial Differential Equations (2012), 22 Seiten. http://hal.archives-ouvertes.fr/hal-00741425

Fechter, S.; Jagle, F.; Boger, M.; Zeiler, C.; Munz, C.-D., Rohde, C.*A discontinuous Galerkin based multiscale method for compressible multiphase flow*,

Proceedings of the 7th ICCFD, Hawaii, 2012.

**2011**

Corli, A.; Rohde, C.

*Singular limits for a parabolic-elliptic regularization of scalar conservation laws*,

preprint, 2011.

Colombo, R. M.; Schleper, V. *Two-phase flows: Non-Smooth well posedness and the compressible to incompressible limit,*

preprint, 2011.

**2010**

Rohde, C*A local and low-order Navier-Stokes-Korteweg system*

Nonlinear partial differential equations and hyperbolic wave phenomena, 315-337, Contemp. Math., 526, Amer. Math. Soc., (2010)

Kissling, F., Rohde, C.*The computation of nonclassical shock waves with a heterogeneous multiscale method*

Netw. Heterog. Media 5 (2010), 661-674