In this paper, we propose and analyze a nodal-continuous and $H^1$-conforming finite element method for the numerical computation of Maxwell's equations, with singular solution in a fractional order Sobolev space $H^r(\Omega)$, where $r$ may take any value in the most interesting interval $(0, 1)$. The key feature of the method is that mass-lumping linear finite element $L^2$ projections act on the curl and divergence partial differential operators so that the singular solution can be sought in a setting of $L^2(\Omega)$ space. We shall use the nodal-continuous linear finite elements, enriched with one element bubble in each element, to approximate the singular and non-$H^1$ solution. Discontinuous and nonhomogeneous media are allowed in the method. Some error estimates are given and a number of numerical experiments for source problems as well as eigenvalue problems are presented to illustrate the superior performance of the proposed method.

In this paper, we develop a nodal discontinuous Galerkin method for solving the time-dependent Maxwell’s equations when metamaterials are involved. Both semi- and fully-discrete schemes are constructed. Numerical stability and error estimate are proved for both schemes. Numerical results are presented to demonstrate that the method is not only efficient but also very effective in solving metamaterial Maxwell’s equatons.

Solutions to conservation laws satisfy the monotonicity property: the number of local extrema is a non-increasing function of time, and local maximum/minimum values decrease/increase monotonically in time. This paper investigates this property from a numerical standpoint. We introduce a class of fully discrete in space and time, high order accurate, difference schemes, called generalized monotone schemes. Convergence toward the entropy solution is proven via a new technique of proof, assuming that the initial data has a finite number of extremum values only, and the flux-function is strictly convex. We define discrete paths of extrema by tracking local extremum values in the approximate solution. In the course of the analysis we establish the pointwise convergence of the trace of the solution along a path of extremum. As a corollary, we obtain a proof of convergence for a MUSCL-type scheme being second order accurate away from sonic points and extrema.

Variational models with \ell_1-norm based regularization, in particular total variation (TV) and its variants, have long been known to offer superior image restoration quality, but processing speed remained a bottleneck, preventing their widespread use in the practice of color image processing. In this paper, by extending the grayscale image deblurring algorithm proposed in [Y. Wang, J. Yang, W. Yin, and Y. Zhang, <i>SIAM J. Imaging Sci.</i>, 1 (2008), pp. 248272], we construct a simple and efficient algorithm for multichannel image deblurring and denoising, applicable to both within-channel and cross-channel blurs in the presence of additive Gaussian noise. The algorithm restores an image by minimizing an energy function consisting of an \ell_1-norm fidelity term and a regularization term that can be either TV, weighted TV, or regularization functions based on higher-order derivatives. Specifically, we use a multichannel

This work is devoted to the numerical computation of complex band structure $\k=\k(\omega)\in\mathbb C^3$ for positive $\omega$ of three dimensional isotropic dispersive or non-dispersive photonic crystals from the perspective of structured quadratic eigenvalue problems (QEPs). Our basic strategy is to fix two degrees of freedom in $\k\in\mathbb C^3$ and to view the remaining one as the eigenvalue of a quadratic operator pencil derived from Maxwell's equations. Then Yee's scheme is employed to discretize $\nabla\times$ and $\k\times$ operators in this quadratic operator pencil. Distinct from the others' works which either ignore or directly exploit the Hamiltonian structure of the spectrum of the resulting QEP, we reformulate this QEP into an equivalent $\top$-palindromic QEP to facilitate the use of superior structure-preserving algorithms. Ultimately we rely on the structured Arnoldi algorithm, namely the G$\top$SHIRA algorithm, to compute eigenvalues of a $\top$-skew-Hamiltonian pair which are near or in $[-2,2]$, a much narrower region than the whole positive real axis in the origin problem. Moreover, to accelerate the inner iterations of the G$\top$SHIRA algorithm, we propose the preconditioning technique, making most of the eigenmatrix, which can essentially be seen as the Kronecker product of three discrete Fourier transformation matrices, of the commutative discretized $\partial_x,\partial_y,\partial_z$ operators. The advantage of our method is discussed in detail and corroborated by several numerical results.