a fluctuating boundary integral method for brownian
play

A fluctuating boundary integral method for Brownian suspensions - PowerPoint PPT Presentation

A fluctuating boundary integral method for Brownian suspensions Yuanxun (Bill) Bao Courant Institute, NYU Collaborators: Aleksandar Donev (Courant) Leslie Greengard (Courant) Eric Keaveny (Imperial) Manas Rachh (Yale) SIAM Computational


  1. A fluctuating boundary integral method for Brownian suspensions Yuanxun (Bill) Bao Courant Institute, NYU Collaborators: Aleksandar Donev (Courant) Leslie Greengard (Courant) Eric Keaveny (Imperial) Manas Rachh (Yale) SIAM Computational Science and Engineering March 2, 2017 Y. Bao (CIMS) Fluctuating BIE

  2. Brownian Dynamics with Hydrodynamic Interactions ⊲ Consider a suspension of N b rigid bodies with configuration Q = { q β , θ β } N b β =1 consisting of positions and orientations im- mersed in a Stokes fluid. ⊲ The Ito stochastic equation of Brownian Dynamics (BD) is d Q 2 W ( t ) + ( k B T ) ∂ Q · N , 1 dt = − N ∂ Q U + (2 k B T N ) where N ( Q ) is the body mobility matrix , U ( Q ) is the potential energy, k B T is the temperature, and W ( t ) is a vector of independent white noise processes. ⊲ Here the stochastic noise amplitude is determined from the � � ∗ 1 1 fluctuation-dissipation balance : N N = N . 2 2 ⊲ The stochastic drift term ∂ Q · N = � j ∂ j N ij is related to the Ito interpretation of the noise. Y. Bao (CIMS) Fluctuating BIE

  3. Brownian Dynamics with Hydrodynamic Interactions ⊲ Consider a suspension of N b rigid bodies with configuration Q = { q β , θ β } N b β =1 consisting of positions and orientations im- mersed in a Stokes fluid. ⊲ The Ito stochastic equation of Brownian Dynamics (BD) is d Q 2 W ( t ) + ( k B T ) ∂ Q · N , 1 dt = − N ∂ Q U + (2 k B T N ) where N ( Q ) is the body mobility matrix , U ( Q ) is the potential energy, k B T is the temperature, and W ( t ) is a vector of independent white noise processes. ⊲ Here the stochastic noise amplitude is determined from the � � ∗ 1 1 fluctuation-dissipation balance : N N = N . 2 2 ⊲ The stochastic drift term ∂ Q · N = � j ∂ j N ij is related to the Ito interpretation of the noise. Y. Bao (CIMS) Fluctuating BIE

  4. Brownian Dynamics with Hydrodynamic Interactions ⊲ Consider a suspension of N b rigid bodies with configuration Q = { q β , θ β } N b β =1 consisting of positions and orientations im- mersed in a Stokes fluid. ⊲ The Ito stochastic equation of Brownian Dynamics (BD) is d Q 2 W ( t ) + ( k B T ) ∂ Q · N , 1 dt = − N ∂ Q U + (2 k B T N ) where N ( Q ) is the body mobility matrix , U ( Q ) is the potential energy, k B T is the temperature, and W ( t ) is a vector of independent white noise processes. ⊲ Here the stochastic noise amplitude is determined from the � � ∗ 1 1 fluctuation-dissipation balance : N N = N . 2 2 ⊲ The stochastic drift term ∂ Q · N = � j ∂ j N ij is related to the Ito interpretation of the noise. Y. Bao (CIMS) Fluctuating BIE

  5. Brownian Dynamics with Hydrodynamic Interactions ⊲ Consider a suspension of N b rigid bodies with configuration Q = { q β , θ β } N b β =1 consisting of positions and orientations im- mersed in a Stokes fluid. ⊲ The Ito stochastic equation of Brownian Dynamics (BD) is d Q 2 W ( t ) + ( k B T ) ∂ Q · N , 1 dt = − N ∂ Q U + (2 k B T N ) where N ( Q ) is the body mobility matrix , U ( Q ) is the potential energy, k B T is the temperature, and W ( t ) is a vector of independent white noise processes. ⊲ Here the stochastic noise amplitude is determined from the � � ∗ 1 1 fluctuation-dissipation balance : N N = N . 2 2 ⊲ The stochastic drift term ∂ Q · N = � j ∂ j N ij is related to the Ito interpretation of the noise. Y. Bao (CIMS) Fluctuating BIE

  6. Hydrodynamic Body Mobility Matrix ⊲ The body mobility matrix N ( Q ) � 0 is a symmetric positive semidefinite (SPD) and it includes hydrodynamic interactions and (periodic) boundary conditions . ⊲ For viscous-dominated flows ( Re → 0), we can assume steady Stokes flow and solve the Stokes mobility problem , U = N F , where U = { u β , ω β } N b β =1 collects the linear and angular velocities, F = { f β , τ β } N b β =1 collects the applied forces and torques. ⊲ At every time step of BD simulation, we need to generate particle velocity in the form of (dropping k B T ), 1 � 2 W . U = N F + N ⊲ This talk: How to accurately and efficiently compute the action of 1 2 ? N and N Y. Bao (CIMS) Fluctuating BIE

  7. Hydrodynamic Body Mobility Matrix ⊲ The body mobility matrix N ( Q ) � 0 is a symmetric positive semidefinite (SPD) and it includes hydrodynamic interactions and (periodic) boundary conditions . ⊲ For viscous-dominated flows ( Re → 0), we can assume steady Stokes flow and solve the Stokes mobility problem , U = N F , where U = { u β , ω β } N b β =1 collects the linear and angular velocities, F = { f β , τ β } N b β =1 collects the applied forces and torques. ⊲ At every time step of BD simulation, we need to generate particle velocity in the form of (dropping k B T ), 1 � 2 W . U = N F + N ⊲ This talk: How to accurately and efficiently compute the action of 1 2 ? N and N Y. Bao (CIMS) Fluctuating BIE

  8. Hydrodynamic Body Mobility Matrix ⊲ The body mobility matrix N ( Q ) � 0 is a symmetric positive semidefinite (SPD) and it includes hydrodynamic interactions and (periodic) boundary conditions . ⊲ For viscous-dominated flows ( Re → 0), we can assume steady Stokes flow and solve the Stokes mobility problem , U = N F , where U = { u β , ω β } N b β =1 collects the linear and angular velocities, F = { f β , τ β } N b β =1 collects the applied forces and torques. ⊲ At every time step of BD simulation, we need to generate particle velocity in the form of (dropping k B T ), 1 � 2 W . U = N F + N ⊲ This talk: How to accurately and efficiently compute the action of 1 2 ? N and N Y. Bao (CIMS) Fluctuating BIE

  9. Hydrodynamic Body Mobility Matrix ⊲ The body mobility matrix N ( Q ) � 0 is a symmetric positive semidefinite (SPD) and it includes hydrodynamic interactions and (periodic) boundary conditions . ⊲ For viscous-dominated flows ( Re → 0), we can assume steady Stokes flow and solve the Stokes mobility problem , U = N F , where U = { u β , ω β } N b β =1 collects the linear and angular velocities, F = { f β , τ β } N b β =1 collects the applied forces and torques. ⊲ At every time step of BD simulation, we need to generate particle velocity in the form of (dropping k B T ), 1 � 2 W . U = N F + N ⊲ This talk: How to accurately and efficiently compute the action of 1 2 ? N and N Y. Bao (CIMS) Fluctuating BIE

  10. First-Kind Boundary Integral Formulation ⊲ Let us first ignore Brownian terms and solve a mobility problem to compute N F . ⊲ For simplicity, consider only a single body Ω. The first-kind bound- ary integral equation for the mobility problem, � G ( q − q ′ ) µ ( q ′ ) d q ′ v ( q ) = u + ω × q = for all q ∈ ∂ Ω , (1) ∂ Ω along with force and torque balance conditions � � µ ( q ) d q = f and q × µ ( q ) d q = τ , (2) ∂ Ω ∂ Ω where µ ( q ∈ ∂ Ω) is the surface traction (single-layer density) and G is the (periodic) Stokeslet. ⊲ Note that one can alternatively use a completed second-kind or a mixed first-second kind formulation for improved conditioning. We only know how to generate Brownian displacements efficiently in the first-kind formulation. Y. Bao (CIMS) Fluctuating BIE

  11. First-Kind Boundary Integral Formulation ⊲ Let us first ignore Brownian terms and solve a mobility problem to compute N F . ⊲ For simplicity, consider only a single body Ω. The first-kind bound- ary integral equation for the mobility problem, � G ( q − q ′ ) µ ( q ′ ) d q ′ v ( q ) = u + ω × q = for all q ∈ ∂ Ω , (1) ∂ Ω along with force and torque balance conditions � � µ ( q ) d q = f and q × µ ( q ) d q = τ , (2) ∂ Ω ∂ Ω where µ ( q ∈ ∂ Ω) is the surface traction (single-layer density) and G is the (periodic) Stokeslet. ⊲ Note that one can alternatively use a completed second-kind or a mixed first-second kind formulation for improved conditioning. We only know how to generate Brownian displacements efficiently in the first-kind formulation. Y. Bao (CIMS) Fluctuating BIE

  12. First-Kind Boundary Integral Formulation ⊲ Let us first ignore Brownian terms and solve a mobility problem to compute N F . ⊲ For simplicity, consider only a single body Ω. The first-kind bound- ary integral equation for the mobility problem, � G ( q − q ′ ) µ ( q ′ ) d q ′ v ( q ) = u + ω × q = for all q ∈ ∂ Ω , (1) ∂ Ω along with force and torque balance conditions � � µ ( q ) d q = f and q × µ ( q ) d q = τ , (2) ∂ Ω ∂ Ω where µ ( q ∈ ∂ Ω) is the surface traction (single-layer density) and G is the (periodic) Stokeslet. ⊲ Note that one can alternatively use a completed second-kind or a mixed first-second kind formulation for improved conditioning. We only know how to generate Brownian displacements efficiently in the first-kind formulation. Y. Bao (CIMS) Fluctuating BIE

Recommend


More recommend