A Robust GPU-Accelerated Kernel Compensation Solver with Novel Discretization for Photonic Crystals in Anisotropic Media
This work provides a robust numerical method for simulating photonic crystals with anisotropic media, which is important for computational electromagnetics and photonics researchers.
The authors developed a GPU-accelerated solver for Maxwell eigenproblems in 3D photonic crystals with anisotropic media, using kernel compensation and a novel discretization for off-diagonal permittivity tensors, achieving robust and accurate results on benchmarks.
This paper develops a robust solver for the Maxwell eigenproblem in 3D photonic crystals with anisotropic media. The solver employs the kernel compensation technique under the framework of Yee's scheme to eliminate null space and enable matrix-free, GPU-accelerated operations via 3D discrete Fourier transform. Furthermore, we propose a novel discretization for permittivity tensor containing off-diagonal entries and prove that the resulting matrix is Hermitian positive definite, which ensures the correctness of the kernel compensation technique. Numerical experiments on several benchmark examples are demonstrated to validate the robustness and accuracy of our scheme.