In this paper, we develop a second-order asymptotic-preserving and positivity-preserving discontinuous Galerkin (DG) scheme for the Kerr-Debye model. By using the approach first introduced by Zhang and Shu in Ref.~\refcite{ZhSh:04} with an energy estimate and Taylor expansion, the asymptotic-preserving property of the semi-discrete DG methods is proved rigorously. In addition, we propose a class of unconditional positivity-preserving implicit-explicit (IMEX) Runge-Kutta methods for the system of ordinary differential equations arising from the semi-discretization of the Kerr-Debye model. The new IMEX Runge-Kutta methods are based on the modification of the strong-stability-preserving (SSP) implicit Runge-Kutta method and have second-order accuracy. The numerical results validate our analysis.

BDDC (Balancing Domain Decomposition by Constraints) and FETI-DP (Dual-Primal Finite Element Tearing and Interconnecting) algorithms with adaptively enriched coarse spaces are developed and analyzed for second order elliptic problems with high contrast and random coefficients. Among many approaches to form adaptive coarse spaces, we consider an approach using eigenvectors of generalized eigenvalues problems defined on each subdomain interface, seeMandel and Sousedk(2007), Galvis and Efendiev(2010), Spillane etal.(2011), Spillane etal.(2013), Klawonn etal.(2015).

In this paper we propose a finite element time-domain method for modeling the optical black holes (OBHs) coupled with the perfectly matched layer (PML) technique. Stability analysis is carried out for the proposed scheme. Simulations of cylindrical, elliptical and square black holes demonstrate that our method is quite effective in modeling OBHs in time domain. To our best knowledge, this is the first OBHs simulation realized by the finite element time-domain method.

METHODS AND APPLICATIONS OF ANALYSIS. 2000 International Press Vol. 7, No. 3, pp. 591-604, September 2000 012

In this work, we study the unsteady axisymmetric mixed convection boundary layer flow and heat transfer of non-Newtonian power law fluid over a cylinder. Different from most classical works, the temperature dependent variable fluid viscosity and thermal conductivity are taken into account in highly coupled velocity and temperature fields. The motion of the fluid can be modeled by a time-dependent nonlinear parabolic system in cylindrical coordinates, which is solved numerically by using Chebyshev spectral method along with the strong stability preserving (SSP) third order Runge-Kutta time discretization. We apply the numerical solver to problems with different power law indices, viscosity parameter, thermal conductivity parameter and Richardson numbers, and compute up to the steady state. The numerical solver is checked by testing the spectral convergence of the numerical approximation to a smooth exact solution of the PDEs with source terms. Moreover, the combined effects of pertinent physical parameters on the flow and heat transfer characteristics are analyzed in detail.