Codes

Here I share some of my codes used in the previous works.

-

Spectual method of solving 2D nonlinear reaction-diffusion equations in 2D circular region

To solve nonlinear reaction-diffusion equations of the form

\[\partial_t f_i=\sum_{k=1}^N\nabla\cdot(h_{ik}\nabla f_k)+g_i(f_1,f_2,\cdots,f_N)\]

in 2D circular region. But note it has unsolved divergence problem on the system boundary.

-

Non-interacting run-and-tumble particles on 2D lattice, using continuous time Gillespie algorithm

I also have an efficient version with Logistic population dynamics on site, and a version with CUDA.

-

Differential dynamic microscopy

Extracting information of particle motion from high-speed timelapse images, based on a numerical inverse Laplace transformation of theoretical model prediction. For details please check my PhD thesis。

-

Molecular dynamics simulation of interacting active particles

-

Simulation of non-interacting Brownian active particles