A new method for two dimensional hyperbolic problems in semi‐unbounded irregular domains

Anewmethod called the full rational differential quadrature method is presented to deal with two dimensional linear and nonlinear hyperbolic problems in semi‐unbounded irregular domains. The spacial and temporal discretizations are both implemented by the rational differential quadrature method (RDQM). The RDQM, proven to be A‐stable (Chen and Tanaka, Comput Mech 28 (2002) 331–338) in the temporal discretization, is much more efficient than the finite difference schemes widely used in earlier works. In addition, the irregular boundary conditions are treated by the direct expansion method(DEM). Numerical experiments show that the present method is of high efficiency and easy to implement. © 2011Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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.     

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.     

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.     

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.     

A Petrov-Galerkin method with quadrature for elliptic boundary value problems

We propose and analyse a fully discrete Petrov–Galerkinmethod with quadrature, for solving second-order, variable coefficient,elliptic boundary value problems on rectangular domains. Inour scheme, the trial space consists of C² splines of degreer 3, the test space consists of C⁰ splines of degree r –2, and we use composite (r – 1)-point Gauss quadrature.We show existence and uniqueness of the approximate solutionand establish optimal order error bounds in H², H¹ and L² norms.     

The versions of the finite element method for problems with boundary layers

We study the uniform approximation of boundary layer functions for , , by the and versions of the finite element method. For the version (with fixed mesh), we prove super-exponential convergence in the range . We also establish, for this version, an overall convergence rate of in the energy norm error which is uniform in , and show that this rate is sharp (up to the term) when robust estimates uniform in are considered. For the version with variable mesh (i.e., the version), we show that exponential convergence, uniform in , is achieved by taking the first element at the boundary layer to be of size . Numerical experiments for a model elliptic singular perturbation problem show good agreement with our convergence estimates, even when few degrees of freedom are used and when is as small as, e.g., . They also illustrate the superiority of the approach over other methods, including a low-order version with optimal ``exponential" mesh refinement. The estimates established in this paper are also applicable in the context of corresponding spectral element methods.

     

Stability of a streamline diffusion finite element method for turning point problems

A one-dimensional singularly perturbed problem with a boundary turning point is considered in this paper. Let V^h be the linear finite element space on a suitable grid . A variant of streamline diffusion finite element method is proved to be almost uniform stable in the sense that the numerical approximation u_h satisfies u-u_h∞C|lnε| inf_vhV^hu-v_h∞, where C is independent with the small diffusion coefficient ε and the mesh . Such stability result is applied to layer-adapted grids to obtain almost ε-uniform second order scheme for turning point problems.     

A numerical method for solving boundary value problems for fractional differential equations

A numerical scheme, based on the Haar wavelet operational matrices of integration for solving linear two-point and multi-point boundary value problems for fractional differential equations is presented. The operational matrices are utilized to reduce the fractional differential equation to system of algebraic equations. Numerical examples are provided to demonstrate the accuracy and efficiency and simplicity of the method.     

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.     

Piecewise reproducing kernel method for singularly perturbed delay initial value problems

A direct application of the reproducing kernel method presented in the previous works cannot yield accurate approximate solutions for singularly perturbed delay differential equations. In this letter, we construct a new numerical method called piecewise reproducing kernel method for singularly perturbed delay initial value problems. Numerical results show that the present method does not share the drawback of standard reproducing kernel method and is an effective method for the considered singularly perturbed delay initial value problems.     

Estimating quadrature errors for an efficient method for quasi-singular boundary integral

The numerical resolution of the boundary integral equations applied to the differential equations of Laplace, Helmholtz and Maxwell requires the handling of quasi-singular integrals with different order of singularity. The numerical approximation of the integral equations of different kinds is made by boundary finite elements. In this paper, we present a complete survey for estimating quadrature errors for the numerical techniques proposed by Huang and Cruse [Q. Huang, T.A. Cruse, Some notes on singular integral techniques in boundary element analysis, Int. J. Numer. Methods Eng. 36 (15) (1993) 2643-2659], to calculate the quasi-singular integrals. To validate the accuracy and efficiency of these techniques and approve our study some numerical examples are presented and discussed.     

Nonlocal boundary value problems for second-order functional differential equations

We study the existence of a solution to a nonlocal boundary value problem for a class of second-order functional differential equations with piecewise constant arguments. The equation studied includes terms depending on the derivative as well as delay arguments in the derivative.     

Numerical solution of the system of second-order boundary value problems using the local radial basis functions based differential quadrature collocation method

In this research, we propose a numerical scheme to solve the system of second-order boundary value problems. In this way, we use the Local Radial Basis Function Differential Quadrature (LRBFDQ) method for approximating the derivative. The LRBFDQ method approximates the derivatives by Radial Basis Functions (RBFs) interpolation using a small set of nodes in the support domain of any node. So the new scheme needs much less computational work than the globally supported RBFs collocation method. We use two techniques presented by Bayona et al. (2011, 2012) [29], [30] to determine the optimal shape parameter. Some examples are presented to demonstrate the accuracy and easy implementation of the new technique. The results of numerical experiments are compared with the analytical solution, finite difference (FD) method and some published methods to confirm the accuracy and efficiency of the new scheme presented in this paper.     

Differential quadrature domain decomposition method for problems on a triangular domain

The differential quadrature method (DQM) has been studied for years and it has been shown by many researchers that the DQM is an attractive numerical method with high efficiency and accuracy. The conventional DQM is mostly effective for one‐dimensional and multidimensional problems with geometrically regular domains. But to deal with problems on a triangular domain, we will meet difficulties. In this article we will study how to solve problems on a triangular domain by using DQM combined with the domain decomposition method (DDM). © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

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     

Asymptotics for mixed Dirichlet–Robin problems in irregular domains

We study the asymptotic behaviour of the solutions of two-dimensional elliptic problems with Robin boundary conditions on the “prefractal” curves approximating the Koch curve type fractals.     

Positive solutions for a class of fractional differential equations with integral boundary conditions

The existence of positive solutions for a class of fractional equations involving the Riemann–Liouville fractional derivative with integral boundary conditions is investigated. By means of the monotone iteration method and some inequalities associated with the Green function, we obtain the existence of a positive solution and establish the iterative sequence for approximating the solution.