From meshing constraints discussed in the previous section it follows that the number of grid points is an even number. Hence we can choose the equidistant mesh symmetric around excluding the point .
We implemented four different discretizations for the derivative : upwinding first and second order and central differences (second and fourth order), see [Bie97]. Special care has to be taken of the boundary points and the points next to the boundary. The discretization at the boundary must be at least exact of order if the internal discretization is to be exact of order . Otherwise the overall order will decrease and be less than .
The discrete formula for the calculation of the current is fixed by the discretization of the free streaming operator. If this is not respected then the calculated current density will not be constant in one dimension and conservation of mass is severely violated. This is further discussed in Section 10.1.
Previous: 8.4 Discretization Up: 8.4 Discretization Next: 8.4.2 Drift Term