24.02.2014 Views

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Cell Adhesion <strong>and</strong> Detachment in Shear Flow 219<br />

Table 1. Simulation parameters.<br />

Parameters Definition simulation value<br />

R cell radius 4.0–5.0 µm<br />

N r receptor number 780<br />

N L lig<strong>and</strong> density 10 6 –10 8 /cm<br />

λ equilibrium bond length 0.2 µm<br />

σ spring constant 0.016 dyne/cm<br />

µ viscosity 0.01–0.014 g/cm·s<br />

ρ fluid density 1.0 g/cm 2<br />

U max shear rate 20–80/s<br />

H c cut-off length 0.4 µm<br />

T temperature 310 K<br />

kf 0 forward reaction rate 100.0/s<br />

kr 0 reverse reaction rate 10.0/s<br />

r 0 reactive compliance 0.02 µm<br />

with the long semi-axis r a equal to 0.5 <strong>and</strong> the short semi-axis r b equal to<br />

0.4. The velocity boundary conditions are as follows: a given constant on the<br />

top boundary, zero on the bottom boundary, <strong>and</strong> periodicity in the horizontal<br />

direction. The fluid <strong>and</strong> cells are at rest <strong>and</strong> the cells are in the contact region<br />

initially (see Fig. 6(a)). We assume that the densities of fluid <strong>and</strong> cells are<br />

1 g/cm 3 . The mesh size h for the flow field is 1/48 <strong>and</strong> the time step △t is<br />

0.001 (unit: 0.1 second). The parameters used in the simulations are given in<br />

Table 1.<br />

We observed the simulations up to t = 100 (10 s), long enough for the flow<br />

to be fully developed. The simulations were conducted at different shear rates<br />

<strong>and</strong> dynamical viscosities, <strong>and</strong> the results are summarized in Table 2. From<br />

the table, we can see, no cells were detached from the wall by the observed<br />

time when the shear rate is 20/s for the dynamical viscosity of 0.01 g/cm-s;<br />

while the detachment percentage increases from 10% to 40% when the shear<br />

rate increases from 30/s to 40/s. All the 20 cells were detached from the wall<br />

when the shear rate is greater than 80/s. Figure 5 shows the effect of shear<br />

rate on cell detachment. This observation qualitatively agrees with the in<br />

vitro experiment [SKE + 99]. We also observed that the detachment percentage<br />

increases from 10% to 35% when the dynamical viscosity is increased from 0.01<br />

to 0.014 (g/cm-s).<br />

Figure 6 shows the snapshots of positions of 20 cells at t =0,5,5.35,<br />

6.06, 9.49, <strong>and</strong> 10 (s), for the simulation with the dynamical viscosity equal<br />

to 0.01 (g/cm-s) <strong>and</strong> the shear rate of 30 (/s). The snapshots quite clearly<br />

depict the process of cell detachment from the wall. All the cells adhered to<br />

the wall at t = 5 s; one cell was about to be detached at t =5.35 s; one cell<br />

was completely detached from the layer at t =6.06 s. We found that during<br />

the early stage of detachment the percentage of the detached cells is highly

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!