Off-plane motion of an oblate capsule in a simple shear flow

This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.We investigate the mechanical equilibrium state of an oblate capsule when its revolution axis is initially off the shear plane. We consider an oblate capsule with an aspect ratio of 0.5 and a strain-hardening membrane. The three-dimensional fluid-structure interaction problem is solved numerically by coupling a finite element method with a boundary integral method. The capsule converges towards the same mechanical equilibrium state whatever the initial orientation. This equilibrium depends on the capillary number Ca, which compares the viscous to the elastic forces and on the viscosity ratio between the internal and external fluids. For = 1, the tumbling and swinging motions, observed when the revolution axis is initially in the shear plane, are mechanically stable until Ca 1; when Ca is further increased, the capsule assumes the rolling motion that is observed when its revolution axis is initially aligned with the vorticity axis. When is increased, the tumbling-to-swinging transition appears for higher Ca and the swinging-to-rolling transition for lower Ca. For 5, the swinging regime completely disappears: depending on Ca, it is then either the tumbling or the rolling motion that is the mechanical equilibrium state

Dynamics of a spherical capsule in a planar hyperbolic flow: influence of bending resistance

International audienceWe consider an initially spherical capsule freely suspended in a planar hyperbolic flow and study the influence of the wall bending resistance on the capsule dynamics. The capsule wall is assumed to be made of a three-dimensional homogeneous elastic material. The fluid-structure interaction between the capsule and the external flow is modeled numerically by coupling a boundary integral method with a shell finite element method. It is found that, for given three-dimensional wall mechanical properties, the capsule deformability is drastically reduced as the bending resistance is increased. But, if one expresses the same results as a function of the two-dimensional mechanical properties of the mid-surface, which is how the capsule wall is modeled in the thin-shell model, the capsule deformed shape is identical to the one predicted for a capsule devoid of bending resistance. The bending rigidity is found to have a negligible influence on the shape and deformation: the capsule main deformation mode is thus solely a function of the elastic stretching of the mid-surface. The wall bending resistance still plays a role locally in the regions where buckling occurs. Its influence is studied in the low flow strength regime, for which wrinkling of the wall is observed to persist at steady state. We show that the wrinkle wavelength only depends on the bending number, which compares the relative importance of bending and shearing phenomena, and provide the correlation law. This result is interesting as it allows bending resistance to be estimated from experiments on capsules in a planar hyperbolic flow at low flow strength

Motion of a spherical capsule in simple shear flow: influence of the bending resistance

National audienceWe simulate the motion of an initially spherical capsule in a simple shear flow in order to determine the influence of the bending resistance on wrinkle formation on the membrane. We use a numerical method coupling a nonlinear shell finite element method for the capsule wall mechanics with a boundary integral method to solve the Stokes equation. The capsule wall is discretized with MITC linear triangular shell finite elements. We find that, at low flow strength, buckling occurs in the central region of the capsule. The number of wrinkles on the membrane decreases with the bending stiffness and above a critical value, wrinkles no longer form. For thickness to radius ratios below 5%, the bending stiffness does not have any significant effect on the overall capsule motion and deformation. The mean capsule shape is identical whether the wall is modeled as a shell or a two-dimensional membrane, which shows that the dynamics of thin capsules is mainly governed by shear elasticity and membrane effects

Coupling boundary integral and shell finite element methods to study the fluid structure interactions of a microcapsule in a simple shear flow

International audienceWe simulate the motion of an initially spherical capsule in a simple shear flow in order to determine the influence of the bending resistance on the formation of wrinkles on the membrane. The fluid structure interactions are obtained numerically coupling a boundary integral method to solve for the Stokes equation with a nonlinear finite element method for the capsule wall mechanics. The capsule wall is discretized with MITC linear triangular shell finite elements. We find that, at low flow strength, buckling occurs in the central region of the capsule. The number of wrinkles on the membrane decreases with the bending stiffness and, above a critical value, wrinkles no longer form. For thickness to radius ratios below 5%, the bending stiffness does not have any significant effect on the overall capsule motion and deformation. The mean capsule shape is identical whether the wall is modeled as a shell or a two-dimensional membrane, which shows that the dynamics of thin capsules is mainly governed by shear elasticity and membrane effects

Online fabrication and characterization of capsule populations with a flow-focusing microfluidic system

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.We have designed a microfluidic system that combines a double flow-focusing setup for calibrated capsule fabrication with a microchannel for the characterization of their mechanical properties. The double flow-focusing system consists of a first Y junction to create the microdroplets and of a second Y junction to introduce the cross-linking agent allowing the membrane formation. The human serum albumin (HSA) aqueous solution for the dispersed solution, hydrophobic phase for the continuous solution and cross-linking agent solution are introduced by means of syringe pumps. A wavy channel after the second junction allows to control the reticulation time. A cylindrical microchannel then enables to deform and characterize the capsules formed. The mechanical properties of the capsule membrane are obtained by inverse analysis (Chu et al. 2011). The results show that the drop size increases with the flow rate ratio between the central and lateral channels and does not change much regardless of the flow rate of the reticulation phase. The mean shear modulus of the capsules fabricated after 23 s of reticulation is of the order of the surface tension of HSA solution with Dragoxat indicating that the reticulation time is too short to form an elastic membrane around the droplet. When the reticulation time is increased to 60 s, the membrane shear modulus is multiplied by a factor of 3 confirming that a solid membrane has formed around the drop

Motion and deformation of capsules flowing in microfluidic channels

Une capsule est une goutte de liquide enveloppée par une membrane fine et déformable. Les propriétés mécaniques de la membrane sont essentielles pour le mouvement de la capsule. L analyse de l écoulement d une suspension de capsules dans un canal microfluidique au moyen d un modèle mécanique est une technique permettant de déterminer les propriétés élastiques de la membrane. Un modèle numérique tridimensionnel a été développé pour résoudre ce problème d interaction fluide-structure en écoulement confiné. Il couple une méthode des intégrales de frontières pour les écoulements des fluides et une méthode éléments finis pour la déformation de la membrane. Le modèle est utilisé pour étudier l écoulement d une capsule initialement sphérique dans des canaux de différentes sections. Dans un canal cylindrique, on montre que l effet de confinement du canal conduit à la compression de la capsule. Cela engendre la formation de plis sur la membrane autour de l axe de l écoulement, phénomène également observé expérimentalement. Dans un canal de section carrée, les effets de la loi constitutive de la membrane, du rapport de taille et du débit d écoulement sur la déformation de la capsule sont systématiquement étudiés. La comparaison entre les résultats expérimentaux et numériques nous permet de déduire les propriétés mécaniques de la membrane d une population de capsules artificielles. Ce travail démontre la faisabilité de la mesure de propriétés mécaniques d une membrane en utilisant une technique microfluidique en canal carré. Il pourrait être étendu par l étude d écoulements instationnaires dans un canal de section variable ou avec bifurcations.A capsule is a liquid droplet enclosed by a thin and deformable membrane. The membrane mechanical properties are critical for the deformation and motion of capsules. The flow of a capsule suspension through a microfluidic channel with dimensions comparable to those of the suspended particles can be used to infer the membrane elastic properties. However a mechanical model of the process is necessary. We present a three-dimensional numerical model to simulate such fluid-structure interaction problem. We use a novel numerical model that couples a boundary integral method for the internal and external fluid flows and a finite element method for the membrane deformation. The model is applied to study the flow of an initially spherical capsule in channels with different cross-sections. In a cylindrical channel with circular cross-section, we show that the confinement effect leads to the compression of the capsule in the hoop direction. The membrane tends to buckle and to fold as observed experimentally. In a microfluidic channel with a square cross-section, the effects of the membrane constitutive law, size ratio and flow strength on the capsule deformation are systematically studied. The comparison between experimental and numerical results allows us to deduce the membrane mechanical properties of a population of artificial capsules. The present work shows that it is possible to measure the membrane mechanical properties by using a microfluidic channel with a square cross-section. It can be extended to unsteady capsule flows in a channel with variable cross-sections or bifurcations.COMPIEGNE-BU (601592101) / SudocSudocFranceF

Optimal design of hydraulic capsule pipelines transporting spherical capsules

Scarcity of fossil fuels is affecting efficiency of established modes of cargo transport within transportation industry. Extensive research is being carried out on improving efficiency of existing modes of cargo transport, as well as to develop alternative means of transporting goods. One such alternative method can be through the use of energy contained within fluid flowing in pipelines in order to transfer goods from one place to another. The present study focuses on the use of advanced numerical modelling tools to simulate the flow within Hydraulic Capsule Pipelines (HCPs) transporting Spherical Capsules with an aim of developing design equations. Hydraulic Capsule Pipeline is the term which refers to the transport of goods in hollow containers, typically of spherical or cylindrical shapes, termed as capsules, being carried along the pipeline by water. HCPs are being used in mineral industries and have potential for use in Oil & Gas sector. A novel modelling technique has been employed to carry out the investigations under various geometric and flow conditions within HCPs. Both qualitative and quantitative flow analysis has been carried out on the flow of spherical shaped capsules in an HCP for both on-shore and off-shore applications. Furthermore, based on Least-Cost Principle, an optimisation methodology has been developed for the design of single stage HCPs. The input to the optimisation model is the solid throughput required from the system, and the outputs are the optimal diameter of the HCPs and the pumping requirements for the capsule transporting system

Accelerated boundary integral method for multiphase flow in non-periodic geometries

An accelerated boundary integral method for Stokes flow of a suspension of deformable particles is presented for an arbitrary domain and implemented for the important case of a planar slit geometry. The computational complexity of the algorithm scales as O(N) or $O(N\log N$), where $N$ is proportional to the product of number of particles and the number of elements employed to discretize the particle. This technique is enabled by the use of an alternative boundary integral formulation in which the velocity field is expressed in terms of a single layer integral alone, even in problems with non-matched viscosities. The density of the single layer integral is obtained from a Fredholm integral equation of the second kind involving the double layer integral. Acceleration in this implementation is provided by the use of General Geometry Ewald-like method (GGEM) for computing the velocity and stress fields driven by a set of point forces in the geometry of interest. For the particular case of the slit geometry, a Fourier-Chebyshev spectral discretization of GGEM is developed. Efficient implementations employing the GGEM methodology are presented for the resulting single and the double layer integrals. The implementation is validated with test problems on the velocity of rigid particles and drops between parallel walls in pressure driven flow, the Taylor deformation parameter of capsules in simple shear flow and the particle trajectory in pair collisions of capsules in shear flow. The computational complexity of the algorithm is verified with results from several large scale multiparticle simulations.Comment: Journal of Computational Physics, to appea

Effective swimming strategies in low Reynolds number flows

The optimal strategy for a microscopic swimmer to migrate across a linear shear flow is discussed. The two cases, in which the swimmer is located at large distance, and in the proximity of a solid wall, are taken into account. It is shown that migration can be achieved by means of a combination of sailing through the flow and swimming, where the swimming strokes are induced by the external flow without need of internal energy sources or external drives. The structural dynamics required for the swimmer to move in the desired direction is discussed and two simple models, based respectively on the presence of an elastic structure, and on an orientation dependent friction, to control the deformations induced by the external flow, are analyzed. In all cases, the deformation sequence is a generalization of the tank-treading motion regimes observed in vesicles in shear flows. Analytic expressions for the migration velocity as a function of the deformation pattern and amplitude are provided. The effects of thermal fluctuations on propulsion have been discussed and the possibility that noise be exploited to overcome the limitations imposed on the microswimmer by the scallop theorem have been discussed.Comment: 14 pages, 5 figure