BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jens Harting

Jens Harting appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dynphysics.comp-phBiological Physicsphysics.chem-phcond-mat.mtrl-scicond-mat.stat-mechphysics.app-phDistributed, Parallel, and Cluster Computingcond-mat.mes-hallphysics.geo-phTissues and Organs

Catalog footprint

What is connected

62works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

62 published item(s)

preprint2022arXiv

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Controlling the dewetting morphologies of thin liquid films by switchable substrates

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Theoretical framework for pairwise microswimmer interactions

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capillary interactions between soft capsules protruding through thin fluid films

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Desorption energy of soft particles from a fluid interface

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Self-Similar Liquid Lens Coalescence

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure and rheology of suspensions of spherical strain-hardening capsules

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Inertial migration of neutrally-buoyant particles in superhydrophobic channels

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

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Predicting optimal hematocrit in silico

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

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Active elastohydrodynamics of vesicles in narrow, blind constrictions

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Controlled Capillary Assembly of Magnetic Janus Particles at Fluid-Fluid Interfaces

Capillary interactions can be used to direct assembly of particles adsorbed at fluid-fluid interfaces. Precisely controlling the magnitude and direction of capillary interactions to assemble particles into favoured structures for materials science purposes is desirable but challenging. In this paper, we investigate capillary interactions between magnetic Janus particles adsorbed at fluid-fluid interfaces. We develop a pair-interaction model that predicts that these particles should arrange into a side--side configuration, and carry out simulations that confirm the predictions of our model. Finally, we investigate the monolayer structures that form when many magnetic Janus particles adsorb at the interface. We find that the particles arrange into long, straight chains exhibiting little curvature, in contrast with capillary interactions between ellipsoidal particles. We further find a regime in which highly ordered, lattice-like monolayer structures form, which can be tuned dynamically using an external magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effect of body deformability on microswimming

In this work we consider the following question: given a mechanical microswimming mechanism, does increased deformability of the swimmer body hinder or promote the motility of the swimmer? To answer this we study a microswimmer model composed of deformable beads connected with springs. We determine the velocity of the swimmer analytically, starting from the forces driving the motion and assuming that the oscillations in the effective radii of the beads are known and are much smaller than the radii themselves. We find that to the lowest order, only the driving frequency mode of the surface oscillations contributes to the swimming velocity, and that this velocity may both rise and fall with the deformability of the beads depending on the spring constant. To test these results, we run immersed boundary lattice Boltzmann simulations of the swimmer, and show that they reproduce both the velocity-promoting and velocity-hindering effects of bead deformability correctly in the predicted parameter ranges. Our results mean that for a general swimmer, its elasticity determines whether passive deformations of the swimmer body, induced by the fluid flow, aid or oppose the motion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inversion of hematocrit partition at microfluidic bifurcations

Partitioning of red blood cells (RBCs) at the level of bifurcations in the microcirculatory system affects many physiological functions yet it remains poorly understood. We address this problem by using T-shaped microfluidic bifurcations as a model. Our computer simulations and in vitro experiments reveal that the hematocrit ($ϕ_0$) partition depends strongly on RBC deformability, as long as $ϕ_0 <20$% (within the normal range in microcirculation), and can even lead to complete deprivation of RBCs in a child branch. Furthermore, we discover a deviation from the Zweifach-Fung effect which states that the child branch with lower flow rate recruits less RBCs than the higher flow rate child branch. At small enough $ϕ_0$, we get the inverse scenario, and the hematocrit in the lower flow rate child branch is even higher than in the parent vessel. We explain this result by an intricate up-stream RBC organization and we highlight the extreme dependence of RBC transport on geometrical and cell mechanical properties. These parameters can lead to unexpected behaviors with consequences on the microcirculatory function and oxygen delivery in healthy and pathological conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiphase lattice Boltzmann simulations for porous media applications -- a review

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Setting the pace of microswimmers: when increasing viscosity speeds up self-propulsion

It has long been known that some microswimmers seem to swim counter-intuitively faster when the viscosity of the surrounding fluid is increased, whereas others slow down. This conflicting dependence of the swimming velocity on the viscosity is poorly understood theoretically. Here we explain that any mechanical microswimmer with an elastic degree of freedom in a simple Newtonian fluid can exhibit both kinds of response to an increase in the fluid viscosity for different viscosity ranges, if the driving is weak. The velocity response is controlled by a single parameter $\varGamma$, the ratio of the relaxation time of the elastic component of the swimmer in the viscous fluid and the swimming stroke period. This defines two velocity-viscosity regimes, which we characterize using the bead-spring microswimmer model and analyzing the different forces acting on the parts of this swimmer. The analytical calculations are supported by lattice-Boltzmann simulations, which accurately reproduce the two velocity regimes for the predicted values of $\varGamma$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Soft particles at a fluid interface

Particles added to a fluid interface can be used as a surface stabilizer in the food, oil and cosmetic industries. As an alternative to rigid particles, it is promising to consider highly deformable particles that can adapt their conformation at the interface. In this study we compute the shapes of soft elastic particles using molecular dynamics simulations of a cross-linked polymer gel, complemented by continuum calculations based on linear elasticity. It is shown that the particle shape is not only affected by the Young's modulus of the particle, but also strongly depends on whether the gel is partially or completely wetting the fluid interface. We find that the molecular simulations for the partially wetting case are very accurately described by the continuum theory. By contrast, when the gel is completely wetting the fluid interface the linear theory breaks down and we reveal that molecular details have a strong influence on the equilibrium shape.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tunable Dipolar Capillary Deformations for Magnetic Janus Particles at Fluid-Fluid Interfaces

Janus particles have attracted significant interest as building blocks for complex materials in recent years. Furthermore, capillary interactions have been identified as a promising tool for directed self-assembly of particles at fluid-fluid interfaces. In this paper, we develop theoretical models describing the behaviour of magnetic Janus particles adsorbed at fluid-fluid interfaces interacting with an external magnetic field. Using numerical simulations, we test the models predictions and show that the magnetic Janus particles deform the interface in a dipolar manner. We suggest how to utilise the resulting dipolar capillary interactions to assemble particles at a fluid-fluid interface, and further demonstrate that the strength of these interactions can be tuned by altering the external field strength, opening up the possibility to create novel, reconfigurable materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Assembling Ellipsoidal Particles at Fluid Interfaces using Switchable Dipolar Capillary Interactions

The fabrication of novel soft materials is an important scientific and technological challenge. We investigate the response of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces to external magnetic fields. By exploiting previously discovered first-order orientation phase transitions we show how to switch on and off dipolar capillary interactions between particles, leading to the formation of distinctive self-assembled structures and allowing dynamic control of the bottom-up fabrication of reconfigurable novel-structured materials

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detachment Energies of Spheroidal Particles from Fluid-Fluid Interfaces

The energy required to detach a single particle from a fluid-fluid interface is an important parameter for designing certain soft materials, for example, emulsions stabilised by colloidal particles, colloidosomes designed for targeted drug delivery, and bio-sensors composed of magnetic particles adsorbed at interfaces. For a fixed particle volume, prolate and oblate spheroids attach more strongly to interfaces because they have larger particle-interface areas. Calculating the detachment energy of spheroids necessitates the difficult measurement of particle-liquid surface tensions, in contrast with spheres, where the contact angle suffices. We develop a simplified detachment energy model for spheroids which depends only on the particle aspect ratio and the height of the particle centre of mass above the fluid-fluid interface. We use lattice Boltzmann simulations to validate the model and provide quantitative evidence that the approach can be applied to simulate particle-stabilized emulsions, and highlight the experimental implications of this validation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Domain and droplet sizes in emulsions stabilized by colloidal particles

Particle-stabilized emulsions are commonly used in various industrial applications. These emulsions can present in different forms, such as Pickering emulsions or bijels, which can be distinguished by their different topologies and rheology. We numerically investigate the effect of the volume fraction and the uniform wettability of the stabilizing spherical particles in mixtures of two fluids. For this, we use the well-established three-dimensional lattice Boltzmann method, extended to allow for the added colloidal particles with non-neutral wetting properties. We obtain data on the domain sizes in the emulsions by using both structure functions and the Hoshen-Kopelman (HK) algorithm, and demonstrate that both methods have their own (dis-)advantages. We confirm an inverse dependence between the concentration of particles and the average radius of the stabilized droplets. Furthermore, we demonstrate the effect of particles detaching from interfaces on the emulsion properties and domain size measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Forced transport of deformable containers through narrow constrictions

We study, numerically and analytically, the forced transport of deformable containers through a narrow constriction. Our central aim is to quantify the competition between the constriction geometry and the active forcing, regulating whether and at which speed a container may pass through the constriction and under what conditions it gets stuck. We focus, in particular, on the interrelation between the force that propels the container and the radius of the channel, as these are the external variables that may be directly controlled in both artificial and physiological settings. We present Lattice-Boltzmann simulations that elucidate in detail the various phases of translocation, and present simplified analytical models that treat two limiting types of these membrane containers: deformational energy dominated by the bending or stretching contribution. In either case we find excellent agreement with the full simulations, and our results reveal that not only the radius but also the length of the constriction determines whether or not the container will pass.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interface deformations affect the orientation transition of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces

Manufacturing new soft materials with specific optical, mechanical and magnetic properties is a significant challenge. Assembling and manipulating colloidal particles at fluid interfaces is a promising way to make such materials. We use lattice-Boltzmann simulations to investigate the response of magnetic ellipsoidal particles adsorbed at liquid-liquid interfaces to external magnetic fields. We provide further evidence for the first-order orientation phase transition predicted by Bresme and Faraudo [Journal of Physics: Condensed Matter 19 (2007), 375110]. We show that capillary interface deformations around the ellipsoidal particle significantly affect the tilt-angle of the particle for a given dipole-field strength, altering the properties of the orientation transition. We propose scaling laws governing this transition, and suggest how to use these deformations to facilitate particle assembly at fluid-fluid interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interplay between microdynamics and macrorheology in vesicle suspensions

The microscopic dynamics of objects suspended in a fluid determines the macroscopic rheology of a suspension. For example, as shown by Danker and Misbah [Phys. Rev. Lett. {\bf 98}, 088104 (2007)], the viscosity of a dilute suspension of fluid-filled vesicles is a non-monotonic function of the viscosity contrast (the ratio between the viscosities of the internal encapsulated and the external suspending fluids) and exhibits a minimum at the critical point of the tank-treading-to-tumbling transition. By performing numerical simulations, we recover this effect and demonstrate that it persists for a wide range of vesicle parameters such as the concentration, membrane deformability, or swelling degree. We also explain why other numerical and experimental studies lead to contradicting results. Furthermore, our simulations show that this effect even persists in non-dilute and confined suspensions, but that it becomes less pronounced at higher concentrations and for more swollen vesicles. For dense suspensions and for spherical (circular in 2D) vesicles, the intrinsic viscosity tends to depend weakly on the viscosity contrast.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interplay of inertia and deformability on rheological properties of a suspension of capsules

The interplay of inertia and deformability has a substantial impact on the transport of soft particles suspended in a fluid. However, to date a thorough understanding of these systems is still missing and only a limited number of experimental and theoretical studies is available. We combine the finite-element, immersed-boundary and lattice-Boltzmann methods to simulate three-dimensional suspensions of soft particles subjected to planar Poiseuille flow at finite Reynolds numbers. Our findings confirm that the particle deformation and inclination increase when inertia is present. We observe that the Segré-Silberberg effect is suppressed with respect to the particle deformability. Depending on the deformability and strength of inertial effects, inward or outward lateral migration of the particles takes place. In particular, for increasing Reynolds numbers and strongly deformable particles, a hitherto unreported distinct flow focusing effect emerges which is accompanied by a non-monotonic behaviour of the apparent suspension viscosity and thickness of the particle-free layer close to the channel walls. This effect can be explained by the behaviour of a single particle and the change of the particle collision mechanism when both deformability and inertia effects are relevant.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice-Boltzmann simulations of the drag force on a sphere approaching a superhydrophobic striped plane

By means of lattice-Boltzmann simulations the drag force on a sphere of radius R approaching a superhydrophobic striped wall has been investigated as a function of arbitrary separation h. Superhydrophobic (perfect-slip vs. no-slip) stripes are characterized by a texture period L and a fraction of the gas area $ϕ$. For very large values of h/R we recover the macroscopic formulae for a sphere moving towards a hydrophilic no-slip plane. For h/R=O(1) and smaller the drag force is smaller than predicted by classical theories for hydrophilic no-slip surfaces, but larger than expected for a sphere interacting with a uniform perfectly slipping wall. At a thinner gap, $h\ll R$ the force reduction compared to a classical result becomes more pronounced, and is maximized by increasing $ϕ$. In the limit of very small separations our simulation data are in quantitative agreement with an asymptotic equation, which relates a correction to a force for superhydrophobic slip to texture parameters. In addition, we examine the flow and pressure field and observe their oscillatory character in the transverse direction in the vicinity of the wall, which reflects the influence of the heterogeneity and anisotropy of the striped texture. Finally, we investigate the lateral force on the sphere, which is detectable in case of very small separations and is maximized by stripes with $ϕ=0.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parallelised Hoshen-Kopelman algorithm for lattice-Boltzmann simulations

We discuss two topics that we have encountered in our lattice-Boltzmann simulations of complex fluids: the sizes of droplets in particle-stabilised emulsions and deformable particles in fluid flow. The common factor in these seemingly disparate subjects is that both represent an opportunity for a novel application of the Hoshen-Kopelman algorithm. This algorithm is based on detecting connected clusters on a lattice and labelling the involved lattice sites such that all sites that are connected share the same label. The assumption of the presence of a lattice makes it a convenient algorithm to use in combination with lattice-Boltzmann simulations. In order to apply this algorithm on the fly during massively parallel simulations, it needs to be parallelised as well. We present our parallel implementation, which is tailored to a common parallelisation scheme for the lattice-Boltzmann method, and compare it to previous work. We then briefly discuss some examples of results obtained using this procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recent advances in the simulation of particle-laden flows

A substantial number of algorithms exists for the simulation of moving particles suspended in fluids. However, finding the best method to address a particular physical problem is often highly non-trivial and depends on the properties of the particles and the involved fluid(s) together. In this report we provide a short overview on a number of existing simulation methods and provide two state of the art examples in more detail. In both cases, the particles are described using a Discrete Element Method (DEM). The DEM solver is usually coupled to a fluid-solver, which can be classified as grid-based or mesh-free (one example for each is given). Fluid solvers feature different resolutions relative to the particle size and separation. First, a multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine resolution) is presented to study the behavior of particle stabilized fluid interfaces and second, a Smoothed Particle Hydrodynamics implementation (mesh-free, meso-scale resolution, similar to the particle size) is introduced to highlight a new player in the field, which is expected to be particularly suited for flows including free surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Timescales of emulsion formation caused by anisotropic particles

Particle stabilized emulsions have received an enormous interest in the recent past, but our understanding of the dynamics of emulsion formation is still limited. For simple spherical particles, the time dependent growth of fluid domains is dominated by the formation of droplets, particle adsorption and coalescence of droplets (Ostwald ripening), which eventually can be almost fully blocked due to the presence of the particles. Ellipsoidal particles are known to be more efficient stabilizers of fluid interfaces than spherical particles and their anisotropic shape and the related additional rotational degrees of freedom have an impact on the dynamics of emulsion formation. In this paper, we investigate this point by means of simple model systems consisting of a single ellipsoidal particle or a particle ensemble at a flat interface as well as a particle ensemble at a spherical interface. By applying combined multicomponent lattice Boltzmann and molecular dynamics simulations we demonstrate that the anisotropic shape of ellipsoidal particles causes two additional timescales to be of relevance in the dynamics of emulsion formation: a relatively short timescale can be attributed to the adsorption of single particles and the involved rotation of particles towards the interface. As soon as the interface is jammed, however, capillary interactions between the particles cause a local reordering on very long timescales leading to a continuous change in the interface configuration and increase of interfacial area. This effect can be utilized to counteract the thermodynamic instability of particle stabilized emulsions and thus offers the possibility to produce emulsions with exceptional stability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accurate lubrication corrections for spherical and non-spherical particles in discretized fluid simulations

Discretized fluid solvers coupled to a Newtonian dynamics method are a popular tool to study suspension flow. As any simulation technique with finite resolution, the lattice Boltzmann method, when coupled to discrete particles using the momentum exchange method, resolves the diverging lubrication interactions between surfaces near contact only insufficiently. For spheres, it is common practice to account for surface-normal lubrication forces by means of an explicit correction term. A method that additionally covers all further singular interactions for spheres is present in the literature as well as a link-based approach that allows for more general shapes but does not capture non-normal interactions correctly. In this paper, lattice-independent lubrication corrections for aspherical particles are outlined, taking into account all leading divergent interaction terms. An efficient implementation for arbitrary spheroids is presented and compared to purely normal and link-based models. Good consistency with Stokesian dynamics simulations of spheres is found. The non-normal interactions affect the viscosity of suspensions of spheres at volume fractions Φ>= 0.3 but already at Φ>= 0.2 for spheroids. Regarding shear-induced diffusion of spheres, a distinct effect is found at 0.1 <= Φ<= 0.5 and even increasing the resolution of the radius to 8 lattice units is no substitute for an accurate modeling of non-normal interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Complex Dynamics of a Bilamellar Vesicle as a Simple Model for Leukocytes

The influence of the internal structure of a biological cell (e.g., a leukocyte) on its dynamics and rheology is not yet fully understood. By using 2D numerical simulations of a bilamellar vesicle (BLV) consisting of two vesicles as a cell model, we find that increasing the size of the inner vesicle (mimicking the nucleus) triggers a tank-treading-to-tumbling transition. A new dynamical state is observed, the undulating motion: the BLV inclination with respect to the imposed flow oscillates while the outer vesicle develops rotating lobes. The BLV exhibits a non-Newtonian behavior with a time-dependent apparent viscosity during its unsteady motion. Depending on its inclination and on its inner vesicle dynamical state, the BLV behaves like a solid or a liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From creeping to inertial flow in porous media: a lattice Boltzmann - Finite Element study

The lattice Boltzmann method has been successfully applied for the simulation of flow through porous media in the creeping regime. Its technical properties, namely discretization, straightforward implementation and parallelization, are responsible for its popularity. However, flow through porous media is not restricted to near zero Reynolds numbers since inertial effects play a role in numerous natural and industrial processes. In this paper we investigate the capability of the lattice Boltzmann method to correctly describe flow in porous media at moderate Reynolds numbers. The selection of the lattice resolution, the collision kernel and the boundary conditions becomes increasingly important and the challenge is to keep artifacts due to compressibility effects at a minimum. The lattice Boltzmann results show an accurate quantitative agreement with Finite Element Method results and evidence the capability of the method to reproduce Darcy's law at low Reynolds numbers and Forchheimer's law at high Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical simulations of complex fluid-fluid interface dynamics

Interfaces between two fluids are ubiquitous and of special importance for industrial applications, e.g., stabilisation of emulsions. The dynamics of fluid-fluid interfaces is difficult to study because these interfaces are usually deformable and their shapes are not known a priori. Since experiments do not provide access to all observables of interest, computer simulations pose attractive alternatives to gain insight into the physics of interfaces. In the present article, we restrict ourselves to systems with dimensions comparable to the lateral interface extensions. We provide a critical discussion of three numerical schemes coupled to the lattice Boltzmann method as a solver for the hydrodynamics of the problem: (a) the immersed boundary method for the simulation of vesicles and capsules, the Shan-Chen pseudopotential approach for multi-component fluids in combination with (b) an additional advection-diffusion component for surfactant modelling and (c) a molecular dynamics algorithm for the simulation of nanoparticles acting as emulsifiers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Towards a continuum model for particle-induced velocity fluctuations in suspension flow through a stenosed geometry

Non-particulate continuum descriptions allow for computationally efficient modeling of suspension flows at scales that are inaccessible to more detailed particulate approaches. It is well known that the presence of particles influences the effective viscosity of a suspension and that this effect has thus to be accounted for in macroscopic continuum models. The present paper aims at developing a non-particulate model that reproduces not only the rheology but also the cell-induced velocity fluctuations, responsible for enhanced diffusivity. The results are obtained from a coarse-grained blood model based on the lattice Boltzmann method. The benchmark system comprises a flow between two parallel plates with one of them featuring a smooth obstacle imitating a stenosis. Appropriate boundary conditions are developed for the particulate model to generate equilibrated cell configurations mimicking an infinite channel in front of the stenosis. The averaged flow field in the bulk of the channel can be described well by a non-particulate simulation with a matched viscosity. We show that our proposed phenomenological model is capable to reproduce many features of the velocity fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of nanoparticles and surfactant on droplets in shear flow

We present three-dimensional numerical simulations, employing the well-established lattice Boltzmann method, and investigate similarities and differences between surfactants and nanoparticles as additives at a fluid-fluid interface. We report on their respective effects on the surface tension of such an interface. Next, we subject a fluid droplet to shear and explore the deformation properties of the droplet, its inclination angle relative to the shear flow, the dynamics of the particles at the interface, and the possibility of breakup. Particles are seen not to affect the surface tension of the interface, although they do change the overall interfacial free energy. The particles do not remain homogeneously distributed over the interface, but form clusters in preferred regions that are stable for as long as the shear is applied. However, although the overall structure remains stable, individual nanoparticles roam the droplet interface, with a frequency of revolution that is highest in the middle of the droplet interface, normal to the shear flow, and increases with capillary number. We recover Taylor's law for small deformation of droplets when surfactant or particles are added to the droplet interface. The effect of surfactant is captured in the capillary number, but the inertia of adsorbed massive particles increases deformation at higher capillary number and eventually leads to easier breakup of the droplet.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flow past superhydrophobic surfaces with cosine variation in local slip length

Anisotropic super-hydrophobic surfaces have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of super-hydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local slip length. We derive approximate formulae for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that are in good agreement with the exact numerical solution and lattice-Bolzmann simulations for any surface slip fraction. The cosine texture can provide a very large effective (forward) slip, but it was found to be less efficient in generating a transverse flow as compared to super-hydrophobic stripes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How does confinement affect the dynamics of viscous vesicles and red blood cells?

Despite its significance in microfluidics, the effect of confinement on the transition from the tank-treading (steady motion) to the tumbling (unsteady motion) dynamical state of deformable micro-particles has not been studied in detail. In this paper, we investigate the dynamics of a single viscous vesicle under confining shear as a general model system for red blood cells, capsules, or viscous droplets. The transition from tank-treading to tumbling motion can be triggered by the ratio between internal and external fluid viscosities. Here, we show that the transition can be induced solely by reducing the confinement, keeping the viscosity contrast constant. The observed dynamics results from the variation of the relative importance of viscous-, pressure-, and lubrication-induced torques exerted upon the vesicle. Our findings are of interest for designing future experiments or microfluidic devices: the possibility to trigger the tumbling-to-tank-treading transition either by geometry or viscosity contrast alone opens attractive possibilities for microrheological measurements as well as the detection and diagnosis of diseased red blood cells in confined flow.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamic interactions in active colloidal crystal microrheology

In dense colloids it is commonly assumed that hydrodynamic interactions do not play a role. However, a found theoretical quantification is often missing. We present computer simulations that are motivated by experiments where a large colloidal particle is dragged through a colloidal crystal. To qualify the influence of long-ranged hydrodynamics, we model the setup by conventional Langevin dynamics simulations and by an improved scheme with limited hydrodynamic interactions. This scheme significantly improves our results and allows to show that hydrodynamics strongly impacts on the development of defects, the crystal regeneration as well as on the jamming behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice Boltzmann simulations of anisotropic particles at liquid interfaces

Complex colloidal fluids, such as emulsions stabilized by complex shaped particles, play an important role in many industrial applications. However, understanding their physics requires a study at sufficiently large length scales while still resolving the microscopic structure of a large number of particles and of the local hydrodynamics. Due to its high degree of locality, the lattice Boltzmann method, when combined with a molecular dynamics solver and parallelized on modern supercomputers, provides a tool that allows such studies. Still, running simulations on hundreds of thousands of cores is not trivial. We report on our practical experiences when employing large fractions of an IBM Blue Gene/P system for our simulations. Then, we extend our model for spherical particles in multicomponent flows to anisotropic ellipsoidal objects rendering the shape of e.g. clay particles. The model is applied to a number of test cases including the adsorption of single particles at fluid interfaces and the formation and stabilization of Pickering emulsions or bijels.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Micro- and nanoscale fluid flow on chemical channels

We study the time evolution and driven motion of thin liquid films lying on top of chemical patterns on a substrate. Lattice-Boltzmann and molecular dynamics methods are used for simulations of the flow of microscopic and nanoscopic films, respectively. Minimization of fluid surface area is used to examine the corresponding equilibrium free energy landscapes. The focus is on motion across patterns containing diverging and converging flow junctions, with an eye towards applications to lab-on-a-chip devices. Both open liquid-vapor systems driven by body forces and confined liquid-liquid systems driven by boundary motion are considered. As in earlier studies of flow on a linear chemical channel, we observe continuous motion of a connected liquid film across repeated copies of the pattern, despite the appearance of pearling instabilities of the interface. Provided that the strength of the driving force and the volume of liquid are not too large, the liquid is confined to the chemical channels and its motion can be directed by small variations in the geometry of the pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tensorial slip of super-hydrophobic channels

We describe a generalization of the tensorial slip boundary condition, originally justified for a thick (compared to texture period) channel, to any channel thickness. The eigenvalues of the effective slip length tensor, however, in general case become dependent on the gap and cannot be viewed as a local property of the surface, being a global characteristic of the channel. To illustrate the use of the tensor formalism we develop a semi-analytical theory of an effective slip in a parallel-plate channel with one super-hydrophobic striped and one hydrophilic surface. Our approach is valid for any local slip at the gas sectors and an arbitrary distance between the plates, ranging from a thick to a thin channel. We then present results of lattice Boltzmann simulations to validate the analysis. Our results may be useful for extracting effective slip tensors from global measurements, such as the permeability of a channel, in experiments or simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

From Bijels to Pickering emulsions: a lattice Boltzmann study

Particle stabilized emulsions are ubiquitous in the food and cosmetics industry, but our understanding of the influence of microscopic fluid-particle and particle-particle interactions on the macroscopic rheology is still limited. In this paper we present a simulation algorithm based on a multicomponent lattice Boltzmann model to describe the solvents combined with a molecular dynamics solver for the description of the solved particles. It is shown that the model allows a wide variation of fluid properties and arbitrary contact angles on the particle surfaces. We demonstrate its applicability by studying the transition from a "bicontinuous interfacially jammed emulsion gel" (bijel) to a "Pickering emulsion" in dependence on the contact angle, the particle concentration, and the ratio of the solvents.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantification of the performance of chaotic micromixers on the basis of finite time Lyapunov exponents

Chaotic micromixers such as the staggered herringbone mixer developed by Stroock et al. allow efficient mixing of fluids even at low Reynolds number by repeated stretching and folding of the fluid interfaces. The ability of the fluid to mix well depends on the rate at which "chaotic advection" occurs in the mixer. An optimization of mixer geometries is a non trivial task which is often performed by time consuming and expensive trial and error experiments. In this paper an algorithm is presented that applies the concept of finite-time Lyapunov exponents to obtain a quantitative measure of the chaotic advection of the flow and hence the performance of micromixers. By performing lattice Boltzmann simulations of the flow inside a mixer geometry, introducing massless and non-interacting tracer particles and following their trajectories the finite time Lyapunov exponents can be calculated. The applicability of the method is demonstrated by a comparison of the improved geometrical structure of the staggered herringbone mixer with available literature data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulations of Blood Flow in Plain Cylindrical and Constricted Vessels with Single Cell Resolution

Understanding the physics of blood is challenging due to its nature as a suspension of soft particles and the fact that typical problems involve different scales. This is valid also for numerical investigations. In fact, many computational studies either neglect the existence of discrete cells or resolve relatively few cells very accurately. The authors recently developed a simple and highly efficient yet still particulate model with the aim to bridge the gap between currently applied methods. The present work focuses on its applicability to confined flows in vessels of diameters up to 100 micrometres. For hematocrit values below 30 percent, a dependence of the apparent viscosity on the vessel diameter in agreement with experimental literature data is found.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two-dimensional Vesicle dynamics under shear flow: effect of confinement

Dynamics of a single vesicle under shear flow between two parallel plates is studied using two-dimensional lattice-Boltzmann simulations. We first present how we adapted the lattice-Boltzmann method to simulate vesicle dynamics, using an approach known from the immersed boundary method. The fluid flow is computed on an Eulerian regular fixed mesh while the location of the vesicle membrane is tracked by a Lagrangian moving mesh. As benchmarking tests, the known vesicle equilibrium shapes in a fluid at rest are found and the dynamical behavior of a vesicle under simple shear flow is being reproduced. Further, we focus on investigating the effect of the confinement on the dynamics, a question that has received little attention so far. In particular, we study how the vesicle steady inclination angle in the tank-treading regime depends on the degree of confinement. The influence of the confinement on the effective viscosity of the composite fluid is also analyzed. At a given reduced volume (the swelling degree) of a vesicle we find that both the inclination angle, and the membrane tank-treading velocity decrease with increasing confinement. At sufficiently large degree of confinement the tank-treading velocity exhibits a non-monotonous dependence on the reduced volume and the effective viscosity shows a nonlinear behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A simplified particulate model for coarse-grained hemodynamics simulations

Human blood flow is a multi-scale problem: in first approximation, blood is a dense suspension of plasma and deformable red cells. Physiological vessel diameters range from about one to thousands of cell radii. Current computational models either involve a homogeneous fluid and cannot track particulate effects or describe a relatively small number of cells with high resolution, but are incapable to reach relevant time and length scales. Our approach is to simplify much further than existing particulate models. We combine well established methods from other areas of physics in order to find the essential ingredients for a minimalist description that still recovers hemorheology. These ingredients are a lattice Boltzmann method describing rigid particle suspensions to account for hydrodynamic long range interactions and---in order to describe the more complex short-range behavior of cells---anisotropic model potentials known from molecular dynamics simulations. Paying detailedness, we achieve an efficient and scalable implementation which is crucial for our ultimate goal: establishing a link between the collective behavior of millions of cells and the macroscopic properties of blood in realistic flow situations. In this paper we present our model and demonstrate its applicability to conditions typical for the microvasculature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Contact angle dependence on the fluid-wall dispersive energy

Vapor-liquid menisci of the truncated and shifted Lennard-Jones fluid between parallel planar walls are investigated by molecular dynamics simulation. Thereby, the characteristic energy of the unlike dispersive interaction between fluid molecules and wall atoms is systematically varied to determine its influence on the contact angle. The temperature is varied as well, covering most of the range between the triple point temperature and the critical temperature of the bulk fluid. The transition between obtuse and acute angles is found to occur at a temperature-independent magnitude of the fluid-wall dispersive interaction energy. On the basis of the present simulation results, fluid-wall interaction potentials can be adjusted to contact angle measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Contact angle determination in multicomponent lattice Boltzmann simulations

Droplets on hydrophobic surfaces are ubiquitous in microfluidic applications and there exists a number of commonly used multicomponent and multiphase lattice Boltzmann schemes to study such systems. In this paper we focus on a popular implementation of a multicomponent model as introduced by Shan and Chen. Here, interactions between different components are implemented as repulsive forces whose strength is determined by model parameters. In this paper we present simulations of a droplet on a hydrophobic surface. We investigate the dependence of the contact angle on the simulation parameters and quantitatively compare different approaches to determine it. Results show that the method is capable of modelling the whole range of contact angles. We find that the a priori determination of the contact angle is depending on the simulation parameters with an uncertainty of 10 to 20%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evaluation of pressure boundary conditions for permeability calculations using the lattice-Boltzmann method

Lattice-Boltzmann (LB) simulations are a common tool to numerically estimate the permeability of porous media. For valuable results, the porous structure has to be well resolved resulting in a large computational effort as well as high memory demands. In order to estimate the permeability of realistic samples, it is of importance to not only implement very efficient codes, but also to choose the most appropriate simulation setup to achieve accurate results. With the focus on accuracy and computational effort, we present a comparison between different methods to apply an effective pressure gradient, efficient boundary conditions, as well as two LB implementations based on pore-matrix and pore-list data structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lattice Boltzmann simulations of liquid flilm drainage between smooth surfaces

Exploring the hydrodynamic boundary of a surface by approaching a colloidal sphere and measuring the occurring drag force is a common experimental technique. However, numerous parameters like the wettability and surface roughness influence the result. In experiments these cannot be separated easily. For a deeper understanding of such surface effects a tool is required that predicts the influence of different surface properties. In this paper we present computer simulations based on the lattice Boltzmann method of a sphere submerged in a Newtonian liquid. We show that our method is able to reproduce the theoretical predictions for flat and noninteracting surfaces. In order to provide high precision simulation results the influence of finite size effects has to be well controlled. Therefore we study the influence of the required system size and resolution of the sphere and demonstrate that already moderate computing resources allow the error to be kept below 1%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantitative analysis of numerical estimates for the permeability of porous media from lattice-Boltzmann simulations

During the last decade, lattice-Boltzmann (LB) simulations have been improved to become an efficient tool for determining the permeability of porous media samples. However, well known improvements of the original algorithm are often not implemented. These include for example multirelaxation time schemes or improved boundary conditions, as well as different possibilities to impose a pressure gradient. This paper shows that a significant difference of the calculated permeabilities can be found unless one uses a carefully selected setup. We present a detailed discussion of possible simulation setups and quantitative studies of the influence of simulation parameters. We illustrate our results by applying the algorithm to a Fontainebleau sandstone and by comparing our benchmark studies to other numerical permeability measurements in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random-roughness hydrodynamic boundary conditions

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rotational behavior of red blood cells in suspension---a mesoscale simulation study

The nature of blood as a suspension of red blood cells makes computational hemodynamics a demanding task. Our coarse-grained blood model, which builds on a lattice Boltzmann method for soft particle suspensions, enables the study of the collective behavior of the order of 10^6 cells in suspension. After demonstrating the viscosity measurement in Kolmogorov flow, we focus on the statistical analysis of the cell orientation and rotation in Couette flow. We quantify the average inclination with respect to the flow and the nematic order as a function of shear rate and hematocrit. We further record the distribution of rotation periods around the vorticity direction and find a pronounced peak in the vicinity of the theoretical value for free model cells even though cell-cell interactions manifest themselves in a substantial width of the distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results allow to assist in the optimization of microchannel flows for large throughput.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Simulation of fluid flow in hydrophobic rough microchannels

Surface effects become important in microfluidic setups because the surface to volume ratio becomes large. In such setups the surface roughness is not any longer small compared to the length scale of the system and the wetting properties of the wall have an important influence on the flow. However, the knowledge about the interplay of surface roughness and hydrophobic fluid-surface interaction is still very limited because these properties cannot be decoupled easily in experiments. We investigate the problem by means of lattice Boltzmann (LB) simulations of rough microchannels with a tunable fluid-wall interaction. We introduce an ``effective no-slip plane'' at an intermediate position between peaks and valleys of the surface and observe how the position of the wall may change due to surface roughness and hydrophobic interactions. We find that the position of the effective wall, in the case of a Gaussian distributed roughness depends linearly on the width of the distribution. Further we are able to show that roughness creates a non-linear effect on the slip length for hydrophobic boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Computer Simulation of Particle Suspensions

Particle suspensions are ubiquitous in our daily life, but are not well understood due to their complexity. During the last twenty years, various simulation methods have been developed in order to model these systems. Due to varying properties of the solved particles and the solvents, one has to choose the simulation method properly in order to use the available compute resources most effectively with resolving the system as well as needed. Various techniques for the simulation of particle suspensions have been implemented at the Institute for Computational Physics allowing us to study the properties of clay-like systems, where Brownian motion is important, more macroscopic particles like glass spheres or fibers solved in liquids, or even the pneumatic transport of powders in pipes. In this paper we will present the various methods we applied and developed and discuss their individual advantages.

0 citations0 reviewsNo code linkedOpen access