No abstract uploaded!

In this paper we develop a leap-frog-type finite element method for modeling the electromagnetic wave propagation in metamaterials. The metamaterial model equations are represented by integrodifferential Maxwell’s equations, which are quite challenging to analyze and solve in that we have to solve a coupled problem with different partial differential equations given in different material regions. Our method is based on a mixed finite element method using edge elements, which can easily handle the tangential continuity of the electric field. Stability analysis and optimal error estimate are carried out for the proposed scheme. The scheme is implemented and confirmed to obey the proved optimal convergence rate by using a smooth analytical solution. Then the scheme is extended to model wave propagation in heterogeneous media composed of metamaterials and free space, and extensive numerical results (using a rectangular edge element, a triangular edge element, and mixed edge elements) demonstrate the effectiveness of our algorithm for modeling the exotic backward wave propagation phenomenon in metamaterials. To the best of our knowledge, this is the first paper with an exhaustive simulation of backward wave propagation in metamaterials using time-domain finite element methods.

In this article, we present a second-order in time implicit-explicit (IMEX) local discontinuous Galerkin (LDG) method for computing the Cahn-Hilliard equation, which describes the phase separation phenomenon. It is well-known that the Cahn-Hilliard equation has a nonlinear stability property, i.e., the free-energy functional decreases with respect to time. The discretized Cahn-Hilliard system modeled by the IMEX LDG method can inherit the nonlinear stability of the continuous model. We apply a stabilization technique and prove unconditional energy stability of our scheme. Numerical experiments are performed to validate the analysis. Computational efficiency can be significantly enhanced by using this IMEX LDG method with a large time step.

In this paper, we introduce a new method to avoid zero reductions in Gröbner basis computation. We call this method LASyz, which stands for Lineal Algebra to compute Syzygies. LASyz uses exhaustively the information of both principal syzygies and non-trivial syzygies to avoid zero reductions. All computation is done using linear algebra techniques. LASyz is easy to understand and implement. The method does not require to compute Gröbner bases of subsequences of generators incrementally and it imposes no restrictions on the reductions allowed. We provide a complete theoretical foundation for the LASyz method and we describe an algorithm to compute Gröbner bases for zero dimensional ideals based on this foundation. A qualitative comparison with similar algorithms is provided and the performance of the algorithm is illustrated with experimental data.

No abstract uploaded!