3D Crank-Nicolson finite difference time domain method for dispersive media

The unconditionally stable Crank-Nicolson finite difference time domain (CN-FDTD) method is extended to incorporate frequency-dependent media in three dimensions. A Gaussian-elimination-based direct sparse solver is used to deal with the large sparse matrix system arising from the formulation. Numerical results validate and confirm that the scheme is unconditionally stable for time steps over the ...