In this paper we present a rigorous derivation of the material parameters for both the cylinder and rectangle cloaking structures. Numerical results using these material parameters are presented to demonstrate the cloaking effect.

In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson–Nernst–Planck (PNP) equations, which have found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.

In this follow-up work, we extend the discontinuous Galerkin (DG) methods previously developed on rectangular meshes \cite{previous} to triangular meshes. The DG schemes in \cite{previous} are both optimally convergent and energy conserving. However, as we shall see in the numerical results section, the DG schemes on triangular meshes only have suboptimal convergence rate. We prove the energy conservation and an error estimate for the semi-discrete schemes. The stability of the fully discrete scheme is proved and its error estimate is stated. We present extensive numerical results with convergence consistent of our error estimate, and simulations of wave propagation in Drude metamaterials to demonstrate the flexibility of triangular meshes.

In this paper, we compute the stationary states of the multicomponent phase-field crystal model by formulating it as a block constrained minimization problem. The original infinite-dimensional non-convex minimization problem is approximated by a finite-dimensional constrained non-convex minimization problem after an appropriate spatial discretization. To efficiently solve the above optimization problem, we propose a so-called adaptive block Bregman proximal gradient (AB-BPG) algorithm that fully exploits the problem's block structure. The proposed method updates each order parameter alternatively, and the update order of blocks can be chosen in a deterministic or random manner. Besides, we choose the step size by developing a practical linear search approach such that the generated sequence either keeps energy dissipation or has a controllable subsequence with energy dissipation. The convergence property of the proposed method is established without the requirement of global Lipschitz continuity of the derivative of the bulk energy part by using the Bregman divergence. The numerical results on computing stationary ordered structures in binary, ternary, and quinary component coupled-mode Swift-Hohenberg models have shown a significant acceleration over many existing methods.

In this paper, a structure-preserving algorithm is developed for the computation of a semi-stabilizing solution of a Generalized Algebraic Riccati Equation (GARE). The semi-stabilizing solution of GAREs has been used to characterize the solvability of the (J, J ′ )-spectral factorization problem in control theory for general rational matrices which may have poles and zeros on the extended imaginary axis. The main difficulty in solving such a GARE lies in the fact that its associated Hamiltonian/skew-Hamiltonian pencil has eigenvalues on the extended imaginary axis. Consequently, it is not clear which eigenspace of the associated Hamiltonian/skew-Hamiltonian pencil can characterize the desired semi-stabilizing solution. That is, it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute. Hence, the well-known generalized eigenspace approach for the classical algebraic Riccati equations cannot be employed directly. The proposed algorithm consists of a structure-preserving doubling algorithm (SDA) and a postprocessing procedure to determine the desired eigenvectors and principal vectors corresponding to the purely imaginary and infinite eigenvalues. Under mild assumptions, linear convergence of rate 1/2 for the SDA is proved. Numerical experiments illustrate that the proposed algorithm performs efficiently and reliably.