2,242 research outputs found
A M\"untz-Collocation spectral method for weakly singular volterra integral equations
In this paper we propose and analyze a fractional Jacobi-collocation spectral
method for the second kind Volterra integral equations (VIEs) with weakly
singular kernel . First we develop a family of fractional
Jacobi polynomials, along with basic approximation results for some weighted
projection and interpolation operators defined in suitable weighted Sobolev
spaces. Then we construct an efficient fractional Jacobi-collocation spectral
method for the VIEs using the zeros of the new developed fractional Jacobi
polynomial. A detailed convergence analysis is carried out to derive error
estimates of the numerical solution in both - and weighted
-norms. The main novelty of the paper is that the proposed method is
highly efficient for typical solutions that VIEs usually possess. Precisely, it
is proved that the exponential convergence rate can be achieved for solutions
which are smooth after the variable change for a
suitable real number . Finally a series of numerical examples are
presented to demonstrate the efficiency of the method
Regularity theory and high order numerical methods for the (1D)-fractional Laplacian
This paper presents regularity results and associated high-order numerical methods for one-dimensional Fractional-Laplacian boundary-value problems. On the basis of a factorization of solutions as a product of a certain edge-singular weight times a ``regular´´ unknown, a characterization of the regularity of solutions is obtained in terms of the smoothness of the corresponding right-hand sides. In particular, for right-hand sides which are analytic in a Bernstein Ellipse, analyticity in the same Bernstein Ellipse is obtained for the ``regular´´ unknown. Moreover, a sharp Sobolev regularity result is presented which completely characterizes the co-domain of the Fractional-Laplacian operator in terms of certain weighted Sobolev spaces introduced in (Babu{s}ka and Guo, SIAM J. Numer. Anal. 2002). The present theoretical treatment relies on a full eigendecomposition for a certain weighted integral operator in terms of the Gegenbauer polynomial basis. The proposed Gegenbauer-based Nystr"om numerical method for the Fractional-Laplacian Dirichlet problem, further, is significantly more accurate and efficient than other algorithms considered previously. The sharp error estimates presented in this paper indicate that the proposed algorithm is spectrally accurate, with convergence rates that only depend on the smoothness of the right-hand side. In particular, convergence is exponentially fast (resp. faster than any power of the mesh-size) for analytic (resp. infinitely smooth) right-hand sides. The properties of the algorithm are illustrated with a variety of numerical results.Fil: Acosta, Gabriel. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Matemática; ArgentinaFil: Borthagaray, Juan Pablo. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Matemática; ArgentinaFil: Bruno, Oscar Ricardo. California Institute Of Technology; Estados UnidosFil: Maas, MartÃn Daniel. Consejo Nacional de Investigaciónes CientÃficas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de AstronomÃa y FÃsica del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de AstronomÃa y FÃsica del Espacio; Argentin
- …