Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2254372

Formats

Article sections

- Abstract
- 1 Introduction
- 2 Theory and Methods
- 3 Results and Discussion
- 4 Conclusions
- Supplementary Material
- References

Authors

Related links

J Biomech. Author manuscript; available in PMC 2009 January 1.

Published in final edited form as:

Published online 2007 September 20. doi: 10.1016/j.jbiomech.2007.07.020

PMCID: PMC2254372

NIHMSID: NIHMS37873

July 23, 2007

In this paper we develop a lattice Boltzmann algorithm to simulate red blood cell (RBC) behavior in shear flows. The immersed boundary method is employed to incorporate the fluid-membrane interaction between the flow field and deformable cells. The cell membrane is treated as a neo-Hookean viscoelastic material and a Morse potential is adopted to model the intercellular interaction. Utilizing the available mechanical properties of RBCs, multiple cells have been studied in shear flows using a two-dimensional approximation. These cells aggregate and form a rouleau under the action of intercellular interaction. The equilibrium configuration is related to the interaction strength. The end cells exhibit concave shapes under weak interaction and convex shapes under strong interaction. In shear flows, such a rouleau-like aggregate will rotate or be separated, depending on the relative strengths of the intercellular interaction and hydrodynamic viscous forces. These behaviors are qualitatively similar to experimental observations and show the potential of this numerical scheme for future studies of blood flow in microvessels.

Red blood cells (RBCs) are an important component in blood because of their large number density (~ 5×10^{6}/mm^{3}) and their crucial role in oxygen transport. Typically, a human RBC has a biconcave shape of ~ 8 *μ*m in diameter and ~ 2 *μ*m in thickness. The interior fluid (cytoplasm) has a viscosity of 6 cP, which is about 5 times of that of the suspending plasma (~ 1.2 cP). The cell membrane is highly deformable so that RBCs can pass through capillaries of as small as 4 *μ*m inner diameter with large deformation (Popel and Johnson, 2005; Mchedlishvili and Maeda, 2001). RBCs can also aggregate and form one-dimensional stacks-of-coins-like rouleaux or three-dimensional aggregates (Popel and Johnson, 2005; Stoltz et al., 1999; Baumler et al., 1999). The process is reversible and the rouleaux and aggregates can be broken by, for example, increasing the flow shear rate. This phenomenon is particularly important in microcirculation, since such rouleaux or aggregates can dramatically influence blood flow in microvessels. However, the underlying mechanism of RBC aggregation is not yet clear. Currently, there exist two theoretical descriptions of this process: the bridging model and the depletion model (Popel and Johnson, 2005; Baumler et al., 1999). The former assumes that macromolecules, such as fibrinogen or dextran, can adhere onto the adjacent RBC surfaces and bridge them together (Merill et al., 1966; Brooks, 1973; Chien and Jan, 1973a). The depletion model attributes the RBC aggregation to a polymer depletion layer between RBC surfaces, which is accompanied by a decrease of the osmotic pressure (Baumler and Donath, 1987; Evans and Needham, 1988). Detailed discussions of these two models can be found elsewhere [see, for example, a review by Baumler et al. (1999)].

The nature of blood flow changes greatly with the vessel diameter. In vessels larger than 200 *μ*m, the blood flow can be accurately modeled as a homogeneous fluid. However, in arterioles and venules smaller than 25 *μ*m, and also in capillaries with diameter of 4–10 *μ*m, the RBCs have to be treated as discrete fluid capsules suspended in the plasma. Significant efforts have been devoted to numerically study the RBC behaviors in various flow situations. For example, Pozrikidis and coworkers (1995; 2001; 2003) have employed the boundary integral method for Stokes flows to investigate RBC deformation and motion in shear and channel flows; Eggleton and Popel (1998) have combined the immersed boundary method (IBM) (Peskin, 1977) with a finite element treatment of the RBC membrane to simulate large three-dimensional RBC deformations in shear flow. Recently, a lattice Boltzmann approach has also been adopted for RBC flows in microvessels, where the RBCs were represented as rigid rods in two dimensions (Migliorini et al., 2002; Sun et al., 2003; Sun and Munn, 2005). Bagchi (2007) has simulated a large RBC population in vessels of size 20-300 *μ*m. However, RBC aggregation was not considered in these studies. Liu et al. (2004) modeled the intercellular interaction through a Morse potential, thus accounting for RBC aggregation. The cell membrane was represented by elastic elements. Bagchi et al. (2005) have extended the IBM approach of Eggleton and Popel (1998) to two-cell systems and introduced the intercellular interaction according to a legand-receptor binding model. Chung et al. (2006) have utilized the theoretical formulation of depletion energy proposed by Neu and Meiselman (2002) to study two rigid elliptical particles in a channel flow. Sun and Munn (2006) have also improved their lattice Boltzmann model by including an interaction force between rigid RBCs.

In this paper, we develop a two-dimensional lattice Boltzmann scheme to simulate multiple deformable RBCs in shear flow. The lattice Boltzmann method (LBM) was chosen to solve the incompressible fluid field for its ability to deal with complex boundary conditions and its advantage for parallel computation (Succi, 2001), both of which could be valuable in our future studies. Different from previous LBM studies of RBC flows (Migliorini et al., 2002; Sun et al., 2003; Sun and Munn, 2005), here we modeled the cells as deformable fluid capsules and hence the membrane mechanics and cytoplasm viscosity could be considered. IBM was also utilized to incorporate the fluid-membrane interaction, and a Morse potential to describe the intercellular interaction. Detailed theory and formulations will be given in the next Section. The algorithm and program have been validated by comparing simulation results with theoretical predictions, and excellent agreement was found. Finally, simulations of multiple deformable RBCs were conducted and the results demonstrated the effects of intercellular interactions and shear rate on the RBC rheological behaviors, which are also in qualitative agreement with experimental observations.

Here we employ the lattice Boltzmann model (LBM) (Zhang and Kwok, 2005; Zhang et al., 2004) to solve the flow field and the immersed boundary method (Peskin, 1977) to incorporate the fluid-membrane interaction. The LBM approach is advantageous for parallel computations due to its local dynamics and can be relatively easily applied to systems with complex boundaries (Succi, 2001). Detailed descriptions of these methods are available elsewhere (also see the Supplemental Materials).

Experiments have shown that the RBC membrane is a neo-Hookean, highly deformable viscoelastic material (Hochmuth and Waugh, 1987). Also it exhibits finite bending resistance that becomes more profound in regions with large curvature (Evans, 1983). For the two-dimensional RBC model in this work, according to Bagchi et al. (2005), the neo-Hookean elastic component of membrane stress can be expressed as

(1)

where *E _{s}* is the membrane shear modulus and

(2)

where *μ _{m}* is the membrane viscosity. In addition, the bending resistance can also be represented by relating the membrane curvature change to the membrane stress with (Pozrikidis, 2001; Bagchi et al., 2005)

(3)

where *E _{b}* is the bending modulus and

(4)

Here **t** and **n** are the local tangential and normal directions on the membrane.

The physiological and pathological importance of RBC aggregation has been realized and extensive experimental investigations have been performed (Popel and Johnson, 2005; Stoltz et al., 1999; Baumler et al., 1999; Kim et al., 2006; Kounov and Petrov, 1999; Chien et al., 1977); however, the underlying mechanisms of the RBC aggregation are still subjects of investigation. Both the bridging and the depletion models can describe certain aggregation phenomena, however, they fail to explain some specific observations (Baumler et al., 1999; Armstrong et al., 1999). In general, it can be assumed that the attractive interaction between RBC surfaces would occur when they are close, and repulsive interaction would occur when the separation distance is sufficiently small. The repulsive interaction includes the steric forces due to the glycocalyx and electrostatic repulsion from the negative charges on RBC surfaces (Liu et al., 2004). In previous studies, to model such intercellular interactions, Bagchi et al. (2005) adopted the formalism of a ligand-receptor dynamics according to the bridging model. However, this description involves a number of parameters whose values are not available from experiments. On the other hand, Chung et al. (2006) employed a mathematical description of the depletion model proposed by Neu and Meiselman (2002). It was noticed that such a description results in a constant (instead of a decaying) attractive force at large separation, which is not physically realistic.

In this study, we follow the approach recently proposed by Liu et al. (2004) to model the intercellular interaction energy *ϕ* as a Morse potential:

(5)

where *r* is the surface separation, *r*_{0} and *D _{e}* are, respectively, the zero force separation and surface energy, and

Before applying the methods described above to RBC flows, several validation simuations have been carried out, and excellent agreement with known theoretical results was observed. All simulations in this paper are performed with a D2Q9 lattice structure (Zhang and Kwok, 2006). The biconcave shape of an RBC is described by the ollowing equation suggested by Evans and Fung (1972):

(6)

where *c*_{0}=0.207, *c*_{1}=2.002, and *c*_{2}=1.122. The non-dimensional coordinates (, ) are normalized by the radius of a human RBC (3.91 *μ*m). The cell membrane is represented by 100 segments and the simulation domain is a 256×128 rectangle in lattice units (50 *μ*m×25 *μ*m, 0.196 *μ*m per lattice unit). For these simulations, we applied periodic boundary condition in the horizontal (*x*) direction and modified bounce-back boundary conditions (Nguyen and Ladd, 2005) on the top and bottom to generate a shear flow field. Other simulation parameters employed in the present work are also listed in Table 1 with references. The intercellular interactions are represented as a Morse potential and evaluated through the Derjaguin approximation (Israelachvili, 1991) for cells with arbitrary shapes. For the parameters in the Morse potential Eq.(5) for the intercellular interaction, experimental values are not available. According to Neu and Meiselman (2002), the zero force distance *r*_{0} is ~14 nm. In our simulations, *r*_{0} is chosen as 2.5 lattice units (0.49 *μ*m) due to the limitation of the grid size and the immersed boundary treatment employed in the computational scheme. This distance is relatively large compared to experimental data (Chien and Jan, 1973b) and the scheme could be improved by adopting a smaller lattice grid and/or the adaptive mesh refinement technology (Toolke et al., 2006). To improve the computational efficiency, as is typical in molecular dynamics simulations, a cut-off distance is adopted. The scaling factor *β* is adjusted so that the attractive force beyond the cut-off distance *r _{c}*=3.52

It has long been recognized that RBCs can aggregate to form stacks-of-coins-like rouleaux in the presence of fibrinogen or other macromolecules (Chien and Jan, 1973b; Skalak and Chien, 1987). In this section, we investigate four RBCs with different interaction strengths *D _{e}* under a no-flow condition. Similar to those employed in previous work (Bagchi et al., 2005; Liu et al., 2004), the cells are initially placed in the domain with a center-to-center distance of 3.52

Rouleaux formed with different intercellular interaction strengths: *D*_{e} = (a) 5.2 × 10^{−8}, (b) 2.1 × 10^{−7}, (c) 2.6 × 10^{−7}, (d) 3.9 × 10^{−7}, (e) 6.5 × 10^{−7}, and (f) 1.3 × **...**

Representative snapshots for the simulation of four biconcave cells in a shear flow *γ* = 20 s^{−1} without intercellular interaction (*D*_{e}=0).

The complex RBC behaviors in flow fields result from the highly coupled interplay of the fluid viscous stresses due to the flow, the membrane forces due to cell deformation, and intercellular interactions. In this section, we will quantitatively investigate these relationships by simulating a system of four biconcave cells in shear flow, which is similar to that observed in hemorheological experiments (Stoltz et al., 1999). The simulation strategy is similar to those employed in previous work (Bagchi et al., 2005; Liu et al., 2004). Rouleaux, formed under intercellular interactions in no-flow condition, as described above, are placed at the center of the simulation domain. A shear flow is applied by specifying the top and bottom boundary velocities, which have the same value but in opposite horizontal directions, in the modified bounce-back boundary conditions (Nguyen and Ladd, 2005). Different shear rates can be obtained by adjusting the boundary velocity.

Experiments show that RBC aggregates can be broken into smaller parts or even individual cells at high shear rate (Popel and Johnson, 2005; Stoltz et al., 1999). Figures Figures22--44 display representative snapshots of the RBC deformation and movement in a shear flow with shear rate *γ*=20 s^{−1}. The slight asymmetry in streamlines is due to the plotting software. The strength of intercellular interaction in the Morse potential Eq.(5) are, respectively, *D _{e}* = 0, 5.2×10

Representative snapshots for the simulation of four biconcave cells in a shear flow *γ* = 20 s^{−1} with a strong intercellular interaction (*D*_{e} = 1.3 × 10^{−6} *μ*J/*μ*m^{2}.

When a weak intercellular interaction is introduced (*D _{e}* = 5.2×10

Representative snapshots for the simulation of four biconcave cells in a shear flow *γ* = 20 s^{−1} with a weak intercellular interaction (*D*_{e} = 5.2 × 10^{−8} *μ*J/*μ*m^{2}.

To examine the hydrodynamic effect on cell behavior, we increase the shear rate to *γ* = 50 s^{−1} for the above cases with weak and strong interaction strengths *D _{e}* = 5.2 × 10

Representative snapshots for the simulation of four biconcave cells in a shear flow *γ* = 50 s^{−1} with a weak intercellular interaction (*D*_{e} = 5.2 × 10^{−8} *μ*J/*μ*m^{2}.

Representative snapshots for the simulation of four biconcave cells in a shear flow *γ* = 50 s^{−1} with a strong intercellular interaction (*D*_{e} = 1.3 × 10^{−6} *μ*J/*μ*m^{2}.

The above simulations show that the dynamic RBC behaviors in shear flows are closely related to the applied shear rate and the intercellular interaction. A larger shear rate produces larger hydrodynamic forces and the RBC rouleaux are less stable and easier split into smaller aggregates or individual cells. On the other hand, a stronger interaction can generate more compact rouleaux and tends to hold the cells together in shear flows. As a result, rouleaux with a stronger interaction can survive a larger shear rate without being broken. These findings are in qualitative agreement with existing experimental observations and therefore demonstrate that our model could be useful in further microscopic hemodynamic and hemorheologic studies. As found by Bagchi et al. (2005), the membrane viscosity appears have no major influences on the cell behaviors in the situations studied here. Effects of other RBC properties, such as deformability and suspending viscosity, will be examined in a future study. We have also observed various cell shapes in these simulations, some of which are similar to previous three-dimensional simulations (Eggleton and Popel, 1998; Pozrikidis, 2001). Goldsmith and Marlow (1979) had also presented RBC deformations in concentrated ghost cell suspensions; however, a direct comparison is difficult due to the different conditions. Liu et al. (2004) had studied a very similar system of four RBCs in shear flow. Our results, in both the flow field and cell movement, are similar to those of Liu et al. (2004) when there is no intercellular interaction involving aggregation forces. However, the aggregation process was not investigated and the cell deformation was less evident in Liu et al. (2004). The small cell deformation could be due to the finite element modeling of the RBC membrane as a thin shell, in which the experimentally measured RBC membrane properties cannot be utilized directly.

We have presented a numerical algorithm to investigate microscopic hemodynamic and hemorheological behaviors of discrete RBCs in shear flow. The immersed boundary treatment has been combined with the lattice Boltzmann method to study deformable fluid capsules in flows. This algorithm has been applied to study the red cell behaviors in shear flow. Crucial factors, including the RBC membrane mechanics, plasma/cytoplasm viscosity difference, intercellular interaction, and the hydrodynamic viscous forces, have been considered in this numerical model. Simulations of four RBCs in shear flows have been performed, and the effects of the intercellular interaction and hydrodynamic forces on the cell aggregating and disaggregating processes have been investigated. Without the intercellular interaction, the cells remain separated and no aggregates are formed. By introducing the interaction, cells are attracted to each other and rouleau-like aggregates can be generated. The aggregate configuration is closely related to the interaction strength. Such an aggregate could rotate as a unit with certain deformation, or be separated under a shear flow, depending on the intercellular interaction and shear strength. These findings are qualitatively consistent with experimental observations, and therefore demonstrate that the proposed algorithm could be useful in microscopic blood flow studies. Extending this algorithm to three-dimensional situations is straightforward, as well as parallelizing it for better computational efficiency. Different intercellular interaction models can also be readily incorporated in the numerical scheme.

It should be mentioned that the Morse potential adopted in this study is only a qualitative description of the still unclear microscopic intercellular interactions. In principle, parameters involved in the Morse potential Eq.(5) could be obtained through comparing simulation results to observable macroscopic behaviors, such as rouleaux configurations, aggregation and dissociation processes, and apparent viscosities at different shear rates. For the sake of simulation efficiency, the zero force distance *r*_{0} in this study was chosen to be relatively large compared to experimental values. Similar assumptions and simplifications were made in previous studies [for example, see Fig. 6 in Liu et al. (2004)]. In addition, since the fluid trapped in a gap between two close RBC surfaces is immobilized, from a viewpoint of hydrodynamics, a difference in the gap sizes might not have significant impact on the global behaviors. Nevertheless, with the present knowledge of the RBC intercellular interactions and with presently used computational technology, the Morse potential could be a reasonable choice to model the RBC aggregation.

(a) Normal and less deformable RBCs flowing in a straight channel and (b) the velocity profiles compared with a pure plasma flow.

Click here to view.^{(537K, pdf)}

Click here to view.^{(319K, avi)}

Click here to view.^{(770K, avi)}

Click here to view.^{(473K, avi)}

Click here to view.^{(318K, avi)}

Click here to view.^{(483K, avi)}

This work was supported by NIH grant HL52684. J. Zhang acknowledges the financial support from the Natural Science and Engineering Research Council of Canada (NSERC) through a postdoctoral fellowship at the Johns Hopkins University. P. C. Johnson is also a Senior Scientist at La Jolla Bioengineering Institute.

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Simulation Movies Available as Supplemental Materials.

- Armstrong JK, Meiselman HJ, Fisher TC. Evidence against macromolecular bridging as the mechanism of red blood cell aggregation induced by nonionic polymers. Biorheology. 1999;36:433–437. [PubMed]
- Bagchi P. Mesoscale simulation of blood flow in small vessels. Biophysical Journal. 2007;92:1858–1877. [PubMed]
- Bagchi P, Johnson PC, Popel AS. Computational fluid dynamic simulation of aggregation of deformable cells in a shear flow. Journal of Biomechanical Engineering. 2005;127:1070–1080. [PubMed]
- Baumler H, Donath E. Does dextran really significantly increase the surface potential of human red blood cells? Studia Biophysica. 1987;120:113–122.
- Baumler H, Neu B, Donath E, Kiesewetter H. Basic phenomena of red blood cell rouleaux formation biorheology. Biorheology. 1999;36:439–442. [PubMed]
- Brooks D. Electrostatic effects in dextran mediated cellular interactions. Journal of Colloid and Interface Science. 1973;43:714–726.
- Chien S, Jan K. Red cell aggregation by macromolecules: Roles of surface adsorption and electrostatic repulsion. Journal of supramolecular structure. 1973a;1:385–409. [PubMed]
- Chien S, Jan K. Ultrastructural basis of the mechanism of rouleaux formation. Microvascular Research. 1973b;5:155–166. [PubMed]
- Chien SS, Simchon S, Abbot RE, Jan KM. Surface adsorption of dextrans on human red cell membrane. Journal of Colloid and Interface Science. 1977;62:461–470.
- Chung B, Johnson PC, Popel AS. Application of chimera grid to modeling cell motion and aggregation in a narrow tube. International Journal of Numerical Methods in Fluids. 2006;53(1):105–128.
- Eggleton CD, Popel AS. Large deformation of red blood cell ghosts in a simple shear flow. Physics of Fluids. 1998;10(8):1834–1845.
- Evans E, Needham D. Attraction between lipid bilayer membranes - comparison of mean field theory and measurements of adhesion energy. Macromolecules. 1988;21:1822–1831.
- Evans EA. Bending elastic modulus of red blood cell membrane derived from buckling instability in micropipette aspiration tests. Biophysical Journal. 1983;43:27–30. [PubMed]
- Evans EA, Fung YC. Improved measurements of the erythrocyte geometry. Microvascular Research. 1972;4:335–347. [PubMed]
- Evans EA, Hochmuth RM. Membrane viscoelasity. Biophysical Journal. 1976;16:1–11. [PubMed]
- Fischer T, Schmid-Schonbein H. Tank treading motion of red blood cell membranes in viscometric flow: Behavior of intracellular and extracellular markers (with film) Blood Cells. 1977;3:351–365.
- Goldsmith HL, Marlow JC. Flow behaviors of erythrocytes ii: particle motions in concentrated suspensions of ghost cells. Journal of Colloid and Interface Science. 1979;71(2):383–407.
- Hochmuth RM, Waugh RE. Erythrocyte membrane elasticity and viscosity. Annual Review of Psychology. 1987;49:209–219. [PubMed]
- Israelachvili J. Intermolecular and Surface Forces. 2nd. edition Academic Press; San Diego, USA: 1991.
- Kim S, Popel AS, Intaglietta M, Johnson PC. Effect of erythrocyte aggregation at normal human levels on functional capillary density in rat spinotrapezius muscle. American Journal of Physiology - Heart and Circulatory Physiology. 2006;290:H941–H947. [PubMed]
- Kounov NB, Petrov VG. Determination of erythrocyte aggregation. Mathematical Biosciences. 1999;157:345–356. [PubMed]
- Liu Y, Zhang L, Wang X, Liu WK. Coupling of Navier-Stokes equations with protein molecular dynamics and its application to hemodynamics. International Journal of Numerical Methods in Fluids. 2004;46:1237–1252.
- Mchedlishvili G, Maeda N. Blood flow structure related to red cell flow: a determinant of blood fluidity in narrow microvessels. Japanese Journal of Physiology. 2001;51(1):19–30. [PubMed]
- Merill E, Gilliland E, Lee T, Salzman E. Blood rheology: Effect of fibrinogen deduced by addition. Circulation Research. 1966;18:437–446. [PubMed]
- Migliorini C, Qian Y, Chen H, Brown EB, Jain RK, Munn LL. Red blood cells augment leukocyte rolling in a virtual blood vessel. Biophysical Journal. 2002;83(4):1834–1841. [PubMed]
- Mills JP, Qie L, Dao M, Lim CT, Suresh S. Nonlinear elastic and viscoelastic deformation of the human red blood cell with optical tweezers. Mechanics and Chemistry of Biosystems. 2004;1:169–180. [PubMed]
- Neu B, Meiselman HJ. Depletion-mediated red blood cell aggregation in polymer solutions. Biophysical Journal. 2002;83:2482–2490. [PubMed]
- Nguyen N-Q, Ladd AJC. Sedimentation of hard-sphere suspensions at low Reynolds number. Journal of Fluid Mechanics. 2005;525:73–104.
- Peskin C. Numerical analysis of blood flow in the heart. Journal of Computational Physics. 1977;25:220–252.
- Popel AS, Johnson PC. Microcirculation and hemorheology. Annual Review of Fluid Mechanics. 2005;37:43–69. [PMC free article] [PubMed]
- Pozrikidis C. Finite deformation of liquid capsules enclosed by elastic membranes in simple shear flow. Journal of Fluid Mechanics. 1995;297:123–152.
- Pozrikidis C. Effect of membrane bending stiffness on the deformation of capsules in simple shear flow. Journal of Fluid Mechanics. 2001;440:269–291.
- Pozrikidis C. Numerical simulation of the flow-induced deformation of red blood cells. Annals of Biomedical Engineering. 2003;31:1194–1205. [PubMed]
- Skalak R, Chien S. Handbook of Bioengineering. McGraw-Hill; New York, NY: 1987.
- Stoltz J, Singh M, Riha P. Hemorheology in Practice. IOS Press; Amsterdam, Netherlands: 1999.
- Succi S. The Lattice Boltzmann Equation. Oxford Univ. Press; Oxford: 2001.
- Sun C, Migliorini C, Munn LL. Red blood cells initiate leukocyte rolling in postcapillary expansions: A lattice Boltzmann analysis. Biophysical Journal. 2003;85(1):208–222. [PubMed]
- Sun C, Munn LL. Particulate nature of blood determines macroscopic rheology: A 2-D lattice Boltzmann analysis. Biophysical Journal. 2005;88(3):1635–1645. [PubMed]
- Sun C, Munn LL. Influence of erythrocyte aggregation on leukocyte margination in postcapillary expansions: A lattice Boltzmann analysis. Physica A. 2006;362:191–196.
- Toolke J, Freudiger S, Krafczyk M. An adaptive scheme using hierarchical grids for lattice boltzmann multi-phase flow simulations. Computers and Fluids. 2006;35(89):820–830.
- Waugh RE, Hochmuth RM. Chapter 60: Mechanics and deformability of hematocytes. In: Bronzino JD, editor. Biomedical Engineering Fundamentals. 3rd edition CRC; Boca Raton, FL: 2006. pp. 60–3.
- Zhang J, Kwok DY. A 2D lattice Boltzmann study on electrohydrodynamic drop deformation with the leaky dielectric theory. Journal of Computational Physics. 2005;206(1):150–161.
- Zhang J, Kwok DY. Contact line and contact angle dynamics in superhydrophobic channels. Langmuir. 2006;22:4998–5004. [PubMed]
- Zhang J, Li B, Kwok DY. Mean-field free-energy approach to the lattice Boltzmann method for liquid-vapor and solid-fluid interfaces. Physical Review E. 2004;69:032602. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |