Analysis of higher order difference method for a pseudo-parabolic equation with delay

WOS: 000504461100008In this paper, the author considers the one dimensional initial-boundary problem for a pseudo-parabolic equation with time delay in second spatial derivative. To solve this problem numerically, the author constructs higher order difference method and obtain the error estimate for its solution. Based on the method of energy estimates the fully discrete scheme is shown to be convergent of order four in space and of order two in time. Some numerical examples illustrate the convergence and effectiveness of the numerical method.Scientific and Technological Research Council of Turkey (TUBITAK)Turkiye Bilimsel ve Teknolojik Arastirma Kurumu (TUBITAK) [1059B191700821]This research was fully supported by Scientific and Technological Research Council of Turkey (TUBITAK) [grant number 1059B191700821]. The author would like to thank the Department of Mathematics of the University of Oklahoma for the hospitality during her work on the project, Nikola Petrov and Murad Ozaydin who were instrumental in arranging her stay in Oklahoma. Ilhame Amirali thanks Nikola Petrov for stimulating discussions

A Petrov-Galerkin Finite Element Method for Fractional Convection-Diffusion Equations

In this work, we develop variational formulations of Petrov-Galerkin type for one-dimensional fractional boundary value problems involving either a Riemann-Liouville or Caputo derivative of order $\alpha\in(3/2, 2)$ in the leading term and both convection and potential terms. They arise in the mathematical modeling of asymmetric super-diffusion processes in heterogeneous media. The well-posedness of the formulations and sharp regularity pickup of the variational solutions are established. A novel finite element method is developed, which employs continuous piecewise linear finite elements and "shifted" fractional powers for the trial and test space, respectively. The new approach has a number of distinct features: It allows deriving optimal error estimates in both $L^2(D)$ and $H^1(D)$ norms; and on a uniform mesh, the stiffness matrix of the leading term is diagonal and the resulting linear system is well conditioned. Further, in the Riemann-Liouville case, an enriched FEM is proposed to improve the convergence. Extensive numerical results are presented to verify the theoretical analysis and robustness of the numerical scheme.Comment: 23 p

A RBF partition of unity collocation method based on finite difference for initial-boundary value problems

Meshfree radial basis function (RBF) methods are popular tools used to numerically solve partial differential equations (PDEs). They take advantage of being flexible with respect to geometry, easy to implement in higher dimensions, and can also provide high order convergence. Since one of the main disadvantages of global RBF-based methods is generally the computational cost associated with the solution of large linear systems, in this paper we focus on a localizing RBF partition of unity method (RBF-PUM) based on a finite difference (FD) scheme. Specifically, we propose a new RBF-PUM-FD collocation method, which can successfully be applied to solve time-dependent PDEs. This approach allows to significantly decrease ill-conditioning of traditional RBF-based methods. Moreover, the RBF-PUM-FD scheme results in a sparse matrix system, reducing the computational effort but maintaining at the same time a high level of accuracy. Numerical experiments show performances of our collocation scheme on two benchmark problems, involving unsteady convection-diffusion and pseudo-parabolic equations

Legendre-Gauss-Lobatto collocation method for solving multi-dimensional systems of mixed Volterra-Fredholm integral equations

Integral equations play a crucial role in many scientific and engineering problems, though solving them is often challenging. This paper addresses the solution of multi-dimensional systems of mixed Volterra-Fredholm integral equations (SMVF-IEs) by means of a Legendre-Gauss-Lobatto collocation method. The one-dimensional case is addressed first. Afterwards, the method is extended to two-dimensional linear and nonlinear SMVF-IEs. Several numerical examples reveal the effectiveness of the approach and show its superiority in comparison to other alternative techniques for treating SMVF-IEs

A Petrov--Galerkin Finite Element Method for Fractional Convection-Diffusion Equations

In this work, we develop variational formulations of Petrov--Galerkin type for one-dimensional fractional boundary value problems involving either a Riemann--Liouville or Caputo derivative of order $\alpha\in(3/2, 2)$ in the leading term and both convection and potential terms. They arise in the mathematical modeling of asymmetric superdiffusion processes in heterogeneous media. The well-posedness of the formulations and sharp regularity pickup of the variational solutions are established. A novel finite element method (FEM) is developed, which employs continuous piecewise linear finite elements and “shifted” fractional powers for the trial and test space, respectively. The new approach has a number of distinct features: it allows the derivation of optimal error estimates in both the $L^2(D)$ and $H^1(D)$ norms; and on a uniform mesh, the stiffness matrix of the leading term is diagonal and the resulting linear system is well conditioned. Further, in the Riemann--Liouville case, an enriched FEM is proposed to improve the convergence. Extensive numerical results are presented to verify the theoretical analysis and robustness of the numerical scheme

Abstract book

Welcome at the International Conference on Differential and Difference Equations & Applications 2015. The main aim of this conference is to promote, encourage, cooperate, and bring together researchers in the fields of differential and difference equations. All areas of differential & difference equations will be represented with special emphasis on applications. It will be mathematically enriching and socially exciting event. List of registered participants consists of 169 persons from 45 countries. The five-day scientific program runs from May 18 (Monday) till May 22, 2015 (Friday). It consists of invited lectures (plenary lectures and invited lectures in sections) and contributed talks in the following areas: Ordinary differential equations, Partial differential equations, Numerical methods and applications, other topics

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described