“…It is convenient to assume that e is located midway between P and E, w midway between W and P, A midway between S and P, and n midway between N and P. The upwind schcme was used to solve the current continuity equation to avoid numerically unstable, whereby space charge density at the grid point outside the control volume is used for faccs for which current leaves the control volume. The direction of current density through an control volume face is determined by the sign of the twopoint diffcrence of the potentials on both sides of the face (Levin and Hoburg, 1990). For example, the component of current dcnsity that enters the control volume from point W to point P is expressed as while the component of current density leaving the control volume from point P to point N is expressed as Therefore, integration of Eq.…”