We introduce an algorithm for solving two-sided space-fractional partial differential equations. The space-fractional derivatives we consider here are left-handed and right-handed Riemann–Liouville fractional derivatives which are expressed by using Hadamard finite-part integrals. We approximate the Hadamard finite-part integrals by using piecewise quadratic interpolation polynomials and obtain a numerical approximation of the space-fractional derivative with convergence order

A new high-order finite difference scheme to approximate the Caputo fractional derivative $\frac{1}{2} \big ( \, _{0}^{C}D^{\alpha}_{t}f(t_{k})+ \, _{0}^{C}D^{\alpha}_{t}f(t_{k-1}) \big ), k=1, 2, \dots, N, $ with the convergence order $O(\Delta t^{4-\alpha}), \, \alpha\in(1,2)$ is obtained when $f^{\prime \prime \prime} (t_{0})=0$, where $\Delta t$ denotes the time step size. Based on this scheme we introduce a finite difference method for solving fractional diffusion wave equation with the convergence order $O(\Delta t^{4-\alpha} + h^2)$, where $h$ denotes the space step size. Numerical examples are given to show that the numerical results are consistent with the theoretical results.

In this paper, a class of discrete wave equations with Dirichlet boundary conditions are obtained by using the center-difference method. For any positive integers m and T, when the existence of time mT-periodic solutions is considered, a strongly indefinite discrete system needs to be established. By using a variant generalized weak linking theorem, a non-resonant superlinear (or superquadratic) result is obtained and the Ambrosetti-Rabinowitz condition is improved. Such a method cannot be used for the corresponding continuous wave equations or the continuous Hamiltonian systems; however, it is valid for some general discrete Hamiltonian systems.

This article discusses the oscillatory behaviour of the differential equation of mixed type.

We consider the optimal convergence rates of the semidiscrete finite element approximations for solving linear space-fractional partial differential equations by using the regularity results for the fractional elliptic problems obtained recently by Jin et al. \cite{jinlazpasrun} and Ervin et al. \cite{ervheuroo}. The error estimates are proved by using two approaches. One approach is to apply the duality argument in Johnson \cite{joh} for the heat equation to consider the error estimates for the linear space-fractional partial differential equations. This argument allows us to obtain the optimal convergence rates under the minimal regularity assumptions for the solution. Another approach is to use the approximate solution operators of the corresponding fractional elliptic problems. This argument can be extended to consider more general linear space-fractional partial differential equations. Numerical examples are given to show that the numerical results are consistent with the theoretical results.

The study analysed unsteady Rivlin-Ericksen fluid and irreversibility of exponentially temperature dependent variable viscosity of hydromagnetic two-step exothermic chemical reactive flow along the channel axis with walls convective cooling. The non-Newtonian Hele-Shaw flow of Rivlin-Erickson fluid is driven by bimolecular chemical kinetic and unvarying pressure gradient. The reactive fluid is induced by periodic changes in magnetic field and time. The Newtons law of cooling is satisfied by the constant heat coolant convection exchange at the wall surfaces with the neighboring regime. The dimensionless non-Newtonian reactive fluid equations are numerically solved using a convergent and consistence semi-implicit finite difference technique which are confirmed stable. The response of the reactive fluid flow to variational increase in the values of some entrenched fluid parameters in the momentum and energy balance equations are obtained. A satisfying equations for the ratio of irreversibility, entropy generation and Bejan number are solved with the results presented graphically and discussed quantitatively. From the study, it was obtained that the thermal criticality conditions with the right combination of thermo-fluid parameters, the thermal runaway can be prevented. Also, the entropy generation can minimize by at low dissipation rate and viscosity.

We consider finite difference methods for solving nonlinear fractional differential equations in the Caputo fractional derivative sense with non-uniform meshes. Under the assumption that the Caputo derivative of the solution of the fractional differential equation is suitably smooth, Li et al. \lq \lq Finite difference methods with non-uniform meshes for nonlinear fractional differential equations\rq\rq, Journal of Computational Physics, 316(2016), 614-631, obtained the error estimates of finite difference methods with non-uniform meshes. However the Caputo derivative of the solution of the fractional differential equation in general has a weak singularity near the initial time. In this paper, we obtain the error estimates of finite difference methods with non-uniform meshes when the Caputo fractional derivative of the solution of the fractional differential equation has lower smoothness. The convergence result shows clearly how the regularity of the Caputo fractional derivative of the solution affect the order of convergence of the finite difference methods. Numerical results are presented that confirm the sharpness of the error analysis.

We introduce a high-order numerical method for solving nonlinear fractional differential equation with non-uniform meshes. We first transform the fractional nonlinear differential equation into the equivalent Volterra integral equation. Then we approximate the integral by using the quadratic interpolation polynomials. On the first subinterval $[t_{0}, t_{1}]$, we approximate the integral with the quadratic interpolation polynomials defined on the nodes $t_{0}, t_{1}, t_{2}$ and in the other subinterval $[t_{j}, t_{j+1}], j=1, 2, \dots N-1$, we approximate the integral with the quadratic interpolation polynomials defined on the nodes $t_{j-1}, t_{j}, t_{j+1}$. A high-order numerical method is obtained. Then we apply this numerical method with the non-uniform meshes with the step size $\tau_{j}= t_{j+1}- t_{j}= (j+1) \mu$ where $\mu= \frac{2T}{N (N+1)}$. Numerical results show that this method with the non-uniform meshes has the higher convergence order than the standard numerical methods obtained by using the rectangle and the trapzoid rules with the same non-uniform meshes.

In this paper, we introduce a numerical method for solving the time-space fractional telegraph equations. The numerical method is based on a quadrature formula approach and a stability condition for the numerical method is obtained. Two numerical examples are given and the stability regions are plotted.

We introduce a modified L1 scheme for solving time fractional partial differential equations and obtain error estimates for smooth and nonsmooth initial data in both homogeneous and inhomogeneous cases. Jin \et (2016, An analysis of the L1 scheme for the subdiffusion equation with nonsmooth data, IMA J. of Numer. Anal., 36, 197-221) established an $O(k)$ convergence rate for the L1 scheme for smooth and nonsmooth initial data for the homogeneous problem, where $k$ denotes the time step size. We show that the modified L1 scheme has convergence rate $O(k^{2-\alpha}), 0< \alpha <1$ for smooth and nonsmooth initial data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the numerical results are consistent with the theoretical results.