We review some recent developments in numerical algorithms to solve the time-dependent Maxwell equations for systems with spatially varying permittivity and permeability. We show that the Suzuki product-formula approach can be used to construct a family of unconditionally stable algorithms, the conventional Yee algorithm, and two new variants of the Yee algorithm that do not require the use of ...

For solving time-dependent convection-dominated partial differential equations (PDEs), which arise frequently in computational physics, high order numerical methods, including finite difference, finite volume, finite element and spectral methods, have been undergoing rapid developments over the past decades. In this article we give a brief survey of two selected classes of high order methods, n...

where x = (x1, . . . , xd) ∈ IR , t > 0. HJ equations arise in many practical areas such as differential games, mathematical finance, image enhancement and front propagation. It is well known that solutions of (1) are Lipschitz continuous but derivatives can become discontinuous even if the initial data is smooth. There is a close relation between HJ equations and hyperbolic conservation laws. ...

In this work, we introduce new finite-difference shock-capturing central schemes on staggered grids. Staggered schemes may have better resolution of the corresponding unstaggered schemes of the same order. They are based on high order non oscillatory reconstruction (ENO or WENO), and a suitable ODE solver for the computation of the integral of the flux. Although they suffer from a more severe s...

SUMMARY Metamaterials are generally defined as a class of artificial effective media which macroscopically exhibit extraordinary electromagnetic properties that may not be found in nature, and are composed of periodically structured dielectric, or magnetic, or metallic materials. This paper reviews recently developed electromagnetic modeling methods of metama-tericals and their inherent basic i...

In the last decades, more or less complex physically-based hydrological models, have been developed that solve the shallow water equations or their approximations using various numerical methods. Model users may not necessarily know the different hypothesis lying behind these development and simplifications, and it might therefore be difficult to judge if a code is well adapted to their objecti...

For many elliptic PDE problems, finite-difference and finite-element methods are the techniques of choice. In a finite-difference approach, a solution uk on a set of discrete gridpoints 1, . . . , k is searched for. The discretized partial differential equation and boundary conditions create linear relationships between the different values of uk. In the finite-element method, the solution is e...

Article history: Received 30 March 2013 Received in revised form 21 July 2013 Accepted 5 February 2014 Available online 14 February 2014

A third-order Energy Stable Weighted Essentially Non–Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter (AIAA 2008–2876, 2008) was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO s...