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.

In this paper, we propose a numerical method to solve isotropic elliptic equations on point cloud by generalizing the point integral method. The idea of the point integral method is to approximate the differential operators by integral operators and discretize the corresponding integral equation on point cloud. The key step is to get the integral approximation. In this paper, with proper kernel function, we get an integral approximation for the elliptic operators with isotropic coefficients. Moreover, the integral approximation has been proved to keep the coercivity of the original elliptic operator. The convergence of the point integral method is also proved.

We show that the gradient mapping of the squared norm of FischerBurmeister function is globally Lipschitz continuous and semismooth, which provides a theoretical basis for solving nonlinear second-order cone complementarity problems via the conjugate gradient method and the semismooth Newtons method.

In this paper, we study the ill posed Perona-Malik equation of image processing [14], and the regularized P-M model i.e. C-model proposed by Catte et al. [4]. The authors present the convex compound of these two models in the form of the system of partial differential equations. The weak solution for the equations is proved in detail. The additive operator splitting (AOS) algorithm for the proposed model is also given. Finally, we show some numeric experiment results on images.

We present a mathematical analysis of the asymptotic preserving scheme proposed in [M. Lemou and L. Mieussens, SIAM J. Sci. Comput., 31, pp. 334-368, 2008] for linear transport equations in kinetic and diffusive regimes. We prove that the scheme is uniformly stable and accurate with respect to the mean free path of the particles. This property is satisfied under an explicitly given CFL condition. This condition tends to a parabolic CFL condition for small mean free paths, and is close to a convection CFL condition for large mean free paths. Ou r analysis is based on very simple energy estimates.