A finite-difference-based numerical simulator for diffusophoretic, self-propelling particles.
Get the source:
git clone https://github.com/NaokiHori/SwimmingParticle
cd SwimmingParticle
Build and run:
make output
make all
./a.out
Check the result:
python visualize.py
There are four parameters which can be changed:
- Péclet number
- Radius of the outer circle
- Number of radial and azimuthal grid points
They are all defined in src/main.c
.
Note that the number of azimuthal grids are to be a power of 2.
Re-compile the source (make all
) when modified.
Governing equation and its numerical treatment are briefly documented here.
- Williamson, J. Comput. Phys. (35), 1980
- Hu et al., Phys. Rev. Lett. (123), 2019
- Wood and Porter, J. Appl. Eng. Math. (6), 2019