This paper is devoted to developing the Il'in-Allen-Southwell (IAS) parameter-uniform difference scheme on uniform meshes for solving strongly coupled systems of singularly perturbed convection-diffusion equations whose solutions may display boundary and/or interior layers, where strong coupling means that the solution components in the system are coupled together mainly through their first derivatives. By decomposing the coefficient matrix of convection term into the Jordan canonical form, we first construct the IAS scheme for 1-D systems and then extend the scheme to 2-D systems by employing an alternating direction technique. The robustness of the developed IAS scheme is illustrated through a series of numerical examples, including the magnetohydrodynamic duct flow problem with a high Hartmann number. Numerical evidence indicates that the IAS scheme is formally second order accurate in the sense that it is second order convergent when the perturbation parameter $\varepsilon$ is not too small and when the $\varepsilon$ is sufficiently small, the scheme is first order convergent in the discrete maximum norm uniformly in $\varepsilon$.

In this paper, we propose and analyze a mixed $H^1$-conforming finite element method for solving Maxwell's equations in terms of electric field and Lagrange multiplier, where the multiplier is introduced accounting for the divergence constraint. We mainly focus on the case that the physical domain is non-convex and its boundary includes reentrant corners or edges, which may lead the solution of Maxwell's equations to be a non-$H^1$ very weak function and thus causes many numerical difficulties. The proposed method is formulated in the stabilized form by adding an additional mesh-dependent stabilization term to the mixed variational formulation. A general framework of stability and error analysis is established. Specifically, a pair of $H^1$-conforming finite elements, namely the $CP_2$-$P_1$ elements, for electric field and multiplier is studied and its stability and error bounds are also derived. Numerical experiments for source problems as well as eigenvalue problems on the $L$-shaped and cracked domains are presented to illustrate the high performance of the proposed mixed $H^1$-conforming finite element method.

In this paper, we derive the CR Reilly'’s formula and its applications to studying of the fi…rst eigenvalue estimate for CR Dirichlet eigenvalue problem and embedded p-minimal hypersurfaces. In particular, we obtain the …first Dirichlet eigenvalue estimate in a compact pseudo-hermitian manifold with boundary and the fi…rst eigenvalue estimate of the tangential sublaplacian on closed oriented embedded p-minimal hypersurfaces in a closed pseudohermitian manifold of vanishing torsion.

In this paper, we consider the heat ‡flow for p-pseudoharmonic maps from a closed Sasakian manifold into a compact Riemannian manifold. We prove global existence and asymptotic convergence of the solution for the p-pseudoharmonic map heat ‡flow provided that the sectional curvature of the target manifold is nonpositive. Moreover, without the curvature assumption on the target manifold, we obtain global existence and asymptotic convergence of the p-pseudoharmonic map heat fl‡ow as well when its initial p-energy is sufficiently small.

In this paper, we generalize the Cao-Yau'’s gradient estimate for the sum of squares of vector …elds up to higher step under ssumption of the generalized curvature-dimension inequality. With its applications, by deriving a curvature-dimension inequality, we are able to obtain the Li-Yau gradient estimate for the CR heat equation in a closed pseudohermitian manifold of nonvanishing torsion tensors. As consequences, we obtain the Harnack inequality and upper bound estimate for the CR heat kernel. 1.