Slide 30
Slide 30 text
PyStokes: Usage
import pystokes , p y f o r c e s
import numpy as np
a , Np = 1 , 100
L , dim = 128 , 3
dt , T = 0.01 , 100
v = np . zeros ( dim∗Np ) ; r = v ; F = v ;
pRbm = pystokes . p e r i o d i c .Rbm(a , Np, L)
f f = p y f o r c e s . f o r c e F i e l d s . Forces (Np)
for t t in range (T) :
f f . sedimentation (F , g=−10)
pRbm. s t o k e s l e t V ( v , r , F)
r = ( r + (F/(0.75∗ a ) + v )∗ dt)%L
12