Configuration for fluid-Poisson model: Element type: Linear Lagrange elements for particle balance equations, quadratic Lagrange elements for Poisson equation, linear discontinuous Lagrange elements for the surface charge density balance equation; Time-stepping: Backward differentiation formula (BDF), with automatically adapted time-step size and the BDF order (in the range from 1 to 5); Solver: Fully coupled constant Newton solver, with the direct linear solver (PARDISO); Tolerance: 0.0001; Mesh details: 1) t: [0, 3.9] microseconds, number of mesh elements: 360194, DoF: 3340893; 2) t: [3.9, 3.94793] microseconds, number of mesh elements: 360194, DoF: 3340893; 3) t: [3.94793, 3.955] microseconds, number of mesh elements: 401562, DoF: 3726070; 4) t: [3.955, 4] microseconds, number of mesh elements: 572190, DoF: 5172073; 5) t: [4, 4.03] microseconds, number of mesh elements: 464810, DoF: 4153330