# Deformation of a Capsule in a Power-Law Shear Flow.

1. IntroductionFlow induced deformation of a capsule consisting of a membrane enclosing an internal medium such as a gel or a liquid is an important problem in fundamental research as well as bioengineering applications. For example, a capsule in shear flow is a fundamental process that is related to erythrocytes (or red blood cells), leukocytes (or white blood cells), and platelets in blood flow [1-6]. Deformation is essential for red blood cells to perform their physiological functions in the circulation of capillaryblood vessels and thus affects the rheology of the blood [6-8]. The deformations of white blood cells and red blood cells can, respectively, affect the immune response and the oxygen load release [9, 10]. The synthetic microcapsules with polymerized interfaces are designed for drug delivery, cosmetic production, and other technical usages [11, 12]. Therefore, great effort has been made to studythis problem (e.g., [1, 4, 6, 8, 10, 12-14]).

Both experimental and numerical methods have been conducted to observe capsule behaviors and the relevant underneath fluid-structure interaction physics. Schmid-Schonbein and Wells [15] and Goldsmith [16] observed that red blood cells tumble like rigid particles at low shear rates while they deform to a steady configuration and direction after which the membrane rotates around the internal liquid (tank-treading movement) at high shear rates. Later, Goldsmith and Marlow [17] and Keller and Skalak [18] found that the viscosity ratio between the liquids inside and outside the cell may also affect the type of behaviors. A higher viscosity inside would cause unsteady tumbling-rotating motion, while a smaller viscosity inside would lead to the tank-treading movement with a stationary shape. These phenomena were captured by Xu et al. [14]. More recently, Dupire et al. [19] reported rolling motion in addition to other behaviors. A hysteresis cycle and two transient dynamics driven by the shear rate (i.e., an intermittent regime during the "tank-treading-to-flipping" transition and a Frisbee-like "spinning" regime during the "rolling-to-tank-treading" transition) were highlighted.

There are several numerical methods that have been used to study capsule dynamics. Examples are the boundary element method (e.g., [20]), arbitrary Lagrangian-Euler method (e.g., [21-23]), immersed finite element method (e.g., [24]), and immersed boundary method (IBM) (e.g., [12-14, 25-34]). Specifically, Zhou and Pozrikidis [20] studied the transient and large deformation of capsules with position-dependent membrane tension. Choi and Kim [21] simulated the motion of red blood cells freely suspended in shear flow to investigate the nature of pairwise interception of red blood cells using a fluid-particle interaction method based on the arbitrary Lagrangian-Eulerian method. Villone et al. [22, 23] studied the effect of the non-Newtonian fluid on flexible particle deformation and migration in shear and channel flows by using the arbitrary Lagrangian-Eulerian method. The Navier-Stokes equations and cell-cell interaction were coupled in the framework of the immersed finite element method and mesh-free method by Y. Liu and W. K. Liu [24] to model complex blood flows with deformable red blood cells within micro and capillary vessels in three dimensions. The transient deformation of a liquid-filled elastic capsule in simple shear flow was studied by Sui et al. [1, 4, 35, 36]. The fluid inertia on the dynamics of deformable particles has been studied by Krnger et al. [32] and Kaoui and Harting [34]. More recently, optical force based separation of particles/capsules was simulated by Chang et al. [37-39]. Still, as far as known to us, the existing numerical simulations seldom consider the non-Newtonian rheology effects on the capsule behaviors, while blood and most fluids involved in biomedical engineering are non-Newtonian fluids [6, 8, 40, 41].

Following the work by Sui et al. [1] and Xu et al. [14], we develop an immersed boundary-lattice Boltzmann method (IB-LBM) to study the non-Newtonian effects on the deformation of a capsule in a shear flow. As a typical rheology, the power-law fluid is used. In this method, the capsuledynamicsand thefluiddynamicsare coupledbyusing the IBM, and the incompressible viscous power-law fluid motion is acquired by solving the lattice Boltzmann equation (LBE).

The rest of this paper is organized as follows. Section 2 briefly introduces the governing equations of the fluid and solid structures and describes the numerical approach. Section 3 presents the numerical results. Final conclusions are given in Section 4.

2. Mathematical Formulation and Numerical Method

2.1. Physical Model and Mathematical Formulation. In this work, a two-dimensional liquid capsule enclosed by an elastic membrane and immersed in an incompressible non-Newtonian fluid is considered, as illustrated in Figure 1 where s is the arch length coordinate, n denotes the surface normal that points into the outer fluid, t denotes the tangent unit vector that points to the increasing arc length, and [U.sub.0] is the velocity applied at both top and bottom walls to form a simple shear flow. The incompressible non-Newtonian fluid dynamics is achieved by using LBM [42, 43]. Great effort has been made in using LBM to solve the complex flows (see several reviews [42-44] for the effort). Many publications have presented the details of LBM; thus we just provide a brief description in this paper and discuss the extension for non-Newtonian fluids. The details of LBM and its applications are referred to the references provided. Using the IB-LBM, the lattice Boltzmann equation (LBE) that governs the viscous flow dynamics and incorporates the traction jump across the interface due to the elastic membrane is written as [1,14, 42, 43, 45, 46]

[mathematical expression not reproducible], (1)

[mathematical expression not reproducible], (2)

[mathematical expression not reproducible], (3)

f (x, t) = [integral] [DELTA]F (s, t) [[delta].sub.D] (x - X (s, t)) ds, (4)

where [g.sub.i] (x, t) is the distribution function for particles with velocity [e.sub.i] at position x and time t, [DELTA]t is the size of the time step, [g.sub.i.sup.eq] (x, t) is the equilibrium distribution function, [[tau].sub.LB] represents the dimensionless relaxation time, [G.sub.i] is the term representing the body force effect on the distribution function, [[omega].sub.i] are the weights, u = (u, v) is the velocity of the fluid, [c.sub.s] is the speed of sound defined by [c.sub.s] = [DELTA]x/[square root of 3[DELTA]t] with [DELTA]x being grid spacing, f is the body force acting on the fluid, [DELTA]F(s, t) is the Lagrangian force density on the fluid by the elastic boundary, X is the position vector on the membrane, and [[delta].sub.D] (x - X (s, t)) is Dirac's delta function.

In the present work, a two-dimensional nine-speed (D2Q9) model is used, as shown in Figure 2. In this model, the nine possible particle velocities are given by

[mathematical expression not reproducible]. (5)

The values of [e.sub.i] ensure that, within one time step, a particle moves to one of the eight neighboring nodes as shown in Figure 2 or stays at its current location. The weights, [[omega].sub.i], are given by [[omega].sub.0] = 4/9 and [[omega].sub.i] = 1/9 for i = 1 to 4 and [[omega].sub.i] = 1/36 for i = 5 to 8. In addition, the relaxation time is related to the kinematic viscosity in the Navier-Stokes equations in terms of

[[tau].sub.LB] = 0.5 + [v/[c.sup.2.sub.s][DELTA]t], (6)

where v = [mu]/[rho] with [mu] being the dynamic viscosity of the ambient fluid and [rho] being the fluid density.

When the particle density distributions are known, the fluid density, velocity, and pressure are then computed from

[mathematical expression not reproducible], (7)

Theoretically the LBM introduced above simulates the compressible viscous flow instead of incompressible viscous one, because the spatial density variation is not zero in LBM simulations. In the applications, the Mach number (Ma = [u.sub.0]/[c.sub.s]) should be low (e.g., Ma [less than or equal to] 0.3) so that the incompressible viscous flow can be correctly simulated. The deduction process from LBE to the incompressible viscous flow governing equations can be found in [47].

The dynamics viscosity is a constant for a Newtonian fluid, while it is dependent on the local shear rate for a non-Newtonian fluid. Without loss of generality, the power-law fluid is taken as a representation of non-Newtonian fluids in the present paper. The rheological equation of state for power-law fluids is defined by [48]

[mu] = [eta][[??].sup.n-1], (8)

[??] = max ([square root of 2[E.sub.ij] [E.sub.ij]], [[??].sub.m]), (9)

where [eta] is the power-law consistency index, n is the power-law fluid behavior index, [??] is the shear rate, and [[??].sub.m] is the minimum shear rate that is applied to avoid the numerical singularity caused by the zero shear rate. The power-law fluids of n < 1, n > 1 and n = 1 are, respectively, the shear-thinning, shear-thickening, and Newtonian fluids. In (9), the Einstein summation convention is applied. In LBM implementation, [E.sub.ij] can be either calculated macroscopicallyby using the finite difference method or locally in mesoscopic scale by using [g.sup.eq] and f [49]. To achieve the non-Newtonian rheology, a shear rate-dependant relaxation time is used which can be obtained by applying the effective viscosity determined by (8) in (6).

Because of the deformation, the membrane develops a transverse shear tension q and a bending moment m. In addition, due to the stretching motion, a tension, [tau], is induced. Consider the force balance of membranes; we acquire

[mathematical expression not reproducible], (11)

Please note that [DELTA]F is the Lagrangian force on the fluid exerted by the elastic body boundary and is opposite to the fluid force on the boundary. To evaluate m and [tau] for the thin membrane, we use Hooke's law which is a relatively simpler constitutive law for modeling small deformation of capsules. Hooke's law states that the tension and the bending moment are linearly related to the stretch and the curvature, respectively. It can be written in the form

[mathematical expression not reproducible], (12)

where [E.sub.B] is the bending coefficient, [E.sub.S] is the stretching coefficient, [s.sub.0] is the initial arch length, [kappa] is the curvature of membrane, and [[kappa].sub.0] is the curvature in the resting configuration. If Es is large so that the stretching deformation is small, Hooke's law works well. Following the work by Sui et al. [1], the capsule membrane is assumed to be infinitely thin so that the bending effect is neglected; that is, [E.sub.B] = 0. Actually, the effect of [E.sub.B] is similar to that of [E.sub.S] when [E.sub.B] is small compared to [E.sub.S] [1, 35]. If [E.sub.B] is large, the capsule may undergo flipping motion [35].

The velocity of a point on the capsule is interpolated from the flow field, and the position of the capsule is updated explicitly; that is,

[mathematical expression not reproducible], (13)

[partial derivative]X (s, t)/[partial derivative]t] = U (s, t), (14)

where U(s, t) is the velocity of the capsule.

In this work, we choose the flow shear rate (e.g., [U.sub.0]/h), density, and the radius of the capsule to nondimension the governing equations and obtain two dimensionless parameters: the Reynolds number Re and dimensionless shear rate G, which are defined by

[mathematical expression not reproducible], (15)

G = [[eta]a/[E.sub.S]] [([U.sub.0]/h]).sup.n], (16)

where a is the radius of the undeformed capsule. G measures the ratio of shear force to the elastic force. For applications where inertia force is important, we can also use Re x G to nondimension the elastic property, which measures the ratio of fluid inertial forces to stretching elastic forces. Please note that the two-dimensional model is used in this work, while red blood cell deformation is a three-dimensional problem; however, the results obtained in this research should show some features common with the three-dimensional simulations, as demonstrated in [1].

2.2. Numerical Method. Similar to [1], the capsule is discretized by [N.sub.f] nodal points which are initially distributed with equal distances. The position of the mth node at time level n is denoted by [X.sup.n.sub.m]. To compute the stretching force at mth node, a finite difference scheme is used; that is,

[mathematical expression not reproducible], (17)

where [DELTA]s is the Lagrangian grid spacing on the membrane and the tension [tau] and tangent vector, t = [partial derivative]X/[partial derivative]s, at the segment center, m + 1/2, are both computed using a secondorder central difference scheme.

The time integration of (14) is calculated according to

[X.sup.n+1] = [X.sup.n] + [DELTA]t[U.sup.n+1]. (18)

In the IBM, a smooth approximation [50] of Dirac's delta function, [[delta].sub.h], is used,

[mathematical expression not reproducible], (19)

In the present simulations, [DELTA]x = [DELTA]y = [DELTA]t (in lattice units) is used.

Now, the computational algorithm can be summarized as follows. Given all values at time step n, the values at time step n+1 can be undated by the following:

(1) Calculate the Lagrangian force density [DELTA][F.sup.n+1] from [X.sup.n] by using (11)-(12).

(2) Spread the Lagrangian force density [DELTA][F.sup.n+1] onto the ambient fluid nodes by using (4), and obtain [f.sup.n+1].

(3) Solve flow field with body force by using the LBM method described by (1)-(3) and (6)-(10).

(4) Update [U.sup.n+1] by using (13).

(5) And finally, update [X.sup.n+1] by using (18).

In the present work, the above-mentioned computational simulation algorithm is implemented in the Fortran 90 programming language.

2.3. Validation. The IB-LBM in this work has been validated and verified in our previous studies (see, e.g., [14, 46]) and has been used to study filament(s) flapping in viscous fluids [51-53], sperm swimming, and cell/particle flows [10, 54]. In the present work, we focus on the validation of non-Newtonian flow by considering a power-law flow in a straight channel which is one of the benchmark problems to validate an in-house computational fluid dynamics solver. As in our previous work [41], we consider a two-dimensional steady laminar developing flow of power-law fluid with a uniform incoming velocity [U.sub.[infinity]] through a rectangular channel of height h and length L, as shown in Figure 3. The physically realistic initial and boundary conditions are given as

[mathematical expression not reproducible], (20)

[mathematical expression not reproducible], (21)

[mathematical expression not reproducible], (22)

[mathematical expression not reproducible], (23)

The computations are performed with the dimensionless domain size (L/h x 1) of 40 x 1 discretized by 2001 x 51 uniform Cartesian nodes. The numerical results in terms of the fully developed velocity profiles are obtained for the Reynolds number (defined by [rho][h.sup.n][U.sup.2-n.sub.[infinity]]/[eta]) of 100 and for three power-law indices; that is, n = 0.6, 1.0, and 1.4. The simulations are performed for sufficiently long time so that the flow in the channel attains a steady state. The fully developed velocity profiles predicted by the numerical simulations are compared in Figure 4 with the corresponding analytical solution for fully developed velocity profile [41,48] for power-law fluid flow in a channel which is given as

[mathematical expression not reproducible]. (24)

From Figure 4, it is found that the present numerical results show a good agreement with the analytical solutions for various values of power-law index, giving us confidence in the reliability and accuracy of the present numerical solution procedure. It is noted from Figure 4 that the shear layer is thinned for n < 1 and thickened for n > 1 compared to the Newtonian fluid case (n = 1).

3. Numerical Results

We first consider the power-law index effect on the deformation of a cylindrical capsule in a shear flow. The Reynolds number is 0.05, which is in the range of normal physiological conditions. The dimensionless shear rate G is 0.04. The computational domain ranges from 0 to 20a in both x-axis and y-axis. The capsule is at the center of the domain, and its membrane is equally discretized into 80 Lagrangian nodes. The grid resolution is [DELTA]x = [DELTA]y = [DELTA]t = 0.1a. The characteristic velocity is set as [U.sub.0] = 0.05 so that the dimensionless relaxation time is 0.5 < [[tau].sub.LB] < 3.0. Such setup is consistent with that used in [1]. To study the power-law index effect, n is set in the range of 0.2 < n < 1.8, covering the shear-thinning, Newtonian, and shear-thickening fluids. In order to quantify the deformation of the capsule, the Taylor shape parameter [D.sub.xy] is introduced [1],

[D.sub.xy] = [L - B/L + B], (25)

where L and B are, respectively, the length and width of a cross-section of the cylindrical capsule.

Figure 5 shows the deformation of the flexible capsule in a shear flow for Reynolds number of Re = 0.05, dimensionless shear rate of G = 0.04, and power-law index of n = 0.2 to 1.8. There are several interesting observations from Figure 5. First, the capsule deforms to a steady shape and then the membrane rotates around the liquid inside (tank-treading motion), which is further indicated by the streamlines in Figure 6. Second, the deformation increases with the powerlaw index. When the fluid is shear-thinning (i.e., n < 1.0), the deformation is smaller compared to the Newtonian fluid case (n = 1.0), while the deformation is larger compared to the Newtonian fluid case for the shear-thickening fluid; that is, n > 1.0. This can be explained by the power-law rheology. When n < 1.0, the effective viscosity near the capsule is smaller compared to the Newtonian fluid, while the effective viscosity near the capsule is higher than that of Newtonian fluid for n> 1.0. Based on the definition of G in (16), the local G is larger for n > 1.0 and smaller for n < 1.0 compared to that of Newtonian fluid. As presented by Sui et al. [1], a larger G corresponds to a larger [G.sub.xy], that is, larger deformation of the capsule. Third, the Taylor shape parameter [G.sub.xy], which is used to quantify the deformation, increases with the power-law index. Finally, it is noted that [G.sub.xy] is approximately linear function of n, as shown in Figure 5(b).

In order to study the Reynolds number effect on the deformation of the capsule, we simulate two additional Reynolds numbers, Re = 0.1 and 0.025. Figure 7 shows the deformation of the flexible capsule in a shear flow for Reynolds number of Re = 0.1 and 0.025, dimensionless shear rate of G = 0.04, and power-law index of 0.2 [less than or equal to] n [less than or equal to] 1.8. It is found that the deformation ([D.sub.xy]) for Re = 0.1 is larger compared to the cases of Re = 0.05 and 0.025. However, the difference is quite small, implying that, in the low Reynolds number regime, for example, Re [less than or equal to] 0.1 in this work, the deformation of the capsule is not significantly affected by the Reynolds numbers used, as the inertial force is ignorable here, and the shear forces and capsule elastic forces are dominant. Therefore, the dimensionless shear rate (G) should significantly affect the deformation of the capsule, which will be further verified by the simulations shown below by varying G.

Finally, we study the shear rate effect on the deformation of the capsule by using G = 0.004 and 0.1 at Re = 0.05. The deformation of the flexible capsule in a shear flow for dimensionless shear rate of G = 0.1 and 0.004, Reynolds number of Re = 0.05, and power-law index of n = 0.2 to 1.8 is shown in Figure 8, from which several interesting observations are obtained. First, the capsule deformation is sensitive to the dimensionless shear rate. This can be explained by the definition of G in (16): G measures the ratio of shear (viscous) forces to the stretching elastic forces, which is the dominant physical process here. A change of this ratio would cause significant difference in the capsule deformation. Second, the power-law index effect is stronger for larger G, as indicated by the slopes of the [G.sub.xy] functions shown in Figure 8(c). This can be explained by the fact that the physical process changes from a shear force dominant to an elastic-force dominant process when G varies from 0.1 to 0.004. For low G, for example, 0.004, the elastic forces are dominant, and thus the shear force change caused by the change of the non-Newtonian rheology is smaller compared to that for larger G, for example, 0.1. Finally, the deformed capsule is obviously biased from elliptical cylinder for large G and n; for example, G = 0.1 and n [greater than or equal to] 1.2. This is caused by the shear-induced torque on the deformed capsule and the decrease of effective bending resistance caused by the shear-induced elongation.

To further discuss the non-Newtonian effect, [chi] = -[mu][omega], which measures the local shear stress, is introduced [40]. Figure 9 shows contours of [chi] for G = 0.04, Re = 0.025, n = 0.6, and 1.4. It finds that [chi] near the long axial ends is larger (i.e., the local shear stress is larger) for n = 1.4 compared to that for n = 0.6. This is a further explanation of larger deformation for larger n.

4. Conclusion

A numerical approach combining the immersed boundary method and the lattice Boltzmann method has been developed for fluid-structure interactions involving non-Newtonian fluids. Without loss of generality, the power-law fluid is taken as a representation of non-Newtonian fluids to present the method. This method couples the flexible structure (e.g., capsule) dynamics and the fluid dynamics by using the immersed boundary method and calculates the incompressible viscous power-law fluid motion by solving the lattice Boltzmann equation. In order to achieve the non-Newtonian rheology, a shear rate-dependant relaxation time is employed.

The non-Newtonian flow solver has been validated by conducting a power-law flow in a straight channel. The power-law index has been varied from 0.6 to 1.4. The present numerical results show a good agreement with the analytical solutions for various values of power-law index, giving us confidence in the reliability and accuracy of the present numerical solution procedure.

To study the non-Newtonian effects on the deformation of a capsule in a power-law shear flow, we have performed simulations by varying the Reynolds number from 0.025 to 0.1, dimensionless shear rate from 0.004 to 0.1, and power-law index from 0.2 to 1.8. It is found that the capsule deformation increases with the power-law index for different Reynolds numbers and nondimensional shear rates. In addition, the Reynolds number does not have significant effect on the capsule deformation in the flow regime considered. Finally, the power-law index effect is stronger for larger dimensionless shear rate compared to smaller values.

http://dx.doi.org/10.1155/2016/7981386

Competing Interests

The author declares that they have no competing interests.

Acknowledgments

Dr. Fang-Bao Tian is the recipient of an Australian Research Council Discovery Early Career Researcher Award (Project no. DE160101098). This research was partly undertaken on the NCI National Facility in Canberra, Australia, which is supported by the Australian Commonwealth Government.

References

[1] Y. Sui, Y. T. Chew, and H. T. Low, "A lattice Boltzmann study on the large deformation of red blood cells in shear flow," International Journal of Modern Physics C, vol. 18, no. 6, pp. 993-1011, 2007.

[2] J. Zhang, P. C. Johnson, and A. S. Popel, "An immersed boundary lattice Boltzmann approach to simulate deformable liquid capsules and its application to microscopic blood flows," Physical Biology, vol. 4, no. 4, pp. 285-295, 2007.

[3] J. O. Barber, J. P. Alberding, J. M. Restrepo, and T. W. Secomb, "Simulated two-dimensional red blood cell motion, deformation, and partitioning in microvessel bifurcations," Annals of Biomedical Engineering, vol. 36, no. 10, pp. 1690-1698, 2008.

[4] Y. Sui, Y. T. Chew, P. Roy, Y. P. Cheng, and H. T. Low, "Dynamic motion of red blood cells in simple shear flow," Physics of Fluids, vol. 20, no. 11, Article ID 112106, 2008.

[5] T. Ye, H. Li, and K. Y. Lam, "Modeling and simulation of microfluid effects on deformation behavior of a red blood cell in a capillary," Microvascular Research, vol. 80, no. 3, pp. 453-463, 2010.

[6] J. B. Freund, "Numerical simulation of flowing blood cells," Annual Review of Fluid Mechanics, vol. 46, pp. 67-95, 2014.

[7] H. Schmid-Schonbein, R. Wells, and J. Goldstone, "Influence of deformability of human red cells upon blood viscosity," Circulation Research, vol. 25, no. 2, pp. 131-143, 1969.

[8] D. Barthes-Biesel, "Motion and deformation of elastic capsules and vesicles in flow," Annual Review of Fluid Mechanics, vol. 48, pp. 25-52, 2016.

[9] F. A. Brayner, H. R. C. Araujo, S. S. Santos et al., "Haemocyte population and ultrastructural changes during the immune response of the mosquito Culex quinquefasciatus to microfilariae of Wuchereria bancrofti," Medical and Veterinary Entomology, vol. 21, no. 1, pp. 112-120, 2007.

[10] X. Yuan-Qing, T. Xiao-Ying, T. Fang-Bao, P. Yu-Hua, X. Yong, and Z. Yan-Jun, "IB-LBM simulation of the haemocyte dynamics in a stenotic capillary," Computer Methods in Biomechanics and Biomedical Engineering, vol. 17, no. 9, pp. 978-985, 2014.

[11] G. O. Fanger, "What good are microcapsules," Chemtech, vol. 7, pp. 397-405, 1974.

[12] W.-X. Huang, C. B. Chang, and H. J. Sung, "Three-dimensional simulation of elastic capsules in shear flow by the penalty immersed boundary method," Journal of Computational Physics, vol. 231, no. 8, pp. 3340-3364, 2012.

[13] Y. Sui, Y. T. Chew, P. Roy, and H. T. Low, "A hybrid method to study flow-induced deformation of three-dimensional capsules," Journal of Computational Physics, vol. 227, no. 12, pp. 6351-6371, 2008.

[14] Y.-Q. Xu, F.-B. Tian, and Y.-L. Deng, "An efficient red blood cell model in the frame of IB-LBM and its application," International Journal of Biomathematics, vol. 6, no. 1, Article ID 1250061, 2013.

[15] H. Schmid-Schonbein and R. Wells, "Fluid drop-like transition of erythrocytes under shear," Science, vol. 165, no. 3890, pp. 288-291, 1969.

[16] H. L. Goldsmith, "Deformation of human red cells in tube flow," Biorheology, vol. 7, no. 4, pp. 235-242, 1971.

[17] H. L. Goldsmith and J. Marlow, "Flow behaviour of erythrocytes. I. Rotation and deformation in dilute suspensions," Proceedings of the Royal Society B: Biological Sciences, vol. 182, no. 1068, pp. 351-384, 1972.

[18] S. R. Keller and R. Skalak, "Motion of a tank-treading ellipsoidal particle in a shear flow," Journal of Fluid Mechanics, vol. 120, pp. 27-47, 1982.

[19] J. Dupire, M. Socol, and A. Viallat, "Full dynamics of a red blood cell in shear flow," Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 51, pp. 20808-20813, 2012.

[20] H. Zhou and C. Pozrikidis, "Deformation of liquid capsules with incompressible interfaces in simple shear flow," Journal of Fluid Mechanics, vol. 283, pp. 175-200, 1995.

[21] C. R. Choi and C. N. Kim, "Fluid-particle interaction simulations ofthe interception of red blood cells in shear flow," Journal of Mechanical Science and Technology, vol. 24, no. 7, pp. 14511463, 2010.

[22] M. M. Villone, F. Greco, M. A. Hulsen, and P. L. Maffettone, "Simulations of an elastic particle in Newtonian and viscoelastic fluids subjected to confined shear flow," Journal of Non-Newtonian Fluid Mechanics, vol. 210, pp. 47-55, 2014.

[23] M. M. Villone, F. Greco, M. A. Hulsen, and P. L. Maffettone, "Numerical simulations of deformable particle lateral migration in tube flow of Newtonian and viscoelastic media," Journal of Non-Newtonian Fluid Mechanics, vol. 234, pp. 105-113, 2016.

[24] Y. Liu and W. K. Liu, "Rheology of red blood cell aggregation by computer simulation," Journal of Computational Physics, vol. 220, no. 1, pp. 139-154, 2006.

[25] C. S. Peskin, Flow patterns around heart valves: a digital computer method for solving the equations of motion [Ph.D. thesis], Yeshiva University, New York, NY, USA, 1972.

[26] M. F. McCracken and C. S. Peskin, "A vortex method for blood flow through heart valves," Journal of Computational Physics, vol. 35, no. 2, pp. 183-205, 1980.

[27] C. S. Peskin and D. M. McQueen, "Fluid dynamics of the heart and its valves," in Case Studies in Mathematical Modeling: Ecology, Physiology, and Cell Biology, H. G. Othmer, F. R. Adler, M. A. Lewis, and J. C. Dallon, Eds., pp. 309-337, Prentice-Hall, Englewood, NJ, USA, 1996.

[28] A. M. Roma, C. S. Peskin, and M. J. Berger, "An adaptive version of the immersed boundary method," Journal of Computational Physics, vol. 153, no. 2, pp. 509-534, 1999.

[29] A. Alexeev, R. Verberg, and A. C. Balazs, "Modeling the motion of microcapsules on compliant polymeric surfaces," Macromolecules, vol. 38, no. 24, pp. 10244-10260, 2005.

[30] Y. Mori and C. S. Peskin, "Implicit second-order immersed boundary methods with boundary mass," Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 25-28, pp. 2049-2067, 2008.

[31] Y. Kim, Y. Seol, M.-C. Lai, and C. S. Peskin, "The immersed boundary method for two-dimensional foam with topological changes," Communications in Computational Physics, vol. 12, no. 2, pp. 479-493, 2012.

[32] T. Kruger, B. Kaoui, and J. Harting, "Interplay of inertia and deformability on rheological properties of a suspension of capsules," Journal of Fluid Mechanics, vol. 751, pp. 725-745, 2014.

[33] H. Zhang, F. X. Trias, A. Oliva et al., "PIBM: particulate immersed boundary method for fluid-particle interaction problems," Powder Technology, vol. 272, pp. 1-13, 2015.

[34] B. Kaoui and J. Harting, "Two-dimensional lattice Boltzmann simulations of vesicles with viscosity contrast," Rheologica Acta, vol. 55, pp. 465-475, 2016.

[35] Y. Sui, Y. T. Chew, P. Roy, X. B. Chen, and H. T. Low, "Transient deformation of elastic capsules in shear flow: effect of membrane bending stiffness," Physical Review E, vol. 75, no. 6, Article ID 066301, 2007.

[36] Y. Sui, H. T. Low, Y. T. Chew, and P. Roy, "Tank-treading, swinging, and tumbling of liquid-filled elastic capsules in shear flow," Physical Review E, vol. 77, no. 1, Article ID 016310, 2008.

[37] C. B. Chang, W.-X. Huang, and H. J. Sung, "Lateral migration of an elastic capsule by optical force in a uniform flow," Physical ReviewE--Statistical, Nonlinear, and Soft Matter Physics, vol. 86, no. 6, Article ID 066306, 2012.

[38] C. B. Chang, W.-X. Huang, K. H. Lee, and H. J. Sung, "Optical separation of ellipsoidal particles in a uniform flow," Physics of Fluids, vol. 26, no. 6, Article ID 062001, 2014.

[39] C. B. Chang, W.-X. Huang, and H. J. Sung, "Cross-type optical separation of elastic oblate capsules in a uniform flow," Journal of Applied Physics, vol. 117, no. 3, Article ID 034701, 2015.

[40] F.-B. Tian, L. Zhu, P.-W. Fok, and X.-Y. Lu, "Simulation of a pulsatile non-Newtonian flow past a stenosed 2D artery with atherosclerosis," Computers in Biology and Medicine, vol. 43, no. 9, pp. 1098-1113, 2013.

[41] F.-B. Tian, R. P. Bharti, and Y.-Q. Xu, "Deforming-spatial-domain/stabilized space-time (DSD/SST) method in computation of non-Newtonian fluid flow and heat transfer with moving boundaries," Computational Mechanics, vol. 53, no. 2, pp. 257-271, 2014.

[42] S. Chen and G. D. Doolen, "Lattice Boltzmann method for fluid flows," Annual Review of Fluid Mechanics, vol. 30, pp. 329-364, 1998.

[43] C. K. Aidun and J. R. Clausen, "Lattice-Boltzmann method for complex flows," Annual Review of Fluid Mechanics, vol. 42, pp. 439-472, 2010.

[44] H.-B. Deng, Y.-Q. Xu, D.-D. Chen, H. Dai, J. Wu, and F.-B. Tian, "On numerical modeling of animal swimming and flight," Computational Mechanics, vol. 52, no. 6, pp. 1221-1242, 2013.

[45] J. Wu and C. Shu, "Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications," Journal of Computational Physics, vol. 228, no. 6, pp. 1963-1979, 2009.

[46] F.-B. Tian, H. Luo, L. Zhu, J. C. Liao, and X.-Y. Lu, "An efficient immersed boundary-lattice Boltzmann method for the hydrodynamic interaction of elastic filaments," Journal of Computational Physics, vol. 230, no. 19, pp. 7266-7283, 2011.

[47] Z. L. Guo, C. G. Zheng, and B. C. Shi, "Discrete lattice effects on the forcing term in the lattice Boltzmann method," Physical Review E, vol. 65, no. 4, Article ID 046308, 2002.

[48] R. P. Chhabra and J. F. Richardson, Non-Newtonian Flow and Applied Rheology, Butterworth-Heinemann, Oxford, UK, 2nd edition, 2008.

[49] A. A. Delouei, M. Nazari, M. H. Kayhani, S. K. Kang, and S. Succi, "Non-Newtonian particulate flow simulation: a directforcing immersed boundary-lattice Boltzmann approach," PhysicaA: Statistical Mechanics and Its Applications, vol. 447, pp. 1-20, 2016.

[50] C. S. Peskin, "The immersed boundary method," Acta Numerica, vol. 11, pp. 479-517, 2002.

[51] F.-B. Tian, H. Luo, L. Zhu, and X.-Y. Lu, "Interaction between a flexible filament and a downstream rigid body," Physical Review E, vol. 82, no. 2, Article ID 026301, 2010.

[52] F.-B. Tian, H. Luo, L. Zhu, and X.-Y. Lu, "Coupling modes of three filaments in side-by-side arrangement," Physics of Fluids, vol. 23, no. 11, Article ID 111903, 2011.

[53] F.-B. Tian, "Role of mass on the stability of flag/flags in uniform flow," Applied Physics Letters, vol. 103, no. 3, Article ID 034101, 2013.

[54] Y. Xu, F. Tian, H. Li, and Y. Deng, "Red blood cell partitioning and blood flux redistribution in microvascular bifurcation," Theoretical and Applied Mechanics Letters, vol. 2, no. 2, Article ID 024001, 2012.

Fang-Bao Tian

School of Engineering and Information Technology, University of New South Wales, Canberra, ACT 2600, Australia

Correspondence should be addressed to Fang-Bao Tian; f.tian@adfa.edu.au

Received 28 April 2016; Revised 10 July 2016; Accepted 11 August 2016

Academic Editor: Giuseppe Pontrelli

Caption: Figure 1: Schematic illustration of a circular liquid capsule immersed in a fluid.

Caption: Figure 2: Nine base vectors representing 9 possible velocity directions in the D2Q9 lattice model.

Caption: Figure 3: Sketch of the power-law fluid flow in a channel.

Caption: Figure 4: Comparison of the present numerical results of the fully developed velocity profiles in a channel with the corresponding analytical profiles (24) for Reynolds number of Re =100 and power-law indices n = 0.6, 1.0, and 1.4.

Caption: Figure 5: Deformation of the flexible capsule in a shear flow for Reynolds number of Re = 0.05, dimensionless shear rate of G = 0.04, and power-law index of n = 0.2 to 1.8: (a) capsule shapes for difference power-law indices (the dashed line is for the Newtonian fluid case where n = 1.0) and (b) Taylor shape parameter as a function of power-law index.

Caption: Figure 6: The streamline pattern inside and outside the capsule at steady state for Re = 0.05, G = 0.04, and n = 1.0.

Caption: Figure 7: Deformation of the flexible capsule in a shear flow for Reynolds number of Re = 0.1 and 0.025, dimensionless shear rate of G = 0.04, and power-law index of n = 0.2 to 1.8: (a) capsule shapes for difference power-law indices at Re = 0.1 (the dashed line is for the Newtonian fluid case where n = 1.0), (b) capsule shapes for difference power-law indices at Re = 0.025, and (c) Taylor shape parameter as a function of power-law index. [G.sub.xy] for Re = 0.05 is shown in (c) for comparison.

Caption: Figure 8: Deformation of the flexible capsule in a shear flow for dimensionless shear rate of G = 0.1 and 0.004, Reynolds number of Re = 0.05, and power-law index of n = 0.2 to 1.8: (a) capsule shapes for difference power-law indices at G = 0.1 (the dashed line is for the Newtonian fluid case where n = 1.0), (b) capsule shapes for difference power-law indices at G = 0.004, and (c) Taylor shape parameter as a function of power-law index. [G.sub.xy] for G = 0.04 is shown in (c) for comparison.

Caption: Figure 9: Contours of [chi] = -[mu][omega] for G = 0.04 and Re = 0.025: (a) n = 0.6 and (b) n = 1.4.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Tian, Fang-Bao |

Publication: | Computational and Mathematical Methods in Medicine |

Date: | Jan 1, 2016 |

Words: | 6072 |

Previous Article: | Analysis of Blood Transfusion Data Using Bivariate Zero-Inflated Poisson Model: A Bayesian Approach. |

Next Article: | Improved CEEMDAN and PSO-SVR Modeling for Near-Infrared Noninvasive Glucose Detection. |

Topics: |