Abstract

We develop a new matrix-based approach to the numerical solution of partial differential equations (PDE) and apply it to the numerical solution of partial fractional differential equations (PFDE). The proposed method is to discretize a given PFDE as a Sylvester Equation, and parameterize the integral surface using matrix algebra. The combination of these two notions results in an algorithm which can solve a general class of PFDE efficiently and accurately by means of an O(n3) algorithm for solving the Sylvester Matrix Equation (over an m×n grid with m∼n). The proposed parametrization of the integral surface allows for the solution with the more general Robin boundary conditions, and allows for high-order approximations to derivative boundary conditions. To achieve our ends, we also develop a new matrix-based approximation to fractional order derivatives. The proposed method is demonstrated by the numerical solution of the fractional diffusion equation with fractional derivatives in both the temporal and spatial directions.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

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.