Search

The traditional element free Galerkin (EFG) approach is constructed on variational weak form that the test and trial functions are shape functions of moving least squares (MLS) approximation. In the current paper, we propose a new version of the EFG method based on the shape functions of reproducing kernel particle method (RKPM). In other words, based on the developed approach in Han and Meng (2001) the fractional Tricomi-type equation will be solved using the new technique. The fractional derivative has been introduced in the Caputo’s sense and is approximated by a finite difference plan of order O ( τ 3 − α ) , 1 < α < 2 . We use the EFG-RKPM to discrete the spatial direction. We illustrate some numerical results on non-rectangular domains. The unconditional stability and convergence of the new technique have been proved. Numerical examples display the theoretical results and the efficiency of the proposed approach. Also, the numerical results are compared with the finite element method (FEM) and EFG-MLS procedure. [ABSTRACT FROM AUTHOR]

The element free Galerkin technique is a meshless method based on the variational weak form in which the test and trial functions are the shape functions of moving least squares approximation. Since the shape functions of moving least squares approximation do not have the Kronecker property thus the Dirichlet boundary condition can not be applied directly and also in this case obtaining an error estimate is not simple. The main aim of the current paper is to propose an error estimate for the extracted numerical scheme from the element free Galerkin method. To this end, we select the fractional cable equation with Dirichlet boundary condition. Firstly, we obtain a time-discrete scheme based on a finite difference formula with convergence order O ( τ 1 + min ⁡ { α , β } ) , then we use the meshless element free Galerkin method, to discrete the space direction and obtain a full-discrete scheme. Also, for calculating the appeared integrals over the boundary and the domain of problem the Gauss–Legendre quadrature rule has been used. In the next, we change the main problem with Dirichlet boundary condition to a new problem with Robin boundary condition. Then, we show that the new technique is unconditionally stable and convergent using the energy method. We show convergence orders of the time discrete scheme and the full discrete scheme are O ( τ 1 + min ⁡ { α , β } ) and O ( r p + 1 + τ 1 + min ⁡ { α , β } ) , respectively. So, we can say that the main aim of this paper is as follows, (1) Transferring the main problem with Dirichlet boundary condition (old problem) to a problem with Robin boundary condition (new problem), (2) Showing that with special condition (when σ → + ∞ ) the solution of the new problem is convergent to the solution of the old problem, (3) Obtaining an error estimate for the new problem. Numerical examples confirm the efficiency and accuracy of the proposed scheme. [ABSTRACT FROM AUTHOR]

Three inverse problems of reconstructing the time-dependent, spacewise- dependent and both initial condition and spacewise-dependent heat source in the one-dimensional heat equation are considered. These problems are reformulated by eliminating the unknown functions using some special assumptions concerning the points in space or time as additional measurements. Then direct techniques are proposed to solve the non-classical boundary value problems. For obtaining the robust and stable approximations, Bernstein multi-scaling and B-spline basis functions in the context of the Ritz–Galerkin method are utilized to immediate passage from differential equations to algebraic equations and afterwards, a Newton-type method is used to produce the admissible solution. The numerical convergence and stability are discussed in the test examples to show that the presented schemes provide accurate and acceptable approximations. [ABSTRACT FROM AUTHOR]

Abstract: In the last three decades, Sinc numerical methods have been extensively used for solving differential equations, not only because of their exponential convergence rate, but also due to their desirable behavior toward problems with singularities. This paper illustrates the application of Sinc-collocation and Sinc-Galerkin methods to the approximate solution of the two-dimensional time dependent Schrödinger equation with nonhomogeneous boundary conditions. Some numerical examples are presented and the proposed methods are compared with each other. [Copyright &y& Elsevier]

• A new meshless method is employed to solve fractional modified distributed-order anomalous sub-diffusion equation. • The mentioned method is based on the shape functions of reproducing kernel particle method and Galekrin idea. • The stability and convergence of the developed method have been studied. • The main purpose of this paper is to propose an error analysis to verify that the solutions of penalty method obtained by applying the essential boundary condition are convergent to the solution of main BVP with essential boundary condition. • Some examples are studied that they show the efficiency of the proposed method. • The numerical results confirm the theoretical concepts. In the Galerkin weak form technique based on various kernels that they do not have δ -Kronecker property, in order to apply the essential boundary condition, there are two straight strategies that one of them is the Lagrange multiplier method and another one is the penalty method. In the penalty method the main boundary value problem (BVP) is converted to a new BVP with Robin boundary condition. So, we obtain a new BVP that it must be solved. The main purpose of this paper is to propose an error analysis to verify that the solutions of penalty method obtained by applying the essential boundary condition are convergent to the solution of main BVP with essential boundary condition. For this aim, we select fractional modified distributed-order anomalous sub-diffusion equation. At the first stage, we propose a second-order difference scheme for the temporal variable. The convergence and stability analysis for the time-discrete scheme are proposed. At the second stage, we derive the full-discrete scheme based on the Galerkin weak form and shape functions of reproducing kernel particle method (RKPM) as the mentioned shape functions do not have the δ -Kronecker property. Furthermore, it is shown that when the penalty parameter goes to infinity then the solutions of BVP with Robin boundary condition are convergent to the solutions of BVP based on the essential boundary condition. The proposed examples verify that the present error estimate is true. [ABSTRACT FROM AUTHOR]