Numerical solver with cip method for Fokker Planck equation of stochastic cooling

A numerical solver for a Fokker-Planck equation in a stochastic momentum cooling process by using a CIP method is developed. The Fokker-Planck equation is separated to a coherent phase and an incoherent phase for a numerical calculation procedure. The incoherent term is written as the diffusion equation, while the coherent term is described by the advection equation. For solving the advection term of the equation, we use the CIP method, and the results are compared with other numerical schemes. It is found that the rational function based CIP method is effective as the numerical solver for the stochastic momentum cooling model.