Squirmer

From testwiki
Jump to navigation Jump to search

Template:Short description

Spherical microswimmer in Stokes flow

The squirmer is a model for a spherical microswimmer swimming in Stokes flow. The squirmer model was introduced by James Lighthill in 1952 and refined and used to model Paramecium by John Blake in 1971.[1] [2] Blake used the squirmer model to describe the flow generated by a carpet of beating short filaments called cilia on the surface of Paramecium. Today, the squirmer is a standard model for the study of self-propelled particles, such as Janus particles, in Stokes flow.[3]

Velocity field in particle frame

Here we give the flow field of a squirmer in the case of a non-deformable axisymmetric spherical squirmer (radius R).[1][2] These expressions are given in a spherical coordinate system.

ur(r,θ)=23(R3r31)B1P1(cosθ)+n=2(Rn+2rn+2Rnrn)BnPn(cosθ),
uθ(r,θ)=23(R32r3+1)B1V1(cosθ)+n=212(nRn+2rn+2+(2n)Rnrn)BnVn(cosθ).

Here Bn are constant coefficients, Pn(cosθ) are Legendre polynomials, and Vn(cosθ)=2n(n+1)θPn(cosθ).
One finds P1(cosθ)=cosθ,P2(cosθ)=12(3cos2θ1),,V1(cosθ)=sinθ,V2(cosθ)=12sin2θ,.
The expressions above are in the frame of the moving particle. At the interface one finds uθ(R,θ)=n=1BnVn and ur(R,θ)=0.

Shaker, β=
Pusher, β=5
Neutral, β=0
Puller, β=5
Shaker, β=
Passive particle
Shaker, β=
Pusher, β=5
Neutral, β=0
Puller, β=5
Shaker, β=
Passive particle
Velocity field of squirmer and passive particle (top row: lab frame, bottom row: swimmer frame, β=B2/|B1| ).

Template:Clear

Swimming speed and lab frame

By using the Lorentz Reciprocal Theorem, one finds the velocity vector of the particle 𝐔=12𝐮(R,θ)sinθdθ=23B1𝐞z. The flow in a fixed lab frame is given by 𝐮L=𝐮+𝐔:

urL(r,θ)=R3r3UP1(cosθ)+n=2(Rn+2rn+2Rnrn)BnPn(cosθ),
uθL(r,θ)=R32r3UV1(cosθ)+n=212(nRn+2rn+2+(2n)Rnrn)BnVn(cosθ).

with swimming speed U=|𝐔|. Note, that limr𝐮L=0 and urL(R,θ)0.

Structure of the flow and squirmer parameter

The series above are often truncated at n=2 in the study of far field flow, rR. Within that approximation, uθ(R,θ)=B1sinθ+12B2sin2θ, with squirmer parameter β=B2/|B1|. The first mode n=1 characterizes a hydrodynamic source dipole with decay 1/r3 (and with that the swimming speed U). The second mode n=2 corresponds to a hydrodynamic stresslet or force dipole with decay 1/r2.[4] Thus, β gives the ratio of both contributions and the direction of the force dipole. β is used to categorize microswimmers into pushers, pullers and neutral swimmers.[5]

Swimmer Type pusher neutral swimmer puller shaker passive particle
Squirmer Parameter β<0 β=0 β>0 β=±
Decay of Velocity Far Field 𝐮1/r2 𝐮1/r3 𝐮1/r2 𝐮1/r2 𝐮1/r
Biological Example E.Coli Paramecium Chlamydomonas reinhardtii

The above figures show the velocity field in the lab frame and in the particle-fixed frame. The hydrodynamic dipole and quadrupole fields of the squirmer model result from surface stresses, due to beating cilia on bacteria, or chemical reactions or thermal non-equilibrium on Janus particles. The squirmer is force-free. On the contrary, the velocity field of the passive particle results from an external force, its far-field corresponds to a "stokeslet" or hydrodynamic monopole. A force-free passive particle doesn't move and doesn't create any flow field.

See also

References

Template:Reflist