We present a hybrid asymptotic/numerical method for the accurate computation of single- and double-layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called “geometrically induced stiffness,” meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.

A variety of problems in computational physics and engineering require the convolu- tion of the heat kernel (a Gaussian) with either discrete sources, densities supported on boundaries, or continuous volume distributions. We present a unified fast Gauss transform for this purpose in two dimensions, making use of an adaptive quad-tree discretization on a unit square which is assumed to contain all sources. Our implementation permits either free-space or periodic boundary conditions to be imposed, and is efficient for any choice of variance in the Gaussian.

In this paper, we propose an unfitted Nitsche's method for computing wave modes in topological materials. The proposed method is based on Nitsche's technique to study the performance-enhanced topological materials which have strongly heterogeneous structures (e.g., the refractive index is piecewise constant with high contrasts). For periodic bulk materials, we use Floquet-Bloch theory and solve an eigenvalue problem on a torus with unfitted meshes. For the materials with a line defect, a sufficiently large domain with zero boundary conditions is used to compute the localized eigenfunctions corresponding to the edge modes. The interfaces are handled by Nitsche's method on an unfitted uniform mesh. We prove the proposed methods converge optimally, and present numerical examples to validate the theoretical results and demonstrate the capability of simulating topological materials.

In this paper, we propose an unfitted Nitsche's method to compute the band structures of phononic crystal with impurities of general geometry. The proposed method does not require the background mesh to fit the interfaces of impurities, and thus avoids the expensive cost of generating body-fitted meshes and simplifies the inclusion of interface conditions in the formulation. The quasi-periodic boundary conditions are handled by the Floquet-Bloch transform, which converts the computation of band structures into an eigenvalue problem with periodic boundary conditions. More importantly, we show the well-posedness of the proposed method using a delicate argument based on the trace inequality, and further prove the convergence by the Babuška-Osborn theory. We achieve the optimal convergence rate at the presence of the impurities of general geometry. We confirm the theoretical results by two numerical examples, and show the capability of the proposed methods for computing the band structures without fitting the interfaces of impurities.

Photonic graphene, a photonic crystal with honeycomb structures, has been intensively studied in both theoretical and applied fields. Similar to graphene which admits Dirac Fermions and topological edge states, photonic graphene supports novel and subtle propagating modes (edge modes) of electromagnetic waves. These modes have wide applications in many optical systems. In this paper, we propose a novel gradient recovery method based on Bloch theory for the computation of topological edge modes in the honeycomb structure. Compared to standard finite element methods, this method provides higher order accuracy with the help of gradient recovery technique. This high order accuracy is highly desired for constructing the propagating electromagnetic modes in applications. We analyze the accuracy and prove the superconvergence of this method. Numerical examples are presented to show the efficiency by computing the edge mode for the P-symmetry and C-symmetry breaking cases in honeycomb structures.