Source author record

Dieter Bothe

Dieter Bothe 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

23works
12topics
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

23 published item(s)

preprint2022arXiv

A Research Software Engineering Workflow for Computational Science and Engineering

University research groups in Computational Science and Engineering (CSE) generally lack dedicated funding and personnel for Research Software Engineering (RSE), which, combined with the pressure to maximize the number of scientific publications, shifts the focus away from sustainable research software development and reproducible results. The neglect of RSE in CSE at University research groups negatively impacts the scientific output: research data - including research software - related to a CSE publication cannot be found, reproduced, or re-used, different ideas are not combined easily into new ideas, and published methods must very often be re-implemented to be investigated further. This slows down CSE research significantly, resulting in considerable losses in time and, consequentially, public funding. We propose a RSE workflow for Computational Science and Engineering (CSE) that addresses these challenges, that improves the quality of research output in CSE. Our workflow applies established software engineering practices adapted for CSE: software testing, result visualization, and periodical cross-linking of software with reports/publications and data, timed by milestones in the scientific publication process. The workflow introduces minimal work overhead, crucial for university research groups, and delivers modular and tested software linked to publications whose results can easily be reproduced. We define research software quality from a perspective of a pragmatic researcher: the ability to quickly find the publication, data, and software related to a published research idea, quickly reproduce results, understand or re-use a CSE method, and finally extend the method with new research ideas.

preprint2022arXiv

On the molecular mechanism behind the bubble rise velocity jump discontinuity in viscoelastic liquids

Bubbles rising in viscoelastic liquids may exhibit a jump discontinuity of the rise velocity as a critical bubble volume is exceeded. The phenomenon has been extensively investigated in the literature, both by means of experiments as well as via numerical simulations. The occurrence of the velocity jump has been associated with a change of the bubble shape under the formation of a pointed tip at the rear end and to the appearance of a so-called negative wake with the liquid velocity behind the bubble, pointing in the opposite direction to that in viscous Newtonian fluids. We revisit this topic, starting with a review of the state of knowledge on the interrelations between the mentioned characteristic features. In search for a convincing explanation of the jump phenomenon, we performed detailed numerical simulations of the transient rise of single bubbles in 3D, allowing for a local polymer molecular conformation tensor analysis. The latter shows that polymer molecules traveling along the upper bubble hemisphere are stretched in the circumferential direction, due to the flow kinematics. Then, depending on the relaxation time scale of the polymer, the stored elastic energy is either unloaded essentially above or below the bubble's equator. In the former case, this slows down the bubble, while the bubble gets accelerated otherwise. In this latter case, the relative velocity of the polymer molecules against the bubble is increased, giving rise to a self-amplification of the effect and thus causing the bubble rise velocity to jump to a higher level. Detailed experimental velocity measurements in the liquid field around the bubble confirmed the conclusion that the ratio of the time scale of the Lagrangian transport of polymer molecules along the bubble contour to the relaxation time scale of the polymer molecules determines the sub- or supercritical state of the bubble motion.

preprint2022arXiv

Sharp-Interface Continuum Thermodynamics of multicomponent fluid systems with interfacial mass

We revisit the sharp-interface continuum thermodynamics of two-phase multicomponent fluid systems with interfacial mass. Since the published work is not fully consistent, we provide a rigorous derivation of the local balance equations and the entropy production rates, including all relevant steps and mathematical tools. Special emphasis is put on an axiomatic form of the entropy principle which, at the same time, allows for an efficient closure process. The latter builds on an appropriate representation of the entropy production rates, which is based on structural information being obtained from a more refined model which includes partial momenta and is given in an appendix. The closure relations obtained for the one-sided bulk-interface species transfer rates are further investigated, which leads to a novel model for mass transfer at fluid interfaces, influenced by the presence of adsorbed surface active agents. We finally compare our results in detail with the existing literature on continuum thermodynamics of two-phase multicomponent fluid systems with interfacial mass.

preprint2020arXiv

A geometry-based model for spreading drops applied to drops on a silicon wafer and a swellable polymer brush film

We investigate the dynamics of spreading in a regime where the shape of the drop is close to a spherical cap. The latter simplification is applicable in the late (viscous) stage of spreading for highly viscous drops with a diameter below the capillary length. Moreover, it applies to the spreading of a drop on a swellable polymer brush, where the complex interaction with the substrate leads to a very slow spreading dynamics. The spherical cap geometry allows to derive a closed ordinary differential equation (ODE) for the spreading if the capillary number is a function of the contact angle as it is the case for empirical contact angle models. The latter approach has been introduced by de Gennes (Reviews of Modern Physics, 1985) for small contact angles. In the present work, we generalize the method to arbitrary contact angles. The method is applied to experimental data of spreading water-glycerol drops on a silicon wafer and spreading water drops on a PNIPAm coated silicon wafer. It is found that the ODE-model is able to describe the spreading kinetics in the case of partial wetting. Moreover, the model can predict the spreading dynamics of spherical cap-shaped droplets if the relationship between the contact angle and the capillary number is universal.

preprint2020arXiv

Boundary conditions for dynamic wetting -- A mathematical analysis

The moving contact line paradox discussed in the famous paper by Huh and Scriven has lead to an extensive scientific discussion about singularities in continuum mechanical models of dynamic wetting in the framework of the two-phase Navier-Stokes equations. Since the no-slip condition introduces a non-integrable and therefore unphysical singularity into the model, various models to relax the singularity have been proposed. Many of the relaxation mechanisms still retain a weak (integrable) singularity, while other approaches look for completely regular solutions with finite curvature and pressure at the moving contact line. In particular, the model introduced recently in (Lukyanov, Pryer, Langmuir 2017) aims for regular solutions through modified boundary conditions. The present work applies the mathematical tool of compatibility analysis to continuum models of dynamic wetting. The basic idea is that the boundary conditions have to be compatible at the contact line in order to allow for regular solutions. Remarkably, the method allows to compute explicit expressions for the pressure and the curvature locally at the moving contact line for regular solutions to the model of Lukyanov and Pryer. It is found that solutions may still be singular for the latter model.

preprint2020arXiv

Contact line advection using the geometrical Volume-of-Fluid method

We consider the interface advection problem by a prescribed velocity field in the special case when the interface intersects the domain boundary, i.e. in the presence of a contact line. This problem emerges from the discretization of continuum models for dynamic wetting. The kinematic evolution equation for the dynamic contact angle (Fricke et al., 2019) expresses the fundamental relationship between the rate of change of the contact angle and the structure of the transporting velocity field. The goal of the present work is to develop an interface advection method that is consistent with the fundamental kinematics and transports the contact angle correctly with respect to a prescribed velocity field. In order to verify the advection method, the kinematic evolution equation is solved numerically and analytically (for special cases). We employ the geometrical Volume-of-Fluid (VOF) method on a structured Cartesian grid to solve the hyperbolic transport equation for the interface in two spatial dimensions. We introduce generalizations of the Youngs and ELVIRA methods to reconstruct the interface close to the domain boundary. Both methods deliver first-order convergent results for the motion of the contact line. However, the Boundary Youngs method shows strong oscillations in the numerical contact angle that do not converge with mesh refinement. In contrast to that, the Boundary ELVIRA method provides linear convergence of the numerical contact angle transport.

preprint2020arXiv

Mass transport in multicomponent compressible fluids: Local and global well-posedness in classes of strong solutions for general class-one models

We consider a system of partial differential equations describing mass transport in a multicomponent isothermal compressible fluid. The diffusion fluxes obey the Fick-Onsager or Maxwell-Stefan closure approach. Mechanical forces result into one single convective mixture velocity, the barycentric one, which obeys the Navier-Stokes equations. The thermodynamic pressure is defined by the Gibbs-Duhem equation. Chemical potentials and pressure are derived from a thermodynamic potential, the Helmholtz free energy, with a bulk density allowed to be a general convex function of the mass densities of the constituents. The resulting PDEs are of mixed parabolic--hyperbolic type. We prove two theoretical results concerning the well-posedness of the model in classes of strong solutions: 1. The solution always exists and is unique for short--times and 2. If the initial data are sufficiently near to an equilibrium solution, the well-posedness is valid on arbitrary large, but finite time intervals. Both results rely on a contraction principle valid for systems of mixed type that behave like the compressible Navier-Stokes equations. The linearised parabolic part of the operator possesses the self map property with respect to some closed ball in the state space, while being contractive in a lower order norm only. In this paper, we implement these ideas by means of precise a priori estimates in spaces of exact regularity.

preprint2020arXiv

On the structure of continuum thermodynamical diffusion fluxes -- A novel closure scheme and its relation to the Maxwell-Stefan and the Fick-Onsager approach

This paper revisits the modeling of multicomponent diffusion within the framework of thermodynamics of irreversible processes. We briefly review the two well-known main approaches, leading to the generalized Fick-Onsager multicomponent diffusion fluxes or to the generalized Maxwell-Stefan equations. The latter approach has the advantage that the resulting fluxes are consistent with non-negativity of the partial mass densities for non-singular and non-degenerate Maxwell-Stefan diffusivities. On the other hand, this approach requires computationally expensive matrix inversions since the fluxes are only implicitly given. We propose and discuss a novel and more direct closure which avoids the inversion of the Maxwell-Stefan equations. It is shown that all three closures are actually equivalent under the natural requirement of positivity for the concentrations, thus revealing the general structure of continuum thermodynamical diffusion fluxes. As a special case, the new closure also gives rise to a core-diagonal diffusion model in which only those cross-effects are present that are necessary to guarantee consistency with total mass conservation, plus a compositional dependence of the diffusivity. This core-diagonal closure turns out to provide a rigorous fundament for recent extensions of the Darken equation from binary mixtures to the general multicomponent case. As an outcome of our investigation, we also address different questions related to the sign of multicomponent thermodynamic or Fickian diffusion coefficients. We show rigorously that in general the second law requires positivity properties for tensors and operators rather than for scalar diffusivities.

preprint2020arXiv

Well-posedness analysis of multicomponent incompressible flow models

In this paper, we extend our study of mass transport in multicomponent isothermal fluids to the incompressible case. For a mixture, incompressibility is defined as the independence of average volume on pressure, and a weighted sum of the partial mass densities stays constant. In this type of models, the velocity field in the Navier-Stokes equations is not solenoidal and, due to different specific volumes of the species, the pressure remains connected to the densities by algebraic formula. By means of a change of variables in the transport problem, we equivalently reformulate the PDE system as to eliminate positivity and incompressibility constraints affecting the density, and prove two type of results: the local-in-time well-posedness in classes of strong solutions, and the global-in-time existence of solutions for initial data sufficiently close to a smooth equilibrium solution.

preprint2019arXiv

A Kinematic Evolution Equation for the Dynamic Contact Angle and some Consequences

We investigate the moving contact line problem for two-phase incompressible flows with a kinematic approach. The key idea is to derive an evolution equation for the contact angle in terms of the transporting velocity field. It turns out that the resulting equation has a simple structure and expresses the time derivative of the contact angle in terms of the velocity gradient at the solid wall. Together with the additionally imposed boundary conditions for the velocity, it yields a more specific form of the contact angle evolution. Thus, the kinematic evolution equation is a tool to analyze the evolution of the contact angle. Since the transporting velocity field is required only on the moving interface, the kinematic evolution equation also applies when the interface moves with its own velocity independent of the fluid velocity. We apply the developed tool to a class of moving contact line models which employ the Navier slip boundary condition. We derive an explicit form of the contact angle evolution for sufficiently regular solutions, showing that such solutions are unphysical. Within the simplest model, this rigorously shows that the contact angle can only relax to equilibrium if some kind of singularity is present at the contact line. Moreover, we analyze more general models including surface tension gradients at the contact line, slip at the fluid-fluid interface and mass transfer across the fluid-fluid interface.

preprint2018arXiv

Highly accurate numerical computation of implicitly defined volumes using the Laplace-Beltrami operator

This paper introduces a novel method for the efficient and accurate computation of the volume of a domain whose boundary is given by an orientable hypersurface which is implicitly given as the iso-contour of a sufficiently smooth level-set function. After spatial discretization, local approximation of the hypersurface and application of the Gaussian divergence theorem, the volume integrals are transformed to surface integrals. Application of the surface divergence theorem allows for a further reduction to line integrals which are advantageous for numerical quadrature. We discuss the theoretical foundations and provide details of the numerical algorithm. Finally, we present numerical results for convex and non-convex hypersurfaces embedded in cuboidal domains, showing both high accuracy and thrid- to fourth-order convergence in space.

preprint2016arXiv

Global wellposedness for a class of reaction-advection-anisotropic-diffusion systems

We prove existence and uniqueness of global solutions for a class of reaction-advection-anisotropic-diffusion systems whose reaction terms have a "triangular structure". We thus extend previous results to the case of time-space dependent anisotropic diffusions and with time-space dependent advection terms. The corresponding models are in particular relevant for transport processes inside porous media and in situations in which additional migration occurs. The proofs are based on optimal $L^p$-maximal regularity results for the general time-dependent linear operator dual to the one involved in the considered systems. As an application, we prove global well-posedness for a prototypical class of chemically reacting systems with mass-action kinetics, involving networks of reactions of the type $C_1+\ldots+C_{P-1} \rightleftharpoons C_{P} $. Finally, we analyze how a classical a priori $L^2$-estimate of the solutions, which holds with this kind of nonlinear reactive terms, extends to our general anisotropic-advection framework. It does extend with the same assumptions for isotropic diffusions and is replaced by an $L^{(N+1)/N}$-estimate in the general situation.

preprint2015arXiv

On the Interface Formation Model for Dynamic Triple Lines

This paper revisits the theory of Y. Shikhmurzaev on forming interfaces as a continuum thermodynamical model for dynamic triple lines. We start with the derivation of the balances for mass, momentum, energy and entropy in a three-phase fluid system with full interfacial physics, including a brief review of the relevant transport theorems on interfaces and triple lines. Employing the entropy principle in the form given in [Bothe & Dreyer, Acta Mechanica, doi:10.1007/s00707-014-1275-1] but extended to this more general case, we arrive at the entropy production and perform a linear closure, except for a nonlinear closure for the sorption processes. Specialized to the isothermal case, we obtain a thermodynamically consistent mathematical model for dynamic triple lines and show that the total available energy is a strict Lyapunov function for this system.

preprint2015arXiv

On the multi-physics of mass-transfer across fluid interfaces

Mass transfer of gaseous components from rising bubbles to the ambient liquid can be described based on continuum mechanical sharp-interface balances of mass, momentum and species mass. In this context, the standard model consists of the two-phase Navier-Stokes equations for incompressible fluids with constant surface tension, complemented by reaction-advection-diffusion equations for all constituents, employing Fick's law. This standard model is inconsistent with the continuity equation, the momentum balance and the second law of thermodynamics. The present paper reports on the details of these severe shortcomings and provides thermodynamically consistent model extensions which are required to capture various phenomena which occur due to the multi-physics of interfacial mass transfer. In particular, we provide a simple derivation of the interface Maxwell-Stefan equations which does not require a time scale separation, while the main contribution is to show how interface concentrations and interface chemical potentials mediate the influence on mass transfer of a transfer component exerted by the change in interface energy due to an adsorbing surfactant.

preprint2014arXiv

A graph-theoretical approach for the computation of connected iso-surfaces based on volumetric data

The existing combinatorial methods for iso-surface computation are efficient for pure visualization purposes, but it is known that the resulting iso-surfaces can have holes, and topological problems like missing or wrong connectivity can appear. To avoid such problems, we introduce a graph-theoretical method for the computation of iso-surfaces on cuboid meshes in $\mathbb{R}^3$. The method for the generation of iso-surfaces employs labeled cuboid graphs $G(V,E,\mathcal{F})$ such that $V$ is the set of vertices of a cuboid $C\subset\mathbb{R}^3$, $E$ is the set of edges of $C$ and $\mathcal{F}\,:\,V\rightarrow [0,1]$. The nodes of $G$ are weighted by the values of $\mathcal{F}$ which represents the volumetric information, e.g.\ from a Volume of Fluid method. Using a given iso-level $c\in (0,1)$, we first obtain all iso-points, i.e.\ points where the value $c$ is attained by the edge-interpolated $\mathcal{F}$-field. The iso-surface is then built from iso-elements which are composed of triangles and are such that their polygonal boundary has only iso-points as vertices. All vertices lie on the faces of a single mesh cell. We give a proof that the generated iso-surface is connected up to the boundary of the domain and it can be decomposed into different oriented components. Two different components may have discrete points or line segments in common. The graph-theoretical method for the computation of iso-surfaces developed in this paper enables to recover local information of the iso-surface that can be used e.g.\ to compute discrete mean curvature and to solve surface PDEs. Concerning the computational effort, the resulting algorithm is as efficient as existing combinatorial methods.

preprint2014arXiv

Continuum thermodynamics of chemically reacting fluid mixtures

We consider viscous, heat conducting mixtures of molecularly miscible chemical species forming a fluid in which the constituents can undergo chemical reactions. Assuming a common temperature for all components, we derive a closed system of partial mass and partial momentum balances plus a mixture balance of internal energy. This is achieved by careful exploitation of the entropy principle and requires appropriate definitions of absolute temperature and chemical potentials, based on an adequate definition of thermal energy excluding diffusive contributions. The resulting interaction forces split into a thermo-mechanical and a chemical part, where the former turns out to be symmetric in case of binary interactions. For chemically reacting systems and as a new result, the chemical interaction force is a contribution being non-symmetric outside of chemical equilibrium. The theory also provides a rigorous derivation of the so-called generalized thermodynamic driving forces, avoiding the use of approximate solutions to the Boltzmann equations. Moreover, using an appropriately extended version of the entropy principle and introducing cross-effects already before closure as entropy invariant couplings between principal dissipative mechanisms, the Onsager symmetry relations become a strict consequence. With a classification of the factors in the binary products of the entropy production according to their parity--instead of the classical partition into so-called fluxes and driving forces--the apparent anti-symmetry of certain couplings is thereby also revealed. If the diffusion velocities are small compared to the speed of sound, the Maxwell-Stefan equations follow in the case without chemistry, thereby neglecting wave phenomena in the diffusive motion. This results in a reduced model with only mass being balanced individually. In the reactive case ...

preprint2014arXiv

Global existence for a class of reaction-diffusion systems with mass action kinetics and concentration-dependent diffusivities

In this work we study the existence of classical solutions for a class of reaction-diffusion systems with quadratic growth naturally arising in mass action chemistry when studying networks of reactions of the type $A_i+A_j \rightleftharpoons A_k$ with Fickian diffusion, where the diffusion coefficients might depend on time, space and on all the concentrations $c_i$ of the chemical species. In the case of one single reaction, we prove global existence for space dimensions $N\leq 5$. In the more restrictive case of diffusion coefficients of the type $d_i(c_i)$, we use an $L^2$-approach to prove global existence for $N\leq 9$. In the general case of networks of such reactions we extend the previous method to get global solutions for general diffusivities if $N\leq 3$ and for diffusion of type $d_i(c_i)$ if $N\leq 5$. In the latter quasi-linear case of $d_i(c_i)$ and for space dimensions $N=2$ and $N=3$, global existence holds for more than quadratic reactions. We can actually allow for more general rate functions including fractional power terms, important in applications. We obtain global existence under appropriate growth restrictions with an explicit dependence on the space dimension $N$.

preprint2014arXiv

Thermodynamically consistent modeling for dissolution/growth of bubbles in an incompressible solvent

We derive mathematical models of the elementary process of dissolution/growth of bubbles in a liquid under pressure control. The modeling starts with a fully compressible version, both for the liquid and the gas phase so that the entropy principle can be easily evaluated. This yields a full PDE system for a compressible two-phase fluid with mass transfer of the gaseous species. Then the passage to an incompressible solvent in the liquid phase is discussed, where a carefully chosen equation of state for the liquid mixture pressure allows for a limit in which the solvent density is constant. We finally provide a simplification of the PDE system in case of a dilute solution.

preprint2013arXiv

Global existence for diffusion-electromigration systems in space dimension three and higher

We prove existence of global weak solutions for the Nernst-Planck-Poisson problem which describes the evolution of concentrations of charged species $X_1, ..., X_P$ subject to Fickian diffusion and chemical reactions in the presence of an electrical field, including in particular the Boltzmann statistics case. In contrast to the existing literature, existence is proved in any dimension. Moreover, we do not need the assumption $P = 2$ nor the assumption of equal diffusivities for all $P$ components. Our approach relies on the intrinsic energy structure and on an adequate nonlinear and curiously more regular approximate problem. The delicate passing to the limit is done in adequate functional spaces which lead to only weak solutions.

preprint2013arXiv

voFoam - A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM

A new parallelized unsplit geometrical Volume of Fluid (VoF) algorithm with support for arbitrary unstructured meshes and dynamic local Adaptive Mesh Refinement (AMR), as well as for two and three dimensional computation is developed. The geometrical VoF algorithm supports arbitrary unstructured meshes in order to enable computations involving flow domains of arbitrary geometrical complexity. The implementation of the method is done within the framework of the OpenFOAM library for Computational Continuum Mechanics (CCM) using the C++ programming language with modern policy based design for high program code modularity. The development of the geometrical VoF algorithm significantly extends the method base of the OpenFOAM library by geometrical volumetric flux computation for two-phase flow simulations. For the volume fraction advection, a novel unsplit geometrical algorithm is developed, which inherently sustains volume conservation utilizing unique Lagrangian discrete trajectories located in the mesh points. This practice completely eliminates the possibility of an overlap between the flux polyhedra and hence significantly increases volume conservation. A new efficient (quadratic convergent) and accurate iterative flux correction algorithm is developed, which avoids topological changes of the flux polyhedra. Our geometrical VoF algorithm is dimension agnostic, providing automatic support for both 2D and 3D computations, following the established practice in OpenFOAM. The geometrical algorithm used for the volume fraction transport has been extended to support dynamic local AMR available in OpenFOAM. Furthermore, the existing dynamic mesh capability of OpenFOAM has been modified to support the geometrical mapping algorithm executed as a part of the dynamic local AMR cycle. The method implementation is fully parallelized using the domain decomposition approach.

preprint2012arXiv

On a Class of Energy Preserving Boundary Conditions for Incompressible Newtonian Flows

We derive a class of energy preserving boundary conditions for incompressible Newtonian flows and prove local-in-time well-posedness of the resulting initial boundary value problems, i.e. the Navier-Stokes equations complemented by one of the derived boundary conditions, in an Lp-setting in domains, which are either bounded or unbounded with almost flat, sufficiently smooth boundary. The results are based on maximal regularity properties of the underlying linearisations, which are also established in the above setting.

preprint2012arXiv

On Two-Phase Flows with Soluble Surfactant

The presence of surfactants has a pronounced effect on the surface tension and, hence, on the stress balance at the phase separating interface of two-phase flows. The transport of momentum induced by the local variations of the capillary forces are known as Marangoni effects. Here we study a model, which assumes the surfactant to be soluble in one of the adjacent bulk phases and which represents a generalization of the two-phase Navier-Stokes equations. Based on maximal Lp-regularity results for suitable linearizations we obtain local well-posedness of this model. We employ recent results from the Lp-theory of two-phase flows without surfactant.

preprint2010arXiv

On the Maxwell-Stefan approach to multicomponent diffusion

We consider the system of Maxwell-Stefan equations which describe multicomponent diffusive fluxes in non-dilute solutions or gas mixtures. We apply the Perron-Frobenius theorem to the irreducible and quasi-positive matrix which governs the flux-force relations and are able to show normal ellipticity of the associated multicomponent diffusion operator. This provides local-in-time wellposedness of the Maxwell-Stefan multicomponent diffusion system in the isobaric, isothermal case.