Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media
Tóm tắt
In this work, we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion, the instantaneous nonlinear cubic Kerr response, and the nonlinear delayed Raman molecular vibrational response. Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al. (J Comput Phys 350: 420–452, 2017) and Lyu et al. (J Sci Comput 89: 1–42, 2021), a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part (i.e., the auxiliary differential equations) modeling the linear and nonlinear dispersion in the material. The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization. A nodal discontinuous Galerkin (DG) method is further applied in space for efficiently handling nonlinear terms at the algebraic level, while preserving the energy stability and achieving high-order accuracy. Indeed with
$$d_E$$
as the number of the components of the electric field, only a
$$d_E\times d_E$$
nonlinear algebraic system needs to be solved at each interpolation node, and more importantly, all these small nonlinear systems are completely decoupled over one time step, rendering very high parallel efficiency. We evaluate the proposed schemes by comparing them with the methods in Bokil et al. (2017) and Lyu et al. (2021) (implemented in nodal form) regarding the accuracy, computational efficiency, and energy stability, by a parallel scalability study, and also through the simulations of the soliton-like wave propagation in one dimension, as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric (TE) mode of the equations.
Từ khóa
Tài liệu tham khảo
Agrawal, G.P.: Nonlinear fiber optics. In: Christiansen, P.L., Sørensen, M.P., Scott, A.C. (eds.) Nonlinear Science at the Dawn of the 21st Century, pp. 195–211. Springer (2000)
Bey, K.S., Oden, J.T., Patra, A.: A parallel hp-adaptive discontinuous Galerkin method for hyperbolic conservation laws. Appl. Numer. Math. 20(4), 321–336 (1996)
Biswas, R., Devine, K.D., Flaherty, J.E.: Parallel, adaptive finite element methods for conservation laws. Appl. Numer. Math. 14(1/2/3), 255–283 (1994)
Bokil, V.A., Cheng, Y., Jiang, Y., Li, F.: Energy stable discontinuous Galerkin methods for Maxwell’s equations in nonlinear optical media. J. Comput. Phys. 350, 420–452 (2017)
Bokil, V.A., Cheng, Y., Jiang, Y., Li, F., Sakkaplangkul, P.: High spatial order energy stable FDTD methods for Maxwell’s equations in nonlinear optical media in one dimension. J. Sci. Comput. 77(1), 330–371 (2018)
Boyd, R.W.: Nonlinear Optics. Academic Press, New York (2003)
Chen, C.-L.: Optical Solitons in Optical Fibers. Wiley, New York (2005)
Cockburn, B., Singler, J.R., Zhang, Y.: Interpolatory HDG method for parabolic semilinear PDEs. J. Sci. Comput. 79(3), 1777–1800 (2019)
Fezoui, L., Lanteri, S.: Discontinuous Galerkin methods for the numerical solution of the nonlinear Maxwell equations in 1d. PhD thesis, INRIA (2015)
Fisher, A., White, D., Rodrigue, G.: An efficient vector finite element method for nonlinear electromagnetic modeling. J. Comput. Phys. 225(2), 1331–1346 (2007)
Gilles, L., Hagness, S.C., Vázquez, L.: Comparison between staggered and unstaggered finite-difference time-domain grids for few-cycle temporal optical soliton propagation. J. Comput. Phys. 161(2), 379–400 (2000)
Greene, J., Taflove, A.: Scattering of spatial optical solitons by subwavelength air holes. Microw. Wirel. Compon. Lett. IEEE 17, 760–762 (2007)
Greene, J.H., Taflove, A.: General vector auxiliary differential equation finite-difference time-domain method for nonlinear optics. Opt. Express 14(18), 8305–8310 (2006)
Hesthaven, J.S., Warburton, T.: Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications. Springer Science & Business Media, Berlin (2007)
Jackson, J.D.: Classical Electrodynamics. Wiley, New York (1999)
Lubin, Z., Greene, J.H., Taflove, A.: FDTD computational study of ultra-narrow TM non-paraxial spatial soliton interactions. IEEE Microw. Wirel. Compon. Lett. 21(5), 228–230 (2011)
Lyu, M., Bokil, V.A., Cheng, Y., Li, F.: Energy stable nodal discontinuous Galerkin methods for nonlinear Maxwell’s equations in multi-dimensions. J. Sci. Comput. 89, 1–42 (2021)
Patel, C.K.N.: Efficient phase-matched harmonic generation in Tellurium with a \({\textrm{CO}}_2\) laser at 10.6 \(\mu\). Phys. Rev. Lett. 15(26), 1027–1030 (1965)
Sørensen, M.P., Webb, G.M., Brio, M., Moloney, J.V.: Kink shape solutions of the Maxwell-Lorentz system. Phys. Rev. E 71(3), 036602 (2005)