| 1. |
- Alghalibi, Dhiya, et al.
(författare)
-
Interface-resolved simulations of particle suspensions in Newtonian, shear thinning and shear thickening carrier fluids
- 2018
-
Ingår i: Journal of Fluid Mechanics. - : Cambridge University Press. - 0022-1120 .- 1469-7645. ; 852, s. 329-357
-
Tidskriftsartikel (refereegranskat)abstract
- We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions (0.1 <= Phi <= 0.4) and particle Reynolds numbers (0.1 <= Rep <= 10). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal Phi. This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter's mean value with the values estimated from the homogenisation theory of Chateau et al. (J. Rheol., vol. 52, 2008, pp. 489-506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the Rep. The wide spectrum of the local shear rate and its dependency on Phi and Rep point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, O(Rep) similar to 10.
|
|
| 2. |
- Banaei, Arash Alizad, et al.
(författare)
-
Numerical simulations of elastic capsules with nucleus in shear flow
- 2017
-
Ingår i: EUROPEAN JOURNAL OF COMPUTATIONAL MECHANICS. - : Taylor & Francis. - 1779-7179 .- 1958-5829. ; 26:1-2, s. 131-153
-
Tidskriftsartikel (refereegranskat)abstract
- The shear-induced deformation of a capsule with a stiff nucleus, a model of eukaryotic cells, is studied numerically. The membrane of the cell and of its nucleus are modelled as a thin elastic material obeying a Neo-Hookean constitutive law. The fluid-structure coupling is obtained using an immersed boundary method. The variations induced by the presence of the nucleus on the cell deformation are investigated when varying the viscosity ratio between the inner and outer fluids, the membrane elasticity and its bending stiffness. The deformation of the eukaryotic cell is smaller than that of the prokaryotic one. The reduction in deformation increases for larger values of the capillary number. The eukaryotic cell remains thicker in itsmiddle part compared to the prokaryotic one, thus making it less flexible to pass through narrow capillaries. For a viscosity ratio of 5, the deformation of the cell is smaller than in the case of uniform viscosity. In addition, for non-zero bending stiffness of the membrane, the deformation decreases and the shape is closer to an ellipsoid. Finally, we compare the results obtained modelling the nucleus as an inner stiffer membrane with those obtained using a rigid particle.
|
|
| 3. |
- Banerjee, Indradumna, 1986-, et al.
(författare)
-
Dynamics of Inertial migration of particles in straight channels
- 2017
-
Konferensbidrag (refereegranskat)abstract
- SUMMARYWe study numerically the entire migration dynamics of spherical and oblate particles in straight rectangular and square cross sectional ducts. The reported results can help in design of straight duct channel based microfluidic systems. KEYWORDS: Inertial microfluidics, Lateral migration, Oblate particles, Straight particles. INTRODUCTIONWe simulate spherical and oblate rigid particles in straight ducts of different aspect ratios using an Immersed Boundary Method. To the best of our knowledge, this is the first time not only the equilibrium position of particles is described, but also the entire migration dynamics of the particle from the initial to final position, including particle trajectory, velocity, rotation and orientation, are investigated. EXPERIMENTAL The fluid is considered incompressible and its motion is governed by the Navier Stokes and Continuity equations. The numerical approach employed is an Immersed Boundary Method (IBM) with two sets of grid points: an equispaced Eulerian mesh for the fluid flow, and Lagrangian grid points uniformly distributed on the surface of the particle. The flow is set up in square and rectangular cross section ducts with no slip and no penetration boundary conditions (Fig.1). RESULTS AND DISCUSSIONWe examine the lateral motion of spherical and oblate particles using the IBM method mentioned above. While simulating three different spheres in a square duct of duct width to sphere diameter ratio H/Ds= [3.5, 5, 10], we find that the particles focus at closest face-cantered equilibrium position from their point of introduction(Fig.2a). We also show the downstream length needed for a sphere to focus, focusing length, as a function of the distance from the vertical duct symmetry line and as a function of Reynolds number(Fig.2b and c respectively). Spherical particles in rectangular duct tend to move laterally toward the longer length wall and then slowly moves towards the equilibrium position at the face-centre along the long wall(fig.3a). We also observe that the focusing length is longer for spherical particles in a rectangular duct, about three times longer than that in square duct (fig. 3b). In case of an oblate particle flowing through a square duct, the lateral motion towards the face centred equilibrium position is similar to that of a sphere (fig.4a), however there is significant tumbling motion of the particle as it tries to reach equilibrium(fig.4b).In a rectangular duct of aspect ratio 2, the oblate particle reaches a steady configuration on the duct symmetry line at the center of the different faces (fig.5a). The focusing length surprisingly is shorter in a rectangular duct for an oblate particle in contrast to its focusing length in a square duct. This is attributed to the higher lateral velocity of the oblate in the second stage of the migration, that with negligible tumbling(fig.5b). The behavior of three oblate particles in a square duct of duct width to longer diameter ratio H/Ds= [3.5, 5, 10] is different compared to a sphere as the largest oblate tend to focus at the duct cross section diagonals compared to the other two which are at face centred equilibrium as in case of a sphere(fig.6a). We attribute this to the rotation rate of the larger particle which is initially increasing and then decreasing(fig.6b).When it comes to focusing lengths, the smaller particles need longer times to reach their final equilibrium(fig.6c). Another interesting behavior we see is the effect of Reynolds number, where it can be seen that the oblate particles show a tilt of 21 degrees when focusing at equilibrium at certain high Reynolds number (fig.7). CONCLUSIONThe results presented employ a highly accurate interface-resolved numerical algorithm, based on the Immersed Boundary Method to study the entire inertial migration of an oblate particle in both square and rectangular ducts and compare it with that of a single sphere. Currently, we apply a volume penalization method and polymeric drag component to the code to solve for viscoelastic effects in circular microcapillaries. ACKNOWLEDGEMENTSThis work was supported by the European Research Council Grant no. ERC-2013-CoG-616186, TRITOS and by the Swedish Research Council Grant no. VR 2014-5001, COST Action MP1305: Flowing matter, and computation time from SNIC. REFERENCES : Lashgari, Iman, et al. Journal of Fluid Mechanics 819 (2017): 540-561.
|
|
| 4. |
- Haque, Simon, et al.
(författare)
-
Stability of fluids with shear-dependent viscosity in the lid-driven cavity
- 2012
-
Ingår i: Journal of Non-Newtonian Fluid Mechanics. - : Elsevier BV. - 0377-0257 .- 1873-2631. ; 173-174, s. 49-61
-
Tidskriftsartikel (refereegranskat)abstract
- The classical problem of the lid-driven cavity extended infinitely in the spanwise direction is considered for non-Newtonian shear-thinning and shear-thickening fluids, where the viscosity is modeled by the Carreau model. Linear stability is used to determine the critical Reynolds number at which the two-dimensional base-flow becomes unstable to three-dimensional spanwise-periodic disturbances. We consider a square cavity, characterized by steady unstable modes, and a shallow cavity of aspect ratio 0.25, where oscillating modes are the first to become unstable for Newtonian fluids. In both cases, the critical Reynolds number first decreases with decreasing power-index n (from shear-thickening to shear-thinning fluids) and then increase again for highly pseudoplastic fluids. In the latter case, this is explained by the thinner boundary layers at the cavity walls and less intense vorticity inside the domain. Interestingly, oscillating modes are found at critical conditions for shear-thickening fluids in a square cavity while the shallow cavity supports a new instability of lower frequency for large enough shear-thinning. Analysis of kinetic energy budgets and structural sensitivity are employed to investigate the physical mechanisms behind the instability.
|
|
| 5. |
- Lashgari, Iman, et al.
(författare)
-
Channel flow of rigid sphere suspensions : Particle dynamics in the inertial regime
- 2016
-
Ingår i: International Journal of Multiphase Flow. - : Elsevier. - 0301-9322 .- 1879-3533. ; 78, s. 12-24
-
Tidskriftsartikel (refereegranskat)abstract
- We consider suspensions of neutrally-buoyant finite-size rigid spherical particles in channel flow and investigate the relation between the particle dynamics and the mean bulk behavior of the mixture for Reynolds numbers 500 ≤ Re ≤ 5000 and particle volume fraction 0 ≤ Φ ≤ 0.3, via fully resolved numerical simulations. Analysis of the momentum balance reveals the existence of three different regimes: laminar, turbulent and inertial shear-thickening depending on which of the stress terms, viscous, Reynolds or particle stress, is the major responsible for the momentum transfer across the channel. We show that both Reynolds and particle stress dominated flows fall into the Bagnoldian inertial regime and that the Bagnold number can predict the bulk behavior although this is due to two distinct physical mechanisms. A turbulent flow is characterized by larger particle dispersion and a more uniform particle distribution, whereas the particulate-dominated flows is associated with a significant particle migration towards the channel center where the flow is smooth laminar-like and dispersion low. Interestingly, the collision kernel shows similar values in the different regimes, although the relative particle velocity and clustering clearly vary with inertia and particle concentration.
|
|
| 6. |
- Lashgari, Iman, et al.
(författare)
-
First instability of the flow of shear-thinning and shear-thickening fluids past a circular cylinder
- 2012
-
Ingår i: Journal of Fluid Mechanics. - : Cambridge University Press (CUP). - 0022-1120 .- 1469-7645. ; 701, s. 201-227
-
Tidskriftsartikel (refereegranskat)abstract
- The first bifurcation and the instability mechanisms of shear-thinning and shear-thickening fluids flowing past a circular cylinder are studied using linear theory and numerical simulations. Structural sensitivity analysis based on the idea of a 'wavemaker' is performed to identify the core of the instability. The shear-dependent viscosity is modelled by the Carreau model where the rheological parameters, i.e. the power-index and the material time constant, are chosen in the range 0.4 <= n <= 1.75 and 0.1 <= lambda <= 100. We show how shear-thinning/shear-thickening effects destabilize/stabilize the flow dramatically when scaling the problem with the reference zero-shear-rate viscosity. These variations are explained by modifications of the steady base flow due to the shear-dependent viscosity; the instability mechanisms are only slightly changed. The characteristics of the base flow, drag coefficient and size of recirculation bubble are presented to assess shear-thinning effects. We demonstrate that at critical conditions the local Reynolds number in the core of the instability is around 50 as for Newtonian fluids. The perturbation kinetic energy budget is also considered to examine the physical mechanism of the instability.
|
|
| 7. |
- Lashgari, Iman, et al.
(författare)
-
Flow regimes of inertial suspensions of finite size particles
- 2015
-
Ingår i: Proceedings - 15th European Turbulence Conference, ETC 2015. - : TU Delft.
-
Konferensbidrag (refereegranskat)abstract
- Inertial regimes in a channel flow of suspension of finite-size neutrally buoyant particles are studied for a wide range of Reynolds numbers: 500 ≤ Re ≤ 5000, and particle volume fractions: 0 ≤ Φ ≤ 0.3. The flow is classified in three different regimes according to the phase-averaged stress budget across the channel [2]. The laminar viscous regime at low Re and Φ where the viscous stress is the dominating term in the budget, the turbulent regime at high Re and relatively low Φ where the momentum is mainly transferred by the action of the Reynolds stress and the inertial shear-thickening regime where the particle stress contributes the most to the significant enhancement of the wall shear stress. Particle distribution and dispersion properties provide additional evidence for the existence of the three different regimes.
|
|
| 8. |
- Lashgari, Iman, 1983-
(författare)
-
Global stability analysis of complex fluids
- 2013
-
Licentiatavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The main focus of this work is on the non-Newtonian effects on the inertial instabilities in shear flows. Both inelastic (Carreau) and elastic models (Oldroyd-B and FENE-P) have been employed to examine the main features of the non-Newtonian fluids; shear-thinning, shear-thickening and elasticity. Several classical configurations have been considered; flow past a circular cylinder, in a lid-driven cavity and in a channel. We have used a wide range of tools for linear stability analysis, modal, non-modal, energy and sensitivity analysis, to determine the instability mechanisms of the non-Newtonian flows and compare them with those of the Newtonian flows. Direct numerical simulations have been also used to prove the results obtained by the linear stability analysis.Significant modifications/alterations in the instability of the different flows have been observed under the action of the non-Newtonian effects. In general, shear-thinning/shear-thickening effects destabilize/stabilize the flow around the cylinder and in a lid driven cavity. Viscoelastic effects both stabilize and destabilize the channel flow depending on the ratio between the viscoelastic and flow time scales. The instability mechanism is just slightly modified in the cylinder flow whereas new instability mechanisms arise in the lid-driven cavity flow. We observe that the non-Newtonian effect can alter the inertial flow at both baseflow and perturbation level (e.g. Carreau fluid past a cylinder or in a lid driven cavity) or it may just affect the perturbations (e.g. Oldroyd-B fluid in channel). In all the flow cases studied, the modifications in the instability dynamics are shown to be strongly connected to the contribution of the different terms in the perturbation kinetic energy budget.
|
|
| 9. |
- Lashgari, Iman, et al.
(författare)
-
Inertial migration of spherical and oblate particles in straight ducts
- 2017
-
Ingår i: Journal of Fluid Mechanics. - : Cambridge University Press (CUP). - 0022-1120 .- 1469-7645.
-
Tidskriftsartikel (refereegranskat)abstract
- We study numerically the inertial migration of a single rigid sphere and an oblate spheroid in straight square and rectangular ducts. A highly accurate interface-resolved numerical algorithm is employed to analyse the entire migration dynamics of the oblate particle and compare it with that of the sphere. Similarly to the inertial focusing of spheres, the oblate particle reaches one of the four face-centred equilibrium positions, however they are vertically aligned with the axis of symmetry in the spanwise direction. In addition, the lateral trajectories of spheres and oblates collapse into an equilibrium manifold before ending at the equilibrium positions, with the equilibrium manifold tangential to lines of constant background shear for both sphere and oblate particles. The differences between the migration of the oblate and sphere are also presented, in particular the oblate may focus on the diagonal symmetry line of the duct cross-section, close to one of the corners, if its diameter is larger than a certain threshold. Moreover, we show that the final orientation and rotation of the oblate exhibit a chaotic behaviour for Reynolds numbers beyond a critical value. Finally, we document that the lateral motion of the oblate particle is less uniform than that of the spherical particle due to its evident tumbling motion throughout the migration. In a square duct, the strong tumbling motion of the oblate in the first stage of the migration results in a lower lateral velocity and consequently longer focusing length with respect to that of the spherical particle. The opposite is true in a rectangular duct where the higher lateral velocity of the oblate in the second stage of the migration, with negligible tumbling, gives rise to shorter focusing lengths.These results can help the design of microfluidic systems for bio-applications.
|
|
| 10. |
- Lashgari, Iman, et al.
(författare)
-
Laminar, Turbulent, and Inertial Shear-Thickening Regimes in Channel Flow of Neutrally Buoyant Particle Suspensions
- 2014
-
Ingår i: Physical Review Letters. - 0031-9007 .- 1079-7114. ; 113:25, s. 254502-
-
Tidskriftsartikel (refereegranskat)abstract
- The aim of this Letter is to characterize the flow regimes of suspensions of finite-size rigid particles in a viscous fluid at finite inertia. We explore the system behavior as a function of the particle volume fraction and the Reynolds number (the ratio of flow and particle inertia to viscous forces). Unlike single-phase flows, where a clear distinction exists between the laminar and the turbulent states, three different regimes can be identified in the presence of a particulate phase, with smooth transitions between them. At low volume fractions, the flow becomes turbulent when increasing the Reynolds number, transitioning from the laminar regime dominated by viscous forces to the turbulent regime characterized by enhanced momentum transport by turbulent eddies. At larger volume fractions, we identify a new regime characterized by an even larger increase of the wall friction. The wall friction increases with the Reynolds number (inertial effects) while the turbulent transport is weakly affected, as in a state of intense inertial shear thickening. This state may prevent the transition to a fully turbulent regime at arbitrary high speed of the flow.
|
|
