A nonlinear and partially implicit finite-difference time-domain scheme with a relaxed stability condition is presented for fast full-wave simulation of electrostatic discharges occurred at a short gap. The nonlinear spark resistance model of Rompe and Weizel is directly incorporated in this scheme. Its stability, accuracy and computational efficiency are assessed in the application to a metallic structure system. Introduction: Full-wave modelling of electrostatic discharges (ESDs) has been developed and extended for different numerical methods [1][2][3][4]. The finite-difference time-domain (FDTD) method is a computationally efficient approach for simulating discharge currents and transient fields caused by ESDs in a system because it is explicit, i.e. requiring no matrix solution. In [1,2], the classical explicit FDTD scheme is extended to include the spark resistance model by Rompe and Weizel (RW) [5]. The extended FDTD scheme is nonlinear but still explicit, and therefore the so-called Courant-Friedrichs-Lewy (CFL) condition must be satisfied to ensure its stability. If a short gap compared with other geometries exists, the maximum stable time step allowed for this scheme can be severely limited by a small mesh width inside the gap. This will make the extended explicit scheme [1, 2] time-consuming.In this Letter, a nonlinear and partially implicit FDTD scheme is presented for fast full-wave simulation of ESD occurred at a short gap. This scheme is based on combining the magnetically mixed Newmark-Leapfrog (MNL) FDTD method [6,7] with the RW model. Its stability condition does not include the mesh width in a spark channel direction, and thus the number of time steps can be reduced by using a larger step size than the CFL condition. As an application, geometry dependency of short-gap ESD caused in a metallic structure system is studied. It is shown that the presented scheme is faster than the extended explicit one, while maintaining the same level of accuracy.