A new augmented method is proposed for elliptic interface problems with a piecewise variable coefficient that has a finite jump across a smooth interface. The main motivation is to get not only a second order accurate solution but also a second order accurate gradient from each side of the interface. Key to the new method is introducing the jump in the normal derivative of the solution as an augmented variable and rewriting the interface problem as a new PDE that consists of a leading Laplacian operator plus lower order derivative terms near the interface. In this way, the leading second order derivative jump relations are independent of the jump in the coefficient that appears only in the lower order terms after the scaling. An upwind type discretization is used for the finite difference discretization at the irregular grid points on or near the interface so that the resulting coefficient matrix is an M-matrix. A multigrid solver is used to solve the linear system of equations, and the GMRES iterative method is used to solve the augmented variable. Second order convergence for the solution and the gradient from each side of the interface is proved in this paper. Numerical examples for general elliptic interface problems confirm the theoretical analysis and efficiency of the new method.

In this work, we investigate wave transmission property through a zero index metamaterial (ZIM) waveguide embedded with rectangular dielectric defects. We show that total reflection and total transmission (cloaking) can be achieved by adjusting the geometric sizes and/or permittivities of the defects. Our work provides another possibility of manipulating wave propagation through ZIM in addition to the widely studied dielectric defects with cylindrical geometries.

No abstract uploaded!

A popular and classical method for finding the best rank one approximation of a real tensor is the higher order power method (HOPM). It is known in the literature that the iterative sequence generated by HOPM converges globally, while the convergence rate can be superlinear, linear or sublinear. In this paper, we examine the convergence rate of HOPM in solving the best rank one approximation problem of real tensors. We first show that the iterative sequence of HOPM always converges globally and provide an explicit eventual sublinear convergence rate. The sublinear convergence rate estimate is in terms of the dimension and the order of the underlying tensor space. Then, we examine the concept of nondegenerate singular vector tuples and show that, if the sequence of HOPM converges to a nondegenerate singular vector tuple, then the global convergence rate is R-linear. We show that, for almost all tensors (in the sense of Lebesgue measure), all the singular vector tuples are nondegenerate, and so, the HOPM “typically” exhibits global R-linear convergence rate. Moreover, without any regularity assumption, we establish that the sequence generated by HOPM always converges globally and R-linearly for orthogonally decomposable tensors with order at least 3. We achieved this by showing that each nonzero singular vector tuple of an orthogonally decomposable tensor with order at least 3 is nondegenerate.

No abstract uploaded!