Next Article in Journal
Effect of NaCl on the Lifetime of Micro- and Nanobubbles
Next Article in Special Issue
Nanomaterials for Biosensing Applications
Previous Article in Journal
Simultaneous Reduction and Functionalization of Graphene Oxide by 4-Hydrazinobenzenesulfonic Acid for Polymer Nanocomposites
Previous Article in Special Issue
A Graphene Oxide-Based Fluorescent Platform for Probing of Phosphatase Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method

1
Department of Mechanical Engineering and Mechanics, Lehigh University, Bethlehem, PA 18015, USA
2
Department of Civil and Environmental Engineering, Lehigh University, Bethlehem, PA 18015, USA
3
School of Mechanics and Engineering, Southwest Jiaotong University, Chengdu 610031, China
4
Bioengineering Program, Lehigh University, Bethlehem, PA 18015, USA
*
Authors to whom correspondence should be addressed.
Nanomaterials 2016, 6(2), 30; https://doi.org/10.3390/nano6020030
Submission received: 16 November 2015 / Revised: 21 January 2016 / Accepted: 25 January 2016 / Published: 5 February 2016
(This article belongs to the Special Issue Nanomaterials for Biosensing Applications)

Abstract

:
Nanodrug-carrier delivery in the blood stream is strongly influenced by nanoparticle (NP) dispersion. This paper presents a numerical study on NP transport and dispersion in red blood cell (RBC) suspensions under shear and channel flow conditions, utilizing an immersed boundary fluid-structure interaction model with a lattice Boltzmann fluid solver, an elastic cell membrane model and a particle motion model driven by both hydrodynamic loading and Brownian dynamics. The model can capture the multiphase features of the blood flow. Simulations were performed to obtain an empirical formula to predict NP dispersion rate for a range of shear rates and cell concentrations. NP dispersion rate predictions from the formula were then compared to observations from previous experimental and numerical studies. The proposed formula is shown to accurately predict the NP dispersion rate. The simulation results also confirm previous findings that the NP dispersion rate is strongly influenced by local disturbances in the flow due to RBC motion and deformation. The proposed formula provides an efficient method for estimating the NP dispersion rate in modeling NP transport in large-scale vascular networks without explicit RBC and NP models.

1. Introduction

Accurately predicting drug delivery is a critical task in drug development research and clinical trials [1,2]. It requires careful consideration of physiological conditions, such as hematocrit level (the volume ratio between red blood cell and the whole blood) [3,4], vessel geometry and flow conditions [5,6,7], drug carrier size and shape [3,8], dissolution rate [9] and external stimuli [10,11]. For small particles in red blood cell (RBC) suspensions, such as nanoparticles (NP) and platelets, recent studies have demonstrated that local flow field disturbances caused by RBC translation and deformation can enhance particle dispersion [3,12,13,14,15,16]. The migration of particles in RBC suspensions under shear has been shown to behave like a random walk process [17,18], with a dispersion rate much larger than thermal diffusion. Therefore, accurate predictions of NP dispersion in RBC suspensions must consider fluid-structure interaction between the immersed solid bodies (particles and cells) and the surrounding fluid. Previously-developed models for predicting NP dispersion in RBC suspensions have relied primarily on empirical data fitting. Aarts et al. experimentally studied shear-induced platelet diffusivity (D), which was fitted with shear rate (η) as a power law D = k η n , where k is a constant and n is a function of hematocrit [19]. However, the model parameters are obtained empirically rather than predicted from the underlying physics. Decuzzi et al. extended the Taylor–Aris theory to calculate an effective NP diffusion rate that considers wall permeability and blood rheology [20,21]. They also reported about a three-fold increase in dispersion rate of 1 µm compared to thermal diffusion [16]. Recently, Fedosov’s group systematically studied micro- and nano-particles in drug delivery, including particle size, shape effect and RBC influence on particle margination and adhesion probabilities [3]. However, there is no analytical formula or quantitative rule to directly predict the NP dispersion rate so far, which is much needed in large-scale drug delivery simulations [20,21,22].
In order to address the deficiencies in previously-developed models for predicting NP dispersion, this paper presents a numerical study on NP dispersion in RBC suspensions that considers the effects of local flow field disturbances due to RBC motion. This study provides insight into the underlying physics driving NP dispersion in these systems and develops simple, yet effective, formulae for predicting dispersion rate as a function of characteristic physiological parameters. These simple predictive formulae will provide an efficient approach for assessing NP dispersion under different flow conditions and hematocrit level, thereby facilitating practical modeling of NP transport and distribution in large-scale vascular systems [22].
The remainder of this paper is organized as follows. The fluid-structure interaction model is introduced in Section 2, including a description of the modeling approach for the fluid environment, RBCs and NPs, as well as the fluid-structure coupling scheme. Section 3 outlines the model setup and test parameters for a parametric study on NP dispersion. Simulation results from the parametric study are presented in Section 4, along with formulae derived from a regression analysis for predicting the NP dispersion rate. Predictions from the formulae are compared to data reported in the literature. Finally, conclusions and recommendations from the study are presented in Section 5.

2. Fluid-Structure Interaction Model

The transport of particles in RBC suspensions is governed by hydrodynamic forces and fluid-structure interaction effects. In order to simulate this two-way coupling between the fluid environment and the immersed soft matter, a numerical fluid-structure interaction (FSI) code was developed utilizing an immersed boundary coupling scheme. Previous research has shown that the immersed boundary method (IBM) is an efficient approach to simulate soft matter and biological cells [23,24,25,26,27,28,29,30], flapping insect wings [31,32,33], harmonic oscillation of thin lamina in fluid [34] and other FSI problems, such as particle settling [35]. The advantage of this approach is two-fold. First, the lattice Boltzmann method (LBM) is very efficient at modeling fluid flow and ideal for parallel computing [36]. Second, the IBM uses an independent solid mesh moving on top of a fixed fluid mesh, thus removing the burden of re-meshing in conventional arbitrary Lagrangian–Eulerian methods [37]. In this study, the fluid was modeled using a lattice Boltzmann scheme, while the RBCs were modeled as elastic membranes. The fluids inside and outside the cell membrane are modeled with the same viscosity, following that in other works [38,39]. This approach greatly simplified the computation and without losing too much generality. It also can capture the multiphase feature of the blood flow. NPs were treated as rigid point elements with motions governed by both hydrodynamic loading and Brownian dynamics. Additional details regarding the fluid-structure interaction model are provided in the following sections.

2.1. Lattice Boltzmann Fluid Model

The lattice Boltzmann method (LBM) is a mesoscale approach to modeling fluid dynamics that has been used extensively in blood flow modeling [24,38,40]. Reviews of the underlying theory for the LBM can be found in the literature [41,42,43,44]. LBM is usually considered as a second order accurate method in space and time [45]. The key concept of the LBM is the transport of particles, characterized by a local density distribution function fi (x,t) in phase space (x, c i ), where t denotes the time and c i denotes the lattice velocity. The dynamics of the LBM involves streaming and collision processes.
f i ( x + Δ t c i , t + Δ t ) f i ( x , t ) Streaming =    Ω    i Collision + F i
where Ω i is the collision operator and F i is the external force term. The simplest and most computationally-efficient collision scheme is the Bhatnagar–Gross–Krook (BGK) scheme [46],
  Ω i = 1 τ ( f i ( x , t ) f i e q ( x , t ) )
where τ is a relaxation constant and f i e q ( x , t ) is the population distribution at equilibrium, which is related to the local macroscale fluid velocity and the speed of sound c s .
f i e q ( u , ρ ) = w i ρ ( 1 + c i u c s 2 + ( c i u ) 2 2 c s 4 u u 2 c s 2 )
where w i is the weight coefficients. An effective fluid viscosity is related to the relaxation time τ and the speed of sound cs,
  ν = c s 2 ( τ 0.5 ) = ( τ 0.5 ) 3
and local macroscale density and velocity are calculated as:
ρ ( x , t ) = i f i ( x , t ) , ρ ( x , t ) u ( x , t ) = i f i ( x , t ) c i + Δ t 2 f
Disturbances in the flow are introduced through a force density term Fi, which can be expressed in terms of an external body force density f and fluid macroscale velocity u ,
F i = ( 1 1 2 τ ) w i ( c i u c s 2 + c i u c s 4 c i ) f
For simplicity, the simulations presented in this paper consider a regular 2D domain with a D2Q9 (2D 9 velocities) fluid lattice, the details of which can be found in [45]. However, it is noted that the modeling approach is readily adaptable to the 3D case (see Figure S9 in the Supplementary Materials) and with more complex geometric configurations. Velocity and non-slip wall boundaries were used. There are several ways to implement velocity boundary conditions, such as Zou/He boundaries [47] and regularized boundary conditions [48]. The basic idea behind Zou/He boundaries is assuming that the bounce back rule applies to the non-equilibrium part of the density population. The regularized boundary condition is to use the moment flux tensor to reconstruct the non-equilibrium part of the density distribution, thus it is more stable, but requires extra calculation. Interested readers are referred to [49] for a summary of different boundary conditions in LBM. In this paper, treatment of velocity boundary conditions followed the recommendations in [47] and is discussed in Section 3.

2.2. Spring Connected Network Cell Membrane Model

Cell models are commonly based on a continuum approach utilizing a strain energy function to define membrane response [50,51,52]. Recently, however, a particle-based cell model governed by molecular dynamics has emerged as a popular alternative due to its simplicity in mathematical description [53,54,55]. These two models have been shown to provide consistent predictions [53]. For the present study, a particle-based RBC model is adopted. The model consists of a network of vertices, Xi, i { 1... N V } , connected with elastic springs, which are effective in axial extension/contraction and flexure, as shown in Figure 1.
Figure 1. (A) Spring connected network cell membrane model. Kinematics for local stretching (B) and bending (C) modes of response.
Figure 1. (A) Spring connected network cell membrane model. Kinematics for local stretching (B) and bending (C) modes of response.
Nanomaterials 06 00030 g001
The potential energy of the membrane is defined as:
V { X i } = V s t r e t c h + V b e n d i n g + V a r e a
The stretch energy mimics the elastic spectrin network and is given by:
V s t r e t c h = 1 2 k s j 1... N s ( L j L j 0 ) 2
where ks is the stretching constant, L j is the length of the spring and L j 0 is the equilibrium spring length. To consider the nonlinear effect of membrane stretching, Fedosov investigated a wormlike chain potential [56]. In a separate study, Wida developed an exponential relationship for mechanical stiffness related to the bond stretch ratio (λ) [57]. In this paper, the same exponential form of the spring constant, k s = k s 0 exp [ 2 ( λ 1 ) ] , where λ is the bond stretch ratio, has been adopted.
The bending energy represents the bending resistance of the lipid bilayer and is defined as:
V b e n d i n g = 1 2 k b j 1... N b ( θ j θ j 0 ) 2
where k b is the bending constant and θ j and θj0 are instantaneous and equilibrium angles between adjacent spring, respectively. In this model, we used the angles when cells are at rest as the equilibrium angles θj0; each angle is also different from the others.
When the cell deforms, the cell surface area and volume remain relatively constant. To maintain a constant area in the 2D model, an artificial restoring force potential V a r e a = k a ( A A 0 ) 2 2 A 0 is introduced to conserve the area where A is the instantaneous area, A 0 is the equilibrium area and k a is a constant.
The equilibrium shape of cells is followed by the analytical curve from Fung [58], as shown in Figure 1. The cell membrane was discretized to 52 nodes and angles. The equilibrium bond length was set to be equal, but the equilibrium angle varies among all of the angles.

2.3. Immersed Boundary Coupling Scheme

The immersed boundary method (IBM) was selected to model the interaction between the fluid and the immersed solids due to the algorithm’s efficiency. The IBM was first proposed by Peskin to study blood flow in the heart [59]. The fluid is solved on a spatially-fixed Eulerian grid, while the immersed solids are modeled using a moving Lagrangian mesh, which is not constrained to the geometric layout of the Eulerian fluid grid. Data are exchanged between the two domains through nodal interpolation. The coupling scheme enforces velocity continuity at the fluid-structure boundary and transfers forces from the structure back into the fluid through an effective force density. This two-way coupling automatically handles immersed body contact and prevents solid penetration through the development of restoring forces in the fluid. The approach has been used for a variety of fluid-structure interaction problems, including the simulation of jelly fish [60], blood flow [4,24,26,27,28,29,30,38] and platelet migration [14]. Comprehensive reviews of the IBM and its applications can be found in [25,61]. It is noticed that the current coupling approach works great for soft deformable solids, an alternative coupling approach, e.g., treating a solid as a moving boundary for the fluid and transferring the force back to the solid, should be used if the solid is relatively rigid, such as a particulate settling study [44,62,63,64].
The immersed structure can be viewed as a parametric surface X(s,t), where the force exerted by the structure on the fluid is interpolated as a source term in the momentum equation using:
f ( x , t ) = F ( s , t ) δ ( x X ( s , t ) ) d s
where F ( s , t ) is the structural force at location s at time of t and ds is the discretized length of the immersed structure. The force is spread to the local fluid nodes through a delta function δ ( r ) ,
δ ( r ) = { 1 4 ( 1 + cos ( π r 2 ) ) 2 r 2 0                      otherwise
Other delta functions can also be used. Peskin listed a few delta functions in [25].
The structure velocity is interpreted from the local fluid field through:
u ( X , t ) = u ( x , t ) δ ( x X ( s , t ) ) d x
The same delta function Equation (11) is used to interpolate cell node velocities from the local fluid flow field.

2.4. Nanoparticle Model

NPs were modeled as rigid bodies with motion governed by hydrodynamic forces and Brownian dynamics [65,66]. Langevin dynamics was used to simulate the motion of particles.
m d u d t = ς u + F c + F r
where ς is the friction coefficient, F c is the conservative force and F r is the random force that satisfies the fluctuation dissipation theorem.
F r ( t ) = 0
F r ( t ) F r ( t ' ) = 2 k B T ς δ ( t t ' ) I
where kBT is the thermal energy, δ(t − t′) is the Dirac delta function and I is the unit-second order tensor.
The solution of the above equation gives u ( t ) = F c + F r ς + C exp ( ς t m ) , where C is a constant. When the time step dt (4.2 × 10−8 s in the present study) is much larger than the relaxation time τ R = m ς (2.2 × 10−9 s), the particle position can be updated with the terminal velocity as u ( t ) = F c + F r ς + u f . A position check was performed to avoid the penetration of particles into the RBCs, e.g., the random displacement induced by thermal fluctuation will be reversed if it will jump inside cells at the next time step. In the diffusion coefficient calculation, periodic boundary conditions were applied to particles when they crossed the channel wall to remove the effect of geometric confinement.
The flow chart of the algorithm for the program is shown in Figure 2, which is a second order accurate Runge–Kutta method based primarily on the midpoint rule [25]. The fluid-structure interaction model was validated for a sphere settling process [67] and a lateral migration process [68] in a viscous fluid (Figures S1–S5). Validation studies for fluid-structure interaction, cell tumbling and tank treading motions (Figure S6, parameters listed in Table S1), as well as NP diffusion (Figures S7–S8) are shown in the Supplementary Materials, which is provided in the supplementary document.
Figure 2. The algorithm flow chart of immersed boundary method.
Figure 2. The algorithm flow chart of immersed boundary method.
Nanomaterials 06 00030 g002

3. Model Setup and Parametric Study

Research has shown that particles in the core region of the vessel migrate toward the cell-free layer regions, where the migration process can be modeled as diffusion [17,18]. This migration is influenced by physical conditions, such as hematocrit level (Ht), cell membrane stiffness (ks), particle size (r), shear rate ( η ), fluid viscosity (ν) and cell size (dc). In this study, two parameters (hematocrit level and shear rate) are considered, while the other parameters are kept constant [13]. NP dispersion is first studied under pure shear flow conditions at different shear rates for a given hematocrit level. This pure shear flow setup was used to explore the NP dispersion rate when cells undergo tumbling and tank treading motion [69,70,71]. Then, the study is extended to investigate NP dispersion in channel flow at different hematocrit levels, which was used to model the blood flow in capillaries with a physiologically-relevant shear rate between 0 and 500 s−1 [72].
In the study, RBCs were assumed to be healthy with typical physical parameters and with a size of 8 μm. NPs were assumed to be spherical with a typical size of 100 nm. For simplicity, surface charges were neglected, so that NPs did not adhere to other NPs or to the RBCs. It is noted that the NP concentration was kept relatively low in order to more readily ascertain the effect of RBC motion on NP dispersion. Dynamic viscosity of the fluid was fixed at 1 × 10−3 Pa·s. Through dimensional analysis, an empirical function between diffusion coefficient, shear rate, cell size and hematocrit was defined as: D d c 2 η = f ( H t ) , where D is the dispersion rate, dc is the cell size, η is the shear rate and H t is hematocrit. The same dimensionless dispersion rate was also suggested in [15]. This formula was validated through simulations presented in later sections.
The test case consisted of a rectangular fluid domain with a 50-µm length, a 25-µm width and a lattice grid size of 0.5 µm, as shown in Figure 3A. Simulation with a finer grid of 0.25 µm showed a similar fluid velocity field. In the shear flow case, the top and bottom surfaces were defined as velocity boundaries, while the left and right edges of the domain were modeled as periodic boundaries. In the channel flow case, a parabolic velocity profile was applied at the left inlet boundary, and the right outlet boundary was modeled as an open condition. Non-slip boundaries were defined along the upper and lower surfaces. Notice that the physical system has to be converted into LB units through reference dimensions, such as length (dx), time (dt) and mass (dm). Here, we select dx = 0.5 µm. The reference mass dm = ρdx3 = 1.28 × 10−16 kg assuming the blood plasma density is 1025 kg/m3. Depending on the Reynolds number, a relaxation time τ is typically recommended between 0.8 and 1.0; see [73,74]. As suggested in [44], a relaxation time τ of 1.0 was used for all simulations for computational efficiency and larger time steps. The time step dt can be determined through normalizing the fluid viscosity τ 0.5 3 = ν d t d x 2 , where ν is the kinematic viscosity of the fluid. Thus, the time step can be calculated as 4.2 × 10−8 s. Other parameters can be normalized, as well, since the basic reference values for length, time and mass have been selected. As all of the simulations performed here have low Reynolds numbers, the compressibility effect of the Mach number can be safely ignored.
RBC membranes were modeled as bi-concave curves with the dimensions shown in Figure 1A. A single RBC was composed of 52 nodes. The cell parameters were selected based on recommended values reported in the literature [53,75], as listed in Table 1. The artificial area constraint ka = 1 was selected so that the area change was within 1%. Periodic boundary conditions (PBC) were applied to the left and right boundaries of the fluid domain for both RBCs and NPs. The area ratio between RBCs and the fluid domain was defined as the hematocrit level.
Table 1. Cell membrane model parameters.
Table 1. Cell membrane model parameters.
ParametersSpecified ValueRecommended RangeReferences
Stretching coefficient ks05 µN/m5~12 µN/m[53,75]
Bending coefficient kb8 × 10−19 J2 × 10−19 to 1 × 10−17 J[53,75]

4. Results and Discussion

4.1. NP Dispersion under Pure Shear Flow

The NP dispersion rate was studied over a range of shear rates for a single layer of three cells. This setup was designed to eliminate the cell-cell interaction between different layers, so that we can focus on the shear rate effect on NP dispersions. Shear rates ranging from 20 to 500 s−1 were selected in order to cover both the RBC tumbling and RBC tank treading regions of the flow regime. For all shear rates investigated in the study, the dimensionless number ηt was held at 25. Therefore, the simulation time is longer for lower shear rate case. Three RBCs and 792 NPs were considered for each simulation. Snapshots of the interaction between NPs and RBCs at shear rates of 40 s−1 and 200 s−1, representative of RBC tumbling and RBC tank treading regions of the flow regime, respectively, are shown in Figure 3 (see the corresponding movies in the Supplementary Materials).
Figure 3. Interaction between nanoparticle (NP) and red blood cell (RBC) at shear rates of (A) 40 s−1 (RBC tumbling) and (B) 200 s−1 (RBC tank treading). The bold red lines outline the RBC membranes, while the green markers denote NPs. Flow streamlines are shown in the background.
Figure 3. Interaction between nanoparticle (NP) and red blood cell (RBC) at shear rates of (A) 40 s−1 (RBC tumbling) and (B) 200 s−1 (RBC tank treading). The bold red lines outline the RBC membranes, while the green markers denote NPs. Flow streamlines are shown in the background.
Nanomaterials 06 00030 g003
The mean square displacement over the y direction at different shear rates was calculated to obtain the NP dispersion rates, as shown in Figure 4. It shows that the dispersion rate is strongly influenced by cell motion. In the RBC tumbling ( η < 40 s−1) and RBC tank treading ( η > 200 s−1) regions of the flow regime, the NP dispersion rate is approximately linear with the shear rate. Between 40 s−1 and 200 s−1, there is a region where RBC motion transits from tumbling to tank treading motion. In this transition region, there is a drop in NP dispersion with increased shear rate. This shows that the cell tumbles first and then it deforms to a shape that cannot sustain tumbling motion, then evolves into the tank treading shape. Thus, the overall average dispersion rate is decreased compared to the pure tumbling regime, which leads to a peak around a shear rate of 40 s−1. For the range of shear rates investigated in the study, the dispersion rate initially increases in the tumbling region, then decreases in the transition region and increases again with the shear rate in the tank treading region. A linear regression model was used to fit both the tumbling (first three data points at a low shear rate) and tank treading data (last three data points at a high shear rate),
{ 7.8 × 10 14 η + 4.7 × 10 12 tumbling       8.5 × 10 15 η + 4.0 × 10 12 tank treading
where η is the shear rate. The formulae indicate that the effect of shear rate on NP dispersion in the tumbling region is roughly an order of magnitude larger than that in the tank treading region. This can be attributed to larger RBC motions in the tumbling region, where RBCs undergo rigid-like rotations that trigger larger local flow disturbances that promote the dispersion of adjacent NP away from the cell. It is also worth noting that the constant terms in the formulae are close to the NP thermal diffusion coefficient. The theoretical diffusion rate for 100-nm particles at a temperature of 300 K is about 4.4 × 10 12 m2/s. This observation agrees with the physical requirement that the dispersion rate should be close to thermal diffusion in the absence of shear flow. Therefore, for a given hematocrit level, and with η < 40 s−1 and η > 200 s−1, the dispersion rate D can be written as:
D = k η + D 0
where D 0 is the thermal diffusion coefficient and k is a constant that depends on the hematocrit level. It is noted that this formula is readily adaptable to different particle sizes, because the constant term D0 already contains the particle size effect. The contribution of RBC motion is represented in the constant k . It is noted that the influence of particle concentration on dispersion rate has been neglected.
Figure 4. NP dispersion rate as a function of shear rate. Error bars indicate the standard variance for three simulations. RBCs undergo tumbling motion at a low shear rate (η< 40 s−1) and tank treading motion at a high shear rate ( η > 200 s−1). In between, there is a transition region. Linear regression lines for the tumbling and tank treading regions are shown, as well.
Figure 4. NP dispersion rate as a function of shear rate. Error bars indicate the standard variance for three simulations. RBCs undergo tumbling motion at a low shear rate (η< 40 s−1) and tank treading motion at a high shear rate ( η > 200 s−1). In between, there is a transition region. Linear regression lines for the tumbling and tank treading regions are shown, as well.
Nanomaterials 06 00030 g004

4.2. NP Dispersion under Channel Flow

It should be noted that the pure shear flow case does not exist in physiological flow. While the shear case provides an understanding of NP dispersion enhancement under different shear flow rates, it is important to simulate NP dispersion in realistic channel flow cases. Figure 5A,B presents snapshots of NP dispersion in a channel flow simulation with a hematocrit of 23.5% and a shear rate of 200 s−1 at 0.26 s and 0.46 s, respectively. The channel width is 25 µm. For these simulations, the specified shear rate was measured as the shear rate at the wall, unless noted otherwise. Due to the increased cell volume, compared to the pure shear flow simulations, the number of NPs was reduced to 378. The NPs were initially positioned in the core region of the channel. Since the shear rate is linearly changing across the channel, RBCs did not exhibit distinctive motions, such as tumbling or tank treading, as shown in the previous pure shear flow. The higher hematocrit and cell-cell interaction also confined the cell motion in the flow. The majority of the cells behaved like a tank treading motion, while some cells in the core region were bended or folded in the channel due to the symmetry of the velocity near the center line of the channel. The RBC motion under other hematocrit levels and shear rates was similar to Figure 5. However, they are not shown here.
Figure 5. Snapshots from a channel flow simulation for a cell-particle mixture with a hematocrit level of 23.5% and a shear rate of 200 s−1 at 0.26 s (A) and 0.46 s (B). Fluid streamlines are shown in the background, while the yellow markers represent 100-nm nanoparticles. For illustration purposes, cells and particles are not shown to scale. A short movie of the simulation is provided in the Supplementary Materials.
Figure 5. Snapshots from a channel flow simulation for a cell-particle mixture with a hematocrit level of 23.5% and a shear rate of 200 s−1 at 0.26 s (A) and 0.46 s (B). Fluid streamlines are shown in the background, while the yellow markers represent 100-nm nanoparticles. For illustration purposes, cells and particles are not shown to scale. A short movie of the simulation is provided in the Supplementary Materials.
Nanomaterials 06 00030 g005
As illustrated in Figure 5A,B, the NPs tend to migrate toward the wall. In order to characterize the NP distribution across the channel, the channel width was divided into bins of 1 µm. The number of NPs within each bin was counted and divided by the total number of NPs to obtain the NP fraction within each bin. The NP fraction across the channel height is plotted in Figure 6A at time points of 0, 0.26 and 0.52 s, for a characteristic shear rate of 200 s−1. Figure 6B presents the NP fraction across the channel height at t = 0.52 s under shear rates of 100 s−1, 200 s−1, 300 s−1 and 500 s−1. The NP fractional values shown in the figure are the average of three sample runs using different random seeds for the NP Brownian motion model. Figure 6B clearly demonstrates that particle migration speed toward the channel walls increases with shear rate.
Figure 6. The NP fraction across the channel height for a hematocrit level of 23.5%. (A) NP fraction at t = 0, 0.26 s and 0.52 s for a shear rate of 200 s−1; and (B) NP fraction at t = 0.52 s for shear rates of 100 s−1, 200 s−1, 300 s−1 and 500 s−1.
Figure 6. The NP fraction across the channel height for a hematocrit level of 23.5%. (A) NP fraction at t = 0, 0.26 s and 0.52 s for a shear rate of 200 s−1; and (B) NP fraction at t = 0.52 s for shear rates of 100 s−1, 200 s−1, 300 s−1 and 500 s−1.
Nanomaterials 06 00030 g006
In order to characterize NP migration speed, the dispersion rate was calculated from the NP mean square displacement. The dispersion rates for different hematocrit levels, and at various shear rates, are shown in Figure 7A. From the pure shear simulation results shown in Equation (14), a modified dimensionless dispersion rate was developed,
D r = D D 0 d c 2 η = f ( H t )
The dimensionless dispersion rate D r is plotted in Figure 7B, where the error bars show standard variance for three sample runs.
Figure 7. Dispersion rate of NPs in blood flow. (A) NP dispersion rate at different hematocrit (Ht) and shear rates; (B) relationship between dimensionless dispersion rate (Dr) and hematocrit (Ht). Error bars show the standard variance from three samples.
Figure 7. Dispersion rate of NPs in blood flow. (A) NP dispersion rate at different hematocrit (Ht) and shear rates; (B) relationship between dimensionless dispersion rate (Dr) and hematocrit (Ht). Error bars show the standard variance from three samples.
Nanomaterials 06 00030 g007
Figure 7 shows that the lateral dispersion of NPs (i.e., migration of NPs toward the vessel walls) is much larger than what is predicted by thermal diffusion alone. This migration is influenced by both the hematocrit level and the shear rate. While the relationship between dispersion rate and shear rate is approximately linear (Figure 7A), the relationship between dispersion and hematocrit is not fully linear (Figure 7B). Nevertheless, a best fit line with reasonable approximation can be written as:
D r = D D 0 d c 2 η = 4 . 643 × 10 3 H t + 5 . 834   × 10 4
Particle dispersion rate predictions from Equation (18) were also compared to dispersion rates published in the literature for platelets [15,76] and 1-µm particles [13]. It is noted that the comparison of Equation (18) with platelets and microparticles was considered due to the lack of the NP dispersion rate in the literature. The predictions and measured dispersion rates are summarized in Table 2. The dimensionless dispersion rate is also plotted in Figure 8. As shown in both the table and the figure, the dispersion rate predictions using Equation (18) are in good agreement with measured rates reported in the literature. Discrepancies between the predicted and measured dispersion rates may be due to the linear correlation assumption for hematocrit and/or the effect of particle concentration, which was very low for this study and was assumed to have a negligible effect on NP dispersion. Another possible reason is that platelets generally have on average an aspect ratio of 1.4, and this may be a confounding factor and contribute to the differences [77,78]. Nevertheless, the order of magnitude of the dispersion rate, as well as its trend with hematocrit are correctly predicted by Equation (18).
Table 2. Comparison of particle dispersion rate predictions from Equation (18) with dispersion rates reported in the literature.
Table 2. Comparison of particle dispersion rate predictions from Equation (18) with dispersion rates reported in the literature.
HtShear (s−1)Dispersion Rate (cm2/s)Prediction (cm2/s)Reference
[0.2, 0.4]400[0.5, 0.68] × 10−6[0.39, 0.63] × 10−6[76]
[0.2, 0.4]1100[1.5, 2.1] × 10−6[1.1, 1.7] × 10−6[76]
[0.1, 0.15, 0.2]44.8[8.2, 11.9, 17.2] × 10−9[31.3, 37.9, 44.6] × 10−9[13]
[0.1, 0.2]804[0.9, 1.4] × 10−7[5.4, 7.8] × 10−7[15]
Figure 8. Comparison of the particle dispersion rate predicted from Equation (18) with the data reported in the literature. The dashed line is the prediction from Equation (18).
Figure 8. Comparison of the particle dispersion rate predicted from Equation (18) with the data reported in the literature. The dashed line is the prediction from Equation (18).
Nanomaterials 06 00030 g008

5. Conclusion and Future Work

This paper presents a numerical study on NP dispersion in blood flow considering the influence of RBC motion and deformation. Analytical formulae, which did not exist from previous modeling results [3,16], were proposed to characterize NP migration to the vessel walls as a function of shear rate at different hematocrit levels. The formula’s predictions agree well with data reported in the literature, thereby facilitating practical modeling of NP transport and distribution in large-scale vascular systems.
In the future, the model will be extended to 3D to explore NP binding and distribution in capillary vessels. The proposed formula for the dispersion rate will also be used to evaluate NP transport and distribution in a large-scale vascular network.

Supplementary Materials

Figures S1–S9, Table S1 and three short movies are available online at https://www.mdpi.com/2079-4991/6/2/30/s1.

Acknowledgments

The study presented in this paper was supported by the National Science Foundation (NSF) Faculty Early Career Development Program (CAREER) Grant CBET-1113040, NSF CBET-1067502, NSF DMS-1516236 and the National Institute of Health (NIH) Grant EB015105.

Author Contributions

J.T. developed the code, performed the simulation and wrote the manuscript. W.K. developed the code, analyzed the data. S.S. contributed to the code and analyzed the data. J.Y. helped design the simulation. Y.L. conceived, designed the simulation and analyzed the data. All revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sanhai, W.R.; Sakamoto, J.H.; Canady, R.; Ferrari, M. Seven Challenges for Nanomedicine. Nat. Nano 2008, 3, 242–244. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, J.; Bradley, R.; Eckmann, D.M.; Ayyaswamy, P.S.; Radhakrishnan, R. Multiscale Modeling of Functionalized Nanocarriers in Targeted Drug Delivery. Curr. Nanosci. 2011, 7, 727–735. [Google Scholar] [CrossRef] [PubMed]
  3. Muller, K.; Fedosov, D.A.; Gompper, G. Margination of micro- and nano-particles in blood flow and its effect on drug delivery. Sci. Rep. 2014, 4. [Google Scholar] [CrossRef] [PubMed]
  4. Tan, J.; Thomas, A.; Liu, Y. Influence of red blood cells on nanoparticle targeted delivery in microcirculation. Soft Matter 2012, 8, 1934–1946. [Google Scholar] [CrossRef] [PubMed]
  5. Doshi, N.; Prabhakarpandian, B.; Rea-Ramsey, A.; Pant, K.; Sundaram, S.; Mitragotri, S. Flow and adhesion of drug carriers in blood vessels depend on their shape: A study using model synthetic microvascular networks. J. Control. Release 2010, 146, 196–200. [Google Scholar] [CrossRef] [PubMed]
  6. Tan, J.; Shah, S.; Thomas, A.; Ou-Yang, H.D.; Liu, Y. The influence of size, shape and vessel geometry on nanoparticle distribution. Microfluid. Nanofluid. 2013, 14, 77–87. [Google Scholar] [CrossRef] [PubMed]
  7. Tan, J.; Wang, S.; Yang, J.; Liu, Y. Coupled particulate and continuum model for nanoparticle targeted delivery. Comput. Struct. 2013, 122, 128–134. [Google Scholar] [CrossRef] [PubMed]
  8. Geng, Y.; Dalhaimer, P.; Cai, S.; Tsai, R.; Tewari, M.; Minko, T.; Discher, D.E. Shape effects of filaments versus spherical particles in flow and drug delivery. Nat. Nanotechnol. 2007, 2, 249–255. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, Y.L.; Shah, S.; Tan, J. Computational Modeling of Nanoparticle Targeted Drug Delivery. Rev. Nanosci. Nanotechnol. 2010, 1, 66–83. [Google Scholar] [CrossRef]
  10. Wang, S.; Zhou, Y.; Tan, J.; Xu, J.; Yang, J.; Liu, Y. Computational modeling of magnetic nanoparticle targeting to stent surface under high gradient field. Comput. Mech. 2014, 53, 403–412. [Google Scholar] [CrossRef] [PubMed]
  11. Guduru, R.; Liang, P.; Runowicz, C.; Nair, M.; Atluri, V.; Khizroev, S. Magneto-electric Nanoparticles to Enable Field-controlled High-Specificity Drug Delivery to Eradicate Ovarian Cancer Cells. Sci. Rep. 2013, 3. [Google Scholar] [CrossRef] [PubMed]
  12. Reasor, D.A.; Mehrabadi, M.; Ku, D.N.; Aidun, C.K. Determination of Critical Parameters in Platelet Margination. Ann. Biomed. Eng. 2013, 41, 238–249. [Google Scholar] [CrossRef] [PubMed]
  13. Saadatmand, M.; Ishikawa, T.; Matsuki, N.; Abdekhodaie, M.J.; Imai, Y.; Ueno, H.; Yamaguchi, T. Fluid particle diffusion through high-hematocrit blood flow within a capillary tube. J. Biomech. 2011, 44, 170–175. [Google Scholar] [CrossRef] [PubMed]
  14. Crowl, L.M.; Fogelson, A.L. Computational model of whole blood exhibiting lateral platelet motion induced by red blood cells. Int. J. Numer. Methods Biomed. Eng. 2010, 26, 471–487. [Google Scholar] [CrossRef] [PubMed]
  15. Zhao, H.; Shaqfeh, E.S.G.; Narsimhan, V. Shear-induced particle migration and margination in a cellular suspension. Phys. Fluids 2012, 24. [Google Scholar] [CrossRef]
  16. Lee, T.-R.; Choi, M.; Kopacz, A.M.; Yun, S.-H.; Liu, W.K.; Decuzzi, P. On the near-wall accumulation of injectable particles in the microcirculation: smaller is not better. Sci. Rep. 2013, 3. [Google Scholar] [CrossRef] [PubMed]
  17. Eckstein, E.C.; Bailey, D.G.; Shapiro, A.H. Self-Diffusion of Particles in Shear-Flow of a Suspension. J. Fluid Mech. 1977, 79, 191–208. [Google Scholar] [CrossRef]
  18. Leighton, D.; Acrivos, A. The Shear-Induced Migration of Particles in Concentrated Suspensions. J. Fluid Mech. 1987, 181, 415–439. [Google Scholar] [CrossRef]
  19. Aarts, P.A.M.M.; Steendijk, P.; Sixma, J.J.; Heethaar, R.M. Fluid shear as a possible mechanism for platelet diffusivity in flowing blood. J. Biomech. 1986, 19, 799–805. [Google Scholar] [CrossRef]
  20. Gentile, F.; Ferrari, M.; Decuzzi, P. The Transport of Nanoparticles in Blood Vessels: The Effect of Vessel Permeability and Blood Rheology. Ann. Biomed. Eng. 2008, 36, 254–261. [Google Scholar] [CrossRef] [PubMed]
  21. Gentile, F.; Decuzzi, P. Time dependent dispersion of nanoparticles in blood vessels. J. Biomed. Sci. Eng. 2010, 3, 517–524. [Google Scholar] [CrossRef]
  22. Hossain, S.S.; Zhang, Y.J.; Liang, X.H.; Hussain, F.; Ferrari, M.; Hughes, T.J.R.; Decuzzi, P. In silico vascular modeling for personalized nanoparticle delivery. Nanomedicine 2013, 8, 343–357. [Google Scholar] [CrossRef] [PubMed]
  23. Kruger, T.; Gross, M.; Raabe, D.; Varnik, F. Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells. Soft Matter 2013, 9, 9008–9015. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, J.F.; Johnson, P.C.; Popel, A.S. Red blood cell aggregation and dissociation in shear flows simulated by lattice Boltzmann method. J. Biomech. 2008, 41, 47–55. [Google Scholar] [CrossRef] [PubMed]
  25. Peskin, C.S. The immersed boundary method. Acta Numer. 2002, 11, 479–517. [Google Scholar] [CrossRef]
  26. Liu, Y.L.; Liu, W.K. Rheology of red blood cell aggregation by computer simulation. J. Comput. Phys. 2006, 220, 139–154. [Google Scholar] [CrossRef]
  27. Liu, Y.; Liu, W.K.; Belytschko, T.; Patankar, N.; To, A.C.; Kopacz, A.; Chung, J.H. Immersed electrokinetic finite element method. Int. J. Numer. Methods Eng. 2007, 71, 379–405. [Google Scholar] [CrossRef]
  28. Liu, Y.; Oh, K.; Bai, J.G.; Chang, C.; Chung, J.; Lee, K.; Liu, W.K. Manipulation of Nanoparticles and Biomolecules by Electric Field and Surface Tension. Comput. Methods Appl. Mech. Eng. 2008, 197, 2156–2172. [Google Scholar] [CrossRef]
  29. Liu, Y.; Zhang, L.T.; Wang, X.; Liu, W.K. Coupling of Navier-Stokes Equations with Protein Molecular Dynamics and Its Application to Hemodynamics. Int. J. Numer. Methods Fluids 2004, 46. [Google Scholar] [CrossRef]
  30. Liu, Y.L.; Chung, J.H.; Liu, W.K.; Ruoff, R.S. Dielectrophoretic assembly of nanowires. J. Phys. Chem. B 2006, 110, 14098–14106. [Google Scholar] [CrossRef] [PubMed]
  31. Wu, J.; Shu, C.; Zhao, N.; Yan, W.W. Fluid Dynamics of Flapping Insect Wing in Ground Effect. J. Bionic Eng. 2014, 11, 52–60. [Google Scholar] [CrossRef]
  32. Wu, J.; Yang, S.C.; Shu, C.; Zhao, N.; Yan, W.W. Ground effect on the power extraction performance of a flapping wing biomimetic energy generator. J. Fluids Struct. 2015, 54, 247–262. [Google Scholar] [CrossRef]
  33. Suzuki, K.; Minami, K.; Inamuro, T. Lift and thrust generation by a butterfly-like flapping wing-body model: Immersed boundary-lattice Boltzmann simulations. J. Fluid Mech. 2015, 767, 659–695. [Google Scholar] [CrossRef]
  34. De Rosis, A.; Leveque, E. Harmonic oscillations of a thin lamina in a quiescent viscous fluid: A numerical investigation within the framework of the lattice Boltzmann method. Comput. Struct. 2015, 157, 209–217. [Google Scholar] [CrossRef]
  35. Wang, Y.; Shu, C.; Teo, C.J.; Wu, J. An immersed boundary-lattice Boltzmann flux solver and its applications to fluid-structure interaction problems. J. Fluids Struct. 2015, 54, 440–465. [Google Scholar] [CrossRef]
  36. Tölke, J. Implementation of a Lattice Boltzmann kernel using the Compute Unified Device Architecture developed by nVIDIA. Comput. Vis. Sci. 2010, 13, 29–39. [Google Scholar] [CrossRef]
  37. Takashi, N.; Hughes, T.J. An arbitrary Lagrangian-Eulerian finite element method for interaction of fluid and a rigid body. Comput. Methods Appl. Mech. Eng. 1992, 95, 115–138. [Google Scholar] [CrossRef]
  38. Kruger, T.; Varnik, F.; Raabe, D. Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method. Comput. Math. Appl. 2011, 61, 3485–3505. [Google Scholar] [CrossRef]
  39. Sui, Y.; Chew, Y.T.; Roy, P.; Low, H.T. A hybrid method to study flow-induced deformation of three-dimensional capsules. J. Comput. Phys. 2008, 227, 6351–6371. [Google Scholar] [CrossRef]
  40. Reasor, D.A.; Clausen, J.R.; Aidun, C.K. Coupling the lattice-Boltzmann and spectrin-link methods for the direct numerical simulation of cellular blood flow. Int. J. Numer. Methods Fluids 2012, 68, 767–781. [Google Scholar] [CrossRef]
  41. Chen, S.; Doolen, G.D. Lattice boltzmann method for fluid flows. Annu. Rev. Fluid Mech. 1998, 30, 329–364. [Google Scholar] [CrossRef]
  42. Luo, L.S. The lattice-gas and lattice boltzmann methods: Past, present, and future. In Proceedings of the International Conference on Applied Computational Fluid Dynamics, Beijing, China, 17–20 October 2000; pp. 52–83.
  43. Succi, S. The Lattice Boltzmann Equation: For Fluid Dynamics and Beyond; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
  44. Ladd, A.J.C. Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann-Equation. 1. Theoretical Foundation. J. Fluid Mech. 1994, 271, 285–309. [Google Scholar] [CrossRef]
  45. Latt, J. Hydrodynamic Limit of Lattice Boltzmann Equations; University of Geneva: Geneva, Switzerland, 2007. [Google Scholar]
  46. Qian, Y.; D’Humières, D.; Lallemand, P. Lattice BGK models for Navier-Stokes equation. EPL 1992, 17. [Google Scholar] [CrossRef]
  47. Zou, Q.S.; He, X.Y. On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys. Fluids 1997, 9, 1591–1598. [Google Scholar] [CrossRef]
  48. Latt, J.; Chopard, B. Lattice Boltzmann method with regularized pre-collision distribution functions. Math. Comput. Simul. 2006, 72, 165–168. [Google Scholar] [CrossRef]
  49. Latt, J.; Chopard, B.; Malaspinas, O.; Deville, M.; Michler, A. Straight velocity boundaries in the lattice Boltzmann method. Phys. Rev. E 2008, 77. [Google Scholar] [CrossRef] [PubMed]
  50. Pozrikidis, C. Effect of membrane bending stiffness on the deformation of capsules in simple shear flow. J. Fluid Mech. 2001, 440, 269–291. [Google Scholar] [CrossRef]
  51. Skalak, R.; Tozeren, A.; Zarda, R.P.; Chien, S. Strain Energy Function of Red Blood-Cell Membranes. Biophys. J. 1973, 13, 245–280. [Google Scholar] [CrossRef]
  52. Kruger, T. Computer Simulation Study of Collective Phenomena in Dense Suspensions of Red Blood Cells under Shear; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  53. Dao, M.; Li, J.; Suresh, S. Molecularly based analysis of deformation of spectrin network and human erythrocyte. Mater. Sci. Eng. C 2006, 26, 1232–1244. [Google Scholar] [CrossRef]
  54. Fedosov, D.A.; Lei, H.; Caswell, B.; Suresh, S.; Karniadakis, G.E. Multiscale Modeling of Red Blood Cell Mechanics and Blood Flow in Malaria. PLoS Comput. Biol. 2011, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Nakamura, M.; Bessho, S.; Wada, S. Spring-network-based model of a red blood cell for simulating mesoscopic blood flow. Int. J. Numer. Methods Biomed. Eng. 2013, 29, 114–128. [Google Scholar] [CrossRef] [PubMed]
  56. Fedosov, D.A.; Caswell, B.; Karniadakis, G.E. A Multiscale Red Blood Cell Model with Accurate Mechanics, Rheology, and Dynamics. Biophys. J. 2010, 98, 2215–2225. [Google Scholar] [CrossRef] [PubMed]
  57. Nakamura, M.; Bessho, S.; Wada, S. Analysis of Red Blood Cell Deformation under Fast Shear Flow for Better Estimation of Hemolysis. Int. J. Numer. Methods Biomed. Eng. 2014, 30, 42–54. [Google Scholar] [CrossRef] [PubMed]
  58. Evans, E.; Fung, Y.-C. Improved measurements of the erythrocyte geometry. Microvasc. Res. 1972, 4, 335–347. [Google Scholar] [CrossRef]
  59. Peskin, C.S.; McQueen, D.M. A three-dimensional computational method for blood flow in the heart I. Immersed elastic fibers in a viscous incompressible fluid. J. Comput. Phys. 1989, 81, 372–405. [Google Scholar] [CrossRef]
  60. Herschlag, G.; Miller, L. Reynolds number limits for jet propulsion: A numerical study of simplified jellyfish. J. Theor. Biol. 2011, 285, 84–95. [Google Scholar] [CrossRef] [PubMed]
  61. Mittal, R.; Iaccarino, G. Immersed Boundary Methods. Annu. Rev. Fluid Mech. 2005, 37, 239–261. [Google Scholar] [CrossRef]
  62. Ladd, A.J. Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 2. Numerical results. J. Fluid Mech. 1994, 271, 311–339. [Google Scholar] [CrossRef]
  63. Aidun, C.K.; Lu, Y.N.; Ding, E.J. Direct analysis of particulate suspensions with inertia using the discrete Boltzmann equation. J. Fluid Mech. 1998, 373, 287–311. [Google Scholar] [CrossRef]
  64. Basagaoglu, H.; Allwein, S.; Succi, S.; Dixon, H.; Carrola, J.T.; Stothoff, S. Two- and three-dimensional lattice Boltzmann simulations of particle migration in microchannels. Microfluid. Nanofluid. 2013, 15, 785–796. [Google Scholar] [CrossRef]
  65. Brańka, A.; Heyes, D. Algorithms for Brownian dynamics simulation. Phys. Rev. E 1998, 58. [Google Scholar] [CrossRef]
  66. Radhakrishnan, R.; Uma, B.; Liu, J.; Ayyaswamy, P.S.; Eckmann, D.M. Temporal multiscale approach for nanocarrier motion with simultaneous adhesion and hydrodynamic interactions in targeted drug delivery. J. Comput. Phys. 2013, 244, 252–263. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, L.; Gerstenberger, A.; Wang, X.D.; Liu, W.K. Immersed finite element method. Comput. Methods Appl. Mech. Eng. 2004, 193, 2051–2067. [Google Scholar] [CrossRef]
  68. Feng, J.; Hu, H.H.; Joseph, D.D. Direct Simulation of Initial-Value Problems for the Motion of Solid Bodies in a Newtonian Fluid. 1. Sedimentation. J. Fluid Mech. 1994, 261, 95–134. [Google Scholar] [CrossRef]
  69. Fischer, T.; Stohr-Lissen, M.; Schmid-Schonbein, H. The red cell as a fluid droplet: tank tread-like motion of the human erythrocyte membrane in shear flow. Science 1978, 202, 894–896. [Google Scholar] [CrossRef] [PubMed]
  70. Sui, Y.; Low, H.T.; Chew, Y.T.; Roy, P. Tank-treading, swinging, and tumbling of liquid-filled elastic capsules in shear flow. Phys. Rev. E 2008, 77. [Google Scholar] [CrossRef] [PubMed]
  71. Pozrikidis, C. Numerical simulation of the flow-induced deformation of red blood cells. Ann. Biomed. Eng. 2003, 31, 1194–1205. [Google Scholar] [CrossRef] [PubMed]
  72. Decuzzi, P.; Ferrari, M. The adhesive strength of non-spherical particles mediated by specific interactions. Biomaterials 2006, 27, 5307–5314. [Google Scholar] [CrossRef] [PubMed]
  73. Kruger, T.; Varnik, F.; Raabe, D. Shear stress in lattice Boltzmann simulations. Phys. Rev. E 2009, 79. [Google Scholar] [CrossRef] [PubMed]
  74. Holdych, D.J.; Noble, D.R.; Georgiadis, J.G.; Buckius, R.O. Truncation error analysis of lattice Boltzmann methods. J. Comput. Phys. 2004, 193, 595–619. [Google Scholar] [CrossRef]
  75. Fedosov, D.A. Multiscale modeling of Blood flow and soft matter. In Applied Mathematics; Brown University: Providence, RI, USA, 2010. [Google Scholar]
  76. Crowl, L.; Fogelson, A.L. Analysis of mechanisms for platelet near-wall excess under arterial blood flow conditions. J. Fluid Mech. 2011, 676, 348–375. [Google Scholar] [CrossRef]
  77. Wang, W.W.; King, M.R. Multiscale Modeling of Platelet Adhesion and Thrombus Growth. Ann. Biomed. Eng. 2012, 40, 2345–2354. [Google Scholar] [CrossRef] [PubMed]
  78. Mody, N.A.; King, M.R. Platelet adhesive dynamics. Part I: Characterization of platelet hydrodynamic collisions and wall effects. Biophys. J. 2008, 95, 2539–2555. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Tan, J.; Keller, W.; Sohrabi, S.; Yang, J.; Liu, Y. Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method. Nanomaterials 2016, 6, 30. https://doi.org/10.3390/nano6020030

AMA Style

Tan J, Keller W, Sohrabi S, Yang J, Liu Y. Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method. Nanomaterials. 2016; 6(2):30. https://doi.org/10.3390/nano6020030

Chicago/Turabian Style

Tan, Jifu, Wesley Keller, Salman Sohrabi, Jie Yang, and Yaling Liu. 2016. "Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method" Nanomaterials 6, no. 2: 30. https://doi.org/10.3390/nano6020030

APA Style

Tan, J., Keller, W., Sohrabi, S., Yang, J., & Liu, Y. (2016). Characterization of Nanoparticle Dispersion in Red Blood Cell Suspension by the Lattice Boltzmann-Immersed Boundary Method. Nanomaterials, 6(2), 30. https://doi.org/10.3390/nano6020030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop