In this paper, based on the basis composed of two sets of splines with distinct local supports, cubic spline quasi-interpolating operators are reviewed on nonuniform type-2 triangulation. The variation diminishing operator is defined by discrete linear functionals based on a fixed number of triangular mesh-points, which can reproduce any polynomial of nearly best degrees. And by means of the modulus of continuity, the estimation of the operator approximating a real sufficiently smooth function is reviewed as well. Moreover, the derivatives of the nearly optimal variation diminishing operator can approximate that of the real sufficiently smooth function uniformly over quasi-uniform type-2 triangulation. And then the convergence results are worked out.

In this paper, a fitted Numerov method is constructed for a class of singularly perturbed one-dimensional parabolic partial differential equations with a small negative shift in the temporal variable.  Similar boundary value problems are associated with a furnace used to process a metal sheet in control theory. Here, the study focuses on the effect of shift on the boundary layer behavior of the solution via finite difference approach.  When the shift parameter is smaller than the perturbation parameter, the shifted term is expanded in Taylor series and an exponentially fitted tridiagonal finite difference scheme is developed. The proposed finite difference scheme is unconditionally stable.  When the shift parameter is larger than the perturbation parameter, a special type of mesh is used for the temporal variable so that the shift lies on the nodal points and an exponentially fitted scheme is developed.  This scheme is also unconditionally stable.  The applicability of the proposed methods is demonstrated by means of two examples.

New numerical techniques are presented for the solution of the two-dimensional time fractional evolution equation in the unit square. In these methods, Galerkin finite element is used for the spatial discretization, and, for the time stepping, new alternating direction implicit (ADI) method based on the backward Euler method combined with the first order convolution quadrature approximating the integral term are considered. The ADI Galerkin finite element method is proved to be convergent in time and in the $L^2$ norm in space. The convergence order is $\mathcal{O}$($k$|ln $k$|+$h^r$), where $k$ is the temporal grid size and $h$ is spatial grid size in the $x$ and $y$ directions, respectively. Numerical results are presented to support our theoretical analysis.

pde2path is a free and easy to use Matlab continuation/bifurcation package for elliptic systems of PDEs with arbitrary many components, on general two dimensional domains, and with rather general boundary conditions. The package is based on the FEM of the Matlab pdetoolbox, and is explained by a number of examples, including Bratu's problem, the Schnakenberg model, Rayleigh-Bénard convection, and von Karman plate equations. These serve as templates to study new problems, for which the user has to provide, via Matlab function files, a description of the geometry, the boundary conditions, the coefficients of the PDE, and a rough initial guess of a solution. The basic algorithm is a one parameter arclength-continuation with optional bifurcation detection and branch-switching. Stability calculations, error control and mesh-handling, and some elementary timeintegration for the associated parabolic problem are also supported. The continuation, branch-switching, plotting etc are performed via Matlab command-line function calls guided by the AUTO style. The software can be downloaded from www.staff.uni-oldenburg.de/hannes.ue ker/pde2path, where also an online documentation of the software is provided such that in this paper we focus more on the mathematics and the example systems.

In this paper, a finite element method is proposed to investigate multiple solutions of the Navier-Stokes equations for an unsteady, laminar, incompressible flow in a porous expanding channel. Dual or triple solutions for the fixed values of the wall suction Reynolds number $R$ and the expansion ratio $α$ are obtained numerically. The computed multiple solutions for the symmetric flow are validated by comparing them with approximate analytic solutions obtained by the similarity transformation and homotopy analysis method. Unlike previous works, our method deals with the Navier-Stokes equations directly and thus has no similarity and other restrictions as in previous works. Finally we use the method to study multiple solutions for three cases of the asymmetric flow (which has not been studied before using the similarity-type techniques).