This article focuses on numerically studying the eigenstructure behavior of generalized eigenvalue problems (GEPs) arising in three dimensional (3D) source-free Maxwell's equations with magnetoelectric coupling effects which model 3D reciprocal chiral media. It is challenging to solve such a large-scale GEP efficiently. We combine the null-space free method with the inexact shift-invert residual Arnoldi method and MINRES linear solver to solve the GEP with a matrix dimension as large as 5,308,416. The eigenstructure is heavily determined by the chirality parameter $\gamma$. We show that all the eigenvalues are real and finite for a small chirality $\gamma$. For a critical value $\gamma = \gamma^*$, the GEP has $2 \times 2$ Jordan blocks at infinity eigenvalues. Numerical results demonstrate that when $\gamma$ increases from $\gamma^*$, the $2 \times 2$ Jordan block will first split into a complex conjugate eigenpair, then rapidly collide \tr{with} the real axis and bifurcate into \tr{positive (resonance) and negative eigenvalues with modulus} smaller than the other existing positive eigenvalues. The resonance band also exhibits an anticrossing interaction. Moreover, the electric and magnetic fields of the resonance modes are localized inside the structure, with only a slight amount of field leaking into the background (dielectric) material.

No abstract uploaded!

Japan’s lunar probe satellite “Kaguya”, China’s lunar probe satellite “Chang’ e-1 lunar probe” left the earth on September 14th, 2007 and October 24th, 2007, began the exploration to the moon. The launching of the China’s lunar satellite is the proud to the global Chinese. However, at the same time as we are proud of the event, we must admit that the pictures published which were sent back by Japanese are clearer and more beautiful than Chinese. While we are saying that the Japan has a better technology, maybe we also need to say that the Japanese can hold better chance than Chinese when taking the pictures. The following two pictures were sent back by the satellite. Picture one is sent back by Japan’s satellite, Picture two is sent back by China’s. It can be seen that there is a long way for China if we want to become stronger at science technology. We also hope that we can have some contribution in the prosperity of China now and in the future.

Wave propagation problems arise in a wide range of applications. The energy conserving property is one of the guiding principles for numerical algorithms, in order to minimize the phase or shape errors after long time integration. In this paper, we develop and analyze a local discontinuous Galerkin (LDG) method for solving the wave equation. We prove optimal error estimates, superconvergence toward a particular projection of the exact solution, and the energy conserving property for the semi-discrete formulation. The analysis is extended to the fully discrete LDG scheme, with the centered second-order time discretization (the leap-frog scheme). Our numerical experiments demonstrate optimal rates of convergence and superconvergence. We also show that the shape of the solution, after long time integration, is well preserved due to the energy conserving property.

No abstract uploaded!