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

Margination of artificially stiffened red blood cells

Authors: Revaz D. Chachanidze, Othmane Aouane, Jens Harting, Christian Wagner, Marc Leonetti

Abstract: Margination, a fundamental process in which leukocytes migrate from the flowing blood to the vessel wall, is well-documented in physiology. However, it is still an open question on how the differences in cell size and stiffness of white and red cells contribute to this phenomenon. To investigate the specific influence of cell stiffness, we conduct experimental and numerical studies on the segregat… ▽ More Margination, a fundamental process in which leukocytes migrate from the flowing blood to the vessel wall, is well-documented in physiology. However, it is still an open question on how the differences in cell size and stiffness of white and red cells contribute to this phenomenon. To investigate the specific influence of cell stiffness, we conduct experimental and numerical studies on the segregation of a binary mixture of artificially stiffened red blood cells within a suspension of healthy cells. The resulting distribution of stiffened cells within the channel is found to depend on the channel geometry, as demonstrated with slit, rectangular, and cylindrical cross-sections. Notably, an unexpected central peak in the distribution of stiffened RBCs, accompanied by fourfold peaks at the corners, emerges in agreement with simulations. Our results unveil a non-monotonic variation in segregation/margination concerning hematocrit and flow rate, challenging the prevailing belief that higher flow rates lead to enhanced margination. △ Less

Submitted 4 September, 2024; originally announced September 2024.

Comments: Accepted for publication in Phys. Rev. Fluids

Towards a universal law for blood flow

Authors: Alexander Farutin, Abdessamad Nait-Ouhra, Gopal Dixit, Mehdi Abbasi, Othmane Aouane, Jens Harting, Chaouqi Misbah

Abstract: Despite decades of research on blood flow, an analogue of Navier-Stokes equations that accurately describe blood flow properties has not been established yet. The reason behind this is that the properties of blood flow seem à priori non universal as they depend on various factors such as global concentration of red blood cells (RBCs) and channel width. Here, we have discovered a universal law when… ▽ More Despite decades of research on blood flow, an analogue of Navier-Stokes equations that accurately describe blood flow properties has not been established yet. The reason behind this is that the properties of blood flow seem à priori non universal as they depend on various factors such as global concentration of red blood cells (RBCs) and channel width. Here, we have discovered a universal law when the stress and strain rate are measured at a given local RBCs concentration. However, the local concentration must be determined in order to close the problem. We propose a non-local diffusion equation of RBCs concentration that agrees with the full simulation. The universal law is exemplified for both shear and pressure driven flows. While the theory is restricted to a simplistic geometry (straight channel) it provides a fundamental basis for future research on blood flow dynamics and could lead to the development of a new theory that accurately describes blood flow properties under various conditions, such as in complex vascular networks. △ Less

Submitted 25 August, 2024; originally announced August 2024.

Phase field simulations of thermal annealing for all-small molecule organic solar cells

Authors: Yasin Ameslon, Olivier J. J. Ronsin, Christina Harreiss, Johannes Will, Stefanie Rechberger Mingjian Wu, Erdmann Spiecker, Jens Harting

Abstract: Interest in organic solar cells (OSCs) is constantly rising in the field of photovoltaic devices. The device performance relies on the bulk heterojunction (BHJ) nanomorphology, which develops during the drying process and additional post-treatment. This work studies the effect of thermal annealing (TA) on an all-small molecule DRCN5T: PC71 BM blend with phase field simulations. The objective is to… ▽ More Interest in organic solar cells (OSCs) is constantly rising in the field of photovoltaic devices. The device performance relies on the bulk heterojunction (BHJ) nanomorphology, which develops during the drying process and additional post-treatment. This work studies the effect of thermal annealing (TA) on an all-small molecule DRCN5T: PC71 BM blend with phase field simulations. The objective is to determine the physical phenomena driving the evolution of the BHJ morphology for a better understanding of the posttreatment/morphology relationship. Phase-field simulation results are used to investigate the impact on the final BHJ morphology of the DRCN5T crystallization-related mechanisms, including nucleation, growth, crystal stability, impingement, grain coarsening, and Ostwald ripening, of the amorphous-amorphous phase separation (AAPS), and of diffusion limitations. The comparison of simulation results with experimental data shows that the morphological evolution of the BHJ under TA is dominated by dissolution of the smallest, unstable DRCN5T crystals and anisotropic growth of the largest crystals. △ Less

Submitted 6 August, 2024; originally announced August 2024.

Comments: 39 pages incl. SI

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

Authors: G. C. Antunes, M. 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 29 July, 2024; v1 submitted 4 July, 2024; originally announced July 2024.

Simulation of perovskite thin layer crystallization with varying evaporation rates

Authors: M. Majewski, S. Qiu, O. Ronsin, L. Lüer, V. M. Le Corre, T. Du, C. J. Brabec, H. -J. Egelhaaf, J. Harting

Abstract: Perovskite solar cells (PSC) are promising potential competitors to established photovoltaic technologies due to their superior efficiency and low-cost solution processability. However, the limited understanding of the crystallization behaviour hinders the technological transition from lab-scale cells to modules. In this work, we perform Phase Field (PF) simulations of the doctor-bladed film forma… ▽ More Perovskite solar cells (PSC) are promising potential competitors to established photovoltaic technologies due to their superior efficiency and low-cost solution processability. However, the limited understanding of the crystallization behaviour hinders the technological transition from lab-scale cells to modules. In this work, we perform Phase Field (PF) simulations of the doctor-bladed film formation to obtain mechanistic and morphological information that is experimentally challenging to access. PF simulations are validated extensively using in- and ex-situ experiments for different solvent evaporation rates. The well-known transition from a film with many pinholes, for a low evaporation rate, to a smooth film, for high evaporation rates, is recovered in simulation and experiment. From the simulation, the transition can be assigned to the change in the ratio of evaporation to crystallization rate because of two distinct mechanisms. Firstly, for larger evaporation rates, nuclei appear at higher concentrations, which favors nucleation as compared to growth. Secondly, the growth of the crystals is confined in a thinner film, which limits their vertical size. Both effects are expected to be valid independent of the specific chemistry of the chosen experimental system, as long as the evaporation time of the solvent is comparable to the crystallization time. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: paper 10 pages, SI 26 pages

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.

Enhancement of bubble transport in porous electrodes and catalysts

Authors: Thomas Scheel, Paolo Malgaretti, Jens Harting

Abstract: We investigate the formation and transport of gas bubbles across a model porous catalyst/electrode using lattice Boltzmann simulations. This approach enables us to systematically examine the influence of a wide range of morphologies, flow velocities, and reaction rates on the efficiency of gas production. By exploring these parameters, we identify critical parameter combinations that significantly… ▽ More We investigate the formation and transport of gas bubbles across a model porous catalyst/electrode using lattice Boltzmann simulations. This approach enables us to systematically examine the influence of a wide range of morphologies, flow velocities, and reaction rates on the efficiency of gas production. By exploring these parameters, we identify critical parameter combinations that significantly contribute to an enhanced yield of gas output. Our simulations reveal the existence of an optimal pore geometry for which the product output is maximized. Intriguingly, we also observe that lower flow velocities improve gas production by leveraging on coalescence-induced bubble detachment from the catalyst. △ Less

Submitted 26 April, 2024; v1 submitted 1 March, 2024; originally announced March 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.

Electroneutrality breakdown for electrolytes embedded in varying-section nanopores

Authors: P. Malgaretti, I. Pagonabarraga, J. Harting

Abstract: We determine the local charge dynamics of a $z-z$ electrolyte embedded in a varying-section channel. By means of an expansion based on the length scale separation between the axial and transverse direction of the channel, we derive closed formulas for the local excess charge for both, dielectric and conducting walls, in $2D$ (planar geometry) as well as in $3D$ (cylindrical geometry). Our results… ▽ More We determine the local charge dynamics of a $z-z$ electrolyte embedded in a varying-section channel. By means of an expansion based on the length scale separation between the axial and transverse direction of the channel, we derive closed formulas for the local excess charge for both, dielectric and conducting walls, in $2D$ (planar geometry) as well as in $3D$ (cylindrical geometry). Our results show that, even at equilibrium, the local charge electroneutrality is broken whenever the section of the channel is not homogeneous for both dielectric and conducting walls as well as for $2D$ and $3D$ channels. Interestingly, even within our expansion, the local excess charge in the fluid can be comparable to the net charge on the walls. We critically discuss the onset of such local electroneutrality breakdown in particular with respect to the correction that it induces on the effective free energy profile experienced by tracer ions. △ Less

Submitted 15 November, 2023; originally announced November 2023.

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.

Closed formula for the transport of micro-nano-particle across model porous media

Authors: P. Malgaretti, J. Harting

Abstract: In the last decade the Fick-Jacobs approximation has been exploited to capture the transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that for fore-aft symmetric channels the flux of non-interacting systems is fully captured by its linear response regime. For this case we derive a very simple fo… ▽ More In the last decade the Fick-Jacobs approximation has been exploited to capture the transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that for fore-aft symmetric channels the flux of non-interacting systems is fully captured by its linear response regime. For this case we derive a very simple formula that captures the correct trends and that can be exploited as a simple tool to design experiments or simulations. Finally, we show that higher order corrections in the flux may appear for non-symmetric channels. △ Less

Submitted 21 January, 2023; originally announced January 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

Modelling diffusive transport of particles interacting with slit nanopore walls: The case of fullerenes in toluene filled alumina pores

Authors: Andreas Baer, Paolo Malgaretti, Malte Kaspereit, Jens Harting, Ana-Suncana Smith

Abstract: Accurate modeling of diffusive transport of nanoparticles across nanopores is a particularly challenging problem. The reason is that for such narrow pores the large surface-to-volume ratio amplifies the relevance of the nanoscopic details and of the effective interactions at the interface with pore walls. Close to the pore wall, there is no clear separation between the length scales associated wit… ▽ More Accurate modeling of diffusive transport of nanoparticles across nanopores is a particularly challenging problem. The reason is that for such narrow pores the large surface-to-volume ratio amplifies the relevance of the nanoscopic details and of the effective interactions at the interface with pore walls. Close to the pore wall, there is no clear separation between the length scales associated with molecular interactions, layering of the solvent at the interface with the pore and the particle size. Therefore, the standard hydrodynamic arguments may not apply and alternative solutions to determining average transport coefficients need to be developed. We here address this problem by offering a multiscale ansatz that uses effective potentials determined from molecular dynamics simulations to parametrise a four state stochastic model for the positional configuration of the particle in the pore. This is in turn combined with diffusivities in the centre of the pore and at the pore wall to calculate the average diffusion constant. We apply this model to the diffusion of fullerenes in a toluene filled slit nanopore and calculate the mean diffusion coefficient as a function of the pore size. We show that the accuracy of our model is affected by the partial slip of the toluene on the pore wall. △ Less

Submitted 14 November, 2022; originally announced November 2022.

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)

Pumping and mixing in active pores

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

Abstract: We show both numerically and analytically that a chemically patterned active pore can act as a micro/nano-pump for fluids, even if it is fore-aft symmetric. This is possible due to a spontaneous symmetry breaking which occurs when advection rather than diffusion is the dominant mechanism of solute transport. We further demonstrate that, for pumping and tuning the flow rate, a combination of geomet… ▽ More We show both numerically and analytically that a chemically patterned active pore can act as a micro/nano-pump for fluids, even if it is fore-aft symmetric. This is possible due to a spontaneous symmetry breaking which occurs when advection rather than diffusion is the dominant mechanism of solute transport. We further demonstrate that, for pumping and tuning the flow rate, a combination of geometrical and chemical inhomogeneities is required. For certain parameter values, the flow is unsteady, and persistent oscillations with a tunable frequency appear. Finally, we find that the flow may lose its axial symmetry and hence promotes mixing in the low Reynolds number regime. △ Less

Submitted 25 October, 2022; v1 submitted 22 March, 2022; originally announced March 2022.

Comments: 16 pages, 8 figures, 1 video ; additional discussion

Journal ref: Phys. Rev. Lett. 129, 188003 (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.

Monolayer Structure of Supramolecular Antagonistic Salt Aggregates

Authors: David Jung, Jens Harting, Marcello Sega

Abstract: The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh$_4$) in water and 3-methylpyridine (3-M… ▽ More The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh$_4$) in water and 3-methylpyridine (3-MP) at the atomistic level. We show that, indeed, the lamellae are formed by a monomolecular layer of the anion, enveloped by 3-MP and hydrated sodium counterions. To understand which thermodynamic forces drive the aggregation, we compare the full-atomistic model with a simplified one for the salt and show that the strong hydrophobic effect granted by the large excluded volume of the anion, together with electrostatic repulsion, suffice to explain the stability of the monomolecular lamellae. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: accepted for publication on J. Phys. Chem. B (2021)

Theoretical framework for pairwise microswimmer interactions

Authors: Sebastian Ziegler, Thomas Scheel, Maxime Hubert, Jens Harting, Ana-Sunčana Smith

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

Submitted 10 February, 2021; originally announced February 2021.

Comments: 29 pages, 6 figures

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.

Equilibrium orientation and adsorption of an ellipsoidal Janus particle at a fluid-fluid interface

Authors: Florian Günther, Qingguang Xie, Jens Harting

Abstract: We investigate the equilibrium orientation and adsorption process of a single, ellipsoidal Janus particle at a fluid-fluid interface. The particle surface comprises equally sized parts that are hydrophobic or hydrophilic. We present free energy models to predict the equilibrium orientation and compare the theoretical predictions with lattice Boltzmann simulations. We find that the deformation of t… ▽ More We investigate the equilibrium orientation and adsorption process of a single, ellipsoidal Janus particle at a fluid-fluid interface. The particle surface comprises equally sized parts that are hydrophobic or hydrophilic. We present free energy models to predict the equilibrium orientation and compare the theoretical predictions with lattice Boltzmann simulations. We find that the deformation of the fluid interface strongly influences the equilibrium orientation of the Janus ellipsoid. The adsorption process of the Janus ellipsoid can lead to different final orientations determined by the interplay of particle aspect ratio and particle wettablity contrast. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 23 pages, 7 figures

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)