a fluctuating immersed boundary method for …donev/flucthydro/rigidibm_dfd...fluid-particle...
TRANSCRIPT
![Page 1: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/1.jpg)
A Fluctuating Immersed Boundary Method forBrownian Suspensions of Rigid Particles
Aleksandar Donev
Courant Institute, New York University
APS DFD MeetingSan Francisco, CA
Nov 23rd 2014
A. Donev (CIMS) RigidIBM 11/2014 1 / 13
![Page 2: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/2.jpg)
Fluid-Particle Coupling
Bent Active Nanorods
Figure: From the Courant Applied Math Lab of Zhang and Shelley [1]
A. Donev (CIMS) RigidIBM 11/2014 3 / 13
![Page 3: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/3.jpg)
Fluid-Particle Coupling
Thermal Fluctuation Flips
QuickTime
A. Donev (CIMS) RigidIBM 11/2014 4 / 13
![Page 4: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/4.jpg)
Fluid-Particle Coupling
Blob/Bead Models
Figure: Blob or “raspberry”models of: a spherical colloid, and a lysozyme [2].
A. Donev (CIMS) RigidIBM 11/2014 5 / 13
![Page 5: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/5.jpg)
Fluid-Particle Coupling
Immersed Rigid Bodies
In the immersed boundary method we extend the fluid velocityeverywhere in the domain,
ρ∂tv + ∇π = η∇2v −∫
Ωλ (q) δ (r − q) dq + ∇ ·
(√2ηkBT W
)∇ · v = 0 everywhere
me u = F +
∫Ωλ (q) dq
Ieω = τ +
∫Ω
[q× λ (q)] dq
v (q, t) = u + q× ω
=
∫v (r, t) δ (r − q) dr for all q ∈ Ω,
where the induced fluid-body force [3] λ (q) is a Lagrangemultiplier enforcing the final no-slip condition (rigidity).
A. Donev (CIMS) RigidIBM 11/2014 6 / 13
![Page 6: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/6.jpg)
Fluid-Particle Coupling
Body Mobility Matrix
Ignoring fluctuations, for viscous-dominated flow we can switch tothe steady Stokes equation.
For a suspension of rigid bodies define the body mobility matrix N ,
[U , Ω]T = N [F , T ]T ,
where the left-hand side collects the linear and angular velocities,and the right hand side collects the applied forces and torques.
The Brownian dynamics of the rigid bodies is given by theoverdamped Langevin equation[
UΩ
]= N
[FT
]+ (2kBTN )
12 ∇ W .
Problem: How to compute N and N12 and simulate the
Brownian motion of the bodies?
A. Donev (CIMS) RigidIBM 11/2014 7 / 13
![Page 7: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/7.jpg)
Fluid-Particle Coupling
Immersed-Boundary Method
Figure: Flow past a rigid cylinder computed using our rigid-bodyimmersed-boundary method at Re = 20. The cylinder is discretized using 121markers/blobs.
A. Donev (CIMS) RigidIBM 11/2014 8 / 13
![Page 8: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/8.jpg)
Fluid-Particle Coupling
Blob Model
The rigid body is discretized through a number of “markers” or“blobs” [4] with positions Q = q1, . . . ,qN.Take an Immersed Boundary Method (IBM) approach and describethe fluid-blob interaction using a localized smooth kernel δa(∆r) withcompact support of size a (regularized delta function).
Our methods:
Work for the steady Stokes regime (Re = 0) as well as finite Reynoldsnumbers because there is no time splitting.Strictly enforce the rigidity constraint.Ensure fluctuation-dissipation balance even in the presence ofnontrivial boundary conditions.Involve no Green’s functions, but rather, use a finite-volumestaggered-grid fluid solver to include hydrodynamics.
A. Donev (CIMS) RigidIBM 11/2014 9 / 13
![Page 9: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/9.jpg)
Fluid-Particle Coupling
Rigid-Body Immersed-Boundary Method
Rigidly-constrained Stokes linear system
∇π − η∇2v = −∑
λiδa (qi − r) +√
2ηkBT ∇ ·W
∇ · v = 0 (Lagrange multiplier is π)∑i
λi = F (Lagrange multiplier is u) (1)∑qi × λi = τ (Lagrange multiplier is ω),∫
δa (qi − r) v (r, t) dr = ui + ωi × qi + slip (Multiplier is λi )
where Λ = λ1, . . . ,λN are the unknown rigidity forces.
1 Specified kinematics (e.g., swimming object): Unknowns are v, π andΛ, while F and τ are outputs (easier).
2 Free bodies (e.g., colloidal suspension): Unknowns are v, π and Λ, uand ω, while F and τ are inputs (harder).
A. Donev (CIMS) RigidIBM 11/2014 10 / 13
![Page 10: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/10.jpg)
Fluid-Particle Coupling
Suspensions of Rigid Bodies
[U , Ω]T = N [F , T ]T ,
The many-body mobility matrix N takes into account higher-orderhydrodynamic interactions,
N =(KM−1K?
)−1,
where the blob mobility matrix M is defined by
Mij = η−1
∫δa(qi − r)G(r, r′)δa(qj − r′) drdr′ (2)
where G is the Green’s function for the Stokes problem (Oseen tensorfor infinite domain), and K is a simple geometric matrix, defined via
K? [U,Ω]T = U + Ω×Q.
A. Donev (CIMS) RigidIBM 11/2014 11 / 13
![Page 11: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/11.jpg)
Fluid-Particle Coupling
Numerical Method
The difficulty is in the numerical method for solving therigidity-constrained Stokes problem: large saddle-point system.
We use an iterative method based on a Schur complement in whichwe approximate the blob mobility matrix analytically relying on neartranslational-invariance of the Peskin IB method [5].
Fast direct solvers (related to FMMs) are required to approximatelycompute the action of M−1.
This works for confined systems, non-spherical particles,finite-Reynolds numbers and even active particles.Can also be extended to semi-rigid structures (e.g., bead-linkpolymer chains).
A. Donev (CIMS) RigidIBM 11/2014 12 / 13
![Page 12: A Fluctuating Immersed Boundary Method for …donev/FluctHydro/RigidIBM_DFD...Fluid-Particle Coupling Immersed-Boundary Method Figure:Flow past a rigid cylinder computed using our](https://reader036.vdocument.in/reader036/viewer/2022070806/5f0426f27e708231d40c91a0/html5/thumbnails/12.jpg)
Fluid-Particle Coupling
References
Daisuke Takagi, Adam B Braunschweig, Jun Zhang, and Michael J Shelley.
Dispersion of self-propelled rods undergoing fluctuation-driven flips.Phys. Rev. Lett., 110(3):038301, 2013.
Jose Garcıa de la Torre, Marıa L Huertas, and Beatriz Carrasco.
Calculation of hydrodynamic properties of globular proteins from their atomic-level structure.Biophysical Journal, 78(2):719–730, 2000.
D. Bedeaux and P. Mazur.
Brownian motion and fluctuating hydrodynamics.Physica, 76(2):247–258, 1974.
F. Balboa Usabiaga, R. Delgado-Buscalioni, B. E. Griffith, and A. Donev.
Inertial Coupling Method for particles in an incompressible fluctuating fluid.Comput. Methods Appl. Mech. Engrg., 269:139–172, 2014.Code available at https://code.google.com/p/fluam.
B. Kallemov, A. Pal Singh Bhalla, B. E. Griffith, and A. Donev.
An immersed boundary method for suspensions of rigid bodies.In preparation, to be submitted to Comput. Methods Appl. Mech. Engrg., 2015.
A. Donev (CIMS) RigidIBM 11/2014 13 / 13