“…( 42) i.e., the p th order moment general concentration partial differential equation under boundary conditions ( 43) and ( 44), a finite difference method based on the Crank Nicholson implicit scheme is employed. Further details of numeric are given in [56,57].…”