Image segmentation is an important task in the domain of computer vision and medical imaging. In natural and medical images, intensity inhomogeneity, i.e. the varying image intensity, occurs often and it poses considerable challenges for image segmentation. In this paper, we propose an efficient variational method for segmenting images with intensity inhomogeneity. The method is inspired by previous works on two-stage segmentation and variational Retinex. Our method consists of two stages. In the first stage, we decouple the image into reflection and illumination parts by solving a convex energy minimization model with either total variation or tight-frame regularisation. In the second stage, we segment the original image by thresholding on the reflection part, and the inhomogeneous intensity is estimated by the smoothly varying illumination part. We adopt a primal dual algorithm to solve the convex model in the first stage, and the convergence is guaranteed. Numerical experiments clearly show that our method is robust and efficient to segment both natural and medical images.

Gibbs phenomenon is the particular manner how a global spectral approximation of a piecewise analytic function behaves at the jump discontinuity. The truncated spectral series has large oscillations near the jump, and the overshoot does not decay as the number of terms in the truncated series increases. There is therefore no convergence in the maximum norm, and convergence in smooth regions away from the discontinuity is also slow. In earlier work, a methodology is proposed to completely overcome this difficulty in the context of spectral collocation methods, resulting in the recovery of exponential accuracy from collocation point values of a piecewise analytic function. In this paper, we extend this methodology to handle spectral collocation methods for functions which are analytic in the open interval but have singularities at end-points. With this extension, we are able to obtain exponential accuracy from collocation point values of such functions. Similar to earlier work, the proof is constructive and uses the Gegenbauer polynomials $C^\lambda_n(x)$. The result implies that the Gibbs phenomenon can be overcome for smooth functions with endpoint singularities.

The entropy solutions of the compressible Euler equations satisfy a minimum principle for the specific entropy. First order schemes such as Godunov-type and Lax-Friedrichs schemes and the second order kinetic schemes also satisfy a discrete minimum entropy principle. In this paper, we show an extension of the positivity-preserving high order schemes for the compressible Euler equations, to enforce the minimum entropy principle for high order finite volume and discontinuous Galerkin (DG) schemes.

The numerical solution of the one-dimensional nonlinear Klein-Gordon equation on an unbounded domain is studied in this paper. Split local absorbing boundary (SLAB) conditions are obtained by the operator splitting method, then the original problem is reduced to an initial boundary value problem on a bounded computational domain, which can be solved by the finite difference method. Several numerical examples are provided to show the advantages and effectiveness of the given method, and some interesting collision behaviors are also observed.

We propose a class of simple and efficient numerical scheme for incompressible fluid equations with coordinate symmetry. By introducing a generalized vorticity-stream formulation, the divergence free constraints are automatically satisfied. In addition, with explicit treatment of the nonlinear terms and local vorticity boundary condition, the Navier–Stokes (MHD, respectively) equation essentially decouples into 2 (4, respectively) scalar equation and thus the scheme is very efficient. Moreover, with proper discretization of the nonlinear terms, the scheme preserves both energy and helicity identities numerically. This is achieved by recasting the nonlinear terms (convection, vorticity stretching, geometric source, Lorentz force and electro-motive force) in terms of Jacobians. This conservative property is valid even in the presence of the pole singularity for axisymmetric flows. The exact conservation of energy and helicity has effectively eliminated excessive numerical viscosity. Numerical examples have demonstrated both accuracy and efficiency of the scheme. Finally, local mesh refinement near the boundary can also be easily incorporated into the scheme without extra cost.