Before revision:  The boundary conditions are enforced weakly through the numerical flux of DG. Numerical results show that the proposed DG-based unit-cell scheme models the device accurately but is significantly faster than the DG scheme that takes into account the whole device  .
After revision:  The boundary conditions are enforced weakly through the numerical flux of DG. Numerical results show that the proposed DG-based unit-cell scheme models the device accurately but is significantly faster than the DG scheme that takes into account the whole device . Additionally, the proposed scheme is used for the first-ever numerical demonstration of optical- and radiation-field screening effects on PCD response. The optical-field screening is found to play a more dominant role in the saturation of PCD output at high levels of optical pump power .