Source author record

Christian Rohde

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

13 published item(s)

preprint2026arXiv

Thermodynamically consistent phase-field modeling and numerical simulation for two-phase fluid-solid dynamics

We introduce a coupled Cahn-Hilliard Navier-Stokes model that governs the two-phase dynamics of a system that consists of a fluid and a solid phase and prove its thermodynamic consistency. Moreover, we present an associated fully-discrete numerical method that relies on a continuous finite element approach and a semi-implicit time-stepping method. As the main theoretical result we show that the fully-discrete method satisfies a discrete analog of the free energy dissipation inequality. Numerical experiments confirm the theoretical findings and show the applicability of the method for realistic settings including an extension to chemically reacting flow. In this context, we provide a preprocessing strategy that enables computing fluid flow in complex geometries given a sharp-interface formulation of the initial phase distribution. Moreover, we briefly introduce different solution strategies for the novel discretization based on the monolithic and partitioned solution paradigms and assess these in a comparative study.

preprint2022arXiv

A Molecular-Continuum Multiscale Model for Inviscid Liquid-Vapor Flow with Sharp Interfaces

The dynamics of compressible liquid-vapor flow depends sensitively on the microscale behavior at the phase boundary. We consider a sharp-interface approach, and propose a multiscale model to describe liquid-vapor flow accurately, without imposing ad-hoc closure relations on the continuum scale. The multiscale model combines the Euler equations on the continuum scale with molecular-scale particle simulations that govern the interface motion. We rely on an interface-preserving moving mesh finite volume method to discretize the continuum-scale sharp-interface flow in a conservative manner. Computational efficiency, while preserving physical properties, is achieved by a surrogate solver for the interface dynamics based on constraint-aware neural networks. The multiscale model is presented in its general form, and applied to regimes of temperature-dependent liquid-vapor flow which have not been accessible before.

preprint2021arXiv

A Finite-Volume Moving-Mesh Method for Two-phase Flow in Fracturing Porous Media

Flow in fractured porous media is modeled frequently by discrete fracture-matrix approaches where fractures are treated as dimensionally reduced manifolds. Generalizing earlier work we focus on two-phase flow in time-dependent fracture geometries including the fracture's aperture. We present the derivation of a reduced model for immiscible two-phase flow in porous media. For the reduced model we present a fully conforming finite-volume discretization coupled with a moving-mesh method. This method permits arbitrary movement of facets of the triangulation while being fully conservative. In numerical examples we show the performance of the scheme and investigate the modeling error of the reduced model.

preprint2020arXiv

On a stochastic Camassa-Holm type equation with higher order nonlinearities

The subject of this paper is a generalized Camassa-Holm equation under random perturbation. We first establish local existence and uniqueness results as well as blow-up criteria for pathwise solutions in the Sobolev spaces $H^s$ with $s>3/2$. Then we analyze how noise affects the dependence of solutions on initial data. Even though the noise has some already known regularization effects, much less is known concerning the dependence on initial data. As a new concept we introduce the notion of stability of exiting times and construct an example showing that multiplicative noise (in Itô sense) cannot improve the stability of the exiting time, and simultaneously improve the continuity of the dependence on initial data. Finally, we obtain global existence theorems and estimate associated probabilities.

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.

preprint2019arXiv

A posteriori error analysis and adaptive non-intrusive numerical schemes for systems of random conservation laws

In this article we consider one-dimensional random systems of hyperbolic conservation laws. We first establish existence and uniqueness of random entropy admissible solutions for initial value problems of conservation laws which involve random initial data and random flux functions. Based on these results we present an a posteriori error analysis for a numerical approximation of the random entropy admissible solution. For the stochastic discretization, we consider a non-intrusive approach, the Stochastic Collocation method. The spatio-temporal discretization relies on the Runge--Kutta Discontinuous Galerkin method. We derive the a posteriori estimator using continuous reconstructions of the discrete solution. Combined with the relative entropy stability framework this yields computable error bounds for the entire space-stochastic discretization error. The estimator admits a splitting into a stochastic and a deterministic (space-time) part, allowing for a novel residual-based space-stochastic adaptive mesh refinement algorithm. We conclude with various numerical examples investigating the scaling properties of the residuals and illustrating the efficiency of the proposed adaptive algorithm.

preprint2019arXiv

Constraint-Aware Neural Networks for Riemann Problems

Neural networks are increasingly used in complex (data-driven) simulations as surrogates or for accelerating the computation of classical surrogates. In many applications physical constraints, such as mass or energy conservation, must be satisfied to obtain reliable results. However, standard machine learning algorithms are generally not tailored to respect such constraints. We propose two different strategies to generate constraint-aware neural networks. We test their performance in the context of front-capturing schemes for strongly nonlinear wave motion in compressible fluid flow. Precisely, in this context so-called Riemann problems have to be solved as surrogates. Their solution describes the local dynamics of the captured wave front in numerical simulations. Three model problems are considered: a cubic flux model problem, an isothermal two-phase flow model, and the Euler equations. We demonstrate that a decrease in the constraint deviation correlates with low discretization errors for all model problems, in addition to the structural advantage of fulfilling the constraint.

preprint2018arXiv

A Particle-based Multiscale Solver for Compressible Liquid-Vapor Flow

To describe complex flow systems accurately, it is in many cases important to account for the properties of fluid flows on a microscopic scale. In this work, we focus on the description of liquid-vapor flow with a sharp interface between the phases. The local phase dynamics at the interface can be interpreted as a Riemann problem for which we develop a multiscale solver in the spirit of the heterogeneous multiscale method, using a particle-based microscale model to augment the macroscopic two-phase flow system. The application of a microscale model makes it possible to use the intrinsic properties of the fluid at the microscale, instead of formulating (ad-hoc) constitutive relations.

preprint2016arXiv

Effect of surfactant concentration and interfacial slip on the flow past a viscous drop at low surface Péclet number

The motion of a viscous drop is investigated when the interface is fully covered with a stagnant layer of surfactant in an arbitrary unsteady Stokes flow for the low surface Péclet number limit. The effect of the interfacial slip coefficient on the behavior of the flow field is also considered. The hydrodynamic problem is solved by the solenoidal decomposition method and the drag force is computed in terms of Faxen's laws using a perturbation ansatz in powers of the surface Péclet number. The analytical expressions for the migration velocity of the drop are also obtained in powers of the surface Péclet number. Further instances corresponding to a given ambient flow as uniform flow, Couette flow, Poiseuille flow are analyzed. Moreover, it is observed that, a surfactant-induced cross-stream migration of the drop occur towards the centre-line in both Couette flow and Poiseuille flow cases. The variation of the drag force and migration velocity is computed for different parameters such as Péclet number, Marangoni number etc.

preprint2016arXiv

On Riemann Solvers and Kinetic Relations for Isothermal Two-Phase Flows with Surface Tension

We consider a sharp-interface approach for the inviscid isothermal dynamics of compressible two-phase flow, that accounts for phase transition and surface tension effects. To fix the mass exchange and entropy dissipation rate across the interface kinetic relations are frequently used. The complete uni-directional dynamics can then be understood by solving generalized two-phase Riemann problems. We present new well-posedness theorems for the Riemann problem and corresponding computable Riemann solvers, that cover quite general equations of state, metastable input data and curvature effects. The new Riemann solver is used to validate different kinetic relations on physically relevant problems including a comparison with experimental data. Riemann solvers are building blocks for many numerical schemes that are used to track interfaces in two-phase flow. It is shown that the new Riemann solver enables reliable and efficient computations for physical situations that could not be treated before.

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.

preprint2015arXiv

Energy-Efficient 5G Outdoor-to-Indoor Communication: SUDAS Over Licensed and Unlicensed Spectrum

In this paper, we study the joint resource allocation algorithm design for downlink and uplink multicarrier transmission assisted by a shared user equipment (UE)-side distributed antenna system (SUDAS). The proposed SUDAS simultaneously utilizes licensed frequency bands and unlicensed frequency bands, (e.g. millimeter wave bands), to enable a spatial multiplexing gain for single-antenna UEs to improve energy efficiency and system throughput of $5$-th generation (5G) outdoor-to-indoor communication. The design of the UE selection, the time allocation to uplink and downlink, and the transceiver processing matrix is formulated as a non-convex optimization problem for the maximization of the end-to-end system energy efficiency (bits/Joule). The proposed problem formulation takes into account minimum data rate requirements for delay sensitive UEs and the circuit power consumption of all transceivers. In order to design a tractable resource allocation algorithm, we first show that the optimal transmitter precoding and receiver post-processing matrices jointly diagonalize the end-to-end communication channel for both downlink and uplink communication via SUDAS. Subsequently, the matrix optimization problem is converted to an equivalent scalar optimization problem for multiple parallel channels, which is solved by an asymptotically globally optimal iterative algorithm. Besides, we propose a suboptimal algorithm which finds a locally optimal solution of the non-convex optimization problem. Simulation results illustrate that the proposed resource allocation algorithms for SUDAS achieve a significant performance gain in terms of system energy efficiency and spectral efficiency compared to conventional baseline systems by offering multiple parallel data streams for single-antenna UEs.

preprint2015arXiv

Resource Allocation for Outdoor-to-Indoor Multicarrier Transmission with Shared UE-side Distributed Antenna Systems

In this paper, we study the resource allocation algorithm design for downlink multicarrier transmission with a shared user equipment (UE)-side distributed antenna system (SUDAS) which utilizes both licensed and unlicensed frequency bands for improving the system throughput. The joint UE selection and transceiver processing matrix design is formulated as a non-convex optimization problem for the maximization of the end-to-end system throughput (bits/s). In order to obtain a tractable resource allocation algorithm, we first show that the optimal transmitter precoding and receiver post-processing matrices jointly diagonalize the end-to-end communication channel. Subsequently, the optimization problem is converted to a scalar optimization problem for multiple parallel channels, which is solved by using an asymptotically optimal iterative algorithm. Simulation results illustrate that the proposed resource allocation algorithm for the SUDAS achieves an excellent system performance and provides a spatial multiplexing gain for single-antenna UEs.