Persistent and carcinogenic polycyclic aromatic hydrocarbons in the north-western coastal marine environment of India

P. U. Singare; J.P. Shirodkar; P. U. Singare; J.P. Shirodkar

doi:10.3934/environsci.2021012

AIMS Environmental Science

2021, Volume 8, Issue 2: 169-189. doi: 10.3934/environsci.2021012

Previous Article Next Article

Research article Special Issues

Persistent and carcinogenic polycyclic aromatic hydrocarbons in the north-western coastal marine environment of India

P. U. Singare ^{1
,
,},
J.P. Shirodkar ²

1.
Department of Chemistry, N.M. Institute of Science, Bhavan's College, Munshi Nagar, Andheri (West), Mumbai 400058, India
2.
Research scholar at Department of Chemistry, Shri Jagdish Prasad Jhabarmal Tibrewala University, Vidyanagari, Jhunjhunu-333 001, (Rajasthan), India

Received: 16 March 2021 Accepted: 28 April 2021 Published: 10 May 2021

The carcinogenic and endocrine-disrupting PAHs were investigated in surface water of the north-western creeks of India. The concentrations of ΣPAHs were found to vary in the range of 114.32-347.04 μg L^-1, (mean 224.78 ± 8.85 μg L^-1), out of which 49.12% contribution is due to ΣC-PAHs. The assessment of toxicity and biological risk arising due to individual C-PAHs was made by calculating their toxic equivalent quantity. The level of individual C-PAHs was reported exceeding the final chronic values, Canadian water quality guideline values and Netherlands maximum permissible concentration values set for the protection of aquatic life. The mean BaP concentration (10.32 ± 2.75 μg L^-1) was above the European Directive 2008/105/EC Environmental Quality Standards (EQS) value; while the sum of BkF + BbF (26.76 μg L^-1) and BghiP + InP (19.59 μg L^-1) were significantly higher than that set by the EQS. The results of the present study will help in understanding the global distribution and fate of PAHs which is required for implementing the necessary steps towards mitigation of the ecotoxicological risk arising due to the existence of such contaminants in the aquatic environment across the world.

Keywords:

Citation: P. U. Singare, J.P. Shirodkar. Persistent and carcinogenic polycyclic aromatic hydrocarbons in the north-western coastal marine environment of India[J]. AIMS Environmental Science, 2021, 8(2): 169-189. doi: 10.3934/environsci.2021012

Related Papers:

[1]	Jin Li, Yongling Cheng . Barycentric rational interpolation method for solving time-dependent fractional convection-diffusion equation. Electronic Research Archive, 2023, 31(7): 4034-4056. doi: 10.3934/era.2023205
[2]	Jin Li, Yongling Cheng . Barycentric rational interpolation method for solving KPP equation. Electronic Research Archive, 2023, 31(5): 3014-3029. doi: 10.3934/era.2023152
[3]	Xiumin Lyu, Jin Li, Wanjun Song . Numerical solution of a coupled Burgers' equation via barycentric interpolation collocation method. Electronic Research Archive, 2025, 33(3): 1490-1509. doi: 10.3934/era.2025070
[4]	Waleed Mohamed Abd-Elhameed, Omar Mazen Alqubori, Abdulrahman Khalid Al-Harbi, Mohammed H. Alharbi, Ahmed Gamal Atta . Generalized third-kind Chebyshev tau approach for treating the time fractional cable problem. Electronic Research Archive, 2024, 32(11): 6200-6224. doi: 10.3934/era.2024288
[5]	Li Tian, Ziqiang Wang, Junying Cao . A high-order numerical scheme for right Caputo fractional differential equations with uniform accuracy. Electronic Research Archive, 2022, 30(10): 3825-3854. doi: 10.3934/era.2022195
[6]	Sahar Albosaily, Wael Mohammed, Mahmoud El-Morshedy . The exact solutions of the fractional-stochastic Fokas-Lenells equation in optical fiber communication. Electronic Research Archive, 2023, 31(6): 3552-3567. doi: 10.3934/era.2023180
[7]	Ziqing Yang, Ruiping Niu, Miaomiao Chen, Hongen Jia, Shengli Li . Adaptive fractional physical information neural network based on PQI scheme for solving time-fractional partial differential equations. Electronic Research Archive, 2024, 32(4): 2699-2727. doi: 10.3934/era.2024122
[8]	Nelson Vieira, M. Manuela Rodrigues, Milton Ferreira . Time-fractional telegraph equation of distributed order in higher dimensions with Hilfer fractional derivatives. Electronic Research Archive, 2022, 30(10): 3595-3631. doi: 10.3934/era.2022184
[9]	Ping Zhou, Hossein Jafari, Roghayeh M. Ganji, Sonali M. Narsale . Numerical study for a class of time fractional diffusion equations using operational matrices based on Hosoya polynomial. Electronic Research Archive, 2023, 31(8): 4530-4548. doi: 10.3934/era.2023231
[10]	Yin Yang, Sujuan Kang, Vasiliy I. Vasil'ev . The Jacobi spectral collocation method for fractional integro-differential equations with non-smooth solutions. Electronic Research Archive, 2020, 28(3): 1161-1189. doi: 10.3934/era.2020064

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]	Kanzari F, Syakti AD, Asia L, et al. (2014) Distributions and sources of persistent organic pollutants (aliphatic hydrocarbons, PAHs, PCBs and pesticides) in surface sediments of an industrialized urban river (Huveaune), France. Sci Total Environ 478: 141-151. doi: 10.1016/j.scitotenv.2014.01.065
[2]	Abdel-Shafy HI, Mansour MSM (2016) A review on polycyclic aromatic hydrocarbons: source, environmental impact, effect on human health and remediation. Egypt J Pet 25: 107-123. doi: 10.1016/j.ejpe.2015.03.011
[3]	Lin Y, Qiu X, Ma Y, et al. (2015) Concentrations and spatial distribution of polycyclic aromatic hydrocarbons (PAHs) and nitrated PAHs (NPAHs) in the atmosphere of North China, and the transformation from PAHs to NPAHs. Environ Pollut 196: 164-170. doi: 10.1016/j.envpol.2014.10.005
[4]	Friedman CL, Pierce JR, Selin NE (2014) Assessing the Influence of Secondary Organic versus Primary Carbonaceous Aerosols on Long-Range Atmospheric Polycyclic Aromatic Hydrocarbon Transport. Environ Sci Technol 48: 3293-3302. doi: 10.1021/es405219r
[5]	Basavaiah N, Mohite RD, Singare PU, et al. (2017) Vertical distribution, composition profiles, sources and toxicity assessment of PAH residues in the reclaimed mudflat sediments from the adjacent Thane Creek of Mumbai. Mar Pollut Bull 118: 112-124. doi: 10.1016/j.marpolbul.2017.02.049
[6]	Singare PU (2015) Studies on Polycyclic Aromatic Hydrocarbons in Sediments of Mithi River of Mumbai, India: Assessment of Sources, Toxicity Risk and Biological Impact. Mar Pollut Bull 101: 232-242.
[7]	Tongo I, Ezemonye L, Akpeh K (2017) Levels, distribution and characterization of polycyclic aromatic hydrocarbons (PAHs) in Ovia river, Southern Nigeria. J Environ Chem Eng 5: 504-512. doi: 10.1016/j.jece.2016.12.035
[8]	Hong WJ, Jia H, Li YF et al. (2016) Polycyclic aromatic hydrocarbons (PAHs) and alkylated PAHs in the coastal seawater, surface sediment and oyster from Dalian, Northeast China. Ecotoxicol Environ Saf 128: 11-20. doi: 10.1016/j.ecoenv.2016.02.003
[9]	Qin XB, Sun HW, Wang CP, et al. (2010) Impacts of crab bioturbation on the fate of polycyclic aromatic hydrocarbons in sediment from the Beitang Estuary of Tianjin, China. Environ Toxicol Chem 29: 1248-1255.
[10]	Gu YG, Lin Q, Lu TT, et al. (2013) Levels, composition profiles and sources of polycyclic aromatic hydrocarbons in surface sediments from Nan'ao Island, a representative mariculture base in South China. Mar Pollut Bull 75: 310-316. doi: 10.1016/j.marpolbul.2013.07.039
[11]	Hawliczek A, Nota B, Cenijn P, et al. (2012) Developmental toxicity and endocrine disrupting potency of 4-azapyrene, benzo[b]fluorene and retene in the zebrafish Danio rerio. Reprod Toxicol 33: 213-223. doi: 10.1016/j.reprotox.2011.11.001
[12]	Liu LY, Wang JZ, Wei GL, et al. (2012) Polycyclic aromatic hydrocarbons (PAHs) in continental shelf sediment of China: implications for anthropogenic influences on coastal marine environment. Environ Pollut 167: 155-162. doi: 10.1016/j.envpol.2012.03.038
[13]	Lewis MA, Russel MJ (2015) Contaminant profiles for surface water, sediment, flora and fauna associated with the mangrove fringe along middle and lower eastern Tampa Bay. Mar Pollut Bull 95: 273-282. doi: 10.1016/j.marpolbul.2015.04.001
[14]	Santana JL, Massone CG, Valdes M, et al. (2015) Occurrence and source appraisal of polycyclic aromatic hydrocarbons (PAHs) in surface waters of the Almendares River, Cuba. Arch Environ Contam Toxicol 69: 143-152. doi: 10.1007/s00244-015-0136-9
[15]	Sarria-Villa R, Ocampo-Duque W, Paez M, et al. (2016) Presence of PAHs in water and sediments of the Colombian Cauca River during heavy rain episodes, and implications for risk assessment. Sci Total Environ 540: 455-465. doi: 10.1016/j.scitotenv.2015.07.020
[16]	Singare PU (2016) Carcinogenic and endocrine disrupting PAHs in the aquatic ecosystem of India. Environ Monit Assess 188: 1-25. doi: 10.1007/s10661-015-4999-z
[17]	Yan J, Liu J, Shi X, et al. (2016) Polycyclic aromatic hydrocarbons (PAHs) in water from three estuaries of China: distribution, seasonal variations and ecological risk assessment. Mar Pollut Bull 109: 471-479. doi: 10.1016/j.marpolbul.2016.05.025
[18]	Manoli E, Samara C, Konstantinou I, et al. (2000) Polycyclic aromatic hydrocarbons in the bulk precipitation and surface waters of Northern Greece. Chemosphere 41: 1845-1855. doi: 10.1016/S0045-6535(00)00134-X
[19]	Mane S, Sundaram S (2014) Studies on some aspects on the biology of green mussel Perna viridis (Linnaeus, 1758) from Versova creek, Mumbai, northwest coast of India. Int Res J Sci Eng 2: 47-50.
[20]	Shirke S, Pinto SM, Kushwaha VK, et al. (2016) Object-based image analysis for the impact of sewage pollution in Malad Creek, Mumbai, India. Environ Monit Assess 188: 95-99. doi: 10.1007/s10661-015-4981-9
[21]	Zeng EY, Yu CC, Tran K (1999) In situ measurements of chlorinated hydrocarbons in the water column off the Palos Verdes Peninsula, California. Environ Sci Technol 33: 392-398. doi: 10.1021/es980561e
[22]	WHO (1998) Polynuclear aromatic hydrocarbons. Guidelines for Drinking-Water Quality, 2nd edition. Addendum to Vol. 2 Health Criteria and Other Supporting Information. World Health Organization, Geneva, pp. 123-152.
[23]	Xiang N, Jiang C, Yang T, et al. (2018) Occurrence and distribution of Polycyclic aromatic hydrocarbons (PAHs) in seawater, sediments and corals from Hainan Island, China. Ecotoxicol Environ Saf 152: 8-15. doi: 10.1016/j.ecoenv.2018.01.006
[24]	Santos E, Souza MRR, Vilela Junior AR, et al. (2018) Polycyclic aromatic hydrocarbons (PAH) in superficial water from a tropical estuarine system: Distribution, seasonal variations, sources and ecological risk assessment. Mar Pollut Bull 127: 352-358. doi: 10.1016/j.marpolbul.2017.12.014
[25]	Nwineewii JD, Marcus AC (2015) Polycyclic Aromatic Hydrocarbons (PAHs) In Surface Water and Their Toxicological Effects in Some Creeks of South East Rivers State (Niger Delta) Nigeria. IOSR J Environ Sci Toxicol Food Technol 9: 27-30.
[26]	Yang D, Qi SH, Zhang Y, et al. (2013) Levels, sources and potential risks of polycyclic aromatic hydrocarbons (PAHs) in multimedia environment along the Jinjiang River mainstream to Quanzhou Bay, China. Mar Pollut Bull 76: 298-306. doi: 10.1016/j.marpolbul.2013.08.016
[27]	Adeniji AO, Okoh OO, Okoh AI (2019) Levels of Polycyclic Aromatic Hydrocarbons in the Water and Sediment of Bufalo River Estuary, South Africa and Their Health Risk Assessment. Arch Environ Contam Toxicol 76: 657-669. doi: 10.1007/s00244-019-00617-w
[28]	Edokpayi JN, Odiyo JO, Popoola OE, et al. (2016) Determination and Distribution of Polycyclic Aromatic Hydrocarbons in Rivers, Sediments and Wastewater Effluents in Vhembe District, South Africa. Int J Environ Res Public Health 13: 387, 1-12. doi: 10.3390/ijerph13040387
[29]	Nekhavhambe TJ., van Ree T, Fatoki OS (2014) Determination and distribution of polycyclic aromatic hydrocarbons in rivers, surface runoff, and sediments in and around Thohoyandou, Limpopo Province, South Africa. Water SA 40: 415-424.
[30]	Eganhouse RP, Simoneit BRT, Kaplan IR (1981) Extractable organic matter in urban stormwater runoff. 2. Molecular characterization. Environ Sci Technol 15: 315-326.
[31]	Hoffman EJ, Mills GL, Latimer JS, et al. (1984) Urban runoff as a source of polycyclic aromatic hydrocarbons to coastal waters. Environ Sci Technol 18: 580-587. doi: 10.1021/es00126a003
[32]	Agarwal T, Khillare P, Shridhar V, et al. (2009) Pattern, sources and toxic potential of PAHs in the agricultural soils of Delhi, India. J Hazard Mater 163: 1033-1039. doi: 10.1016/j.jhazmat.2008.07.058
[33]	Xing XL, Qi S, Zhang J, et al. (2011) Spatial distribution and source diagnosis of polycyclic aromatic hydrocarbons in soils from Chengdu Economic Region, Sichuan Province, western China. J Geochem Explor 110: 146-154. doi: 10.1016/j.gexplo.2011.05.001
[34]	Sprovieri M, Feo ML, Prevedello L, et al. (2007) Heavy metals, polycyclic aromatic hydrocarbons and polychlorinated biphenyls in surface sediments of the Naples harbor (southern Italy). Chemosphere 67: 998-1009. doi: 10.1016/j.chemosphere.2006.10.055
[35]	Ravindra K, Wauters E, Grieken RV (2008) Variation in particulate PAHs levels and their relation with the transboundary movement of the air masses. Sci Total Environ 396: 100-110. doi: 10.1016/j.scitotenv.2008.02.018
[36]	Tobiszewski M, Namiesnik J (2012) PAH diagnostic ratios for the identification of pollution emission sources. Environ Pollut 162: 110-119. doi: 10.1016/j.envpol.2011.10.025
[37]	Cao ZH, Wang YQ, Ma YM, et al. (2005) Occurrence and distribution of polycyclic aromatic hydrocarbons in reclaimed water and surface water of Tianjin, China. J Hazard Mater 122: 51-59. doi: 10.1016/j.jhazmat.2005.04.003
[38]	Boonyatumanond R, Wattayakorn G, Togo A, et al. (2006) Distribution and origins of polycyclic aromatic hydrocarbons (PAHs) in riverine, estuarine, and marine sediments in Thailand. Mar Pollut Bull52: 942-956. doi: 10.1016/j.marpolbul.2005.12.015
[39]	Mostert MMR., Ayoko GA, Kokot S (2010) Application of chemometrics to analysis of soil pollutants. Trends Anal Chem 29: 430-435.
[40]	Mai BX, Qi SH, Zeng EY, et al. (2003) Distribution of polycyclic aromatic hydrocarbons in the coastal region off Macao, China: assessment of input sources and transport pathways using compositional analysis. Environ Sci Technol 37: 4855-4863. doi: 10.1021/es034514k
[41]	Rocher V, Azimi S, Moilleron R, et al. (2004) Hydrocarbons and heavy metals in the different sewer deposits in the Le Marais' catchment (Paris, France): stocks, distributions and origins. Sci Total Environ 323: 107-122. doi: 10.1016/j.scitotenv.2003.10.010
[42]	Wang XC, Sun S, Ma HQ, et al. (2006) Sources and distribution of aliphatic and polyaromatic hydrocarbons in sediments of Jiaozhou Bay, Qingdao, China. Mar Pollut Bull 52: 129-138. doi: 10.1016/j.marpolbul.2005.08.010
[43]	Montuori P, Aurino S, Garzonio F, et al. (2016) Distribution, sources and ecological risk assessment of polycyclic aromatic hydrocarbons in water and sediments from Tiber River and estuary, Italy. Sci Total Environ 566-567: 1254-1267. doi: 10.1016/j.scitotenv.2016.05.183
[44]	Zhang W, Zhang S, Wan C, et al. (2008) Source diagnostics of polycyclic aromatic hydrocarbons in urban road runoff, dust, rain and canopy throughfall. Environ Pollut 153: 594-601. doi: 10.1016/j.envpol.2007.09.004
[45]	Chung MK, Hu R, Cheung KC, et al. (2007) Pollutants in Hongkong soils: polycyclicaromatic hydrocarbons. Chemosphere 67: 464-473. doi: 10.1016/j.chemosphere.2006.09.062
[46]	Li G, Xia X, Yang Z, et al. (2006) Distribution and sources of polycyclic aromatic hydrocarbons in the middle and lower reaches of the Yellow River, China. Environ Pollut 144: 985-993. doi: 10.1016/j.envpol.2006.01.047
[47]	De La Torre-Roche RJ, Lee WY, Campos-Diaz SI (2009) Soil-borne polycyclic aromatic hydrocarbons in El Paso, Texas: analysis of a potential problem in the United States/Mexico border region. J Hazard Mater 163: 946-958. doi: 10.1016/j.jhazmat.2008.07.089
[48]	Akyuz M, Cabuk H (2010) Gaseparticle partitioning and seasonal variation of polycyclic aromatic hydrocarbons in the atmosphere of Zonguldak, Turkey. Sci Total Environ 408: 5550-5558. doi: 10.1016/j.scitotenv.2010.07.063
[49]	Dhananjayan V, Muralidharan S, Peter VR (2012) Occurrence and distribution of polycyclic aromatic hydrocarbons in water and sediment collected along the Harbour Line, Mumbai, India. Int J Oceanogr Article ID 403615, 7.
[50]	Katsoyiannis A, Sweetman AJ, Jones KC (2011) PAH molecular diagnostic ratios applied to atmospheric sources: a critical evaluation using two Decades of source Inventory and air concentration data from the UK. Environ Sci Technol 45: 8897-8906. doi: 10.1021/es202277u
[51]	Pozo K, Perra G, Menchi V, et al. (2011) Levels and spatial distribution of polycyclic aromatic hydrocarbons (PAHs) in sediments from Lenga Estuary, central Chile. Mar Pollut Bull 62: 1572-1576. doi: 10.1016/j.marpolbul.2011.04.037
[52]	Law RJ, Dawes VJ, Woodhead RJ, et al. (1997) Polycyclic aromatic hydrocarbons (PAH) in seawater around England and Wales. Mar Pollut Bull 34: 306-322.
[53]	Barron MG, Podrabsky T, Ogle S, et al. (1999) Are aromatic hydrocarbons the primary determinant of petroleum toxicity to aquatic organisms? Aquat Toxicol 46: 253-268. doi: 10.1016/S0166-445X(98)00127-1
[54]	Agroudy NA, Soliman YA, Hamed MA, et al. (2017) Distribution of PAHs in Water, Sediments Samples of Suez Canal During 2011. J Aquat Pol Toxicol 1: 1-10.
[55]	Pohl A, Kostecki M, Jureczko I, et al. (2018) Polycyclic aromatic hydrocarbons in water and bottom sediments of a shallow, lowland dammed reservoir (on the example of the reservoir Blachownia, South Poland). Arch Environ Prot 44: 10-23.
[56]	US Environmental Protection Agency (USEPA), (2012) Regional screening levels for chemical contaminants at superfund sites. Regional screening table. User's guide. (Access date: November 2012). < http://www.epa.gov/reg3hwmd/risk/human/rb-concentration_table/usersguide.htm > .
[57]	Di Toro DM, McGrath JA, Hansen DJ (2000) Technical basis for narcotic chemicals and polycyclic aromatic hydrocarbon criteria. I. Water and tissue. Environ Toxicol Chem 19: 1951-1970.
[58]	Kalf DF, Crommentuijn T, van de Plassche EJ (1997) Environmental quality objectives for 10 polycyclic aromatic hydrocarbons (PAHs). Ecotoxicol Environ Saf 36: 89-97. doi: 10.1006/eesa.1996.1495
[59]	Canadian Council of Ministers of the Environment (CCME) (2010) Canadian Soil Quality Guidelines, Carcinogenic and Other Polycyclic Aromatic Hydrocarbons (PAHs)-Environmental and Human Health Effects. ISBN 978-1-896997-94-0 PDF.
[60]	Yang B, Xue N, Zhou L, et al. (2012) Risk assessment and sources of polycyclic aromatic hydrocarbons in agricultural soils of Huanghuai plain, China. Ecotoxicol Environ Saf 84: 304-310. doi: 10.1016/j.ecoenv.2012.07.027
[61]	Omayma EA, Sawsan AM, El Nady MM (2016) Application of polycyclic aromatic hydrocarbons in identification of organic pollution in seawater around Alexandria coastal area, Egypt. J Environ Life Sci 1: 39-55.
[62]	Daisey JM, Leyko MA, Kneip TJ (1979) Source identification and allocation of polynuclear aromatic hydrocarbon compounds in the New York City aerosol: methods and applications. In: Jones, P.W., Leber, P. (Eds.), Polynuclear Aromatic Hydrocarbons. Ann Arbor Science, Ann Arbor, pp. 201-215.
[63]	Harrison RM, Smith DJT, Luhana L (1996) Source apportionment of atmospheric polycyclic aromatic hydrocarbons collected from an urban location in Birmingham, UK. Environ Sci Technol 30: 825-832. doi: 10.1021/es950252d
[64]	Rogge WF, Hildemann LM, Mazurek MA, et al. (1993) Source of fine organic aerosol 2. Noncatalyst and catalyst-equipped automobiles and heavy-duty diesel trucks. Environ Sci Technol 27: 636-651.
[65]	Gschwend PM, Hites RA (1981) Fluxes of polycyclic aromatic hydrocarbons to marine and lacustrine sediments in the northeastern United States. Geochimica et Cosmochimica Acta 45: 2359-2367. doi: 10.1016/0016-7037(81)90089-2
[66]	Mitra S, Bianchi TS, Mckee BA, et al. (2002) Black carbon from the Mississippi River: quantities, sources and potential implications for the global carbon cycle. Environ Sci Technol 36: 2296-2302. doi: 10.1021/es015834b
[67]	Masclet P, Bresson MA, Mouvier G (1987) Polycyclic aromatic hydrocarbons emitted by power station, and influence of combustion conditions. Fuel 66: 556-562. doi: 10.1016/0016-2361(87)90163-3

This article has been cited by:

1.	Jin Li, Yongling Cheng, Barycentric rational interpolation method for solving 3 dimensional convection–diffusion equation, 2024, 304, 00219045, 106106, 10.1016/j.jat.2024.106106
2.	Jin Li, Yongling Cheng, Spectral collocation method for convection-diffusion equation, 2024, 57, 2391-4661, 10.1515/dema-2023-0110

Reader Comments

Your name:*

Email:*
© 2021 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

AIMS Environmental Science

1.6 2.9

Metrics

Article views(2331) PDF downloads(134) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Environmental Science

Persistent and carcinogenic polycyclic aromatic hydrocarbons in the north-western coastal marine environment of India

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

Figures and Tables

Other Articles By Authors

Catalog

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

AIMS Environmental Science

Persistent and carcinogenic polycyclic aromatic hydrocarbons in the north-western coastal marine environment of India

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

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

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