Showing 1–50 of 110 results for author: Harting, J

The interplay of shape and catalyst distribution in the yield of compressible flow microreactors

Authors: G. C. Antunes, J. Jiménez-Sanchéz, P. Malgaretti, J. Bachmann, J. Harting

Abstract: We develop a semi-analytical model for transport in structured catalytic microreactors, where both reactant and product are compressible fluids. Making use of the lubrication and Fick-Jacobs approximations, we reduce the three-dimensional governing equations to an effective one-dimensional set of equations. Our model captures the effect of compressibility, of corrugations in the shape of the react… ▽ More We develop a semi-analytical model for transport in structured catalytic microreactors, where both reactant and product are compressible fluids. Making use of the lubrication and Fick-Jacobs approximations, we reduce the three-dimensional governing equations to an effective one-dimensional set of equations. Our model captures the effect of compressibility, of corrugations in the shape of the reactor, as well as of an inhomogeneous catalytic coating of the reactor walls. We show that in the weakly compressible limit (e.g., liquid-phase reactors), the distribution of catalyst does not influence the reactor yield, which we verify experimentally. Beyond this limit, we show that introducing inhomogeneities in the catalytic coating and corrugations to the reactor walls can improve the yield. △ Less

Submitted 8 July, 2024; v1 submitted 4 July, 2024; originally announced July 2024.

A lattice Boltzmann approach for acoustic manipulation

Authors: E. Castro-Avila, P. Malgaretti, J. Harting, J. D. Muñoz

Abstract: We employ a lattice Boltzmann method to compute the acoustic radiation force produced by standing waves on a compressible object. Instead of simulating the fluid mechanics equations directly, the proposed method uses a lattice Boltzmann model that reproduces the wave equation, together with a kernel interpolation scheme, to compute the first order perturbations of the pressure and velocity fields… ▽ More We employ a lattice Boltzmann method to compute the acoustic radiation force produced by standing waves on a compressible object. Instead of simulating the fluid mechanics equations directly, the proposed method uses a lattice Boltzmann model that reproduces the wave equation, together with a kernel interpolation scheme, to compute the first order perturbations of the pressure and velocity fields on the object's surface and, from them, the acoustic radiation force. The procedure reproduces with excellent accuracy the theoretical expressions by Gor'kov and Wei for the sphere and the disk, respectively, even with a modest number of lattice Boltzmann cells. The proposed method shows to be a promising tool for simulating phenomena where the acoustic radiation force plays a relevant role, like acoustic tweezers or the acoustic manipulation of microswimmers, with applications in medicine and engineering. △ Less

Submitted 13 July, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Chemically reactive thin films: dynamics and stability

Authors: Tilman Richter, Paolo Malgaretti, Thomas M. Koller, Jens Harting

Abstract: Catalyst particles or complexes suspended in liquid films can trigger chemical reactions leading to inhomogeneous concentrations of reactants and products in the film. We demonstrate that the sensitivity of the liquid film's gas-liquid surface tension to these inhomogeneous concentrations strongly impacts the film stability. Using linear stability analysis, we identify novel scenarios in which the… ▽ More Catalyst particles or complexes suspended in liquid films can trigger chemical reactions leading to inhomogeneous concentrations of reactants and products in the film. We demonstrate that the sensitivity of the liquid film's gas-liquid surface tension to these inhomogeneous concentrations strongly impacts the film stability. Using linear stability analysis, we identify novel scenarios in which the film can be either stabilized or destabilized by the reactions. Furthermore, we find so far unrevealed rupture mechanisms which are absent in the chemically inactive case. The linear stability predictions are confirmed by numerical simulations, which also demonstrate that the shape of chemically active droplets can depart from the spherical cap and that unsteady states such as traveling and standing waves might appear. Finally, we critically discuss the relevance of our predictions by showing that the range of our selected parameters is well accessible by typical experiments. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Crystalline Morphology Formation in Phase-Field Simulations of Binary Mixtures

Authors: Maxime Siber, Olivier J. J. Ronsin, Jens Harting

Abstract: Understanding the morphology formation process of solution-cast photoactive layers (PALs) is crucial to derive design rules for optimized and reliable third generation solar cell fabrication. For this purpose, a Phase-Field (PF) computational framework dedicated to the simulation of PAL processing has recently been developed. In this study focused on non-evaporating, crystallizing binary mixtures,… ▽ More Understanding the morphology formation process of solution-cast photoactive layers (PALs) is crucial to derive design rules for optimized and reliable third generation solar cell fabrication. For this purpose, a Phase-Field (PF) computational framework dedicated to the simulation of PAL processing has recently been developed. In this study focused on non-evaporating, crystallizing binary mixtures, distinct crystalline morphology formation pathways are characterized by a systematic exploration of the model's parameter space. It is identified how, depending on material properties, regular, dilution-enhanced, diffusion-limited and demixing-assisted crystallization can take place, and which associated structures then arise. A comprehensive description of the thermodynamic and kinetic mechanisms that respectively drive these separate crystallization modes is provided. Finally, comparisons with experimental results reported in the literature highlight the promising potential of PF simulations to support the determination of process-structure relationships for improved PAL production. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Turning catalytically active pores into active pumps

Authors: G. C. Antunes, P. Malgaretti, J. Harting

Abstract: We develop a semi-analytical model of self-diffusioosmotic transport in active pores, which includes advective transport and the inverse chemical reaction which consumes solute. In previous work (Phys. Rev. Lett. 129, 188003, 2022), we have demonstrated the existence of a spontaneous symmetry breaking in fore-aft symmetric pores that enables them to function as a micropump. We now show that this p… ▽ More We develop a semi-analytical model of self-diffusioosmotic transport in active pores, which includes advective transport and the inverse chemical reaction which consumes solute. In previous work (Phys. Rev. Lett. 129, 188003, 2022), we have demonstrated the existence of a spontaneous symmetry breaking in fore-aft symmetric pores that enables them to function as a micropump. We now show that this pumping transition is controlled by three timescales. Two timescales characterize advective and diffusive transport. The third timescale corresponds to how long a solute molecule resides in the pore before being consumed. Introducing asymmetry to the pore (either via the shape or the catalytic coating) reveals a second type of advection-enabled transitions. In asymmetric pores, the flow rate exhibits discontinuous jumps and hysteresis loops upon tuning the parameters that control the asymmetry. This work demonstrates the interconnected roles of shape and catalytic patterning in the dynamics of active pores, and shows how to design a pump for optimum performance. △ Less

Submitted 5 October, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Journal ref: J. Chem. Phys. 159, 134903 (2023)

A sharp interface approach for wetting dynamics of coated droplets and soft particles

Authors: Francesca Pelusi, Fabio Guglietta, Marcello Sega, Othmane Aouane, Jens Harting

Abstract: The wetting dynamics of liquid particles, from coated droplets to soft capsules, holds significant technological interest. Motivated by the need to simulate liquid metal droplet with an oxidize surface layer, in this work we introduce a computational scheme that allows to simulate droplet dynamics with general surface properties and model different levels of interface stiffness, describing also ca… ▽ More The wetting dynamics of liquid particles, from coated droplets to soft capsules, holds significant technological interest. Motivated by the need to simulate liquid metal droplet with an oxidize surface layer, in this work we introduce a computational scheme that allows to simulate droplet dynamics with general surface properties and model different levels of interface stiffness, describing also cases that are intermediate between pure droplets and capsules. Our approach is based on a combination of the immersed boundary (IB) and the lattice Boltzmann (LB) methods. Here, we validate our approach against the theoretical predictions in the context of shear flow and static wetting properties and we show its effectiveness in accessing the wetting dynamics, exploring the ability of the scheme to address a broad phenomenology. △ Less

Submitted 22 July, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 14 pages, 6 figures

Journal ref: Physics of Fluids 35, 082126 (2023)

Cell-free layer development and spatial organization of healthy and rigid red blood cells in a microfluidic bifurcation

Authors: Yazdan Rashidi, Othmane Aouane, Alexis C. Darras, Thomas John, Jens Harting, Christian Wagner, Steffen M. Recktenwald

Abstract: Bifurcations and branches in the microcirculation dramatically affect blood flow as they determine the spatiotemporal organization of red blood cells (RBCs). Such changes in vessel geometries can further influence the formation of a cell-free layer (CFL) close to the vessel walls. Biophysical cell properties, such as their deformability, which is impaired in various diseases, are often thought to… ▽ More Bifurcations and branches in the microcirculation dramatically affect blood flow as they determine the spatiotemporal organization of red blood cells (RBCs). Such changes in vessel geometries can further influence the formation of a cell-free layer (CFL) close to the vessel walls. Biophysical cell properties, such as their deformability, which is impaired in various diseases, are often thought to impact blood flow and affect the distribution of flowing RBCs. This study investigates the flow behavior of healthy and artificially hardened RBCs in a bifurcating microfluidic T-junction. We determine the RBC distribution across the channel width at multiple positions before and after the bifurcation. Thus, we reveal distinct focusing profiles in the feeding mother channel for rigid and healthy RBCs that dramatically impact the cell organization in the successive daughter channels. Moreover, we experimentally show how the characteristic asymmetric CFLs in the daughter vessels develop along their flow direction. Complimentary numerical simulations indicate that the buildup of the CFL is faster for healthy than for rigid RBCs. Our results provide fundamental knowledge to understand the partitioning of rigid RBC as a model of cells with pathologically impaired deformability in complex in vitro networks. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Suspensions of viscoelastic capsules: effect of membrane viscosity on transient dynamics

Authors: Fabio Guglietta, Francesca Pelusi, Marcello Sega, Othmane Aouane, Jens Harting

Abstract: Membrane viscosity is known to play a central role in the transient dynamics of isolated viscoelastic capsules by decreasing their deformation, inducing shape oscillations and reducing the loading time, that is, the time required to reach the steady-state deformation. However, for dense suspensions of capsules, our understanding of the influence of the membrane viscosity is minimal. In this work,… ▽ More Membrane viscosity is known to play a central role in the transient dynamics of isolated viscoelastic capsules by decreasing their deformation, inducing shape oscillations and reducing the loading time, that is, the time required to reach the steady-state deformation. However, for dense suspensions of capsules, our understanding of the influence of the membrane viscosity is minimal. In this work, we perform a systematic numerical investigation based on coupled immersed boundary -- lattice Boltzmann (IB-LB) simulations of viscoelastic spherical capsule suspensions in the non-inertial regime. We show the effect of the membrane viscosity on the transient dynamics as a function of volume fraction and capillary number. Our results indicate that the influence of membrane viscosity on both deformation and loading time strongly depends on the volume fraction in a non-trivial manner: dense suspensions with large surface viscosity are more resistant to deformation but attain loading times that are characteristic of capsules with no surface viscosity, thus opening the possibility to obtain richer combinations of mechanical features. △ Less

Submitted 3 July, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Viscous to inertial coalescence of liquid lenses: a lattice Boltzmann investigation

Authors: Thomas Scheel, Qingguang Xie, Marcello Sega, Jens Harting

Abstract: Liquid lens coalescence is an important mechanism involved in many industrial and scientific applications. It has been investigated both theoretically and experimentally, yet it is numerically very challenging to obtain consistent results over the wide ranges of surface tension and viscosity values that are necessary to capture the asymptotic temporal behavior in the viscous and inertial limits. W… ▽ More Liquid lens coalescence is an important mechanism involved in many industrial and scientific applications. It has been investigated both theoretically and experimentally, yet it is numerically very challenging to obtain consistent results over the wide ranges of surface tension and viscosity values that are necessary to capture the asymptotic temporal behavior in the viscous and inertial limits. We report results of massively parallel simulations based on the color gradient lattice Boltzmann method, which overcome these limitations, and investigate the scaling laws of both regimes. For the two-dimensional case we find good agreement with the similarity solution of the thin-sheet equation, where in the viscous regime the connecting bridge grows linearly with time and in the inertial regime proportionally to $t^{2/3}$. In three dimensions, the viscous growth of the bridge also exhibits a linear time dependence, while in the inertial regime the growth of both the bridge height and the bridge width is proportional to $t^{1/2}$. △ Less

Submitted 14 June, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Comments: 11 pages, 8 figures

Inertial focusing of a dilute suspension in pipe flow

Authors: Othmane Aouane, Marcello Sega, Bastian Bäuerlein, Kerstin Avila, Jens Harting

Abstract: The dynamics of rigid particle suspensions in a wall-bounded laminar flow present several non-trivial and intriguing features, including particle ordering, lateral transport, and the appearance of stable, preferential locations like the Segré-Silberberg annulus. The formation of more than one annulus is a particularly puzzling phenomenon that is still not fully explained. Here, we present numerica… ▽ More The dynamics of rigid particle suspensions in a wall-bounded laminar flow present several non-trivial and intriguing features, including particle ordering, lateral transport, and the appearance of stable, preferential locations like the Segré-Silberberg annulus. The formation of more than one annulus is a particularly puzzling phenomenon that is still not fully explained. Here, we present numerical simulation results of a dilute suspension of particles in (periodic) pipe flow based on the lattice Boltzmann and the discrete element methods (DEM). Our simulations provide access to the full radial position history of the particles while traveling downstream. This allows to accurately quantify the transient and steady states. We observe the formation of the secondary, inner annulus and show that its position invariably shifts toward the Segré-Silberberg one if the channel is sufficiently long, proving that it is, in fact, a transient feature for Reynolds numbers (Re) up to 600. We quantify the variation of the channel focusing length ($L_s/2R$) with Re. Interestingly and unlike the theoretical prediction for a point-like particle, we observe that $L_s/2R$ increases with Re for both the single particle and the suspension. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Formation of crystalline bulk heterojunctions in organic solar cells: insights from phase-field simulations

Authors: Olivier Ronsin, Jens Harting

Abstract: The performance of organic solar cells strongly depends on the bulk heterojunction (BHJ) morphology of the photoactive layer. This BHJ forms during the drying of the wet-deposited solution, because of physical processes such as crystallization and/or liquid liquid phase separation (LLPS). However, the process-structure relationship remains insufficiently understood. In this work, a recently develo… ▽ More The performance of organic solar cells strongly depends on the bulk heterojunction (BHJ) morphology of the photoactive layer. This BHJ forms during the drying of the wet-deposited solution, because of physical processes such as crystallization and/or liquid liquid phase separation (LLPS). However, the process-structure relationship remains insufficiently understood. In this work, a recently developed, coupled phase field fluid mechanics framework is used to simulate the BHJ formation upon drying. For the first time, this allows to investigate the interplay between all the relevant physical processes (evaporation, crystal nucleation and growth, liquid demixing, composition-dependent kinetic properties), within a single coherent theoretical framework. Simulations for the model system P3HT-PCBM are presented. The comparison with previously reported in-situ characterization of the drying structure is very convincing: the morphology formation pathways, crystallization kinetics, and final morphology are in line with experimental results. The final BHJ morphology is a subtle mixture of pure crystalline donor and acceptor phases, pure and mixed amorphous domains, which depends on the process parameters and material properties. The expected benefit of such an approach is to identify physical design rules for ink formulation and processing conditions to optimize the cell performance. It could be applied to recent organic material systems in the future. △ Less

Submitted 28 September, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

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

Authors: Olivier J. J. Ronsin, Jens Harting

Abstract: 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 m… ▽ More 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. △ Less

Submitted 27 June, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Liquid film rupture beyond the thin-film equation: a multi-component lattice Boltzmann study

Authors: Francesca Pelusi, Marcello Sega, Jens Harting

Abstract: Under the condition of partial surface wettability, thin liquid films can be destabilized by small perturbations and rupture into droplets. As successfully predicted by the thin film equation (TFE), the rupture dynamics are dictated by the liquid-solid interaction. The theory describes the latter using the disjoining pressure or, equivalently, the contact angle. The introduction of a secondary flu… ▽ More Under the condition of partial surface wettability, thin liquid films can be destabilized by small perturbations and rupture into droplets. As successfully predicted by the thin film equation (TFE), the rupture dynamics are dictated by the liquid-solid interaction. The theory describes the latter using the disjoining pressure or, equivalently, the contact angle. The introduction of a secondary fluid can lead to a richer phenomenology thanks to the presence of different fluid/surface interaction energies but has so far not been investigated. In this work, we study the rupture of liquid films with different heights immersed in a secondary fluid using a multi-component lattice Boltzmann (LB) approach. We investigate a wide range of surface interaction energies, equilibrium contact angles, and film thicknesses. We found that the rupture time can differ by about one order of magnitude for identical equilibrium contact angles but different surface free energies. Interestingly, the TFE describes the observed breakup dynamics qualitatively well, up to equilibrium contact angles as large as 130$^\circ$. A small film thickness is a much stricter requirement for the validity of the TFE, and agreement with LB results is found only for ratios $ε=h/L$ of the film height $h$ and lateral system size $L$ such as $ε\lesssim\times10^{-3}$. △ Less

Submitted 16 June, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Journal ref: Physics of Fluids 34, 062109 (2022)

Lattice Boltzmann simulations of two linear microswimmers using the immersed boundary method

Authors: D. Geyer, S. Ziegler, A. Sukhov, M. Hubert, A. -S. Smith, O. Aouane, P. Malgaretti, J. Harting

Abstract: The performance of a single or the collection of microswimmers strongly depends on the hydrodynamic coupling among their constituents and themselves. We present a numerical study for a single and a pair of microswimmers based on lattice Boltzmann method (LBM) simulations. Our numerical algorithm consists of two separable parts. Lagrange polynomials provide a discretization of the microswimmers and… ▽ More The performance of a single or the collection of microswimmers strongly depends on the hydrodynamic coupling among their constituents and themselves. We present a numerical study for a single and a pair of microswimmers based on lattice Boltzmann method (LBM) simulations. Our numerical algorithm consists of two separable parts. Lagrange polynomials provide a discretization of the microswimmers and the lattice Boltzmann method captures the dynamics of the surrounding fluid. The two components couple via an immersed boundary method. We present data for a single swimmer system and our data also show the onset of collective effects and, in particular, an overall velocity increment of clusters of swimmers. △ Less

Submitted 14 March, 2023; v1 submitted 24 February, 2022; originally announced February 2022.

Journal ref: Commun. Comput. Phys., 33 (2023), pp. 39-56

Controlling the dewetting morphologies of thin liquid films by switchable substrates

Authors: Stefan Zitz, Andrea Scagliarini, Jens Harting

Abstract: 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, wher… ▽ More 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. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Fluids 8 (2023) L122001

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

Authors: Moritz Lehmann, Mathias J. Krause, Giorgio Amati, Marcello Sega, Jens Harting, Stephan Gekle

Abstract: 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 car… ▽ More 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. △ Less

Submitted 31 January, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 30 pages, 20 figures, 4 tables, 2 code listings

Squeezing multiple soft particles into a constriction: transition to clogging

Authors: Clément Bielinski, Othmane Aouane, Jens Harting, Badr Kaoui

Abstract: We study numerically how multiple deformable capsules squeeze into a constriction. This situation is largely encountered in microfluidic chips designed to manipulate living cells, which are soft entities. We use fully three-dimensional simulations based on the lattice Boltzmann method to compute the flow of the suspending fluid, and on the immersed boundary method to achieve the two-way fluid-stru… ▽ More We study numerically how multiple deformable capsules squeeze into a constriction. This situation is largely encountered in microfluidic chips designed to manipulate living cells, which are soft entities. We use fully three-dimensional simulations based on the lattice Boltzmann method to compute the flow of the suspending fluid, and on the immersed boundary method to achieve the two-way fluid-structure interaction. The mechanics of the capsule membrane elasticity is computed with the finite element method. We obtain two main states: continuous passage of the particles, and their blockage that leads to clogging the constriction. The transition from one state to another is dictated by the ratio between the size of the capsules and the constriction width, and by the capsule membrane deformability. This latter is found to enhance particle passage through narrower constrictions, where rigid particles with similar diameter are blocked and lead to clogging. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 8 pages, 5 figures

Capillary interactions, aggregate formation and the rheology of particle-laden flows: a lattice Boltzmann study

Authors: Lei Yang, Marcello Sega, Steffen Leimbach, Sebastian Kolb, Jürgen Karl, Jens Harting

Abstract: The agglomeration of particles caused by the formation of capillary bridges has a decisive impact on the transport properties of a variety of at a first sight very different systems such as capillary suspensions, fluidized beds in chemical reactors, or even sand castles. Here, we study the connection between the microstructure of the agglomerates and the rheology of fluidized suspensions using a c… ▽ More The agglomeration of particles caused by the formation of capillary bridges has a decisive impact on the transport properties of a variety of at a first sight very different systems such as capillary suspensions, fluidized beds in chemical reactors, or even sand castles. Here, we study the connection between the microstructure of the agglomerates and the rheology of fluidized suspensions using a coupled lattice Boltzmann and discrete element method approach. We address the influence of the shear rate, the secondary fluid surface tension, and the suspending liquid viscosity. The presence of capillary interactions promotes the formation of either filaments or globular clusters, leading to an increased suspension viscosity. Unexpectedly, filaments have the opposite effect on the viscosity as compared to globular clusters, decreasing the suspension viscosity at larger capillary interaction strengths. In addition, we show that the suspending fluid viscosity also has a non-trivial influence on the effective viscosity of the suspension, a fact usually not taken into account by empirical models. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Instability of particle inertial migration in shear flow

Authors: Evgeny S. Asmolov, Tatiana V. Nizkaya, Jens Harting, Olga I. Vinogradova

Abstract: In a shear flow particles migrate to their equilibrium positions in the microchannel. Here we demonstrate theoretically that if particles are inertial, this equilibrium can become unstable due to the Saffman lift force. We derive an expression for the critical Stokes number that determines the onset of instable equilibrium. We also present results of lattice Boltzmann simulations for spherical par… ▽ More In a shear flow particles migrate to their equilibrium positions in the microchannel. Here we demonstrate theoretically that if particles are inertial, this equilibrium can become unstable due to the Saffman lift force. We derive an expression for the critical Stokes number that determines the onset of instable equilibrium. We also present results of lattice Boltzmann simulations for spherical particles and prolate spheroids to validate the analysis. Our work provides a simple explanation of several unusual phenomena observed in earlier experiments and computer simulations, but never interpreted before in terms of the unstable equilibrium. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 9 pages, 8 figures

Two-dimensional Cahn-Hilliard simulations for coarsening kinetics of spinodal decomposition in binary mixtures

Authors: Björn König, Olivier J. J. Ronsin, Jens Harting

Abstract: The evolution of the microstructure due to spinodal decomposition in phase separated mixtures has a strong impact on the final material properties. In the late stage of coarsening, the system is characterized by the growth of a single characteristic length scale $L\sim C t^α$. To understand the structure-property relationship, the knowledge of the coarsening exponent $α$ and the coarsening rate co… ▽ More The evolution of the microstructure due to spinodal decomposition in phase separated mixtures has a strong impact on the final material properties. In the late stage of coarsening, the system is characterized by the growth of a single characteristic length scale $L\sim C t^α$. To understand the structure-property relationship, the knowledge of the coarsening exponent $α$ and the coarsening rate constant $C$ is mandatory. Since the existing literature is not entirely consistent, we perform phase field simulations based on the Cahn-Hilliard equation. We restrict ourselves to binary mixtures using a symmetric Flory-Huggins free energy and a constant mobility term and show that the coarsening for off-critical mixtures is slower than the expected $t^{1/3}$-growth. Instead, we find $α$ to be dependent on the mixture composition and thus from the morphology. Finally, we propose a model to describe the complete coarsening kinetics including the rate constant $C$. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Journal ref: Phys. Chem. Chem. Phys., 2021,23, 24823-24833

Phase-field simulation of liquid-vapor equilibrium and evaporation of fluid mixtures

Authors: Olivier J. J. Ronsin, DongJu Jang, Hans-Joachim Egelhaaf, Christoph J. Brabec, Jens Harting

Abstract: In solution-processing of thin films, the material layer is deposited from a solution composed of several solutes and solvents. The final morphology and hence the properties of the film often depend on the time needed for the evaporation of the solvents. This is typically the case for organic photoactive or electronic layers. Therefore, it is important to be able to predict the evaporation kinetic… ▽ More In solution-processing of thin films, the material layer is deposited from a solution composed of several solutes and solvents. The final morphology and hence the properties of the film often depend on the time needed for the evaporation of the solvents. This is typically the case for organic photoactive or electronic layers. Therefore, it is important to be able to predict the evaporation kinetics of such mixtures. We propose here a new phase-field model for the simulation of evaporating fluid mixtures and simulate their evaporation kinetics. Similar to the Hertz-Knudsen theory, the local liquid-vapor equilibrium is assumed to be reached at the film surface and evaporation is driven by diffusion away from this gas layer. In the situation where the evaporation is purely driven by the liquid-vapor equilibrium, the simulations match the behavior expected theoretically from the free energy: for evaporation of pure solvents, the evaporation rate is constant and proportional to the vaporpressure. For mixtures, the evaporation rate is in general strongly time-dependent because of the changing composition of the film. Nevertheless, for highly non-ideal mixtures, such as poorly compatible fluids or polymer solutions, the evaporation rate becomes almost constant in the limit of low Biot numbers. The results of the simulation have been successfully compared to experiments on a polystyrene-toluene mixture. The model allows to take into account deformations of the liquid-vapor interface and therefore to simulate film roughness or dewetting. △ Less

Submitted 9 November, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Phoretic colloids close to and trapped at fluid interfaces

Authors: Paolo Malgaretti, Jens Harting

Abstract: The active motion of phoretic colloids leads them to accumulate at boundaries and interfaces. Such an excess accumulation, with respect to their passive counterparts, makes the dynamics of phoretic colloids particularly sensitive to the presence of boundaries and pave new routes to externally control their single particle as well as collective behavior. Here we review some recent theoretical resul… ▽ More The active motion of phoretic colloids leads them to accumulate at boundaries and interfaces. Such an excess accumulation, with respect to their passive counterparts, makes the dynamics of phoretic colloids particularly sensitive to the presence of boundaries and pave new routes to externally control their single particle as well as collective behavior. Here we review some recent theoretical results about the dynamics of phoretic colloids close to and adsorbed at fluid interfaces in particular highlighting similarities and differences with respect to solid-fluid interfaces. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Lattice Boltzmann simulations of stochastic thin film dewetting

Authors: Stefan Zitz, Andrea Scagliarini, Jens Harting

Abstract: We study numerically the effect of thermal fluctuations and of variable fluid-substrate interactions on the spontaneous dewetting of thin liquid films. To this aim, we use a recently developed lattice Boltzmann method for thin liquid film flows, equipped with a properly devised stochastic term. While it is known that thermal fluctuations yield shorter rupture times, we show that this is a general… ▽ More We study numerically the effect of thermal fluctuations and of variable fluid-substrate interactions on the spontaneous dewetting of thin liquid films. To this aim, we use a recently developed lattice Boltzmann method for thin liquid film flows, equipped with a properly devised stochastic term. While it is known that thermal fluctuations yield shorter rupture times, we show that this is a general feature of hydrophilic substrates, irrespective of the contact angle. The ratio between deterministic and stochastic rupture times, though, decreases with $θ$. Finally, we discuss the case of fluctuating thin film dewetting on chemically patterned substrates and its dependence on the form of the wettability gradients. △ Less

Submitted 11 August, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 11 pages, 11 figures

Journal ref: Phys. Rev. E 104, 034801 (2021)

Capillary-bridge Forces Between Solid Particles: Insights from Lattice Boltzmann Simulations

Authors: Lei Yang, Marcello Sega, Jens Harting

Abstract: Liquid capillary-bridge formation between solid particles has a critical influence on the rheological properties of granular materials and, in particular, on the efficiency of fluidized bed reactors. The available analytical and semi-analytical methods have inherent limitations, and often do not cover important aspects, like the presence of non-axisymmetric bridges. Here, we conduct numerical simu… ▽ More Liquid capillary-bridge formation between solid particles has a critical influence on the rheological properties of granular materials and, in particular, on the efficiency of fluidized bed reactors. The available analytical and semi-analytical methods have inherent limitations, and often do not cover important aspects, like the presence of non-axisymmetric bridges. Here, we conduct numerical simulations of the capillary bridge formation between equally and unequally-sized solid particles using the lattice Boltzmann method, and provide an assessment of the accuracy of different families of analytical models. We find that some of the models taken into account are shown to perform better than others. However, all of them fail to predict the capillary force for contact angles larger than $π/2$, where a repulsive capillary force attempts to push the solid particle outwards to minimize the surface energy, especially at a small separation distance. △ Less

Submitted 16 June, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 15

Journal ref: AIChE Journal 67, e17350 (2021)

Regimes of motion of magnetocapillary swimmers

Authors: Alexander Sukhov, Maxime Hubert, Galien Grosjean, Oleg Trosman, Sebastian Ziegler, Ylona Collard, Nicolas Vandewalle, Ana-Suncana Smith, Jens Harting

Abstract: The dynamics of a triangular magnetocapillary swimmer is studied using the lattice Boltzmann method. Performing extensive numerical simulations taking into account the coupled dynamics of the fluid-fluid interface and of magnetic particles floating on it and driven by external magnetic fields we identify several regimes of the swimmer motion. In the regime of high frequencies the swimmer's maximum… ▽ More The dynamics of a triangular magnetocapillary swimmer is studied using the lattice Boltzmann method. Performing extensive numerical simulations taking into account the coupled dynamics of the fluid-fluid interface and of magnetic particles floating on it and driven by external magnetic fields we identify several regimes of the swimmer motion. In the regime of high frequencies the swimmer's maximum velocity is centered around the particle's inverse coasting time. Modifying the ratio of surface tension and magnetic forces allows to study the swimmer propagation in the regime of significantly lower frequencies mainly defined by the strength of the magnetocapillary potential. Finally, introducing a constant magnetic contribution in each of the particles in addition to their magnetic moment induced by external fields leads to another regime characterised by strong in-plane swimmer reorientations that resemble experimental observations. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: 12 pages, 12 figures

Journal ref: Eur. Phys. J. E 44, 59 (2021)

Transport of neutral and charged nanorods across varying-section channels

Authors: Paolo Malgaretti, Jens Harting

Abstract: We study the dynamics of neutral and charged rods embedded in varying-section channels. By means of systematic approximations, we derive the dependence of the local diffusion coefficient on both the geometry and charge of the rods. This microscopic insight allows us to provide predictions for the permeability of varying-section channels to rods with diverse lengths, aspect ratios and charge. Our a… ▽ More We study the dynamics of neutral and charged rods embedded in varying-section channels. By means of systematic approximations, we derive the dependence of the local diffusion coefficient on both the geometry and charge of the rods. This microscopic insight allows us to provide predictions for the permeability of varying-section channels to rods with diverse lengths, aspect ratios and charge. Our analysis shows that the dynamics of charged rods is sensitive to the geometry of the channel and that their transport can be controlled by tuning both the shape of the confining walls and the charge of the rod. Interestingly, we find that the channel permeability does not depend monotonically on the charge of the rod. This opens the possibility of a novel mechanism to separate charged rods. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Lattice Boltzmann simulations of drying suspensions of soft particles

Authors: Maarten Wouters, Othmane Aouane, Marcello Sega, Jens Harting

Abstract: The ordering of particles in the drying process of a colloidal suspension is crucial in determining the properties of the resulting film. For example, microscopic inhomogeneities can lead to the formation of cracks and defects that can deteriorate the quality of the film considerably. This type of problem is inherently multiscale and here we study it numerically, using our recently developed metho… ▽ More The ordering of particles in the drying process of a colloidal suspension is crucial in determining the properties of the resulting film. For example, microscopic inhomogeneities can lead to the formation of cracks and defects that can deteriorate the quality of the film considerably. This type of problem is inherently multiscale and here we study it numerically, using our recently developed method for the simulation of soft polymeric capsules in multicomponent fluids. We focus on the effect of the particle softness on the film microstructure during the drying phase and how it relates to the formation of defects. We quantify the order of the particles by measuring both the Voronoi entropy and the isotropic order parameter. Surprisingly, both observables exhibit a non-monotonic behaviour when the softness of the particles is increased. We further investigate the correlation between the interparticle interaction and the change in the microstructure during the evaporation phase. We observe that the rigid particles form chain-like structures that tend to scatter into small clusters when the particle softness is increased. △ Less

Submitted 9 June, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Controllable Capillary Assembly of Magnetic Ellipsoidal Janus Particles into Tunable Rings, Chains and Hexagonal Lattices

Authors: Qingguang Xie, Jens Harting

Abstract: Colloidal assembly at fluid interfaces has a great potential for the bottom-up fabrication of novel structured materials. However, challenges remain in realizing controllable and tunable assembly of particles into diverse structures. Herein, we report the capillary assembly of magnetic ellipsoidal Janus particles at a fluid-fluid interface. Depending on their tilt angle, i.e. the angle the particl… ▽ More Colloidal assembly at fluid interfaces has a great potential for the bottom-up fabrication of novel structured materials. However, challenges remain in realizing controllable and tunable assembly of particles into diverse structures. Herein, we report the capillary assembly of magnetic ellipsoidal Janus particles at a fluid-fluid interface. Depending on their tilt angle, i.e. the angle the particle main axis forms with the fluid interface, these particles deform the interface and generate capillary dipoles or hexapoles. Driven by capillary interactions, multiple particles thus assemble into chain-, hexagonal lattice- and ring-like structures, which can be actively controlled by applying an external magnetic field. We predict a field-strength phase diagram in which various structures are present as stable states. Owing to the diversity, controllability, and tunability of assembled structures, magnetic ellipsoidal Janus particles at fluid interfaces could therefore serve as versatile building blocks for novel materials. △ Less

Submitted 27 December, 2020; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 6 pages, 4 figures, + supplementary material

Inertial migration of oblate spheroids in a plane channel

Authors: Tatiana V. Nizkaya, Anna S. Gekova, Jens Harting, Evgeny S. Asmolov, Olga I. Vinogradova

Abstract: We discuss an inertial migration of oblate spheroids in a plane channel, where steady laminar flow is generated by a pressure gradient. Our lattice Boltzmann simulations show that spheroids orient in the flow, so that their minor axis coincides with the vorticity direction (a log-rolling motion). Interestingly, for spheroids of moderate aspect ratios, the equilibrium positions relative to the chan… ▽ More We discuss an inertial migration of oblate spheroids in a plane channel, where steady laminar flow is generated by a pressure gradient. Our lattice Boltzmann simulations show that spheroids orient in the flow, so that their minor axis coincides with the vorticity direction (a log-rolling motion). Interestingly, for spheroids of moderate aspect ratios, the equilibrium positions relative to the channel walls depend only on their equatorial radius $a$. By analysing the inertial lift force we argue that this force is proportional to $a^3b$, where $b$ is the polar radius, and conclude that the dimensionless lift coefficient of the oblate spheroid does not depend on $b$, and is equal to that of the sphere of radius $a$. △ Less

Submitted 3 November, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 7 pages, 8 figures

The Scallop Theorem and Swimming at the Mesoscale

Authors: Maxime Hubert, Oleg Trosman, Ylona Collard, Alexander Sukhov, Jens Harting, Nicolas Vandewalle, Ana-Suncana Smith

Abstract: By synergistically combining modeling, simulation and experiments, we show that there exists a regime of self-propulsion in which the inertia in the fluid dynamics can be separated from that of the swimmer. This is demonstrated by the motion of an asymmetric dumbbell that, despite deforming in a reciprocal fashion, self-propagates in a fluid due to a non-reciprocal Stokesian flow field. The latter… ▽ More By synergistically combining modeling, simulation and experiments, we show that there exists a regime of self-propulsion in which the inertia in the fluid dynamics can be separated from that of the swimmer. This is demonstrated by the motion of an asymmetric dumbbell that, despite deforming in a reciprocal fashion, self-propagates in a fluid due to a non-reciprocal Stokesian flow field. The latter arises from the difference in the coasting times of the two constitutive beads. This asymmetry acts as a second degree of freedom, recovering the scallop theorem at the mesoscopic scale. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: Main article: 4 pages, 3 figures. SI (for methods and calculations): 7 pages, no figure, 1 movie

MSC Class: 76D07

Journal ref: Phys. Rev. Lett. 126, 224501 (2021)

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions

Authors: David Jung, Maximilian Johannes Uttinger, Paolo Malgaretti, Wolfgang Peukert, Johannes Walter, Jens Harting

Abstract: 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… ▽ More 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$. △ Less

Submitted 25 February, 2022; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: 11 pages, 5 figures

Journal ref: Soft Matter 2022

Capillary interactions between soft capsules protruding through thin fluid films

Authors: Maarten Wouters, Othmane Aouane, Marcello Sega, Jens Harting

Abstract: 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… ▽ More 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. △ Less

Submitted 19 October, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Journal ref: Soft Matter 18, 2157-2167 (2022)

Desorption energy of soft particles from a fluid interface

Authors: Hadi Mehrabian, Jacco H. Snoeijer, Jens Harting

Abstract: 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… ▽ More 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. △ Less

Submitted 18 August, 2020; v1 submitted 19 June, 2020; originally announced June 2020.

Comments: 11 pages, 10 figures

Bio-inspired acousto-magnetic microswarm robots with upstream motility

Authors: Daniel Ahmed, David Hauri, Alexander Sukhov, Dubon Rodrigue, Maranta Gian, Jens Harting, Bradley Nelson

Abstract: The ability to propel against flows, i.e., to perform positive rheotaxis, can provide exciting opportunities for applications in targeted therapeutics and non-invasive surgery. To date, no biocompatible technologies exist for navigating microparticles upstream when they are in a background fluid flow. Inspired by many naturally occurring microswimmers such as bacteria, spermatozoa, and plankton th… ▽ More The ability to propel against flows, i.e., to perform positive rheotaxis, can provide exciting opportunities for applications in targeted therapeutics and non-invasive surgery. To date, no biocompatible technologies exist for navigating microparticles upstream when they are in a background fluid flow. Inspired by many naturally occurring microswimmers such as bacteria, spermatozoa, and plankton that utilize the non-slip boundary conditions of the wall to exhibit upstream propulsion, here, we report on the design and characterization of self-assembled microswarms that can execute upstream motility in a combination of external acoustic and magnetic fields. Both acoustic and magnetic fields are safe to humans, non-invasive, can penetrate deeply into the human body, and are well-developed in clinical settings. The combination of both fields can overcome the limitations encountered by single actuation methods. The design criteria of the acoustically-induced reaction force of the microswarms, which is needed to perform rolling-type motion, are discussed. We show quantitative agreement between experimental data and our model that captures the rolling behaviour. The upstream capability provides a design strategy for delivering small drug molecules to hard-to-reach sites and represents a fundamental step toward the realization of micro- and nanosystem-navigation against the blood flow. △ Less

Submitted 26 March, 2021; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 22 pages

Journal ref: Nat. Mach Intel. 3, 116 (2021)

Structure and rheology of suspensions of spherical strain-hardening capsules

Authors: Othmane Aouane, Andrea Scagliarini, Jens Harting

Abstract: 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 su… ▽ More 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. △ Less

Submitted 6 November, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

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

Authors: Olivier J. J. Ronsin, DongJu Jang, Hans-Joachim Egelhaaf, Christoph J. Brabec, Jens Harting

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2020; v1 submitted 15 January, 2020; originally announced January 2020.

Journal ref: Phys. Chem. Chem. Phys., 2020, 22, 6638--6652

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

Authors: Olivier J. J. Ronsin, Jens Harting

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 6 February, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Journal ref: Energy Technol. 2020, 1901468

Self-Similar Liquid Lens Coalescence

Authors: Michiel A. Hack, Walter Tewes, Qingguang Xie, Charu Datt, Kirsten Harth, Jens Harting, Jacco H. Snoeijer

Abstract: 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… ▽ More 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. △ Less

Submitted 9 March, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Journal ref: Phys. Rev. Lett. 124, 194502 (2020)

Numerical simulations of self-diffusiophoretic colloids at fluid interfaces

Authors: T. Peter, P. Malgaretti, N. Rivas, A. Scagliarini, J. Harting, S. Dietrich

Abstract: The dynamics of active colloids is very sensitive to the presence of boundaries and interfaces which therefore can be used to control their motion. Here we analyze the dynamics of active colloids adsorbed at a fluid-fluid interface. By using a mesoscopic numerical approach which relies on an approximated numerical solution of the Navier-Stokes equation, we show that when adsorbed at a fluid interf… ▽ More The dynamics of active colloids is very sensitive to the presence of boundaries and interfaces which therefore can be used to control their motion. Here we analyze the dynamics of active colloids adsorbed at a fluid-fluid interface. By using a mesoscopic numerical approach which relies on an approximated numerical solution of the Navier-Stokes equation, we show that when adsorbed at a fluid interface, an active colloid experiences a net torque even in the absence of a viscosity contrast between the two adjacent fluids. In particular, we study the dependence of this torque on the contact angle of the colloid with the fluid-fluid interface and on its surface properties. We rationalize our results via an approximate approach which accounts for the appearance of a local friction coefficient. By providing insight into the dynamics of active colloids adsorbed at fluid interfaces, our results are relevant for two-dimensional self assembly and emulsion stabilization by means of active colloids. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Inertial migration of neutrally-buoyant particles in superhydrophobic channels

Authors: Tatiana V. Nizkaya, Evgeny S. Asmolov, Jens Harting, Olga I. Vinogradova

Abstract: 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 superhydr… ▽ More 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. △ Less

Submitted 25 December, 2019; v1 submitted 17 August, 2019; originally announced August 2019.

Comments: 10 pages, 10 figures

Journal ref: Phys. Rev. Fluids 5, 014201 (2020)

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

Authors: Sebastian Ziegler, Maxime Hubert, Nicolas Vandewalle, Jens Harting, Ana-Sunčana-Smith

Abstract: 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 s… ▽ More 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. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Comments: 19 pages, 6 figures

Capillary assemblies in a rotating magnetic field

Authors: Galien Grosjean, Maxime Hubert, Ylona Collard, Alexander Sukhov, Jens Harting, Ana-Suncana Smith, Nicolas Vandewalle

Abstract: Small objects floating on a fluid have a tendency to aggregate due to capillary forces. This effect has been used, with the help of a magnetic induction field, to assemble submillimeter metallic spheres into a variety of structures, whose shape and size can be tuned. Under time-varying fields, these assemblies can propel themselves due to a breaking of time reversal symmetry in their adopted shape… ▽ More Small objects floating on a fluid have a tendency to aggregate due to capillary forces. This effect has been used, with the help of a magnetic induction field, to assemble submillimeter metallic spheres into a variety of structures, whose shape and size can be tuned. Under time-varying fields, these assemblies can propel themselves due to a breaking of time reversal symmetry in their adopted shapes. In this article, we study the influence of an in-plane rotation of the magnetic field on these structures. Various rotational modes have been observed with different underlying mechanisms. The magnetic properties of the particles cause them to rotate individually. Dipole-dipole interactions in the assembly can cause the whole structure to align with the field. Finally, non-reciprocal deformations can power the rotation of the assembly. Symmetry plays an important role in the dynamics, as well as the frequency and amplitude of the applied field. Understanding the interplay of these effects is essential, both to explain previous observations and to develop new functions for these assemblies. △ Less

Submitted 13 November, 2019; v1 submitted 16 July, 2019; originally announced July 2019.

Comments: 10 pages, 9 figures

Journal ref: Soft Matter, 2019, 15, 9093 - 9103

Strict equivalence between Maxwell-Stefan and fast-mode theory for multicomponent polymer mixtures

Authors: Olivier J. J. Ronsin, Jens Harting

Abstract: The applicability of theories describing the kinetic evolution of fluid mixtures depends on the underlying physical assumptions. The Maxwell-Stefan equations, widely used for miscible fluids, express forces depending on coupled fluxes. They need to be inverted to recover a Fickian form which is generally impossible analytically. Moreover, the concentration dependence of the diffusivities has to be… ▽ More The applicability of theories describing the kinetic evolution of fluid mixtures depends on the underlying physical assumptions. The Maxwell-Stefan equations, widely used for miscible fluids, express forces depending on coupled fluxes. They need to be inverted to recover a Fickian form which is generally impossible analytically. Moreover, the concentration dependence of the diffusivities has to be modelled, e.g. by the multicomponent Darken equation. Cahn-Hilliard type equations are preferred for immiscible mixtures, whereby different assumptions on the coupling of fluxes lead to the slow-mode and fast-mode theories. For two components, these were derived from the Maxwell-Stefan theory in the past. Here, we prove that the fast-mode theory and the generalized Maxwell-Stefan theory together with the multicomponent Darken equation are strictly equivalent even for multicomponent systems with very different molecular sizes. Our findings allow to reduce the choice of a suitable theory to the most efficient algorithm for solving the underlying equations. △ Less

Submitted 19 July, 2019; v1 submitted 14 June, 2019; originally announced June 2019.

Journal ref: Macromolecules 52, 6035-6044 (2019)

Structural Characterization of an Ionic Liquid in bulk and in nano-confined environment from MD simulations

Authors: Natasa Vucemilovic-Alagic, Radha D. Banhatti, Robert Stepic, Christian R. Wick, Daniel Berger, Mario Gaimann, Andreas Bear, Jens Harting, David M. Smith, Ana-Suncana Smith

Abstract: This article contains data on structural characterization of the [C2Mim][NTf2] in bulk and in nano-confined environment obtained using MD simulations. These data supplement those presented in the paper Insights from Molecular Dynamics Simulations on Structural Organization and Diffusive Dynamics of an Ionic Liquid at Solid and Vacuum Interfaces, where force fields with three different charge metho… ▽ More This article contains data on structural characterization of the [C2Mim][NTf2] in bulk and in nano-confined environment obtained using MD simulations. These data supplement those presented in the paper Insights from Molecular Dynamics Simulations on Structural Organization and Diffusive Dynamics of an Ionic Liquid at Solid and Vacuum Interfaces, where force fields with three different charge methods and three charge scaling factors were used for the analysis of the IL in the bulk, at the interface with the vacuum and the IL film in the contact with a hydroxylated alumina surface. Here, we present details on the construction of the model systems in an extended detailed methods section. Furthermore, for best parametrization, structural and dynamic properties of IL in different environment are studied with certain features presented herein. △ Less

Submitted 15 May, 2019; originally announced May 2019.

Comments: 13 pages, 12 figures. arXiv admin note: substantial text overlap with arXiv:1903.09450

Mesoscale simulation of soft particles with tunable contact angle in multi-component fluids

Authors: Maarten Wouters, Othmane Aouane, Timm Krueger, Jens Harting

Abstract: Soft particles at fluid interfaces play an important role in many aspects of our daily life, such as the food industry, paints and coatings, and medical applications. Analytical methods are not capable of describing the emergent effects of the complex dynamics of suspensions of many soft particles, whereas experiments typically either only capture bulk properties or require invasive methods. Compu… ▽ More Soft particles at fluid interfaces play an important role in many aspects of our daily life, such as the food industry, paints and coatings, and medical applications. Analytical methods are not capable of describing the emergent effects of the complex dynamics of suspensions of many soft particles, whereas experiments typically either only capture bulk properties or require invasive methods. Computational methods are therefore a great tool to complement experimental work. However, an efficient and versatile numerical method is needed to model dense suspensions of many soft particles. In this article we propose a method to simulate soft particles in a multi-component fluid, both at and near fluid-fluid interfaces, based on the lattice Boltzmann method, and characterize the error stemming from the fluid-structure coupling for the particle equilibrium shape when adsorbed onto a fluid-fluid interface. Furthermore, we characterize the influence of the preferential contact angle of the particle surface and the particle softness on the vertical displacement of the center of mass relative to the fluid interface. Finally, we demonstrate the capability of our model by simulating a soft capsule adsorbing onto a fluid-fluid interface with a shear flow parallel to the interface, and the covering of a droplet suspended in another fluid by soft particles with different wettability. △ Less

Submitted 30 August, 2019; v1 submitted 25 April, 2019; originally announced April 2019.

Comments: 16 pages, 14 figures

Journal ref: Phys. Rev. E 100, 033309 (2019)

Optimal motion of triangular magnetocapillary swimmers

Authors: Alexander Sukhov, Sebastian Ziegler, Qingguang Xie, Oleg Trosman, Jayant Pande, Galien Grosjean, Maxime Hubert, Nicolas Vandewalle, Ana-Suncana Smith, Jens Harting

Abstract: A system of ferromagnetic particles trapped at a liquid-liquid interface and subjected to a set of magnetic fields (magnetocapillary swimmers) is studied numerically using a hybrid method combining the pseudopotential lattice Boltzmann method and the discrete element method. After investigating the equilibrium properties of a single, two and three particles at the interface, we demonstrate a contr… ▽ More A system of ferromagnetic particles trapped at a liquid-liquid interface and subjected to a set of magnetic fields (magnetocapillary swimmers) is studied numerically using a hybrid method combining the pseudopotential lattice Boltzmann method and the discrete element method. After investigating the equilibrium properties of a single, two and three particles at the interface, we demonstrate a controlled motion of the swimmer formed by three particles. It shows a sharp dependence of the average center-of-mass speed on the frequency of the time-dependent external magnetic field. Inspired by experiments on magnetocapillary microswimmers, we interpret the obtained maxima of the swimmer speed by the optimal frequency centered around the characteristic relaxation time of a spherical particle. It is also shown that the frequency corresponding to the maximum speed grows and the maximum average speed decreases with increasing inter-particle distances at moderate swimmer sizes. The findings of our lattice Boltzmann simulations are supported by bead-spring model calculations. △ Less

Submitted 28 June, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: 10 pages, 11 figures

Journal ref: J. Chem. Phys. 151, 124707 (2019)

The effect of the liquid layer thickness on the dissolution of immersed surface droplets

Authors: Qingguang Xie, Jens Harting

Abstract: Droplets on a liquid-immersed solid surface are key elements in many applications, such as high-throughput chemical analysis and droplet-templated porous materials. Such surface droplets dissolve when the surrounding liquid is undersaturated and the dissolution process is usually treated analogous to a sessile droplet evaporating in air. Typically, theoretical models predict the mass loss rate of… ▽ More Droplets on a liquid-immersed solid surface are key elements in many applications, such as high-throughput chemical analysis and droplet-templated porous materials. Such surface droplets dissolve when the surrounding liquid is undersaturated and the dissolution process is usually treated analogous to a sessile droplet evaporating in air. Typically, theoretical models predict the mass loss rate of dissolving droplets as a function of droplet geometrical factors (radius, constant angle), and droplet material properties (diffusion constant and densities), where the thickness of the surrounding liquid layer is neglected. Here, we investigate, both numerically and theoretically, the effect of the thickness of the surrounding liquid layer on the lifetime of dissolving surface droplets. We perform $3D$ lattice Boltzmann simulations and find that the lifetime is proportional to the thickness of the liquid layer. We propose a theoretical model based on a quasistatic diffusion equation which confirms the numerical results. Our results offer insight to the fundamental understanding of dissolving surface droplets and can provide valuable guidelines for the design of devices where the droplet lifetime is of importance. △ Less

Submitted 5 July, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 8 pages,7 figures

Journal ref: Soft Matter, 2019,15, 6461-6468

How antagonistic salts cause nematic ordering and behave like diblock copolymers

Authors: David Jung, Nicolas Rivas, Jens Harting

Abstract: We present simulation results and an explanatory theory on how antagonistic salts affect the spinodal decomposition of binary fluid mixtures. We find that spinodal decomposition is arrested and complex structures form only when electrostatic ion-ion interactions are small. In this case fluid and ion concentrations couple and the charge field can be approximated as a polynomial function of the rela… ▽ More We present simulation results and an explanatory theory on how antagonistic salts affect the spinodal decomposition of binary fluid mixtures. We find that spinodal decomposition is arrested and complex structures form only when electrostatic ion-ion interactions are small. In this case fluid and ion concentrations couple and the charge field can be approximated as a polynomial function of the relative fluid concentrations alone. When the solvation energy associated with transfering an ion from one fluid phase to the other is of the order of a few k_BT, the coupled fluid and charge fields evolve according to the Ohta-Kawasaki free energy functional. This allows us to accurately predict structure sizes and reduce the parameter space to two dimensionless numbers. The lamellar structures induced by the presence of antagonistic salt in our simulations exhibit a high degree of nematic ordering and the growth of ordered domains over time seems to follow a power law. This power law carries a time exponent proportional to the salt concentration. We reproduce and interpret neutron scattering data from previous experiments of similar systems. The dissolution of structures at high salt concentrations observed in these experiments agrees with our simulations and we explain it as the result of a vanishing surface tension due to electrostatic contributions. We conclude by presenting preliminary 3D results showing the same morphologies as predicted by the Ohta-Kawasaki model as a function of volume fraction and suggesting that our findings from 2D systems remain valid in 3D. △ Less

Submitted 19 February, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 13 pages, 7 figures

Journal ref: Journal of Chemical Physics 150, 064912 (2019)

Optimal cell transport in straight channels and networks

Authors: Alexander Farutin, Zaiyi Shen, Gael Prado, Vassanti Audemar, Hamid Ez-Zahraouy, Abdelilah Benyoussef, Benoit Polack, Jens Harting, Petia M. Vlahovska, Thomas Podgorski, Gwennou Coupier, Chaouqi Misbah

Abstract: Flux of rigid or soft particles (such as drops, vesicles, red blood cells, etc.) in a channel is a complex function of particle concentration, which depends on the details of induced dissipation and suspension structure due to hydrodynamic interactions with walls or between neighboring particles. Through two-dimensional and three-dimensional simulations and a simple model that reveals the contribu… ▽ More Flux of rigid or soft particles (such as drops, vesicles, red blood cells, etc.) in a channel is a complex function of particle concentration, which depends on the details of induced dissipation and suspension structure due to hydrodynamic interactions with walls or between neighboring particles. Through two-dimensional and three-dimensional simulations and a simple model that reveals the contribution of the main characteristics of the flowing suspension, we discuss the existence of an optimal volume fraction for cell transport and its dependence on the cell mechanical properties. The example of blood is explored in detail, by adopting the commonly used modeling of red blood cells dynamics. We highlight the complexity of optimization at the level of a network, due to the antagonist evolution of local volume fraction and optimal volume fraction with the channels diameter. In the case of the blood network, the most recent results on the size evolution of vessels along the circulatory network of healthy organs suggest that the red blood cell volume fraction (hematocrit) of healthy subjects is close to optimality, as far as transport only is concerned. However, the hematocrit value of patients suffering from diverse red blood cel pathologies may strongly deviate from optimality. △ Less

Submitted 18 October, 2018; originally announced October 2018.

Comments: arXiv admin note: substantial text overlap with arXiv:1802.05353

Journal ref: Phys. Rev. Fluids 3, 103603 (2018)

A lattice Boltzmann method for thin liquid film hydrodynamics

Authors: Stefan Zitz, Andrea Scagliarini, Suryanarayana Maddu Kondaiah, Anton A. Darhuber, Jens Harting

Abstract: We propose a novel approach to the numerical simulation of thin film flows, based on the lattice Boltzmann method. We outline the basic features of the method, show in which limits the expected thin film equations are recovered and perform validation tests. The numerical scheme is applied to the viscous Rayleigh-Taylor instability of a thin film and to the spreading of a sessile drop towards its e… ▽ More We propose a novel approach to the numerical simulation of thin film flows, based on the lattice Boltzmann method. We outline the basic features of the method, show in which limits the expected thin film equations are recovered and perform validation tests. The numerical scheme is applied to the viscous Rayleigh-Taylor instability of a thin film and to the spreading of a sessile drop towards its equilibrium contact angle configuration. We show that the Cox-Voinov law is satisfied, and that the effect of a tunable slip length on the substrate is correctly captured. We address, then, the problem of a droplet sliding on an inclined plane, finding that the Capillary number scales linearly with the Bond number, in agreement with experimental results. At last, we demonstrate the ability of the method to handle heterogenous and complex systems by showcasing the controlled dewetting of a thin film on a chemically structured substrate. △ Less

Submitted 30 August, 2019; v1 submitted 17 October, 2018; originally announced October 2018.

Comments: 10 pages, 9 figures

Journal ref: Phys. Rev. E 100, 033313 (2019)