OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 18, Iss. 13 — Jun. 21, 2010
  • pp: 13679–13692

Acceleration of FDTD mode solver by high-performance computing techniques

Lin Han, Yanping Xi, and Wei-Ping Huang  »View Author Affiliations

Optics Express, Vol. 18, Issue 13, pp. 13679-13692 (2010)

View Full Text Article

Enhanced HTML    Acrobat PDF (1136 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



A two-dimensional (2D) compact finite-difference time-domain (FDTD) mode solver is developed based on wave equation formalism in combination with the matrix pencil method (MPM). The method is validated for calculation of both real guided and complex leaky modes of typical optical waveguides against the bench-mark finite-difference (FD) eigen mode solver. By taking advantage of the inherent parallel nature of the FDTD algorithm, the mode solver is implemented on graphics processing units (GPUs) using the compute unified device architecture (CUDA). It is demonstrated that the high-performance computing technique leads to significant acceleration of the FDTD mode solver with more than 30 times improvement in computational efficiency in comparison with the conventional FDTD mode solver running on CPU of a standard desktop computer. The computational efficiency of the accelerated FDTD method is in the same order of magnitude of the standard finite-difference eigen mode solver and yet require much less memory (e.g., less than 10%). Therefore, the new method may serve as an efficient, accurate and robust tool for mode calculation of optical waveguides even when the conventional eigen value mode solvers are no longer applicable due to memory limitation.

© 2010 OSA

OCIS Codes
(030.4070) Coherence and statistical optics : Modes
(200.4960) Optics in computing : Parallel processing
(230.0230) Optical devices : Optical devices
(230.7370) Optical devices : Waveguides

ToC Category:
Physical Optics

Original Manuscript: April 26, 2010
Revised Manuscript: June 4, 2010
Manuscript Accepted: June 5, 2010
Published: June 10, 2010

Lin Han, Yanping Xi, and Wei-Ping Huang, "Acceleration of FDTD mode solver by high-performance computing techniques," Opt. Express 18, 13679-13692 (2010)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. E. Schweig and W. Bridges, “Computer analysis of dielectric waveguides: A finite-difference method,” IEEE Trans. Microw. Theory Tech. 32(5), 531–541 (1984). [CrossRef]
  2. M. Stern, “Semivectorial polarised finite difference method for opticalwaveguides with arbitrary index profiles,” IEE Proc. Optoelectron. 135, 56–63 (1988). [CrossRef]
  3. W. P. Huang, S. T. Chu, A. Goss, and S. K. Chaudhuri, “A scalar finite-difference time-domain approach to guided-wave optics,” IEEE Photon. Technol Lett. 3, 524–526 (1991).
  4. A. Fallahkhair, K. Li, and T. Murphy, “Vector finite difference modesolver for anisotropic dielectric waveguides,” J. Lightwave Technol. 26(11), 1423–1431 (2008). [CrossRef]
  5. C. Xu, W. Huang, M. Stern, and S. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc. Optoelectron. 141(5), 281–286 (1994). [CrossRef]
  6. B. Rahman and J. Davies, “Finite-element analysis of optical and microwave waveguide problems,” IEEE Trans. Microw. Theory Tech. 32(1), 20–28 (1984). [CrossRef]
  7. J. Lee, D. Sun, and Z. Cendes, “Full-wave analysis of dielectric waveguides using tangential vector finite elements,” IEEE Trans. Microw. Theory Tech. 39(8), 1262–1271 (1991). [CrossRef]
  8. M. D. Feit and J. A. Fleck., “Computation of mode properties in optical fiber waveguides by a propagating beam method,” Appl. Opt. 19(7), 1154–1164 (1980). [CrossRef] [PubMed]
  9. C. Xu, W. Huang, and S. Chaudhuri, ““Efficient and accurate vector mode calculations by beam propagationmethod,” Lightwave Technology,” Journalism 11, 1209–1215 (1993).
  10. Y. Tsuji and M. Koshiba, “Guided-mode and leaky-mode analysis by imaginary distance beam propagation method based on finite element scheme,” J. Lightwave Technol. 18(4), 618–623 (2000). [CrossRef]
  11. A. Zhao, J. Juntunen, and A. Räisänen, “Analysis of hybrid modes in channel multilayer optical waveguides with the compact 2-D FDTD method,” Microw. Opt. Technol. Lett. 15(6), 398–403 (1997). [CrossRef]
  12. S. Xiao, R. Vahldieck, and H. Jin, “Full-wave analysis of guided wave structures using a novel 2-D FDTD,” IEEE Microw. Guid. Wave Lett. 2(5), 165–167 (1992). [CrossRef]
  13. G. Zhou and X. Li, “Wave equation-based semivectorial compact 2-D-FDTD method for optical waveguide modal analysis,” J. Lightwave Technol. 22(2), 677–683 (2004). [CrossRef]
  14. T. Sarkar and O. Pereira, “Using the matrix pencil method to estimate the parameters of a sum of complex exponentials,” IEEE Antennas Propag. Mag. 37(1), 48–55 (1995). [CrossRef]
  15. Z. Sacks, D. Kingsland, R. Lee, and J. Lee, “A perfectly matched anisotropic absorber for use as an absorbingboundary condition,” IEEE Trans. Antenn. Propag. 43(12), 1460–1463 (1995). [CrossRef]
  16. S. Gedney, “An anisotropic perfectly matched layer-absorbing medium for thetruncation of FDTD lattices,” IEEE Trans. Antenn. Propag. 44(12), 1630–1639 (1996). [CrossRef]
  17. A. Taflove and S. Hagness, Computational electrodynamics: the finite-difference time-domain method (Artech House Norwood, MA, 1995).
  18. S. Ryoo, C. Rodrigues, S. Baghsorkhi, S. Stone, D. Kirk, and W. Hwu, “Optimization principles and application performance evaluation of a multithreaded GPU using CUDA,” in (ACM, 2008), 73–82.
  19. S. Krakiwsky, L. Turner, and M. Okoniewski, “Acceleration of finite-difference time-domain (FDTD) using graphics processor units (GPU),” in Microwave Symposium Digest, IEEE MTT-S International, 2004), 1033- 1036.
  20. A. Balevic, L. Rockstroh, A. Tausendfreund, S. Patzelt, G. Goch, and S. Simon, “Accelerating simulations of light scattering based on finite-difference time-domain method with general purpose gpus,” in 2008), 327–334.
  21. M. Okoniewski, “Vector wave equation 2-D-FDTD method for guided wave problems,” IEEE Microw. Guid. Wave Lett. 3(9), 307–309 (1993). [CrossRef]
  22. J. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys. 114(2), 185–200 (1994). [CrossRef]
  23. J. Goell, “A circular-harmonic computer analysis of rectangular dielectric waveguides,” Bell Syst. Tech. J. 48, 2133–2160 (1969).
  24. J. Heaton, M. Bourke, S. Jones, B. Smith, K. Hilton, G. Smith, J. Birbeck, G. Berry, S. Dewar, and D. Wight, “Optimization of deep-etched, single-mode GaAs/AlGaAs optical waveguides using controlled leakage into the substrate,” J. Lightwave Technol. 17(2), 267–281 (1999). [CrossRef]
  25. W. Huang and C. Xu, “Simulation of three-dimensional optical waveguides by a full-vectorbeam propagation method,” IEEE J. Quantum Electron. 29(10), 2639–2649 (1993). [CrossRef]
  26. W. Chew and W. Weedon, “A 3-D perfectly matched medium from modified Maxwell's equations with stretched coordinates,” Microw. Opt. Technol. Lett. 7(13), 599–604 (1994). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited