Simulating Growth Kinetics in a Data-Parallel 3D Lattice Photobioreactor

Algorithm 3

Kawasaki site exchange computations.

choose a random neighbour site 𝑗

introduce low probability bias 𝑏 if upward neighbour is

chosen

compute energy change if 𝑖,𝑗 exchanged

if energy falls then

 accept change and do exchange

else

 compute Metropolis probability 𝑃

 obtain random probability 𝑟₁

 accept change conditionally on 𝑟₁<𝑃+𝑏

end if