BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Dieter Bothe

Dieter Bothe contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynmath-phmath.APmath.MPmath.NANumerical Analysisphysics.comp-phcond-mat.softphysics.chem-phSoftware Engineering

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

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

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

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

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

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access