A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.

A Novel Accurate Method forMulti-Term Time-Fractional Nonlinear Diffusion Equations in Arbitrary Domains

Anovel accuratemethod is proposed to solve a broad variety of linear and nonlinear(1+1)-dimensional and(2+1)-dimensional multi-term time-fractional partial differential equations with spatial operators of anisotropic diffusivity.For(1+1)-dimensional problems,analytical solutions that satisfy the boundary requirements are derived.Such solutions are numerically calculated using the trigonometric basis approximation for(2+1)-dimensional problems.With the aid of these analytical or numerical approximations,the original problems can be converted into the fractional ordinary differential equations,and solutions to the fractional ordinary differential equations are approximated by modified radial basis functions with time-dependent coefficients.An efficient backward substitution strategy that was previously provided for a single fractional ordinary differential equation is then used to solve the corresponding systems.The straightforward quasilinearization technique is applied to handle nonlinear issues.Numerical experiments demonstrate the suggested algorithm’s superior accuracy and efficiency.

A Novel Method for Solving Time-Dependent 2D Advection-Diffusion-Reaction Equations to Model Transfer in Nonlinear Anisotropic Media

This paper presents a new numerical technique for solving initial and bound-ary value problems with unsteady strongly nonlinear advection diffusion reaction(ADR)equations.The method is based on the use of the radial basis functions(RBF)for the approximation space of the solution.The Crank-Nicolson scheme is used for approximation in time.This results in a sequence of stationary nonlinear ADR equations.The equations are solved sequentially at each time step using the proposed semi-analytical technique based on the RBFs.The approximate solution is sought in the form of the analytical expansion over basis functions and contains free parameters.The basis functions are constructed in such a way that the expansion satisfies the boundary conditions of the problem for any choice of the free parameters.The free parameters are determined by substitution of the expansion in the equation and collocation in the solution domain.In the case of a nonlinear equation,we use the well-known procedure of quasilinearization.This transforms the original equation into a sequence of the linear ones on each time layer.The numerical examples confirm the high accuracy and robustness of the proposed numerical scheme.

A Numerical-Analytical Method for Time-Fractional Dual-Phase-Lag Models of Heat Transfer

The aim of this paper is to present the backward substitution method for solving a class of fractional dual-phase-lag models of heat transfer.The proposed method is based on the Fourier series expansion along the spatial coordinate over the orthonormal basis formed by the eigenfunctions of the corresponding Sturm-Liouville problem.This Fourier expansion of the solution transforms the original fractional par-tial differential equation into a sequence of multi-term fractional ordinary differential equations.These fractional equations are solved by the use of the backward substi-tution method.The numerical examples with temperature-jump boundary condition and parameters of the tissue confirm the high accuracy and efficiency of the proposed numerical scheme.