A hybrid electromagnetics (EM)-circuit simulation method employing the discontinuous Galerkin finite element time domain method (DGFETD) is developed to model single lumped port networks comprised of both linear and non-linear elements. The whole computational domain is split into two subsystems. One is the EM subsystem that is analyzed by the DGFETDwhile another is the circuit circuit subsystem that is modeled by the Modified Nodal Analysis method (MNA) to generate a circuit subsystem. The coupling between the EM and circuit subsystems is achieved through a lumped port. Due to the local properties of DGFETD operations, only small coupling matrix equation systems are involved. To solve non-linear devices, the standard Newton-Raphson method is applied to solve the established non-linear system equations. Numerical examples are presented to validate the proposed algorithm.