Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary

S. Bandyopadhyay; M. Chhetri; B. B. Delgado; N. Mavinga; R. Pardo; S. Bandyopadhyay; M. Chhetri; B. B. Delgado; N. Mavinga; R. Pardo

doi:10.3934/era.2022107

Electronic Research Archive

2022, Volume 30, Issue 6: 2121-2137. doi: 10.3934/era.2022107

Previous Article Next Article

Research article Special Issues

Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary

1.
UNC Greensboro, Greensboro, NC 27402, USA
2.
Universidad Autónoma de Aguascalientes, 20131 Aguascalientes, Mexico
3.
Swarthmore College, Swarthmore, PA 19081, USA
4.
Universidad Complutense de Madrid, 28040 Madrid, Spain

Received: 14 October 2021 Revised: 29 December 2021 Accepted: 16 January 2022 Published: 14 April 2022

This paper deals with the existence of weak solutions for semilinear elliptic equation with nonlinearity on the boundary. We establish the existence of a maximal and a minimal weak solution between an ordered pair of sub- and supersolution for both monotone and nonmonotone nonlinearities. We use iteration argument when the nonlinearity is monotone. For the nonmonotone case, we utilize the surjectivity of a pseudomonotone and coercive operator, Zorn's lemma and a version of Kato's inequality.

Keywords:

Citation: S. Bandyopadhyay, M. Chhetri, B. B. Delgado, N. Mavinga, R. Pardo. Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary[J]. Electronic Research Archive, 2022, 30(6): 2121-2137. doi: 10.3934/era.2022107

Related Papers:

[1]	Yang Song, Beiyan Yang, Jimin Wang . Stability analysis and security control of nonlinear singular semi-Markov jump systems. Electronic Research Archive, 2025, 33(1): 1-25. doi: 10.3934/era.2025001
[2]	Dandan Zuo, Wansheng Wang, Lulu Zhang, Jing Han, Ling Chen . Non-fragile sampled-data control for synchronizing Markov jump Lur'e systems with time-variant delay. Electronic Research Archive, 2024, 32(7): 4632-4658. doi: 10.3934/era.2024211
[3]	Hui Zhao, Hanhong Cui, Yuan Zhao, Xuewu Dai . A neural network-based adaptive fault-tolerant cooperation control for multiple trains with unknown parameters. Electronic Research Archive, 2025, 33(6): 3931-3949. doi: 10.3934/era.2025174
[4]	Ramalingam Sakthivel, Palanisamy Selvaraj, Oh-Min Kwon, Seong-Gon Choi, Rathinasamy Sakthivel . Robust memory control design for semi-Markovian jump systems with cyber attacks. Electronic Research Archive, 2023, 31(12): 7496-7510. doi: 10.3934/era.2023378
[5]	Chao Ma, Hang Gao, Wei Wu . Adaptive learning nonsynchronous control of nonlinear hidden Markov jump systems with limited mode information. Electronic Research Archive, 2023, 31(11): 6746-6762. doi: 10.3934/era.2023340
[6]	Bin Wang . Random periodic sequence of globally mean-square exponentially stable discrete-time stochastic genetic regulatory networks with discrete spatial diffusions. Electronic Research Archive, 2023, 31(6): 3097-3122. doi: 10.3934/era.2023157
[7]	Meng Gao, Xiaohui Ai . A stochastic Gilpin-Ayala nonautonomous competition model driven by mean-reverting OU process with finite Markov chain and Lévy jumps. Electronic Research Archive, 2024, 32(3): 1873-1900. doi: 10.3934/era.2024086
[8]	Xinling Li, Xueli Qin, Zhiwei Wan, Weipeng Tai . Chaos synchronization of stochastic time-delay Lur'e systems: An asynchronous and adaptive event-triggered control approach. Electronic Research Archive, 2023, 31(9): 5589-5608. doi: 10.3934/era.2023284
[9]	Yueying Liu, Mengping Sun, Zhen Wang, Xiangyun Lin, Cuihua Zhang . Nash equilibrium strategies for non-zero-sum differential games of SDEs with time-varying coefficient and infinite Markov jumps. Electronic Research Archive, 2025, 33(4): 2525-2542. doi: 10.3934/era.2025112
[10]	Shuang Liu, Tianwei Xu, Qingyun Wang . Effect analysis of pinning and impulsive selection for finite-time synchronization of delayed complex-valued neural networks. Electronic Research Archive, 2025, 33(3): 1792-1811. doi: 10.3934/era.2025081

Abstract

1. Introduction

Lots of physical phenomena can be expressed by the FC equation, including, inter alia, dissipative and dispersive partial differential equations (PDEs). In this paper, we consider the FC equation

$\begin{eqnarray} \frac{\partial \phi(t,s)}{\partial t} = -\mu_{0}\ _{0}C_{t}^{1-\alpha_{1}}\phi(t,s)+\ _{0}C_{t}^{1-\alpha_{2}}\frac{\partial^2 \phi(t,s)}{\partial s^2} +f(t,s) ,0 \le s\le 1,0 \le t \le T, \end{eqnarray}$

(1.1)

$\begin{eqnarray} \phi(0,s) = 0, \phi(1,s) = 0, s\in [0,T], \end{eqnarray}$

(1.2)

$\begin{eqnarray} \phi(t,0) = \varphi(t), t\in [0,1] \end{eqnarray}$

(1.3)

where $\mu_{0}\in R, 0 \ < \alpha_{1}, \alpha_{2} < 1$ are constants. There are some definitions of fractional derivatives, such as the Caputo type, Riemann-Liouville type and so on. In the following, we adopt the Caputo type time fractional-order partial derivative as

$\begin{eqnarray} _{0}C_{t}^{\alpha}\phi(t) = \frac{1}{\Gamma(1-\alpha)}\int_{0}^{t}\frac{\phi'(t)}{(t-\tau)^{\alpha}}d\tau, \end{eqnarray}$

(1.4)

and $\Gamma(\alpha)$ is the $\Gamma$ function.

In ^[1], a scheme combining the finite difference method in the time direction and a spectral method in the space direction was proposed. In ^[2], two implicit compact difference schemes for the FC equation were studied, this scheme was proved to be stable, and the convergence order $O(\tau+h^4)$ was given. In ^[3], a two-dimensional FC equation was solved by orthogonal spline collocation (OSC) methods for space discretization and finite difference method for time, which was proved to be unconditionally stable. In^[4], the FC equation with two time Riemann-Liouville derivatives was solved by an explicit numerical method; and the accuracy, stability and convergence of this method were studied. In ^[5], FC equation with two fractional time derivatives were considered, and two new implicit numerical methods for the FC equation were proposed, respectively. The stability and convergence of these methods were also investigated. In ^[6], nonlinear FC equation was solved by a two-grid algorithm with the finite element (FE) method. A time second-order fully discrete two-grid FE scheme and the space direction were approximated. In ^[7], the discrete Crank-Nicolson (CN) finite element method was obtained by the finite difference in time and the finite element in space to approximate the FC equation, the stability and error estimate were analyzed in detail and the optimal convergence rate was obtained. In ^[8], the FC equation involving two integro-differential operators was solved by semi-discrete finite difference approximation, and the scheme was proved unconditionally stable. In reference ^[9], numerical integration with the reproducing kernel gradient smoothing integration are constructed. In reference ^[10], recursive moving least squares (MLS) approximation was constructed.

Like the above methods to solve the FC equation by finite difference approach or finite element method, the time direction and space direction were solved separatively. In the following, we presented the BRIM to solve the time direction and space direction of FC equation at the same time. Lagrange interpolation has been presented by mathematician Lagrange to fitting data to be a certain function. When the number $n$ increases, there are Runge phenomenon that the interpolation result deviates from the original function. In order to avoid the Runge phenomenon, among them, barycentric interpolation was developed in 1960s to overcome it. In recent years, linear rational interpolation (LRI) was proposed by Floater ^[14,15,16] and error of linear rational interpolation ^[11,12,13] is also proved. The barycentric interpolation collocation method (BICM) has developed by Wang et al.^[25,26] and the algorithm of BICM has used for linear/non-linear problems ^[27,28]. In recent research, Volterra integro-differential equation (VIDE) ^[17,21], heat equation (HE) ^[18], biharmonic equation (BE) ^[19], telegraph equation (TE) ^[20], generalized Poisson equations ^[22], fractional reaction-diffusion equation ^[23] and KPP equation ^[24] have been studied by the linear barycentric rational interpolation method (LBRIM) and their convergence rate are also proved.

In this paper, BRIM has been used to solve the FC equation. As the fractional derivative is the nonlocal operator, the spectral method is developed to solve the FC equation and the coefficient matrix is the full matrix. The fractional derivative of the FC equation is changed to nonsingular integral by the order of density function plus one. New Gauss formula is constructed to compute it simply and matrix equation of discrete FC equation is obtained by the unknown function replaced by barycentric rational interpolation basis function. Then, the convergence rate of BRIM is proved.

2. Matrix equation of FC equation

As there is singularity in Eq (1.1), the numerical methods cannot get high accuracy, by fractional integration to second part of (1.1) to overcome the difficulty of singularity. We get

$\begin{align} & _{0}C_{t}^{\alpha}\phi(t,s) \\ & = \frac{1}{\Gamma(\xi-\alpha)}\int_{0}^{t}\frac{\partial^{\xi} \phi(\tau,s)}{\partial \tau^{\xi}}\frac{d\tau}{(\tau-t)^{\alpha+1-\xi}} \\ & = \frac{1}{(\xi-\alpha)\Gamma(\xi-\alpha)}\left[\frac{\partial^{\xi} \phi(0,s)}{\partial t^{\xi}}t^{\xi-\alpha}+\int_{0}^{t}\frac{\partial^{\xi+1} \phi(\tau,s)}{\partial \tau^{\xi+1}}\frac{d\tau}{(t-\tau)^{\alpha-\xi}}\right] \\ & = \Gamma_{\alpha}^{\xi}\left[\frac{\partial^{\xi} \phi(0,s)}{\partial t^{\xi}}t^{\xi-\alpha}+\int_{0}^{t}\frac{\partial^{\xi+1} \phi(\tau,s)}{\partial \tau^{\xi+1}}\frac{d\tau}{(t-\tau)^{\alpha-\xi}}\right], \end{align}$

(2.1)

where $\Gamma_{\alpha}^{\xi} = \frac{1}{(\xi-\alpha)\Gamma(\xi-\alpha)}$ .

Combining (2.1) and (1.1), we have

$\begin{eqnarray} & \frac{\partial \phi}{\partial t}+\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\left[\frac{\partial^{\xi} \phi(0,s)}{\partial t^{\xi}}t^{\xi-\alpha_{1}}+\int_{0}^{t}\frac{\partial^{\xi+1} \phi(\tau,s)}{\partial \tau^{\xi+1}}\frac{d\tau}{(t-\tau)^{\alpha_{1}-\xi}}\right] \\& = \Gamma_{\alpha_{2}}^{\xi}\left[\frac{\partial^{\xi+2} \phi(0,s)}{\partial t^{\xi}\partial s^2 }s^{\xi-\alpha_{2}}+\int_{0}^{t}\frac{\partial^{\xi+3} \phi(\tau,s)}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(s-\tau)^{\alpha_{2}-\xi}}\right]+f(t,s). \end{eqnarray}$

(2.2)

In the following, we give the discrete formula of FC equation and to get the matrix equation from BRIM.

Let

$\begin{eqnarray} \phi(t,s) = \sum\limits_{j = 1}^{m}R_{j}(t)\phi_{j}(s) \end{eqnarray}$

(2.3)

where

$\phi(t_{i},s) = \phi_{i}(s),i = 1,2,\cdots,m$

and

$\begin{equation} R_{j}(t) = \frac{ \frac{\lambda_{j}}{t-t_{j}}}{ \sum\limits_{k = 1}^{n} \frac{\lambda_{k}}{t-t_{k}}} \end{equation}$

(2.4)

where

$\lambda_{k} = \sum\limits_{j\in J_{k}}(-1)^j\prod\limits_{i = j,j \neq k}^{j+d_{t}}\frac{1}{t_{k}-t_{i}},\ \ \ J_{k} = \{j\in \{ 0,1,\cdots, l-d_{t}\}:k-d_{t}\leq j \leq k\}$

is the basis function ^[18]. Taking (2.3) into Eq (2.2),

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}R'_{j}(t)\phi_{j}(s) +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi_{j}(s)t^{\xi-\alpha_{1}}+\int_{0}^{t} \phi_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t-\tau)^{\alpha_{1}-\xi}}\right] \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi^{(2)}_{j}(s)t^{\xi-\alpha_{2}} +\int_{0}^{t} \phi^{(2)}_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t-\tau)^{\alpha_{2}-\xi}}\right]+f(t,s). \end{eqnarray}$

(2.5)

By taking $0 = t_{1} < t_{2} < \cdots < t_{m} = T, a = s_{1} < s_{2} < \cdots < s_{n} = b$ with $h_{t} = T/m, h_{s} = (b-s)/n$ or uninform as Chebychev point $s = \cos((0:m)'\pi/m), t = \cos((0:n)'\pi/n)$ , we get

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}R'_{j}(t_{i})\phi_{j}(s) +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi_{j}(s)t_{i}^{\xi-\alpha_{1}}+\int_{0}^{t_{i}} \phi_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}}\right] \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi^{(2)}_{j}(s)t_{i}^{\xi-\alpha_{2}} +\int_{0}^{t_{i}} \phi^{(2)}_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}}\right]+f(t_{i},s), \end{eqnarray}$

(2.6)

by noting the notation, $R_{j}(t_{i}) = \delta_{ij}, R'_{j}(t_{i}) = R_{ij}^{(1, 0)}$ , where $R_{ij}^{(1, 0)}$ is the first order derivative of barycentric matrix. Equation (2.6) can be written as

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}R_{ij}^{(1,0)}\phi_{j}(s) +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi_{j}(s)t_{i}^{\xi-\alpha_{1}}+\int_{0}^{t_{i}} \phi_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}}\right] \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\left[R^{(\xi)}_{j}(0) \phi^{(2)}_{j}(s)t_{i}^{\xi-\alpha_{2}} +\int_{0}^{t_{i}} \phi^{(2)}_{j}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}}\right]+f(t_{i},s). \end{eqnarray}$

(2.7)

Similarly as the discrete $t$ for $s$ , we get

$\begin{eqnarray} \phi_{j}(s) = \sum\limits_{k = 1}^{n}R_{k}(s)\phi_{ik} \end{eqnarray}$

(2.8)

where $\phi_{i}(s_{j}) = \phi(t_{i}, s_{j}) = \phi_{ij}, i = 1, \cdots, m; j = 1, \cdots, n$ and

$\begin{equation} R_{i}(s) = \frac{ \frac{w_{i}}{s-s_{i}}}{ \sum\limits_{k = 1}^{m} \frac{w_{k}}{s-s_{k}}} \end{equation}$

(2.9)

where

$w_{i} = \sum\limits_{j\in J_{i}}(-1)^j\prod\limits_{k = j,j \neq i}^{j+d_{s}}\frac{1}{s_{i}-s_{k}},\ J_{i} = \{j\in \{ 0,1,\cdots, m-d_{s}\}:i-d_{s}\leq j \leq i\},$

is the basis function ^[18].

Taking (2.8) into Eq (2.7), we get

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}R_{ij}^{(1,0)}R_{k}(s)\phi_{ik} +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{k}(s)t_{i}^{\xi-\alpha_{1}}+\int_{0}^{t_{i}} R_{k}(s)\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}}\right]\phi_{ik} \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{k}^{(2)}(s)t_{i}^{\xi-\alpha_{2}} +\int_{0}^{t_{i}} R^{(2)}_{k}(s) \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}}\right]\phi_{ik}+f(t_{i},s). \end{eqnarray}$

(2.10)

By taking $s_{1}, s_{2}, \cdots, s_{n}$ at the mesh-point, we get

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}R_{ij}^{(1,0)}R_{k}(s_{l})\phi_{ik} +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{k}(s_{l})t_{i}^{\xi-\alpha_{1}}+\int_{0}^{t_{i}} R_{k}(s_{l})\frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}}\right]\phi_{ik} \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{k}^{(2)}(s_{l})t_{i}^{\xi-\alpha_{2}} +\int_{0}^{t_{i}} R^{(2)}_{k}(s_{l}) \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}}\right]\phi_{ik}+f(t_{i},s_{l}). \end{eqnarray}$

(2.11)

By noting the notation, $R_{k}(s_{l}) = \delta_{kl}, R''_{k}(s_{l}) = R_{ij}^{(0, 2)}$ , where $R_{ij}^{(0, 2)}$ is the second order derivative of barycentrix matrix.

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}R_{ij}^{(1,0)}\delta_{kl}\phi_{ik} +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) \delta_{kl}t_{i}^{\xi-\alpha_{1}}+\delta_{kl}\int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau) d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}}\right]\phi_{ik} \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{ij}^{(0,2)}t_{i}^{\xi-\alpha_{2}} +R_{ij}^{(0,2)}\int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}}\right]\phi_{ik}+f(t_{i},s_{l}), \end{eqnarray}$

(2.12)

where

$\begin{array}{ll}{ R_{k}(\tau) = \frac{ \frac{\lambda_{k}}{\tau-\tau_{k}}}{ \sum\limits_{k = 0}^{n} \frac{\lambda_{k}}{\tau-\tau_{k}}}} \end{array}$

and

$\left\{ \begin{aligned} R_{i}^{'}(\tau)& = R_{i}(\tau)\left[- \frac{1}{\tau-\tau_{k}}+ \frac{\sum\limits_{s = 0}^{l} \frac{\lambda_{k}}{(\tau-\tau_{k})^{2}}}{\sum\limits_{s = 0}^{l} \frac{\lambda_{k}}{\tau-\tau_{k}}}\right], \\{\vdots}\\ R_{i}^{(\xi+1)}(\tau)& = [R_{i}^{(\xi)}(\tau)]^{'},\xi\in\mathbf{N^{+}}. \end{aligned}\right.$

The integral term of (2.12) can be written as

$\begin{equation} \int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}} = Q_{j}^{\alpha_{1}}(t_{i}) = Q^{\alpha_{1}}_{ji}, \end{equation}$

(2.13)

$\begin{equation} \int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}} = Q_{j}^{\alpha_{2}}(t_{i}) = Q^{\alpha_{2}}_{ji}, \end{equation}$

(2.14)

then we get

$\begin{eqnarray} & \sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}R_{ij}^{(1,0)}\delta_{kl}\phi_{ik} +\mu_{0}\Gamma_{\alpha_{1}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) \delta_{kl}t_{i}^{\xi-\alpha_{1}}+\delta_{kl}Q_{j}^{\alpha_{2}}(t_{i})\right]\phi_{ik} \\& = \Gamma_{\alpha_{2}}^{\xi}\sum\limits_{j = 1}^{m}\sum\limits_{k = 1}^{n}\left[R^{(\xi)}_{j}(0) R_{ij}^{(0,2)}t_{i}^{\xi-\alpha_{2}} +R_{ij}^{(0,2)}Q_{j}^{\alpha_{1}}(t_{i})\right]\phi_{ik}+f(t_{i},s_{l}). \end{eqnarray}$

(2.15)

The integral (2.12) is calculated by

$\begin{eqnarray} Q_{j}^{\alpha_{1}}(t_{i}) = \int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{1}-\xi}} : = \sum\limits_{i = 1}^{g}R^{(\xi+1)}_{i}(\tau_{i}^{\theta,\alpha_{1}})G_{i}^{\theta,\alpha_{1}}, \end{eqnarray}$

(2.16)

and

$\begin{eqnarray} Q_{j}^{\alpha_{2}}(t_{i}) = \int_{0}^{t_{i}} \frac{R^{(\xi+1)}_{j}(\tau)d\tau}{(t_{i}-\tau)^{\alpha_{2}-\xi}} : = \sum\limits_{i = 1}^{g}R^{(\xi+1)}_{i}(\tau_{i}^{\theta,\alpha_{2}})G_{i}^{\theta,\alpha_{2}}, \end{eqnarray}$

(2.17)

where $G_{i}^{\theta, \alpha_{1}}, G_{i}^{\theta, \alpha_{2}}$ are Gauss weights and $\tau_{i}^{\theta, \alpha_{1}}, \tau_{i}^{\theta, \alpha_{2}}$ are Gauss points with weights $(t_{i}-\tau)^{\xi-\alpha_{1}}, (t_{i}-\tau)^{\xi-\alpha_{2}}$ , see reference ^[22].

Equation systems (2.15) can be written as

$\begin{eqnarray} & \left[{R}^{(01)}\otimes{I}_{n}+\Gamma_{\alpha_{2}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{2}}}\right)\right]\left[\begin{array}{cc}{{\phi}_{11}}\\ {\vdots}\\ {{\phi}_{1n} } \\ {{\phi}_{n1}} \\{\vdots} \\{{\phi}_{mn}}\end{array}\right] \\& - \left[\mu_{0}\Gamma_{\alpha_{1}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{1}}}\right)\right]\left[\begin{array}{cc}{{\phi}_{11}}\\ {\vdots}\\ {{\phi}_{1n} } \\ {{\phi}_{n1}} \\{\vdots} \\{{\phi}_{mn}}\end{array}\right] = \left[\begin{array}{cc}{{f}_{11}}\\ {\vdots}\\ {{f}_{1n}} \\ {{f}_{n1}} \\{\vdots} \\{{f}_{mn}}\end{array}\right], \end{eqnarray}$

(2.18)

$I_{m}$ and $I_{n}$ are identity matrices, $\otimes$ is Kronecker product.

Then Eq (2.18) can be noted as

$\begin{equation} \left[{R}^{(01)}\otimes{I}_{n}+\Gamma_{\alpha_{2}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{2}}}\right)-\mu_{0}\Gamma_{\alpha_{1}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{1}}}\right)\right]\Phi = F \end{equation}$

(2.19)

and

$\begin{equation} R\Phi = F, \end{equation}$

(2.20)

with $R = {R}^{(01)}\otimes{I}_{n}+\Gamma_{\alpha_{2}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{2}}}\right)-\mu_{0}\Gamma_{\alpha_{1}}^{\xi} \left({M}_{1}^{(\xi0)}\otimes {I}_{n}+ {I}_{m}\otimes {Q^{\alpha_{1}}}\right)$ and $\Phi = [{{\phi}_{11}} {\dots} {{\phi}_{1n} } {\dots} {{\phi}_{n1}} {\dots} {{\phi}_{mn}}]^{T}, F = [{{f}_{11}} {\dots} {{f}_{1n} } {\dots} {{f}_{n1}} {\dots} {{f}_{mn}}]^{T}$ .

The boundary condition can be solved by substitution method, additional method or elimination method, see ^[26]. We adopt substitution method and additional method to deal with boundary condition.

3. Convergence rate of FC equation

In this part, error estimate of the FC equation is given with $r_n(s) = \sum\limits_{i = 1}^{n} r_{i}(s) \phi_{i}$ to replace $\phi(s)$ , where $r_{i}(s)$ is defined as (2.9) and $\phi_{i} = \phi(s_{i})$ . We also define

$\begin{eqnarray} e(s): = \phi(s)-r_n(s) = (s-s_{i}) \cdots(s-s_{i+d})\phi\left[s_{i}, s_{i+1}, \ldots, s_{i+d}, s\right], \end{eqnarray}$

(3.1)

see reference ^[18].

Then we have

Lemma 1. For $e(s)$ be defined by (3.1) and $\phi(s)\in C^{d+2}[a, b], d = 1, 2, \cdots$ , there

$\begin{equation} \left|e^{(k)}(s)\right| \leq C h^{d-k+1},k = 0,1,\cdots. \end{equation}$

(3.2)

For the FC equation, rational interpolation function of $\phi(t, s)$ is defined as $r_{mn}(t, s)$

$\begin{equation} r_{mn}(t,s) = \frac{ \sum\limits_{i = 1}^{m+d_{s}} \sum\limits_{j = 1}^{n+d_{t}} \frac{w_{i, j}}{\left(s-s_{i}\right)\left(t-t_{j}\right)} {\phi}_{i, j}}{ \sum\limits_{i = 1}^{m+d_{s}} \sum\limits_{j = 1}^{n+d_{t}} \frac{w_{i, j}}{\left(s-s_{i}\right)\left(t-t_{j}\right)}} \end{equation}$

(3.3)

where

$\begin{equation} w_{i, j} = (-1)^{i-d_{s}+j-d_{t}} \sum\limits_{k_{1} \in J_{i}}\prod\limits_{h_{1} = k_{1}, h_{1} \neq j}^{k_{1}+d_{s}}\frac{1}{\left|s_{i}-s_{h_{1}}\right|} \sum\limits_{k_{2} \in J_{i}} \prod\limits_{h_{2} = k_{2}, h_{2} \neq j}^{k_{2}+d_{t}} \frac{1}{\left|t_{j}-t_{h_{2}}\right|}. \end{equation}$

(3.4)

We define $e(t, s)$ be the error of $\phi(t, s)$ as

$\begin{eqnarray} e(t,s) :& = & \phi(t,s)-r_{mn}(t,s) \\ & = & \left(s-s_{i}\right) \cdots\left(s-s_{i+d_{s}}\right)\phi\left[s_{i}, s_{i+1}, \ldots, s_{i+d_1}, s;t\right] \\ &+&\left(t-t_{j}\right) \cdots\left(t-t_{j+d_{t}}\right)\phi\left[s; t_{j}, t_{j+1}, \ldots, t_{j+d_2}, t\right] \\ &-&\left(s-s_{i}\right) \cdots\left(s-s_{i+d_{s}}\right)\left(t-t_{j}\right) \cdots\left(t-t_{j+d_{t}}\right)\phi\left[s_{i}, s_{i+1}, \ldots, s_{i+d_1}, s; t_{j}, t_{j+1}, \ldots, t_{j+d_2}, t\right]. \end{eqnarray}$

(3.5)

With similar analysis of Lemma 1, we have

Theorem 1. For $e(t, s)$ defined as (3.5) and $\phi(t, s) \in C^{d_{s}+2}[a, b]\times C^{d_{t}+2}[0, T]$ , then we have

$\begin{equation} \left|e^{(k_{1},k_{2})}\left(s,t\right)\right| \leq C (h_{s}^{d_{s}-k_{1}+1}+h_{t}^{d_{t}-k_{2}+1}), k_{1},k_{2} = 0,1,\cdots. \end{equation}$

(3.6)

Let $\phi(s_{m}, t_{n})$ be the approximate function of $\phi(t, s)$ and $L$ to be bounded operator, there holds

$\begin{eqnarray} L \phi(t_{m},s_{n}) = f(t_{m},s_{n}) \end{eqnarray}$

(3.7)

and

$\begin{eqnarray} \lim\limits_{m,n\rightarrow \infty}L \phi(t_{m},s_{n}) = \phi(t,s). \end{eqnarray}$

(3.8)

Then we get

Theorem 2. For $\phi(t_{m}, s_{n}):L\phi(t_{m}, s_{n}) = \phi(t, s)$ and $L$ defined as (3.7), there

$|\phi(t,s)-\phi(t_{m},s_{n})|\le C(h^{d_{s}-1}+\tau^{d_{t}-1}).$

Proof. By

$\begin{eqnarray} \begin{array}{ll} &L\phi(t,s)-L\phi(t_{m},s_{n}) \\& = \frac{\partial \phi(t,s)}{\partial t}- \ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi(t,s)}{\partial s^2} +\mu_{0}\ _{0}C_{t}^{1-\alpha_{2}}\phi(t,s)-f(t,s) \\&- \left[\frac{\partial \phi(t_{m},s_{n})}{\partial t}-\ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi(t_{m},s_{n})}{\partial s^2} +\mu_{0}\ _{0}C_{t}^{1-\alpha_{2}}\phi(t_{m},s_{n})-f(t_{m},s_{n}) \right] \\& = \frac{\partial \phi}{\partial t}-\frac{\partial \phi}{\partial t}(t_{m},s_{n}) -\left[_{0}C_{t}^{1-\alpha_{1}} \frac{\partial^2 \phi}{\partial s^2}-\ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi}{\partial s^2}(s_{m},t_{n})\right] \\&+\mu_{0} \left[_{0}C_{t}^{1-\alpha_{2}}\phi(t,s)-\ _{0}C_{t}^{1-\alpha_{2}}(t_{m},s_{n})) \right] \\&-[f(t,s)-f(t_{m},s_{n})] \\&: = E_{1}(t,s)+E_{2}(t,s)+E_{3}(t,s)+E_{4}(t,s), \end{array} \end{eqnarray}$

(3.9)

here

$E_{1}(t,s) = \frac{\partial \phi}{\partial t}-\frac{\partial \phi}{\partial t}(t_{m},s_{n}),$

$E_{2}(t,s) = _{0}C_{t}^{1-\alpha_{1}} \frac{\partial^2 \phi}{\partial s^2}-\ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi}{\partial s^2}(t_{m},s_{n}),$

$E_{3}(t,s) = \mu_{0} \left[_{0}C_{t}^{1-\alpha_{2}}\phi(t,s)-\ _{0}C_{t}^{1-\alpha_{2}}(t_{m},s_{n})) \right],$

$E_{4}(t,s) = f(t,s)-f(t_{m},s_{n}).$

As for $E_{1}(t, s)$ , we get

$\begin{equation} \begin{array}{ll} E_{1}(t,s) = \left|\frac{\partial \phi}{\partial t}(t,s)-\frac{\partial \phi }{\partial t}(t_{m},s_{n})\right| \\ = \left|\frac{\partial \phi }{\partial t}(t,s)-\frac{\partial \phi }{\partial t}(t_{m},s)+\frac{\partial \phi }{\partial t}(t_{m},s)-\frac{\partial \phi }{\partial t}(t_{m},s_{n})\right| \\ \le \left|\frac{\partial \phi }{\partial t}(t,s)-\frac{\partial \phi }{\partial t}(t_{m},s)\right|+\left|\frac{\partial \phi }{\partial t}(t_{m},s)-\frac{\partial \phi }{\partial t}(t_{m},s_{n})\right| \\ \\ = \left| \frac{ \sum\limits_{i = 1}^{m-d_{s}} (-1)^{i}\frac{\partial \phi }{\partial t}[s_{i}, s_{i+1}, \ldots, s_{i+d_1}, s_{n},t]} { \sum\limits_{i = 1}^{m-d_{s}} \lambda_{i}(s)}\right| \\\\ \nonumber+ \left|\frac{ \sum\limits_{j = 1}^{n-d_{t}}(-1)^{j}\frac{\partial \phi }{\partial t}[ t_{j}, t_{j+1}, \ldots, t_{j+d_2}, s_{n},t_{m}]} { \sum\limits_{j = 1}^{n-d_{t}}\lambda_{j}(t)}\right| \\\\ = \left|\frac{\partial e }{\partial t}(t_{m},s)\right|+\left|\frac{\partial e}{\partial t}(t_{m},s_{n})\right|, \end{array} \end{equation}$

we get

$\begin{equation} |E_{1}(t,s)| \le C (h^{d_{s}}+\tau^{d_{t}}). \end{equation}$

(3.10)

As $E_{2}(t, s)$ , we have

$\begin{equation} \begin{array}{ll} E_{2}(t,s) = _{0}C_{t}^{1-\alpha_{1}} \frac{\partial^2 \phi}{\partial s^2}- \ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi}{\partial s^2}(t_{m},s_{n}) \\ = \Gamma_{\alpha_{2}}^{\xi}\left[\frac{\partial^{\xi+2} \phi(0,s)}{\partial t^{\xi}\partial s^2 }s^{\xi-\alpha_{2}}+\int_{0}^{t}\frac{\partial^{\xi+3} \phi(\tau,s)}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t-\tau)^{\alpha_{2}-\xi}}\right] \\ - \Gamma_{\alpha_{2}}^{\xi}\left[\frac{\partial^{\xi+2} \phi(0,s_{n})}{\partial t^{\xi}\partial s^2 }s_{n}^{\xi-\alpha_{2}}+\int_{0}^{t_{m}}\frac{\partial^{\xi+3}\phi(\tau,s_{n})}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t_{m}-\tau)^{\alpha_{2}-\xi}}\right] \\ = \Gamma_{\alpha_{2}}^{\xi}\left[\frac{\partial^{\xi+2} \phi(0,s)}{\partial t^{\xi}\partial s^2 }s^{\xi-\alpha_{2}}-\frac{\partial^{\xi+2} \phi(0,s_{n})}{\partial t^{\xi}\partial s^2 }s_{n}^{\xi-\alpha_{2}}\right] \\ + \Gamma_{\alpha_{2}}^{\xi}\left[\int_{0}^{t}\frac{\partial^{\xi+3} \phi(\tau,s)}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t-\tau)^{\alpha_{2}-\xi}}-\int_{0}^{t_{m}}\frac{\partial^{\xi+3}\phi(\tau,s_{n})}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t_{m}-\tau)^{\alpha_{2}-\xi}}\right] \end{array} \end{equation}$

(3.11)

and

$\begin{equation} \begin{array}{ll} |E_{2}(t,s)| \\ \le \left|\Gamma_{\alpha_{2}}^{\xi}\left[\frac{\partial^{\xi+2} \phi(0,s)}{\partial t^{\xi}\partial s^2 }s^{\xi-\alpha_{2}}-\frac{\partial^{\xi+2} \phi(0,s_{n})}{\partial t^{\xi}\partial s^2 }s_{n}^{\xi-\alpha_{2}}\right]\right| \\ + \left|\Gamma_{\alpha_{2}}^{\xi}\left[\int_{0}^{t}\frac{\partial^{\xi+3} \phi(\tau,s)}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t-\tau)^{\alpha_{2}-\xi}}-\int_{0}^{t_{m}}\frac{\partial^{\xi+3}\phi(\tau,s_{n})}{\partial \tau^{\xi+1} \partial s^{2}}\frac{d\tau}{(t_{m}-\tau)^{\alpha_{2}-\xi}}\right]\right| \\ \le |\Gamma_{\alpha_{2}}^{\xi}| \left|\frac{\partial^{\xi+2} \phi}{\partial t^{\xi} \partial s^{2}}(0,s)-\frac{\partial^{\xi+2} \phi }{\partial t^{\xi} \partial s^{2}}(0,s_{n})\right|+ |\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t,s)-\frac{\partial^{\xi+3} \phi }{\partial t^{\xi+1} \partial s^{2}}(t_{m},s_{n})\right| \\ : = E_{21}(t,s)+E_{22}(t,s) \end{array} \end{equation}$

(3.12)

where

$\begin{equation} \begin{array}{ll} E_{21}(t,s) = |\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+2} \phi}{\partial t^{\xi} \partial s^{2}}(0,s)-\frac{\partial^{\xi+2} \phi }{\partial t^{\xi} \partial s^{2}}(0,s_{n})\right|, \\ E_{22}(t,s) = |\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t,s)-\frac{\partial^{\xi+3} \phi }{\partial t^{\xi+1} \partial s^{2}}(t_{m},s_{n})\right|. \end{array} \end{equation}$

(3.13)

Now we estimate $E_{21}(t, s)$ and $E_{22}(t, s)$ part by part, for the second part we have

$\begin{equation} \begin{array}{ll} E_{22}(t,s) = |\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t,s)-\frac{\partial^{\xi+3} \phi }{\partial t^{\xi+1} \partial s^{2}}(t_{m},s_{n})\right| \\ = |\Gamma_{\alpha_{2}}^{\xi}| \left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t,s)-\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t_{m},s)+\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t_{m},s)-\frac{\partial^{\xi+3} \phi }{\partial t^{\xi+1} \partial s^{2}}(t_{m},s_{n})\right| \\ \le |\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t,s)-\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t_{m},s)\right|+|\Gamma_{\alpha_{2}}^{\xi}|\left|\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1} \partial s^{2}}(t_{m},s)-\frac{\partial^{\xi+3} \phi }{\partial t^{\xi+1} \partial s^{2}}(t_{m},s_{n})\right| \\ \\ = |\Gamma_{\alpha_{2}}^{\xi}| \left| \frac{ \sum\limits_{i = 1}^{m-d_{s}} (-1)^{i}\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1}\partial s^2}[s_{i}, s_{i+1}, \ldots, s_{i+d_1}, s_{n},t]} { \sum\limits_{i = 1}^{m-d_{s}} \lambda_{i}(s)}\right| \\\\ \nonumber+ |\Gamma_{\alpha_{2}}^{\xi}| \left|\frac{ \sum\limits_{j = 1}^{n-d_{t}}(-1)^{j}\frac{\partial^{\xi+3} \phi}{\partial t^{\xi+1}\partial s^2}[ t_{j}, t_{j+1}, \ldots, t_{j+d_2}, s_{n},t_{m}]} { \sum\limits_{j = 1}^{n-d_{t}}\lambda_{j}(t)}\right| \\\\ = |\Gamma_{\alpha_{2}}^{\xi}| \left|\frac{\partial^{\xi+3} e}{\partial t^{\xi+1}\partial s^2}(t_{m},s)\right|+ |\Gamma_{\alpha_{2}}^{\xi}| \left|\frac{\partial^{\xi+3} e}{\partial t^{\xi+1}\partial s^2}(t_{m},s_{n})\right|, \end{array} \end{equation}$

then we have

$\begin{equation} |E_{22}(t,s)|\le \left|\frac{\partial^{\xi+3} e}{\partial t^{\xi+1}\partial s^2}(t_{m},s)\right|+\left|\frac{\partial^{\xi+3} e}{\partial t^{\xi+1}\partial s^2}(t_{m},s_{n})\right|\le C (h^{d_{s}-\xi}+\tau^{d_{t}-1}). \end{equation}$

(3.14)

For $E_{21}(t, s)$ , we get

$\begin{equation} |E_{21}(t,s)|\le C (h^{d_{s}+1-\xi}+\tau^{d_{t}-1}). \end{equation}$

(3.15)

Similarly as $E_{2}(t, s)$ , for $E_{3}(t, s)$ we have

$\begin{equation} |E_{3}(t,s)| \le C (h^{d_{s}}+\tau^{d_{t}}). \end{equation}$

(3.16)

Combining (3.9), (3.14), (3.16) together, proof of Theorem 2 is completed.

4. Numerical examples

In this part, one example is presented to test the theorem. The nonuniform partition in this experiment defined as second kind of Chybechev point $s = \cos((0:m)'\pi/m), t = \cos((0:n)'\pi/n)$ .

Example 1. Consider the FC equation

$\frac{\partial \phi}{\partial t} = \ _{0}C_{t}^{1-\alpha_{1}}\frac{\partial^2 \phi}{\partial s^2}\phi(t,s) -\mu_{0}\ _{0}C_{t}^{1-\alpha_{2}}\phi(t,s)+f(t,s) ,0 \le s\le 1,0 \le t \le T,$

with the analysis solutions is

$\phi(t,s) = t^2\sin(\pi s),$

with the initial condition

$\phi(s, 0) = 0,$

and boundary condition

$\phi(0, t) = \phi(1, t) = 0,$

and

$f(t,s) = 2\left(t+\frac{\pi^2 t^{1+\alpha_{1}}}{\Gamma(2+\alpha_{1})} +\frac{ t^{1+\alpha_{2}}}{\Gamma(2+\alpha_{2})}\right)\sin(\pi s).$

In and , errors of $m = n = 10$ , $[a, b] = [0, 1]$ and $m = n = 10, d_{t} = d_{s} = 7$ , $[a, b] = [0, 1]$ in Example 1. (a) uniform; (b) nonuniform for FC equation by rational interpolation collocation methods are presented, respectively. From the figure, we know that the precision can reach to $10^{-6}$ for both uniform and nonuniform partition.

Figure 1. Errors of $m = n = 10$ , $[a, b] = [0, 1]$ in Example 1 (a) uniform; (b) nonuniform.

DownLoad: Full-Size Img PowerPoint

Figure 2. Errors of $m = n = 10, d_{t} = d_{s} = 7$ , $[a, b] = [0, 1]$ in Example 1 (a) uniform; (b) nonuniform.

DownLoad: Full-Size Img PowerPoint

In , errors of the FC equation with $m = n = 10, \alpha_1 = \alpha_2 = 0.2$ for substitution methods and additional methods are presented, there are nearly no difference for the two methods. Additional method is more simple than substitution methods to add the boundary condition. In the following, we choosing the substitution method to deal with the boundary condition.

Table 1. Errors of FC equation with $m = n = 10, \alpha_1 = \alpha_2 = 0.2$ .

	method of substitution		method of additional
	uniform	nonuniform	uniform	nonuniform
Larange	1.4662e-06	2.1919e-08	2.7900e-07	1.4310e-07
Rational	1.3038e-05	2.4541e-07	4.9788e-06	1.4310e-07

| Show Table

DownLoad: CSV

Errors of the FC equation for $\alpha_1 = 0.4, \alpha_{2} = 0.6, d_{t} = d_{s} = 5$ with $t = 0.1, 0.9, 1, 5, 10, 15$ are presented under the uniform and nonuniform in . As the time variable become from $0.5$ to $15$ , there are high accuracy for our methods. We can improve the accuracy by increasing $m, n$ or choosing the parameter $d_{t}, d_{s}$ approximately which means our methods is useful.

Table 2. Errors of FC equation for $\alpha_1 = 0.4, \alpha_{2} = 0.6, d_{t} = d_{s} = 5$ .

	uniform	nonuniform	uniform	nonuniform
$t$	$(12, 12)$	$(12, 12)$	$(12, 12) d_{t}=d_{s}=5$	$(12, 12) d_{t}=d_{s}=5$
0.5	2.1021e-11	3.8250e-09	6.8506e-06	1.6436e-06
1	9.0394e-13	4.4206e-10	4.6667e-06	7.8141e-07
5	6.1833e-12	5.6655e-08	2.3777e-04	4.2230e-05
10	1.0094e-12	8.5622e-07	1.9813e-04	1.5634e-05
15	3.5397e-12	1.8827e-05	8.5498e-04	8.2551e-05

| Show Table

DownLoad: CSV

In , errors of $\alpha_1 = 0.01, 0.1, 0.3, 0.5, 0.9, 0.99$ under uniform with $m = n = 10, d_{t} = 5, d_{s} = 5$ with $\alpha_2 = 0.1, 0.4, 0.6, 0.8, 0.99$ are presented under the uniform partition. From the table, we know that for different $\alpha_1, \alpha_2$ our methods have high accuracy with little number $m$ and $n$ . In the following table, numerical results are presented to test our theorem. From and , error of uniform for $\alpha_1 = \alpha_2 = 0.2, d_{s} = 5$ with different $d_{t}$ are given, the convergence rate is $O(h^{d_{t}})$ . From , with space variable uniform for $\alpha_1 = \alpha_2 = 0.2, d_{t} = 5$ , the convergence rate is $O(h^{7})$ , we will investigate in future paper. For and , the errors of Chebyshev partition for $s$ and $t$ are presented. For $d_{t} = 5$ , the convergence rate is $O(h^{d_{s}})$ in , while in , the convergence rate is $O(h^{d_{t}})$ which agrees with our theorem.

Table 3. Errors of $\alpha_1$ under uniform with $m = n = 10, d_{t} = 5, d_{s} = 5$ .

$\alpha_1$	$\alpha_2=0.1$	$\alpha_2=0.4$	$\alpha_2=0.6$	$\alpha_2=0.8$	$\alpha_2=0.99$
0.01	1.0153e-04	1.0246e-04	1.0300e-04	1.0346e-04	1.0384e-04
0.1	1.2753e-05	1.2865e-05	1.2930e-05	1.2987e-05	1.3033e-05
0.3	2.7464e-05	2.7704e-05	2.7845e-05	2.7971e-05	2.8074e-05
0.5	4.5746e-06	4.6152e-06	4.6399e-06	4.6609e-06	4.6794e-06
0.9	9.0295e-06	9.1193e-06	9.1240e-06	9.2142e-06	9.2479e-06
0.99	1.8981e-06	1.8247e-06	1.5293e-06	1.9193e-06	2.0670e-06

| Show Table

DownLoad: CSV

Table 4. Errors of uniform for $\alpha_1 = \alpha_2 = 0.2, d_{s} = 5$ .

$m, n$	$d_{t}=2$		$d_{t}=3$		$d_{t}=4$		$d_{t}=5$
8	1.3626e-02		6.9619e-03		2.0708e-03		9.8232e-04
10	9.6780e-03	1.5332	3.4354e-03	3.1653	6.9542e-04	4.8900	3.2829e-04	4.9117
12	7.0485e-03	1.7389	1.9408e-03	3.1320	2.9186e-04	4.7621	1.3132e-04	5.0255
14	5.4466e-03	1.6725	1.2017e-03	3.1097	1.4211e-04	4.6686	6.0148e-05	5.0654

| Show Table

DownLoad: CSV

Table 5. Errors of uniform for $\alpha_1 = \alpha_2 = 0.2, d_{t} = 5$ .

$m, n$	$d_{s}=2$		$d_{s}=3$		$d_{s}=4$
8	4.9495e-04		4.9492e-04		4.9486e-04
10	1.0051e-04	7.1443	1.0053e-04	7.1431	1.0053e-04	7.1426
12	2.7700e-05	7.0690	2.7711e-05	7.0679	2.7714e-05	7.0673
14	9.4272e-06	6.9921	9.4315e-06	6.9917	9.4314e-06	6.9925

| Show Table

DownLoad: CSV

Table 6. Errors of non-uniform partition with $\alpha_1 = \alpha_2 = 0.2, d_{t} = 5$ .

$m, n$	$d_{s}=2$		$d_{s}=3$		$d_{s}=4$
8	2.8113e-05		2.8110e-05		2.8108e-05
10	2.1197e-05	1.2654	2.1196e-05	1.2652	2.1195e-05	1.2651
12	6.6990e-06	6.3180	6.6989e-06	6.3178	6.6988e-06	6.3176
14	1.6712e-06	9.0069	1.6712e-06	9.0068	1.6712e-06	9.0067

| Show Table

DownLoad: CSV

Table 7. Errors of non-uniform partition $\alpha_1 = \alpha_2 = 0.2, d_{s} = 5$ .

$m, n$	$d_{t}=2$		$d_{t}=3$		$d_{t}=4$		$d_{t}=5$
8	3.1539e-02		8.7995e-03		2.1930e-03		3.3004e-04
10	2.4329e-02	1.1632	4.0288e-03	3.5010	2.7133e-04	9.3648	2.2278e-04	1.7613
12	1.5223e-02	2.5716	1.9127e-03	4.0859	9.5194e-05	5.7449	5.1702e-05	8.0116
14	1.1407e-02	1.8721	1.1143e-03	3.5049	3.5772e-05	6.3493	1.1369e-05	9.8255

| Show Table

DownLoad: CSV

In the following table, $\alpha_{1} = 0.4, \alpha_{2} = 0.6$ is chosen to present numerical results. From and , error of uniform partition $d_{t} = 5$ with different $d_{s}$ are given, the convergence rate is $O(h^{7})$ . From , with space variable $s, d_{s} = 5$ , the convergence rate is $O(h^{d_{t}})$ which agrees with our theorem.

Table 8. Errors of uniform with $\alpha_{1} = 0.4, \alpha_{2} = 0.6, d_{t} = 5$ .

$m, n$	$d_{s}=2$		$d_{s}=3$		$d_{s}=4$
8	4.9427e-04		4.9426e-04		4.9414e-04
10	1.0035e-04	7.1455	1.0041e-04	7.1427	1.0041e-04	7.1413
12	2.7639e-05	7.0720	2.7674e-05	7.0684	2.7684e-05	7.0669
14	9.3984e-06	6.9977	9.4153e-06	6.9942	9.4254e-06	6.9895

| Show Table

DownLoad: CSV

Table 9. Errors of uniform with $\alpha_{1} = 0.4, \alpha_{2} = 0.6, d_{s} = 5$ .

$m, n$	$d_{t}=1$		$d_{t}=2$		$d_{t}=3$		$d_{t}=4$
8	1.3587e-02		6.9513e-03		2.0677e-03		9.8084e-04
10	9.6497e-03	1.5334	3.4314e-03	3.1637	6.9462e-04	4.8884	3.2791e-04	4.9102
12	7.0259e-03	1.7404	1.9389e-03	3.1311	2.9157e-04	4.7613	1.3118e-04	5.0249
14	5.4269e-03	1.6752	1.2005e-03	3.1096	1.4198e-04	4.6682	6.0090e-05	5.0648

| Show Table

DownLoad: CSV

For and , the errors of Chebyshev partition for non-uniform with $\alpha_{1} = 0.4, \alpha_{2} = 0.6$ are presented. For $d_{t} = 5$ , the convergence rate is $O(h^{7})$ in , while in , the convergence rate is $O(h^{d_{t}})$ which agrees with our theorem.

Table 10. Errors of non-uniform with $\alpha_{1} = 0.4, \alpha_{2} = 0.6, d_{s} = 5$ .

$m, n$	$d_{t}=1$		$d_{t}=2$		$d_{t}=3$		$d_{t}=4$
8	3.1481e-02		8.7825e-03		2.1876e-03		3.2930e-04
10	2.4263e-02	1.1671	4.0219e-03	3.5000	2.7124e-04	9.3553	2.2231e-04	1.7606
12	1.5185e-02	2.5704	1.9076e-03	4.0912	9.5106e-05	5.7481	5.1649e-05	8.0057
14	1.1373e-02	1.8751	1.1117e-03	3.5026	3.5733e-05	6.3504	1.1365e-05	9.8211

| Show Table

DownLoad: CSV

Table 11. Errors of non-uniform with $\alpha_{1} = 0.4, \alpha_{2} = 0.6, d_{t} = 5$ .

$m, n$	$d_{s}=2$		$d_{s}=3$		$d_{s}=4$
8	2.8065e-05		2.8059e-05		2.8056e-05
10	2.1156e-05	1.2665	2.1154e-05	1.2660	2.1153e-05	1.2656
12	6.6875e-06	6.3168	6.6874e-06	6.3164	6.6873e-06	6.3161
14	1.6693e-06	9.0033	1.6693e-06	9.0031	1.6693e-06	9.0030

| Show Table

DownLoad: CSV

5. Concluding remarks

In this paper, BRIM was used to solve the (1+1) dimensional FC equation that is presented. For fractional-order PDEs, the convergence order is seriously affected by the orders of fractional derivatives. By fractional integration, the singularity of the fractional derivative of the FC equation can be changed to nonsingular integral, with adding one order to the derivatives of density function. So there are no effects on the orders of fractional derivatives. The singularity of fractional derivative is overcome by the integral to density function from the singular kernel. For the arbitrary fractional derivative, the new Gauss formula is constructed to calculated it simply. For the Diriclet boundary condition, the FC equation is changed to the discrete FC equation and the matrix equation of it is given. In the future, the FC equation with Nuemann condition can be solved by BRIM, and high dimensional FC equation can also be studied by our methods.

Acknowledgments

The work of Jin Li was supported by Natural Science Foundation of Shandong Province (Grant No. ZR2022MA003).

Conflicts of interest

The authors declare that they have no conflicts of interest.

References

[1]	M. Cuesta Leon, Existence results for quasilinear problems via ordered sub- and supersolutions, Ann. Fac. Sci. Toulouse Math., 6 (1997), 591–608. http://www.numdam.org/item?id=AFST_1997_6_6_4_591_0
[2]	J. Schoenenberger-Deuel, P. Hess, A criterion for the existence of solutions of non-linear elliptic boundary value problems, Proc. Roy. Soc. Edinburgh Sect. A, 74 (1976), 49–54. https://doi-org.proxy.swarthmore.edu/10.1017/s030821050001653x doi: 10.1017/s030821050001653x
[3]	J. M. Arrieta, R. Pardo, A. Rodríguez-Bernal, Infinite resonant solutions and turning points in a problem with unbounded bifurcation, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 20 (2010), 2885–2896. https://doi.org/10.1142/S021812741002743X doi: 10.1142/S021812741002743X
[4]	J. M. Arrieta, R. Pardo, A. Rodríguez-Bernal, Bifurcation and stability of equilibria with asymptotically linear boundary conditions at infinity, Proc. Roy. Soc. Edinburgh Sect. A, 137 (2007), 225–252.
[5]	J. M. Arrieta, R. Pardo, A. Rodríguez-Bernal, Equilibria and global dynamics of a problem with bifurcation from infinity, J. Differ. Equ., 246 (2009), 2055–2080. https://doi.org/10.1016/j.jde.2008.09.002 doi: 10.1016/j.jde.2008.09.002
[6]	P. Liu, J. Shi, Bifurcation of positive solutions to scalar reaction-diffusion equations with nonlinear boundary condition, J. Differ. Equ., 264 (2018), 425–454. https://doi.org/10.1016/j.jde.2017.09.014 doi: 10.1016/j.jde.2017.09.014
[7]	N. Mavinga, Generalized eigenproblem and nonlinear elliptic equations with nonlinear boundary conditions, Proc. Roy. Soc. Edinburgh Sect. A, 142 (2012), 137–153. https://doi.org/10.1017/S0308210510000065 doi: 10.1017/S0308210510000065
[8]	N. Mavinga, R. Pardo, Bifurcation from infinity for reaction-diffusion equations under nonlinear boundary conditions, Proc. Roy. Soc. Edinburgh Sect. A, 147 (2017), 649–671. https://doi.org/10.1017/S0308210516000251 doi: 10.1017/S0308210516000251
[9]	C. Morales-Rodrigo, A. Suárez, Some elliptic problems with nonlinear boundary conditions, in Spectral theory and nonlinear analysis with applications to spatial ecology, World Sci. Publ., Hackensack, NJ, 2005,175–199. https://doi.org/10.1142/9789812701589_0009
[10]	J. M. Arrieta, A. N. Carvalho, A. Rodríguez-Bernal, Parabolic problems with nonlinear boundary conditions and critical nonlinearities, J. Differ. Equ., 156 (1999), 376–406. https://doi.org/10.1006/jdeq.1998.3612 doi: 10.1006/jdeq.1998.3612
[11]	R. S. Cantrell, C. Cosner, On the effects of nonlinear boundary conditions in diffusive logistic equations on bounded domains, J. Differ. Equ., 231 (2006), 768–804. https://doi.org/10.1016/j.jde.2006.08.018 doi: 10.1016/j.jde.2006.08.018
[12]	A. A. Lacey, J. R. Ockendon, J. Sabina, Multidimensional reaction diffusion equations with nonlinear boundary conditions, SIAM J. Appl. Math., 58 (1998), 1622–1647. https://doi.org/10.1137/S0036139996308121 doi: 10.1137/S0036139996308121
[13]	C. V. Pao, Nonlinear parabolic and elliptic equations, Plenum Press, New York, 1992.
[14]	K. Akô, On the Dirichlet problem for quasi-linear elliptic differential equations of the second order, J. Math. Soc. Japan, 13 (1961), 45–62, https://doi.org/10.2969/jmsj/01310045 doi: 10.2969/jmsj/01310045
[15]	H. Amann, On the existence of positive solutions of nonlinear elliptic boundary value problems, Indiana Univ. Math. J., 21 (1971), 125–146, https://doi.org/10.1512/iumj.1971.21.21012 doi: 10.1512/iumj.1971.21.21012
[16]	D. H. Sattinger, Monotone methods in nonlinear elliptic and parabolic boundary value problems, Indiana Univ. Math. J., 21 (1971), 979–1000. https://doi.org/10.1512/iumj.1972.21.21079 doi: 10.1512/iumj.1972.21.21079
[17]	H. Amann, M. G. Crandall, On some existence theorems for semi-linear elliptic equations, Indiana Univ. Math. J., 27 (1978), 779–790. https://doi.org/10.1512/iumj.1978.27.27050 doi: 10.1512/iumj.1978.27.27050
[18]	E. N. Dancer, G. Sweers, On the existence of a maximal weak solution for a semilinear elliptic equation, Differ. Integral Equ., 2 (1989), 533–540.
[19]	D. Motreanu, A. Sciammetta, E. Tornatore, A sub-supersolution approach for {N}eumann boundary value problems with gradient dependence, Nonlinear Anal. Real World Appl., 54 (2020), 103096. https://doi.org/10.1016/j.nonrwa.2020.103096 doi: 10.1016/j.nonrwa.2020.103096
[20]	H. Amann, Nonlinear elliptic equations with nonlinear boundary conditions, in New developments in differential equations (Proc. 2nd Scheveningen Conf., Scheveningen, 1975), 1976, 43–63. North–Holland Math. Studies, Vol. 21.
[21]	R. A. Adams, J. J. F. Fournier, Sobolev spaces, vol. 140 of Pure and Applied Mathematics (Amsterdam), 2nd edition, Elsevier/Academic Press, Amsterdam, 2003.
[22]	S. Carl, V. K. Le, D. Motreanu, Nonsmooth variational problems and their inequalities, Springer Monographs in Mathematics, Springer, New York, 2007. https://doi.org/10.1007/978-0-387-46252-3 Comparison principles and applications.
[23]	A. Kufner, O. John, S. Fučík, Function spaces, Noordhoff International Publishing, Leyden; Academia, Prague, 1977, Monographs and Textbooks on Mechanics of Solids and Fluids; Mechanics: Analysis.
[24]	P. Hess, On the solvability of nonlinear elliptic boundary value problems, Indiana Univ. Math. J., 25 (1976), 461–466. https://doi.org/10.1512/iumj.1976.25.25036 doi: 10.1512/iumj.1976.25.25036
[25]	J.-L. Lions, Quelques méthodes de résolution des problèmes aux limites non linéaires, Dunod; Gauthier-Villars, Paris, 1969.
[26]	H. Brezis, A. C. Ponce, Kato's inequality up to the boundary, Commun. Contemp. Math., 10 (2008), 1217–1241. https://doi.org/10.1142/S0219199708003241 doi: 10.1142/S0219199708003241
[27]	H. Brezis, Functional analysis, Sobolev spaces and partial differential equations, Universitext, Springer, New York, 2011.

This article has been cited by:

1.	Xueying Ding, Jianquan Lu, Haitao Li, Yang Liu, Recent developments of Boolean networks with switching and constraints, 2023, 54, 0020-7721, 2765, 10.1080/00207721.2023.2250490
2.	Arnon Ploymukda, Pattrawut Chansangiam, Metric geometric means with arbitrary weights of positive definite matrices involving semi-tensor products, 2023, 8, 2473-6988, 26153, 10.3934/math.20231333
3.	Chenyang Bian, Leihao Du, Zhipeng Zhang, 2024, Synchronization of Probabilistic Boolean Networks under State-Flipped Control, 979-8-3503-7922-8, 1874, 10.1109/YAC63405.2024.10598672
4.	Guodong Zhao, Haitao Li, Ting Hou, Survey of semi-tensor product method in robustness analysis on finite systems, 2023, 20, 1551-0018, 11464, 10.3934/mbe.2023508
5.	Shao Shao, Linying Xiang, Multi-sensor fusion particle filtering for Boolean networks with multi-step randomly-delayed measurements, 2023, 547, 09252312, 126386, 10.1016/j.neucom.2023.126386
6.	Jingjing Dong, Xiaofeng Ma, Lanlan He, Xin Huang, Jianping Zhou, Energy-to-peak control for switched systems with PDT switching, 2023, 31, 2688-1594, 5267, 10.3934/era.2023268
7.	Dandan Zuo, Wansheng Wang, Lulu Zhang, Jing Han, Ling Chen, Non-fragile sampled-data control for synchronizing Markov jump Lur'e systems with time-variant delay, 2024, 32, 2688-1594, 4632, 10.3934/era.2024211
8.	Liqing Wang, Zheng-Guang Wu, Ying Shen, Estimation of State and Mode for Boolean Networks With Markov Jump Parameters, 2024, 54, 2168-2216, 4188, 10.1109/TSMC.2024.3378334
9.	Guodong Zhao, Wenxin Zhang, Haitao Li, Robustness analysis for switched cyber‐physical systems via logical networks, 2025, 1561-8625, 10.1002/asjc.3583

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Electronic Research Archive

1.1 1.7

Metrics

Article views(1831) PDF downloads(116) Cited by(1)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Electronic Research Archive

Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary

Related Papers:

Abstract

1. Introduction

2. Matrix equation of FC equation

3. Convergence rate of FC equation

4. Numerical examples

5. Concluding remarks

Acknowledgments

Conflicts of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Catalog

Electronic Research Archive

Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary

Related Papers:

Abstract

1. Introduction

2. Matrix equation of FC equation

3. Convergence rate of FC equation

4. Numerical examples

5. Concluding remarks

Acknowledgments

Conflicts of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog