Abstract
A discontinuous Galerkin (DG) framework is developed for multiphysics simulation of plasmonic photoconductive devices. The nonequilibrium steady state is modeled by the coupled system of Poisson and drift-diffusion (DD) equations and this system is solved using a nonlinear DG scheme. The optical-to-terahertz (THz) conversion process is modeled by the coupled system of Maxwell and time dependent DD equations and this system is solved using a time domain DG scheme. Numerical experiments show that the proposed scheme provides accurate results.