method
NeumannY.
prepare_rhs
Compute the distribution function at the equilibrium with the value on the border.
simulation class