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

Jens Harting

Jens Harting contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softphysics.flu-dynphysics.comp-phphysics.app-phphysics.chem-phBiological Physicscond-mat.mtrl-scicond-mat.stat-mechDistributed, Parallel, and Cluster Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
9topics
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

18 published item(s)

preprint2022arXiv

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions

We present guidelines to estimate the effect of electrostatic repulsion in sedimenting dilute particle suspensions. Our results are based on combined Langevin dynamics and lattice Boltzmann simulations for a range of particle radii, Debye lengths and particle concentrations. They show a simple relationship between the slope $K$ of the concentration-dependent sedimentation velocity and the range $χ$ of the electrostatic repulsion normalized by the average particle-particle distance. When $χ\to 0$, the particles are too far away from each other to interact electrostatically and $K=6.55$ as predicted by the theory of Batchelor. As $χ$ increases, $K$ likewise increases as if the particle radius increased in proportion to $χ$ up to a maximum around $χ=0.4$. Over the range $χ=0.4-1$, $K$ relaxes exponentially to a concentration-dependent constant consistent with known results for ordered particle distributions. Meanwhile the radial distribution function transitions from a disordered gas-like to a liquid-like form. Power law fits to the concentration-dependent sedimentation velocity similarly yield a simple master curve for the exponent as a function of $χ$, with a step-like transition from 1 to 1/3 centered around $χ= 0.6$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit and novel 16-bit number formats

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop novel 16-bit formats - based on a modified IEEE-754 and on a modified Posit standard - that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method) and finally the impact of a raindrop (based on a Volume-of-Fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Phase-field simulations of the morphology formation in evaporating crystalline multicomponent films

In numerous solution-processed thin films, a complex morphology resulting from liquid-liquid phase separation (LLPS) or from polycrystallization arises during the drying or subsequent processing steps. The morphology has a strong influence on the performance of the final device but unfortunately the process-structure relationship is often poorly and only qualitatively understood. This is because many different physical mechanisms (miscibility, evaporation, crystallization, diffusion, advection) are active at potentially different time scales, and because the kinetics plays a crucial role: the morphology develops until it is kinetically quenched far from equilibrium. In order to unravel the various possible structure formation pathways, we propose a unified theoretical framework that takes into account all these physical phenomena. This phase-field simulation tool is based on the Cahn-Hilliard equations for diffusion and the Allen-Cahn equation for crystallization and evaporation, which are coupled to the equations for the dynamics of the fluid. We discuss and verify the behavior of the coupled model based on simple test cases. Furthermore, we illustrate how this framework allows to investigate the morphology formation in a drying film undergoing evaporation-induced LLPS and crystallization, which is typically a situation encountered, e.g., in organic photovoltaics applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Controlling the dewetting morphologies of thin liquid films by switchable substrates

Switchable and adaptive substrates emerged as valuable tools for the control of wetting and actuation of droplet motion. Here we report a computational study of the dynamics of an unstable thin liquid film deposited on a switchable substrate, modelled with a space and time varying contact angle. With a static pattern, all the fluid is drained into droplets located around contact angle minima, whereas for a sufficiently large rate of wettability variation a state consisting of metastable rivulets is observed. A criterion discriminating whether rivulets can be observed or not is identified in terms of a single dimensionless parameter. Finally, we show and explain theoretically how the film rupture times, droplet shape and rivulet life time depend on the pattern wavelength and speed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Theoretical framework for pairwise microswimmer interactions

Hydrodynamic interactions are crucial for determining the cooperative behavior of microswimmers at low Reynolds numbers. Here we provide a comprehensive analysis of the scaling and strength of the interactions in the case of a pair of three-sphere swimmers with intrinsic elasticity. Both stroke-based and force-based microswimmers are analyzed using an analytic perturbative approach. Following a detailed analysis of the passive interactions, as well as active translations and rotations, we find that the mapping between the stroke-based and force-based swimmers is only possible in a low driving frequency regime where the characteristic time scale is smaller than the viscous one. Furthermore, we find that for swimmers separated by up to hundreds of swimmer lengths, swimming in pairs speeds up the self propulsion, due to the dominant quadrupolar hydrodynamic interactions. Finally, we find that the long term behavior of the swimmers, while sensitive to initial relative positioning, does not depend on the pusher or puller nature of the swimmer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A phase-field model for the evaporation of thin film mixtures

The performance of solution-processed solar cells strongly depends on the geometrical structure and roughness of the photovoltaic layers formed during film drying. During the drying process, the interplay of crystallization and liquid-liquid demixing leads to the structure formation on the nano- and microscale and to the final rough film. In order to better understand how the film structure can be improved by process engineering, we aim at theoretically investigating these systems by means of phase-field simulations. We introduce an evaporation model based on the Cahn-Hilliard equation for the evolution of the fluid concentrations coupled to the Allen-Cahn equation for the liquid-vapour phase transformation. We demonstrate its ability to match the experimentally measured drying kinetics and study the impact of the parameters of our model. Furthermore, the evaporation of solvent blends and solvent-vapour annealing are investigated. The dry film roughness emerges naturally from our set of equations, as illustrated through preliminary simulations of spinodal decomposition and film drying on structured substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capillary interactions between soft capsules protruding through thin fluid films

When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying. Here, we use our recently-introduced numerical method that couples a fluid described using the lattice Boltzmann approach to a finite element description of deformable objects to investigate the drying process of a film with suspended soft particles. Our measured menisci deformations and lateral capillary forces, which agree well with previous theoretical and experimental works in case of rigid particles, show the deformations become smaller with increasing particles softness, resulting in weaker lateral interaction forces. At large interparticle distances, the force approaches that of rigid particles. Finally, we investigate the time dependent formation of particle clusters at the late stages of the film drying.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Catalytic flow with a coupled Finite Difference -- Lattice Boltzmann scheme

Many catalyst devices employ flow through porous structures, which leads to a complex macroscopic mass and heat transport. To unravel the detailed dynamics of the reactive gas flow, we present an all-encompassing model, consisting of thermal lattice Boltzmann model by Kang et al., used to solve the heat and mass transport in the gas domain, coupled to a finite differences solver for the heat equation in the solid via thermal reactive boundary conditions for a consistent treatment of the reaction enthalpy. The chemical surface reactions are incorporated in a flexible fashion through flux boundary conditions at the gas-solid interface. We scrutinize the thermal FD-LBM by benchmarking the macroscopic transport in the gas domain as well as conservation of the enthalpy across the solid-gas interface. We exemplify the applicability of our model by simulating the reactive gas flow through a microporous material catalysing the so-called water-gas-shift reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Desorption energy of soft particles from a fluid interface

The efficiency of soft particles to stabilize emulsions is examined by measuring their desorption free energy, i.e., the mechanical work required to detach the particle from a fluid interface. Here, we consider rubber-like elastic as well as microgel particles, using coarse-grained molecular dynamics simulations. The energy of desorption is computed for two and three-dimensional configurations by means of the mean thermodynamic integration method. It is shown that the softness affects the particle-interface binding in two opposing directions as compared to rigid particles. On the one hand, a soft particle spreads at the interface and thereby removes a larger unfavorable liquid-liquid contact area compared to rigid particles. On the other hand, softness provides the particle with an additional degree of freedom to get reshaped instead of deforming the interface, resulting in a smaller restoring force during the detachment. It is shown that the first effect prevails so that a soft spherical particle attaches to the fluid interface more strongly than rigid spheres. Finally, we consider microgel particles both in the swollen and in the collapsed state. Surprisingly, we find that the latter has a larger binding energy. All results are rationalised using thermodynamic arguments and thereby offer detailed insights into the desorption energy of soft particles from fluid interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Role of the interplay between spinodal decomposition and crystal growth in the morphological evolution of crystalline bulk heterojunctions

The stability of organic solar cells is strongly affected by the morphology of the photoactive layers, whose separated crystalline and/or amorphous phases are kinetically quenched far from their thermodynamic equilibrium during the production process. The evolution of these structures during the lifetime of the cell remains poorly understood. In this paper, a phase-field simulation framework is proposed, handling liquid-liquid demixing and polycrystalline growth at the same time in order to investigate the evolution of crystalline immiscible binary systems. We find that initially, the nuclei trigger the spinodal decomposition, while the growing crystals quench the phase coarsening in the amorphous mixture. Conversely, the separated liquid phases guide the crystal growth along the domains of high concentration. It is also demonstrated that with a higher crystallization rate, in the final morphology, single crystals are more structured and form percolating pathways for each material with smaller lateral dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Self-Similar Liquid Lens Coalescence

A basic feature of liquid drops is that they can merge upon contact to form a larger drop. In spite of its importance to various applications, drop coalescence on pre-wetted substrates has received little attention. Here, we experimentally and theoretically reveal the dynamics of drop coalescence on a thick layer of a low-viscosity liquid. It is shown that these so-called "liquid lenses" merge by the self-similar vertical growth of a bridge connecting the two lenses. Using a slender analysis, we derive similarity solutions corresponding to the viscous and inertial limits. Excellent agreement is found with the experiments without any adjustable parameters, capturing both the spatial and temporal structure of the flow during coalescence. Finally, we consider the crossover between the two regimes and show that all data of different lens viscosities collapse on a single curve capturing the full range of the coalescence dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure and rheology of suspensions of spherical strain-hardening capsules

We investigate the rheology of strain-hardening spherical capsules, from the dilute to the concentrated regime under a confined shear flow using three-dimensional numerical simulations. We consider the effect of capillary number, volume fraction and membrane inextensibility on the particle deformation and on the effective suspension viscosity and normal stress differences of the suspension. The suspension displays a shear-thinning behaviour which is a characteristic of soft particles such as emulsion droplets, vesicles, strain-softening capsules, and red blood cells. We find that the membrane inextensibility plays a significant role on the rheology and can almost suppress the shear-thinning. For concentrated suspensions a non-monotonic dependence of the normal stress differences on the membrane inextensibility is observed, reflecting a similar behaviour in the particle shape. The effective suspension viscosity, instead, grows and eventually saturates, for very large inextensibilities, approaching the solid particle limit. In essence, our results reveal that strain-hardening capsules share rheological features with both soft and solid particles depending on the ratio of the area dilatation to shear elastic modulus. Furthermore, the suspension viscosity exhibits a universal behaviour for the parameter space defined by the capillary number and the membrane inextensibility, when introducing the particle geometrical changes at the steady-state in the definition of the volume fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A general perturbative approach for bead-based microswimmers reveals rich self-propulsion phenomena

Studies of model microswimmers have significantly contributed to the understanding of the principles of self-propulsion we have today. However, only a small number of microswimmer types have been amenable to analytic modeling, and further development of such approaches is necessary to identify the key features of these active systems. Here we present a general perturbative calculation scheme for swimmers composed of beads interacting by harmonic potentials, driven by an arbitrary force protocol. The hydrodynamic interactions are treated using the Oseen and Rotne-Pragner approximations. We validate our approach by using 3 bead assemblies and comparing the results with the numerically obtained full-solutions of the governing equations of motion, as well as with the existing analytic models for a linear and a triangular swimmer geometries. While recovering the relation between the force and swimming velocity, our detailed analysis and the controlled level of approximation allow us to find qualitative differences already in the far field flow of the devices. Consequently, we are able to identify a behavior of the swimmer that is richer than predicted in previous models. Given its generality, the framework can be applied to any swimmer geometry, driving protocol and beads interactions, as well as in many swimmers problems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Inertial migration of neutrally-buoyant particles in superhydrophobic channels

At finite Reynolds numbers particles migrate across flow streamlines to their equilibrium positions in microchannels. Such a migration is attributed to an inertial lift force, and it is well-known that the equilibrium location of neutrally-buoyant particles is determined only by their size and the Reynolds number. Here we demonstrate that the decoration of a bottom wall of the channel by superhydrophobic grooves provides additional possibilities for manipulation of neutrally-buoyant particles. It is shown that the effective anisotropic hydrodynamic slip of such a bottom wall can be readily used to alter the equilibrium positions of particles and to generate their motion transverse to the pressure gradient. These results may guide the design of novel inertial microfluidic devices for efficient sorting of neutrally-buoyant microparticles by their size.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Predicting optimal hematocrit in silico

Optimal hematocrit $H_o$ maximizes oxygen transport. In healthy humans, the average hematocrit $H$ is in the range of 40-45$\%$, but it can significantly change in blood pathologies such as severe anemia (low $H$) and polycythemia (high $H$). Whether the hematocrit level in humans corresponds to the optimal one is a long standing physiological question. Here, using numerical simulations with the Lattice Boltzmann method and two mechanical models of the red blood cell (RBC) we predict the optimal hematocrit, and explore how altering the mechanical properties of RBCs affects $H_o$. We develop a simplified analytical theory that accounts for results obtained from numerical simulations and provides insight into the physical mechanisms determining $H_o$. Our numerical and analytical models can easily be modified to incorporate a wide range of mechanical properties of RBCs as well as other soft particles thereby providing means for the rational design of blood substitutes. Our work lays the foundations for systematic theoretical study of the optimal hematocrit and its link with pathological RBCs associated with various diseases (e.g. sickle cell anemia, diabetes mellitus, malaria, elliptocytosis).

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Active elastohydrodynamics of vesicles in narrow, blind constrictions

Fluid-resistance limited transport of vesicles through narrow constrictions is a recurring theme in many biological and engineering applications. Inspired by the motor-driven movement of soft membrane-bound vesicles into closed neuronal dendritic spines, here we study this problem using a combination of passive three-dimensional simulations and a simplified semi-analytical theory for active transport of vesicles that are forced through such constrictions by molecular motors. We show that the motion of these objects is characterized by two dimensionless quantities related to the geometry and the strength of forcing relative to the vesicle elasticity. We use numerical simulations to characterize the transit time for a vesicle forced by fluid pressure through a constriction in a channel, and find that relative to an open channel, transport into a blind end leads to the formation of an effective lubrication layer that strongly impedes motion. When the fluid pressure forcing is complemented by forces due to molecular motors that are responsible for vesicle trafficking into dendritic spines, we find that the competition between motor forcing and fluid drag results in multistable dynamics reminiscent of the real system. Our study highlights the role of non-local hydrodynamic effects in determining the kinetics of vesicular transport in constricted geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Agglomeration and filtration of colloidal suspensions with DVLO interactions in simulation and experiment

Cake filtration is a widely used solid-liquid separation process. However, the high flow resistance of the nanoporous filter cake lowers the efficiency of the process significantly. The structure and thus the permeability of the filter cakes depend on the compressive load acting on the particles, the particles size, and the agglomeration of the particles. The latter is determined by the particle charge and the ionic strength of the suspension, as described by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. In this paper, we propose a combined stochastic rotation dynamics (SRD) and molecular dynamics (MD) methodology to simulate the cake formation. The simulations give further insight into the dependency of the filter cakes' structure on the agglomeration of the particles, which cannot be accessed experimentally. The permeability, as investigated with lattice Boltzmann (LB) simulations of flow through the discretized cake, depends on the particle size and porosity, and thus on the agglomeration of the particles. Our results agree qualitatively with experimental data obtained from colloidal boehmite suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random-roughness hydrodynamic boundary conditions

We report results of lattice Boltzmann simulations of a high-speed drainage of liquid films squeezed between a smooth sphere and a randomly rough plane. A significant decrease in the hydrodynamic resistance force as compared with that predicted for two smooth surfaces is observed. However, this force reduction does not represent slippage. The computed force is exactly the same as that between equivalent smooth surfaces obeying no-slip boundary conditions, but located at an intermediate position between peaks and valleys of asperities. The shift in hydrodynamic thickness is shown to depend on the height and density of roughness elements. Our results do not support some previous experimental conclusions on very large and shear-dependent boundary slip for similar systems.

0 citations0 reviewsNo code linkedOpen access