Reproducing kernel method for singularly perturbed turning point problems having twin boundary layers

A numerical method is proposed for solving singularly perturbed turning point problems exhibiting twin boundary layers based on the reproducing kernel method (RKM). The original problem is reduced to two boundary layers problems and a regular domain problem. The regular domain problem is solved by using the RKM. Two boundary layers problems are treated by combining the method of stretching variable and the RKM. The boundary conditions at transition points are obtained by using the continuity of the approximate solution and its first derivatives at these points. Two numerical examples are provided to illustrate the effectiveness of the present method. The results compared with other methods show that the present method can provide very accurate approximate solutions.     

On discontinuous Galerkin finite element method for singularly perturbed delay differential equations

A nonsymmetric discontinuous Galerkin FEM with interior penalties has been applied to one-dimensional singularly perturbed problem with a constant negative shift. Using higher order polynomials on Shishkin-type layer-adapted meshes, a robust convergence has been proved in the corresponding energy norm. Numerical experiments support theoretical findings.     

The characteristic finite difference streamline diffusion method for convection-dominated diffusion problems

In this paper, we consider the characteristic finite difference streamline diffusion method for two-dimensional convection-dominated diffusion problems. The scheme is combined the method of characteristics with the finite difference streamline diffusion (FDSD) method to create the characteristic FDSD (C-FDSD) procedures. Stability analysis and error estimate of the C-FDSD method are deduced. The scheme not only realizes the purpose of lowering the time-truncation error, using larger time step for solving the convection-dominated diffusion problems, but also keeps the favorable stability and high precision of the FDSD method. Finally, numerical experiments are presented to illustrate the availability of the scheme.     

Convergence analysis of an upwind mixed element method for advection diffusion problems

We consider a upwinding mixed element method for a system of first order partial differential equations resulting from the mixed formulation of a general advection diffusion problem. The system can be used to model the transport of a contaminant carried by a flow. We use the lowest order Raviart-Thomas mixed finite element space. We show the first order convergence both for concentration and concentration flux in L²(Ω).     

The Lie-group shooting method for solving nonlinear singularly perturbed boundary value problems

A new computational method for solving the second-order nonlinear singularly perturbed boundary value problems (SPBVPs) is provided in this paper. In order to overcome a highly singular behavior very near to the boundary as being not easy to treat by numerical method, we adopt a coordinate transformation from an x-domain to a t-domain via a rescaling technique, which can reduce the singularity within the boundary layer. Then, we construct a Lie-group shooting method (LGSM) to search a missing initial condition through the finding of a suitable value of a parameter r ∈ [0, 1]. Moreover, we can derive a closed-form formula to express the initial condition in terms of r, which can be determined properly by an accurate matching to the right-boundary condition. Numerical examples are examined, showing that the present approach is highly efficient and accurate.     

The hp finite element method for singularly perturbed problems in nonsmooth domains

We consider the numerical approximation of singularly perturbed elliptic boundary value problems over nonsmooth domains. We use a decomposition of the solution that contains a smooth part, a corner layer part and a boundary layer part. Explicit guidelines for choosing mesh‐degree combinations are given that yield finite element spaces with robust approximation properties. In particular, we construct an hp finite element space that approximates all components uniformly, at a near exponential rate. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 63–89, 1999     

Piecewise shooting reproducing kernel method for linear singularly perturbed boundary value problems

In this letter, a new numerical method is proposed for solving second order linear singularly perturbed boundary value problems with left layers. Firstly a piecewise reproducing kernel method is proposed for second order linear singularly perturbed initial value problems. By combining the method and the shooting method, an effective numerical method is then proposed for solving second order linear singularly perturbed boundary value problems. Two numerical examples are used to show the effectiveness of the present method.     

Global superconvergence of finite element methods for parabolic inverse problems

In this article we transform a large class of parabolic inverse problems into a nonclassical parabolic equation whose coefficients consist of trace type functionals of the solution and its derivatives subject to some initial and boundary conditions. For this nonclassical problem, we study finite element methods and present an immediate analysis for global superconvergence for these problems, on basis of which we obtain a posteriori error estimators. This research was supported in part by the Shahid Beheshti University, the National Basic Research Program of China (2007CB814906), the National Natural Science Foundation of China (10471103 and 10771158), Social Science Foundation of the Ministry of Education of China (Numerical methods for convertible bonds, 06JA630047), Tianjin Natural Science Foundation (07JCYBJC14300).     

A mixed finite element method for fourth order eigenvalue problems

A mixed finite element method for approximating eigenpairs of IV order elliptic eigenvalue problems with Dirichlet boundary conditions has been given. The method can be applied to the vibration analysis of anisotropic/orthotropic/isotropic/biharmonic plates. Computer implementation procedures for this mixed method are given along with the results of numerical experiments.     

一类具有转向点的三阶微分方程的边值问题的解的高阶渐近展开

利用上下解方法,研究如下一类具有转向点的三阶微分方程的边值问题{ε~2y″′=f(t,y,ε)y″ g(t,y,ε)y′ h(t,y,ε),a相似文献   

On a rational differential quadrature method in irregular domains for problems with boundary layers

A rational differential quadrature method in irregular domains (RDQMID) is investigated to deal with a kind of singularly perturbed problems with boundary layers. Through a transformation, the boundary layer, which may be not straight, is transformed into a segment of a line parallel to one of the Cartesian axes. The rational differential quadrature method (RDQM) is applied to discretize the governing equation. Finally, a direct expansion method of the boundary conditions (DEMBC) is raised to deal with the boundary conditions. Numerical experiments show that RDQMID is of high accuracy, efficiency and easy to programme.     

Convergence analysis of a finite volume method via a new nonconforming finite element method

The article is devoted to the study of convergence properties of a Finite Volume Method (FVM) using Voronoi boxes for discretization. The approach is based on the construction of a new nonconforming Finite Element Method (FEM), such that the system of linear equations coincides completely with that for the FVM. Thus, by proving convergence properties of the FEM, we obtain similar ones of the FVM. In this article, the investigations are restricted to the Poisson equation. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14:213–231, 1998     

A mixed finite element method for the unilateral contact problem in elasticity

In this paper, we provide a new mixed finite element approximation of the varia-tional inequality resulting from the unilateral contact problem in elasticity. We use the continuous piecewise P2-P1 finite element to approximate the displacement field and the normal stress component on the contact region. Optimal convergence rates are obtained under the reasonable regularity hypotheses. Numerical example verifies our results.     

Preconditioners for the p-version of the Galerkin method for a coupled finite element/boundary element system

We propose and analyze efficient preconditioners for solving systems of equations arising from the p-version for the finite element/boundary element coupling. The first preconditioner amounts to a block Jacobi method, whereas the second one is partly given by diagonal scaling. We use the generalized minimum residual method for the solution of the linear system. For our first preconditioner, the number of iterations of the GMRES necessary to obtain a given accuracy grows like log² p, where p is the polynomial degree of the ansatz functions. The second preconditioner, which is more easily implemented, leads to a number of iterations that behave like p log³ p. Computational results are presented to support this theory. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 47–61, 1998     

The Euler-Galerkin finite element method for nonlocal diffusion problems with a p-Laplace-type operator

This paper is devoted to the study of the finite element method for a class of non-linear nonlocal diffusion problems associated with p-Laplace-type operator. Using the Euler–Galerkin finite element method, the convergence and a priori error estimates for the semi-discrete as well as fully-discrete formulations are established.     

A finite difference analysis of a streamline diffusion method on a Shishkin mesh

We consider streamline diffusion finite element methods applied to a singularly perturbed convection–diffusion two‐point boundary value problem whose solution has a single boundary layer. To analyse the convergence of these methods, we rewrite them as finite difference schemes. We first consider arbitrary meshes, then, in analysing the scheme on a Shishkin mesh, we consider two formulations on the fine part of the mesh: the usual streamline diffusion upwinding and the standard Galerkin method. The error estimates are given in the discrete L _∞ norm; in particular we give the first analysis that shows precisely how the error depends on the user-chosen parameter _τ0 specifying the mesh. When _τ0 is too small, the error becomes O(1), but for _τ0 above a certain threshold value, the error is small and increases either linearly or quadratically as a function of . Numerical tests support our theoretical results. This revised version was published online in August 2006 with corrections to the Cover Date.     

A novel method for solving a class of singularly perturbed boundary value problems based on reproducing kernel method

In this paper, a novel method is presented for solving a class of singularly perturbed boundary value problems. Firstly the original problem is reformulated as a new boundary value problem whose solution does not change rapidly via a proper transformation; then the reproducing kernel method is employed to solve the boundary value new problem. Numerical results show that the present method can provide very accurate analytical approximate solutions.