In loving memory of my beloved miniature dachshund Maddie (16 March 2002 – 16 March 2020). We consider nonlocal differential equations with convolution coefficients of the form

\[{-}M\Big(\big(a*(g\circ |u|)\big)(1)\Big)u''(t)=\lambda f\big(t,u(t)\big),\quad t\in(0,1), \]

$g$

$M$

$M$

$f$

Following the seminal paper by Bourgain, Brezis, and Mironescu, we focus on the asymptotic behaviour of some nonlocal functionals that, for each $u\in L^2(\mathbb {R}^N)$, are defined as the double integrals of weighted, squared difference quotients of $u$. Given a family of weights $\{\rho _{\varepsilon} \}$, $\varepsilon \in (0,\,1)$, we devise sufficient and necessary conditions on $\{\rho _{\varepsilon} \}$ for the associated nonlocal functionals to converge as $\varepsilon \to 0$ to a variant of the Dirichlet integral. Finally, some comparison between our result and the existing literature is provided.

Introducing a pair-parameter matrix Mittag–Leffler function, we study the uniqueness and Hyers–Ulam stability to a new fractional nonlinear partial integro-differential equation with variable coefficients and a mixed boundary condition using Banach’s contractive principle as well as Babenko’s approach in a Banach space. These investigations have serious applications since uniqueness and stability analysis are essential topics in various research fields. The techniques used also work for different types of differential equations with initial or boundary conditions, as well as integral equations. Moreover, we present a Python code to compute approximate values of our newly established pair-parameter matrix Mittag–Leffler functions, which extend the multivariate Mittag–Leffler function. A few examples are given to show applications of the key results obtained.

In this work, we study the existence of solutions of nonlinear fractional coupled system of $\varphi $-Hilfer type in the frame of Banach spaces. We improve a property of a measure of noncompactness in a suitably selected Banach space. Darbo’s fixed point theorem is applied to obtain a new existence result. Finally, the validity of our result is illustrated through an example.

On all Bergman–Besov Hilbert spaces on the unit disk, we find self-adjoint weighted shift operators that are differential operators of half-order whose commutators are the identity, thereby obtaining uncertainty relations in these spaces. We also obtain joint average uncertainty relations for pairs of commuting tuples of operators on the same spaces defined on the unit ball. We further identify functions that yield equality in some uncertainty inequalities.

We apply capacities to explore the space–time fractional dissipative equation: (0.1)

$$ \begin{align} \left\{\begin{aligned} &\partial^{\beta}_{t}u(t,x)=-\nu(-\Delta)^{\alpha/2}u(t,x)+f(t,x),\quad (t,x)\in\mathbb R^{1+n}_{+},\\ &u(0,x)=\varphi(x),\ x\in\mathbb R^{n}, \end{aligned}\right. \end{align} $$

where $\alpha>n$ and $\beta \in (0,1)$. In this paper, we focus on the regularity and the blow-up set of mild solutions to (0.1). First, we establish the Strichartz-type estimates for the homogeneous term $R_{\alpha ,\beta }(\varphi )$ and inhomogeneous term $G_{\alpha ,\beta }(g)$, respectively. Second, we obtain some space–time estimates for $G_{\alpha ,\beta }(g).$ Based on these estimates, we prove that the continuity of $R_{\alpha ,\beta }(\varphi )(t,x)$ and the Hölder continuity of $G_{\alpha ,\beta }(g)(t,x)$ on $\mathbb {R}^{1+n}_+,$ which implies a Moser–Trudinger-type estimate for $G_{\alpha ,\beta }.$ Then, for a newly introduced $L^{q}_{t}L^p_{x}$-capacity related to the space–time fractional dissipative operator $\partial ^{\beta }_{t}+(-\Delta )^{\alpha /2},$ we perform the geometric-measure-theoretic analysis and establish its basic properties. Especially, we estimate the capacity of fractional parabolic balls in $\mathbb {R}^{1+n}_+$ by using the Strichartz estimates and the Moser–Trudinger-type estimate for $G_{\alpha ,\beta }.$ A strong-type estimate of the $L^{q}_{t}L^p_{x}$-capacity and an embedding of Lorentz spaces are also derived. Based on these results, especially the Strichartz-type estimates and the $L^{q}_{t}L^p_{x}$-capacity of fractional parabolic balls, we deduce the size, i.e., the Hausdorff dimension, of the blow-up set of solutions to (0.1).

In this paper, we concern with a backward problem for a nonlinear time fractional wave equation in a bounded domain. By applying the properties of Mittag-Leffler functions and the method of eigenvalue expansion, we establish some results about the existence and uniqueness of the mild solutions of the proposed problem based on the compact technique. Due to the ill-posedness of backward problem in the sense of Hadamard, a general filter regularization method is utilized to approximate the solution and further we prove the convergence rate for the regularized solutions.

In this study, we investigate the intial value problem (IVP) for a time-fractional fourth-order equation with nonlinear source terms. More specifically, we consider the time-fractional biharmonic with exponential nonlinearity and the time-fractional Cahn–Hilliard equation. By using the Fourier transform concept, the generalized formula for the mild solution as well as the smoothing effects of resolvent operators are proved. For the IVP associated with the first one, by using the Orlicz space with the function $\Xi (z)={\textrm {e}}^{|z|^{p}}-1$ and some embeddings between it and the usual Lebesgue spaces, we prove that the solution is a global-in-time solution or it shall blow up in a finite time if the initial value is regular. In the case of singular initial data, the local-in-time/global-in-time existence and uniqueness are derived. Also, the regularity of the mild solution is investigated. For the IVP associated with the second one, some modifications to the generalized formula are made to deal with the nonlinear term. We also establish some important estimates for the derivatives of resolvent operators, they are the basis for using the Picard sequence to prove the local-in-time existence of the solution.

We consider nonlocal equations of the general form

\begin{equation} \left(a*u''\right)(\cdot)+\lambda f\big(\cdot,u(\cdot)\big)=0.\nonumber \end{equation}

The fractional derivatives include nonlocal information and thus their calculation requires huge storage and computational cost for long time simulations. We present an efficient and high-order accurate numerical formula to speed up the evaluation of the Caputo fractional derivative based on the L2-1σ formula proposed in [A. Alikhanov, J. Comput. Phys., 280 (2015), pp. 424-438], and employing the sum-of-exponentials approximation to the kernel function appeared in the Caputo fractional derivative. Both theoretically and numerically, we prove that while applied to solving time fractional diffusion equations, our scheme not only has unconditional stability and high accuracy but also reduces the storage and computational cost.

We propose a hybrid spectral element method for fractional two-point boundary value problem (FBVPs) involving both Caputo and Riemann-Liouville (RL) fractional derivatives. We first formulate these FBVPs as a second kind Volterra integral equation (VIEs) with weakly singular kernel, following a similar procedure in [16]. We then design a hybrid spectral element method with generalized Jacobi functions and Legendre polynomials as basis functions. The use of generalized Jacobi functions allow us to deal with the usual singularity of solutions at t = 0. We establish the existence and uniqueness of the numerical solution, and derive a hptype error estimates under L2(I)-norm for the transformed VIEs. Numerical results are provided to show the effectiveness of the proposed methods.

The computational work and storage of numerically solving the time fractional PDEs are generally huge for the traditional direct methods since they require total memory and work, where NT and NS represent the total number of time steps and grid points in space, respectively. To overcome this difficulty, we present an efficient algorithm for the evaluation of the Caputo fractional derivative of order α∈(0,1). The algorithm is based on an efficient sum-of-exponentials (SOE) approximation for the kernel t–1–α on the interval [Δt, T] with a uniform absolute error ε. We give the theoretical analysis to show that the number of exponentials Nexp needed is of order for T≫1 or for TH1 for fixed accuracy ε. The resulting algorithm requires only storage and work when numerically solving the time fractional PDEs. Furthermore, we also give the stability and error analysis of the new scheme, and present several numerical examples to demonstrate the performance of our scheme.

In this paper, we first apply cosine radial basis function neural networks to solve the fractional differential equations with initial value problems or boundary value problems. In the examples, we successfully obtained the numerical solutions for the fractional Riccati equations and fractional Langevin equations. The computer graphics and numerical solutions show that this method is very effective.

The second order weighted and shifted Grünwald difference (WSGD) operators are developed in [Tian, Zhou and Deng, Math. Comput., 84 (2015), pp. 1703–1727] to solve space fractional partial differential equations. Along this direction, we further design a new family of second order WSGD operators; by properly choosing the weighted parameters, they can be effectively used to discretize space (Riemann-Liouville) fractional derivatives. Based on the new second order WSGD operators, we derive a family of difference schemes for the space fractional advection diffusion equation. By von Neumann stability analysis, it is proved that the obtained schemes are unconditionally stable. Finally, extensive numerical experiments are performed to demonstrate the performance of the schemes and confirm the convergence orders.

In this paper the fractional Euler-Lagrange equation is considered. The fractional equation with the left and right Caputo derivatives of order α ∈ (0,1] is transformed into its corresponding integral form. Next, we present a numerical solution of the integral form of the considered equation. On the basis of numerical results, the convergence of the proposed method is determined. Examples of numerical solutions of this equation are shown in the final part of this paper.

An efficient high order numerical method is presented to solve the mobile-immobile advection-dispersion model with the Coimbra time variable-order fractional derivative, which is used to simulate solute transport in watershed catchments and rivers. On establishing an efficient recursive algorithm based on the properties of Jacobi polynomials to approximate the Coimbra variable-order fractional derivative operator, we use spectral collocation method with both temporal and spatial discretisation to solve the time variable-order fractional mobile-immobile advection-dispersion model. Numerical examples then illustrate the effectiveness and high order convergence of our approach.

In this paper we present the basic tools of a fractional function theory in higher dimensions by means of a fractional correspondence to the Weyl relations via fractional Riemann–Liouville derivatives. A Fischer decomposition, Almansi decomposition, fractional Euler and Gamma operators, monogenic projection, and basic fractional homogeneous powers are constructed. Moreover, we establish the fractional Cauchy–Kovalevskaya extension (FCK extension) theorem for fractional monogenic functions defined on ℝd . Based on this extension principle, fractional Fueter polynomials, forming a basis of the space of fractional spherical monogenics, i.e. fractional homogeneous polynomials, are introduced. We study the connection between the FCK extension of functions of the form xPl and the classical Gegenbauer polynomials. Finally, we present an example of an FCK extension.

In this paper, we introduce the definition of a convex real valued function $f$ defined on the set of integers, $\mathbb{Z}$ . We prove that $f$ is convex on $\mathbb{Z}$ if and only if ${{\Delta }^{2}}f\,\ge \,0$ on $\mathbb{Z}$ . As a first application of this new concept, we state and prove discrete Hermite–Hadamard inequality using the basics of discrete calculus (i.e., the calculus on $\mathbb{Z}$ ). Second, we state and prove the discrete fractional Hermite–Hadamard inequality using the basics of discrete fractional calculus. We close the paper by defining the convexity of a real valued function on any time scale.

High-dimensional two-sided space fractional diffusion equations with variable diffusion coefficients are discussed. The problems can be solved by an implicit finite difference scheme that is proven to be uniquely solvable, unconditionally stable and first-order convergent in the infinity norm. A nonsingular multilevel circulant pre-conditoner is proposed to accelerate the convergence rate of the Krylov subspace linear system solver efficiently. The preconditoned matrix for fast convergence is a sum of the identity matrix, a matrix with small norm, and a matrix with low rank under certain conditions. Moreover, the preconditioner is practical, with an O(N logN) operation cost and O(N) memory requirement. Illustrative numerical examples are also presented.

A finite difference method which is second-order accurate in time and in space is proposed for two-dimensional fractional percolation equations. Using the Fourier transform, a general approximation for the mixed fractional derivatives is analyzed. An approach based on the classical Crank-Nicolson scheme combined with the Richardson extrapolation is used to obtain temporally and spatially second-order accurate numerical estimates. Consistency, stability and convergence of the method are established. Numerical experiments illustrating the effectiveness of the theoretical analysis are provided.