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

Timm Krüger

Timm Krüger 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-dynBiological Physicsphysics.comp-phComputational Engineering, Finance, and Sciencecond-mat.othercond-mat.stat-mechphysics.med-phQuantitative Methods

Catalog footprint

What is connected

21works
9topics
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

21 published item(s)

preprint2020arXiv

Spatio-temporal dynamics of dilute red blood cell suspensions in a microchannel flow at low Reynolds number

Microfluidic technologies are commonly used for the manipulation of red blood cell (RBC) suspensions and analyses of flow-mediated biomechanics. To enhance the performance of microfluidic devices, understanding the dynamics of the suspensions processed within is crucial. We report novel aspects of the spatio-temporal dynamics of RBC suspensions flowing through a typical microchannel at low Reynolds number. Through experiments with dilute RBC suspensions, we find an off-centre two-peak (OCTP) profile of cells contrary to the centralised distribution commonly reported for low-inertia flows. This is reminiscent of the well-known "tubular pinch effect" which arises from inertial effects. However, given the conditions of negligible inertia in our experiments, an alternative explanation is needed for this OCTP profile. Our massively-parallel simulations of RBC flow in real-size microfluidic dimensions using the immersed-boundary-lattice-Boltzmann method (IB-LBM) confirm the experimental findings and elucidate the underlying mechanism for the counterintuitive RBC pattern. By analysing the RBC migration and cell-free layer (CFL) development within a high-aspect-ratio channel, we show that such a distribution is co-determined by the spatial decay of hydrodynamic lift and the global deficiency of cell dispersion in dilute suspensions. We find a CFL development length greater than 46 and 28 hydraulic diameters in the experiment and simulation, respectively, exceeding typical lengths of microfluidic designs. Our work highlights the key role of transient cell distribution in dilute suspensions, which may negatively affect the reliability of experimental results if not taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Modelling ternary fluids in contact with elastic membranes

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modelling approach and taking into account the thermodynamics laws, we derive the equations governing the ternary fluid flow and dynamics of the membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method, employed for resolving the evolution equations of the ternary fluid in an Eulerian description, coupled to the immersed boundary method, allowing for the membrane equations of motion to be solved in a Lagrangian system. The configuration of an elastic capsule placed at a fluid-fluid interface is considered for validation purposes. Systematic simulations are performed for a detailed comparison with reference numerical results obtained by Surface Evolver, and the Galilean invariance of the proposed model is also proven. The proposed approach is versatile, and a wide range of geometries can be simulated. To demonstrate this, the problem of a capillary bridge formed between two deformable capsules is investigated here.

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

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
preprint2016arXiv

Ternary Free Energy Lattice Boltzmann Model with Tunable Surface Tensions and Contact Angles

We present a new ternary free energy lattice Boltzmann model. The distinguishing feature of our model is that we are able to analytically derive and independently vary all fluid-fluid surface tensions and the solid surface contact angles. We carry out a number of benchmark tests: (i) double emulsions and liquid lenses to validate the surface tensions, (ii) ternary fluids in contact with a square well to compare the contact angles against analytical predictions, and (iii) ternary phase separation to verify that the multicomponent fluid dynamics is accurately captured. Additionally we also describe how the model here presented here can be extended to include an arbitrary number of fluid components.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Breakdown of deterministic lateral displacement efficiency for non-dilute suspensions: a numerical study

We investigate the effect of particle volume fraction on the efficiency of deterministic lateral displacement (DLD) devices. DLD is a popular passive sorting technique for microfluidic applications. Yet, it has been designed for treating dilute suspensions, and its efficiency for denser samples is not well known. We perform 3D simulations based on the immersed-boundary, lattice-Boltzmann and finite-element methods to model the flow of red blood cells (RBCs) in different DLD devices. We quantify the DLD efficiency in terms of appropriate "failure" probabilities and RBC counts in designated device outlets. Our main result is that the displacement mode breaks down upon an increase of RBC volume fraction, while the zigzag mode remains relatively robust. This suggests that the separation of larger particles (such as white blood cells) from a dense RBC background is simpler than separating smaller particles (such as platelets) from the same background. The observed breakdown stems from non-deterministic particle collisions interfering with the designed deterministic nature of DLD devices. Therefore, we postulate that dense suspension effects generally hamper efficient particle separation in devices based on deterministic principles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effect of tube diameter and capillary number on platelet margination and near-wall dynamics

The effect of tube diameter $D$ and capillary number $Ca$ on platelet margination in blood flow at $\approx 37\%$ tube haematocrit is investigated. The system is modelled as three-dimensional suspension of deformable red blood cells and nearly rigid platelets using a combination of the lattice-Boltzmann, immersed boundary and finite element methods. Results show that margination is facilitated by a non-diffusive radial platelet transport. This effect is important near the edge of the cell-free layer, but it is only observed for $Ca > 0.2$, when red blood cells are tank-treading rather than tumbling. It is also shown that platelet trapping in the cell-free layer is reversible for $Ca \leq 0.2$. Only for the smallest investigated tube ($D = 10 μ\text{m}$) margination is essentially independent of $Ca$. Once platelets have reached the cell-free layer, they tend to slide rather than tumble. The tumbling rate is essentially independent of $Ca$ but increases with $D$. Tumbling is suppressed by the strong confinement due to the relatively small cell-free layer thickness at $\approx 37\%$ tube haematocrit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuations and diffusion in sheared athermal suspensions of deformable particles

We analyze fluctuations of particle displacements and stresses in a sheared athermal suspension of elastic capsules (red blood cells). Upon variation of the volume fraction from the dilute up to the highly concentrated regime, our numerical simulations reveal different characteristic power-law regimes of the fluctuation variances and relaxation times. In the jammed phase and at high shear rates, anomalous scaling exponents are found that deviate from pure dimensional predictions. The observed behavior is rationalized via kinetic arguments and a dissipation balance model that takes into account the local fluid flows between the particles. Our findings support the view that the rheology of dense suspensions is essentially governed by the non-affine displacements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computer simulations reveal complex distribution of haemodynamic forces in a mouse retina model of angiogenesis

There is currently limited understanding of the role played by haemodynamic forces on the processes governing vascular development. One of many obstacles to be overcome is being able to measure those forces, at the required resolution level, on vessels only a few micrometres thick. In the current paper, we present an in silico method for the computation of the haemodynamic forces experienced by murine retinal vasculature (a widely used vascular development animal model) beyond what is measurable experimentally. Our results show that it is possible to reconstruct high-resolution three-dimensional geometrical models directly from samples of retinal vasculature and that the lattice-Boltzmann algorithm can be used to obtain accurate estimates of the haemodynamics in these domains. We generate flow models from samples obtained at postnatal days (P) 5 and 6. Our simulations show important differences between the flow patterns recovered in both cases, including observations of regression occurring in areas where wall shear stress gradients exist. We propose two possible mechanisms to account for the observed increase in velocity and wall shear stress between P5 and P6: i) the measured reduction in typical vessel diameter between both time points, ii) the reduction in network density triggered by the pruning process. The methodology developed herein is applicable to other biomedical domains where microvasculature can be imaged but experimental flow measurements are unavailable or difficult to obtain.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deformability-based red blood cell separation in deterministic lateral displacement devices - a simulation study

We show, via three-dimensional immersed-boundary-finite-element-lattice-Boltzmann simulations, that deformability-based red blood cell (RBC) separation in deterministic lateral displacement (DLD) devices is possible. This is due to the deformability-dependent lateral extension of RBCs and enables us to predict a priori which RBCs will be displaced in a given DLD geometry. Several diseases affect the deformability of human cells. Malaria-infected RBCs or sickle cells, for example, tend to become stiffer than their healthy counterparts. It is therefore desirable to design microfluidic devices which can detect those diseases based on the cells' deformability fingerprint, rather than preparing samples using expensive and time-consuming biochemical preparation steps. Our findings should be helpful in the development of new methods for sorting cells and particles by deformability.

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

Impact of blood rheology on wall shear stress in a model of the middle cerebral artery

Perturbations to the homeostatic distribution of mechanical forces exerted by blood on the endothelial layer have been correlated with vascular pathologies including intracranial aneurysms and atherosclerosis. Recent computational work suggests that in order to correctly characterise such forces, the shear-thinning properties of blood must be taken into account. To the best of our knowledge, these findings have never been compared against experimentally observed pathological thresholds. In the current work, we apply the three-band diagram (TBD) analysis due to Gizzi et al. to assess the impact of the choice of blood rheology model on a computational model of the right middle cerebral artery. Our results show that, in the model under study, the differences between the wall shear stress predicted by a Newtonian model and the well known Carreau-Yasuda generalized Newtonian model are only significant if the vascular pathology under study is associated with a pathological threshold in the range 0.94 Pa to 1.56 Pa, where the results of the TBD analysis of the rheology models considered differs. Otherwise, we observe no significant differences.

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
preprint2013arXiv

Choice of boundary condition for lattice-Boltzmann simulation of moderate Reynolds number flow in complex domains

Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlet/outlets at arbitrary orientations with respect to the lattice, intermediate Reynolds number, and unsteady flow. Simple bounce-back is one of the most commonly used, simplest, and most computationally efficient boundary conditions, but many others have been proposed. We implement three other methods applicable to complex geometries (Guo, Zheng and Shi, Phys Fluids (2002); Bouzdi, Firdaouss and Lallemand, Phys. Fluids (2001); Junk and Yang Phys. Rev. E (2005)) in our open-source application \HemeLB{}. We use these to simulate Poiseuille and Womersley flows in a cylindrical pipe with an arbitrary orientation at physiologically relevant Reynolds (1--300) and Womersley (4--12) numbers and steady flow in a curved pipe at relevant Dean number (100--200) and compare the accuracy to analytical solutions. We find that both the Bouzidi-Firdaouss-Lallemand and Guo-Zheng-Shi methods give second-order convergence in space while simple bounce-back degrades to first order. The BFL method appears to perform better than GZS in unsteady flows and is significantly less computationally expensive. The Junk-Yang method shows poor stability at larger Reynolds number and so cannot be recommended here. The choice of collision operator (lattice Bhatnagar-Gross-Krook vs.\ multiple relaxation time) and velocity set (D3Q15 vs.\ D3Q19 vs.\ D3Q27) does not significantly affect the accuracy in the problems studied.

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

Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells

Via computer simulations, we provide evidence that the shear rate induced red blood cell tumbling-to-tank-treading transition also occurs at quite high volume fractions, where collective effects are important. The transition takes place as the ratio of suspension stress to the characteristic cell membrane stress exceeds a certain value, independent of volume fraction and cell deformability. This value coincides with that for a transition from an orientationally less ordered to a highly ordered phase. The average cell deformation does not show any signature of the transition, but rather follows a simple scaling law independent of volume fraction.

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
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
preprint2010arXiv

Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method

The deformation of an initially spherical capsule, freely suspended in simple shear flow, can be computed analytically in the limit of small deformations [D. Barthes-Biesel, J. M. Rallison, The Time-Dependent Deformation of a Capsule Freely Suspended in a Linear Shear Flow, J. Fluid Mech. 113 (1981) 251-267]. Those analytic approximations are used to study the influence of the mesh tessellation method, the spatial resolution, and the discrete delta function of the immersed boundary method on the numerical results obtained by a coupled immersed boundary lattice Boltzmann finite element method. For the description of the capsule membrane, a finite element method and the Skalak constitutive model [R. Skalak et al., Strain Energy Function of Red Blood Cell Membranes, Biophys. J. 13 (1973) 245-264] have been employed. Our primary goal is the investigation of the presented model for small resolutions to provide a sound basis for efficient but accurate simulations of multiple deformable particles immersed in a fluid. We come to the conclusion that details of the membrane mesh, as tessellation method and resolution, play only a minor role. The hydrodynamic resolution, i.e., the width of the discrete delta function, can significantly influence the accuracy of the simulations. The discretization of the delta function introduces an artificial length scale, which effectively changes the radius and the deformability of the capsule. We discuss possibilities of reducing the computing time of simulations of deformable objects immersed in a fluid while maintaining high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inertia effects and stress accumulation in a constricted duct: A combined experimental and lattice Boltzmann study

We experimentally and numerically investigate the flow of a Newtonian fluid through a constricted geometry for Reynolds numbers in the range $0.1 - 100$. The major aim is to study non-linear inertia effects at larger Reynolds numbers (>10) on the shear stress evolution in the fluid. This is of particular importance for blood flow as some biophysical processes in blood are sensitive to shear stresses, e.g., the initialization of blood clotting. We employ the lattice Boltzmann method for the simulations. The conclusion of the predictions is that the peak value of shear stress in the constriction grows disproportionally fast with the Reynolds number which leads to a non-linear shear stress accumulation. As a consequence, the combination of constricted blood vessel geometries and large Reynolds numbers may increase the risk of undesired blood clotting.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Shear stress in lattice Boltzmann simulations

A thorough study of shear stress within the lattice Boltzmann method is provided. Via standard multiscale Chapman-Enskog expansion we investigate the dependence of the error in shear stress on grid resolution showing that the shear stress obtained by the lattice Boltzmann method is second order accurate. This convergence, however, is usually spoiled by the boundary conditions. It is also investigated which value of the relaxation parameter minimizes the error. Furthermore, for simulations using velocity boundary conditions, an artificial mass increase is often observed. This is a consequence of the compressibility of the lattice Boltzmann fluid. We investigate this issue and derive an analytic expression for the time-dependence of the fluid density in terms of the Reynolds number, Mach number and a geometric factor for the case of a Poiseuille flow through a rectangular channel in three dimensions. Comparison of the analytic expression with results of lattice Boltzmann simulations shows excellent agreement.

0 citations0 reviewsNo code linkedOpen access