Abstract
We consider the discretization in time if a fractional order diffusion equation. The approximation is based on a further development of the approach of using Laplace transformation to represent the solution as a contour integral, evaluated to high accuracy by quadrature. This technique reduces the problem to a finite set of elliptic equations with complex coefficients, which may be solved in parallel. Three different methods, using 2N + 1 quadrature points, are discussed. The first has an error of order O (e(-cN)) away from t = 0, whereas the second and third methods are uniformly accurate of order O (e(-C root N)). Unlike the first and second methods, the third method does not use the Laplace transform of the forcing term. The basic analysis of the time discretizaiton takes place in a Banach space setting and uses a resolvent esitmate for the associated elliptic operator. The methods are combined with finite element discretization in the spatial variable to yield fully discrete methods.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.