Resolving the Gibbs phenomenon via a discontinuous basis in a mode solver for open optical systems

Parry Y. Chen, Yonatan Sivan

    Research output: Contribution to journalArticlepeer-review

    10 Scopus citations

    Abstract

    Partial differential equations are frequently solved using a global basis, such as the Fourier series, due to excellent convergence. However, convergence becomes impaired when discontinuities are present due to the Gibbs phenomenon, negatively impacting simulation speed and possibly generating spurious solutions. We resolve this by supplementing the smooth global basis with an inherently discontinuous basis, incorporating knowledge of the location of the discontinuities. The solution's discontinuities are reproduced with exponential convergence, expediting simulations. The highly constrained discontinuous basis also eliminates the freedom to generate spurious solutions. We employ the combined smooth and discontinuous bases to construct a solver for the modes of a resonator in an open electromagnetic system. These modes can then expand any scattering problem for any source configuration or incidence condition without further numerics, enabling ready access and physical insight into the spatial variation of Green's tensor. Solving for the modes is the most numerically intensive and difficult step of modal expansion methods, so our mode solver overcomes the last major impediment to the use of modal expansion for open systems.

    Original languageEnglish
    Article number110004
    JournalJournal of Computational Physics
    Volume429
    DOIs
    StatePublished - 15 Mar 2021

    Keywords

    • Discontinuous basis
    • Exponential convergence
    • Gibbs phenomenon
    • Modal expansion
    • Open systems
    • Step discontinuity

    ASJC Scopus subject areas

    • Numerical Analysis
    • Modeling and Simulation
    • Physics and Astronomy (miscellaneous)
    • General Physics and Astronomy
    • Computer Science Applications
    • Computational Mathematics
    • Applied Mathematics

    Fingerprint

    Dive into the research topics of 'Resolving the Gibbs phenomenon via a discontinuous basis in a mode solver for open optical systems'. Together they form a unique fingerprint.

    Cite this