previous up next contents Previous: 8.4 Discretization Up: 8.4 Discretization Next: 8.4.2 Drift Term

8.4.1 Free Term

In the stationary case without scattering Equation 8.1 degenerates for $ k = 0$ as the Wigner potential $ V_{\mathrm{w}}$ is zero for $ k = 0$. If also the potential is zero, i.e., the free case, then left ($ k > 0$) and right ($ k < 0$) going modes do not couple at all. Accordingly, the solution might become discontinuous for $ k = 0$. Discontinuity of the Wigner function for $ k = 0$ is considered as a feature of ballistic transport and is unavoidable under certain types of inflow boundary conditions. A possibility to get an equation for $ k = 0$ is to divide Equation 8.1 by $ k$ and to consider the limit equation for $ k \to 0$.

From meshing constraints discussed in the previous section it follows that the number of grid points $ N_k$ is an even number. Hence we can choose the equidistant mesh symmetric around $ k = 0$ excluding the point $ k = 0$.

We implemented four different discretizations for the derivative $ \partial_x$: 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 $ m-1$ if the internal discretization is to be exact of order $ m$. Otherwise the overall order will decrease and be less than $ m$.

The discrete formula for the calculation of the current $ j$ 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 up next contents Previous: 8.4 Discretization Up: 8.4 Discretization Next: 8.4.2 Drift Term


R. Kosik: Numerical Challenges on the Road to NanoTCAD