Multi-domain spectral approach to rational-order fractional derivatives

C. Klein,N. Stoilov
2024-01-09
Abstract:We propose a method to numerically compute fractional derivatives (or the fractional Laplacian) on the whole real line via Riesz fractional integrals. The compactified real line is divided into a number of intervals, thus amounting to a multi-domain approach; after transformations in accordance with the underlying $Z_{q}$ curve ensuring analyticity of the respective integrands, the integrals over the different domains are computed with a Clenshaw-Curtis algorithm. As an example, we consider solitary waves for fractional Korteweg-de Vries equations and compare these to results obtained with a discrete Fourier transform.
Numerical Analysis,Analysis of PDEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to efficiently calculate fractional derivatives (or fractional Laplacian operators) on the real number line. Specifically, the author proposes a multi - domain spectral method based on Riesz fractional integrals for the numerical calculation of fractional derivatives. This method achieves efficient calculation by dividing the compactified real number line into multiple intervals and using the Clenshaw - Curtis algorithm to calculate the integrals on different intervals. ### Background and Motivation Fractional derivatives have attracted great interest in many applications, such as dealing with flow in heterogeneous media, modeling viscoelasticity in biological and medical applications, anomalous diffusion, and financial flows. Unlike classical derivatives, fractional derivatives are non - local, which enables them to incorporate memory effects and non - local effects. Especially in the context of dispersive partial differential equations (PDEs), non - local effects are particularly important in water wave theory. For example, the Whitham equation is a one - dimensional non - local equation with the dispersion relation of the exact water wave equation. ### Specific Problems The specific problems in the paper are: 1. **How to efficiently calculate fractional derivatives on the real number line**: - The author proposes a multi - domain spectral method. By dividing the real number line into multiple intervals and applying coordinate transformations on each interval, the integrand function is made smooth, and then the Clenshaw - Curtis algorithm is used for numerical integration. - This method combines the high - precision and locality advantages of spectral methods and can allocate higher resolution where needed, such as in regions where the function is highly oscillatory. 2. **Verify the effectiveness of the method**: - The author tests this method through known examples and applies it to the solitary wave solutions of the fractional Korteweg - de Vries (KdV) equation. - They use the Newton - Krylov method to construct these solitary wave solutions and compare them with the results obtained previously using the discrete Fourier transform (DFT) method. ### Mathematical Expressions The fractional derivative \( D^\alpha u(x) \) can be defined through the Riesz fractional integral as: \[ D^\alpha u(x) = \frac{1}{2\Gamma(1-\alpha) \sin\left(\frac{(1-\alpha)\pi}{2}\right)} \frac{\partial}{\partial x} \int_{-\infty}^{\infty} \frac{\text{sgn}(x - y) u(y)}{|x - y|^\alpha} dy \] where \( 0 < \alpha < 1 \), \( \Gamma \) is the gamma function, and \( \text{sgn} \) is the sign function. ### Application Examples The paper uses a specific example to verify the effectiveness of the method, that is, the solitary wave solutions of the fractional KdV equation: \[ u_t - D^\alpha u_x + uu_x = 0 \] The form of the solitary wave solution is \( u(x, t) = Q_c(x - ct) \), where \( c > 0 \) is a real constant. These solitary waves satisfy the equation: \[ cQ_c + D^\alpha Q_c - \frac{1}{2} Q_c^2 = 0 \] ### Conclusion The method proposed in the paper shows high precision and efficiency in calculating fractional derivatives, especially when dealing with non - local effects and highly oscillatory regions. By comparing with the DFT method, the effectiveness of this method is verified.