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 a novel iterative convolution-thresholding method (ICTM) that is applicable to a range of variational models for image segmentation. A variational model usually minimizes an energy functional consisting of a fidelity term and a regularization term. In the ICTM, the interface between two different segment domains is implicitly represented by their characteristic functions. The fidelity term is usually written as a linear functional of the characteristic functions and the regularized term is approximated by a functional of characteristic functions in terms of heat kernel convolution. This allows us to design an iterative convolution-thresholding method to minimize the approximate energy. The method is simple, efficient and enjoys the energy-decaying property. Numerical experiments show that the method is easy to implement, robust and applicable to various image segmentation models.

No abstract uploaded!

We show that the classical fourth order accurate compact finite difference scheme with high order strong stability preserving time discretizations for convection diffusion problems satisfies a weak monotonicity property, which implies that a simple limiter can enforce the bound- preserving property without losing conservation and high order accuracy. Higher order accurate compact finite difference schemes satisfying the weak monotonicity will also be discussed.

In this paper, we develop a three-dimensional parallel solver using the fifth order high-resolution weighted essentially non-oscillatory (WENO) finite difference scheme to perform extensive simulation for three-dimensional gaseousdetonations. A careful study is conducted for the propagation modes of three-dimensional gaseous detonation wave-front structures in a long square duct. The numerical results indicate that, the instability of detonation, overdrive factor, transverse dimension and initial perturbation are the main factors that influence the appearance of spinning detonation.