This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer. 

This paper deals with numerical solution to the multi-term time fractional diffusion equation in a finite domain. An implicit finite difference scheme is established based on Caputo's definition to the fractional derivatives, and the upper and lower bounds to the spectral radius of the coefficient matrix of the difference scheme are estimated, with which the unconditional stability and convergence are proved. The numerical results demonstrate the effectiveness of the theoretical analysis, and the method and technique can also be applied to other kinds of time/space fractional diffusion equations.

The anisotropic perfectly matched layer (APML) defines a continuous vector field outside a rectangle domain and performs the complex coordinate stretching along the vector field. Inspired by [Z. Chen et al., Inverse Probl. Imag., 7, (2013):663--678] and based on the idea of the shortest distance, we propose a new approach to construct the vector field which still allows us to prove the exponential decay of the stretched Green function without the constraint on the thickness of the PML layer. Moreover, by using the reflection argument, we prove the stability of the PML problem in the PML layer and the convergence of the PML method. Numerical experiments are also included.

Based on polyhedral splines, some multivariate splines of different orders with given supports over arbitrary topological meshes are developed. Schemes for choosing suitable families of multivariate splines based on pre-given meshes are discussed. Those multivariate splines with inner knots and boundary knots from the related meshes are used to generate rational spline shapes with related control points. Steps for up to $C^2$-surfaces over the meshes are designed. The relationship among the meshes and their knots, the splines and control points is analyzed. To avoid any unexpected discontinuities and get higher smoothness, a heart-repairing technique to adjust inner knots in the multivariate splines is designed.
With the theory above, bivariate $C^1$-quadratic splines over rectangular meshes are developed. Those bivariate splines are used to generate rational $C^1$-quadratic surfaces over the meshes with related control points and weights. The properties of the surfaces are analyzed. The boundary curves and the corner points and tangent planes, and smooth connecting conditions of different patches are presented. The $C^1$−continuous connection schemes between two patches of the surfaces are presented.

This paper presents a computational method for solving stochastic Ito-Volterra integral equations. First,  Haar wavelets and their properties are employed to derive a general procedure for forming the stochastic operational matrix of Haar wavelets. Then, application of this stochastic operational matrix for solving stochastic Ito-Volterra integral equations is explained. The convergence and error analysis of the proposed method are investigated. Finally, the efficiency of the presented method is confirmed by some examples.

Adaptive numerical methods for solving partial differential equations (PDEs) that control the movement of grid points are called moving mesh methods. In this paper, these methods are examined in the case where a separate PDE, that depends on a monitor function, controls the behavior of the mesh. This results in a system of PDEs: one controlling the mesh and another solving the physical problem that is of interest.  For a class of monitor functions resembling the arc length monitor, a trade off between computational efficiency in solving the moving mesh system and the accuracy level of the solution to the physical PDE is demonstrated. This accuracy is measured in the density of mesh points in the desired portion of the domain where the function has steep gradient. The balance of computational efficiency versus accuracy is illustrated numerically with both the arc length monitor and a monitor that minimizes certain interpolation errors.  Physical solutions with steep gradients in small portions of their domain are considered for both the analysis and the computations.

Image segmentation is a fundamental problem in both image processing and computer vision with numerous applications. In this paper, we propose a two-stage image segmentation scheme based on inexact alternating direction method. Specifically, we first solve the convex variant of the Mumford-Shah model to get the smooth solution, and the segmentation is then obtained  by applying the K-means clustering method to the solution. Some numerical comparisons are arranged to show the effectiveness of our proposed schemes by segmenting many kinds of images such as artificial images, natural images, and brain MRI images.

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multi-material system and explain the exchange of energy among electrons, ions and photons. Their highly nonlinear, strong discontinuous and tightly coupled phenomena always make the numerical solution of such equations extremely challenging. In this paper, we construct two finite volume element schemes both satisfying the discrete conservation property. One of them can well preserve the positivity of analytical solutions, while the other one does not satisfy this property. To fix this defect, two as repair techniques are designed. In addition, as the numerical simulation of 2-D 3-T equations is very time consuming, we also devise a mesh adaptation algorithm to reduce the cost. Numerical results show that these new methods are practical and efficient in solving this kind of problems.