We consider the problem of reconstructing unknown inclusions inside a thermal conductor from boundary measurements, which arises from active thermography and is formulated as an inverse boundary value problem for the heat equation. In our previous works, we proposed a sampling-type method for reconstructing the boundary of the unknown inclusion and gave its rigorous mathematical justification. In this paper, we continue our previous works and provide a further investigation of the reconstruction method from both the theoretical and numerical points of view. First, we analyze the solvability of the Neumann-to-Dirichlet map gap equation and establish a relation of its solution to the Green function of an interior transmission problem for the inclusion. This naturally provides a way of computing this Green function from the Neumann-to-Dirichlet map. Our new findings reveal the essence of the reconstruction method. A convergence result for noisy measurement data is also proved. Second, based on the heat layer potential argument, we perform a numerical implementation of the reconstruction method for the homogeneous inclusion case. Numerical results are presented to show the efficiency and stability of the proposed method.

We propose second order accurate discontinuous Galerkin (DG) schemes which satisfy a strict maximum principle for general nonlinear convection-diffusion equations on unstructured triangular meshes. Motivated by genuinely high order maximum-principle-satisfying DG schemes for hyperbolic conservation laws, we prove that under suitable time step restriction for forward Euler time stepping, for general nonlinear convection-diffusion equations, the same scaling limiter coupled with second order DG methods preserves the physical bounds indicated by the initial condition while maintaining uniform second order accuracy. Similar to the purely convection cases, the limiters are mass conservative and easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. Following the idea in earlier work, we extend the schemes to two-dimensional convection-diffusion equations on triangular meshes. There are no geometric constraints on the mesh such as angle acuteness. Numerical results including incompressible Navier-Stokes equations are presented to validate and demonstrate the effectiveness of the numerical methods.

Solitary waves bifurcated from edges of Bloch bands in two-dimensional periodic media are determined both analytically and numerically in the context of a two-dimensional nonlinear Schrödinger equation with a periodic potential. Using multiscale perturbation methods, the envelope equations of solitary waves near Bloch bands are analytically derived. These envelope equations reveal that solitary waves can bifurcate from edges of Bloch bands under either focusing or defocusing nonlinearity, depending on the signs of the second-order dispersion coefficients at the edge points. Interestingly, at edge points with two linearly independent Bloch modes, the envelope equations lead to a host of solitary wave structures, including reduced-symmetry solitons, dipole-array solitons, vortex-cell solitons, and so on—many of which have not been reported before to our knowledge. It is also shown analytically that the centers of envelope solutions can be positioned at only four possible locations at or between potential peaks. Numerically, families of these solitary waves are directly computed both near and far away from the band edges. Near the band edges, the numerical solutions spread over many lattice sites, and they fully agree with the analytical solutions obtained from the envelope equations. Far away from the band edges, solitary waves are strongly localized, with intensity and phase profiles characteristic of individual families.

In this paper, we discuss high order finite difference weighted essentially non-oscillatory (WENO) schemes, coupled with total variation diminishing (TVD) Runge-Kutta (RK) temporal integration, for solving the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology. Since the solutions to this system are non-negative, we discuss a positivity-preserving limiter without compromising accuracy. Analysis is performed to justify the maintanance of third order spatial / temporal accuracy when the limiters are applied to a third order finite difference scheme and third order TVD-RK time discretization for solving this model. Numerical results are also provided to demonstrate these methods up to fifth order accuracy.

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.