An approach based on the pseudospectral method for fractional telegraph equations

Haifa Bin Jebreen; Beatriz Hernández-Jiménez; Haifa Bin Jebreen; Beatriz Hernández-Jiménez

doi:10.3934/math.20231496

AIMS Mathematics

2023, Volume 8, Issue 12: 29221-29238. doi: 10.3934/math.20231496

Previous Article Next Article

Research article Special Issues

An approach based on the pseudospectral method for fractional telegraph equations

Haifa Bin Jebreen ^{1
,
,},
Beatriz Hernández-Jiménez ²

1.
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
2.
Departamento de Economía, Métodos Cuantitativos e Historia Económica, Universidad Pablo de Olavide, 41013 Sevilla, Spain

Received: 16 August 2023 Revised: 08 October 2023 Accepted: 12 October 2023 Published: 26 October 2023
MSC : 54A25, 65M70, 65Bxx, 35L20

We aim to implement the pseudospectral method on fractional Telegraph equation. To implement this method, Chebyshev cardinal functions (CCFs) are considered bases. Introducing a matrix representation of the Caputo fractional derivative (CFD) via an indirect method and applying it via the pseudospectral method helps to reduce the desired problem to a system of algebraic equations. The proposed method is an effective and accurate numerical method such that its implementation is easy. Some examples are provided to confirm convergence analysis, effectiveness and accuracy.

Keywords:

Citation: Haifa Bin Jebreen, Beatriz Hernández-Jiménez. An approach based on the pseudospectral method for fractional telegraph equations[J]. AIMS Mathematics, 2023, 8(12): 29221-29238. doi: 10.3934/math.20231496

Related Papers:

[1]	Batirkhan Turmetov, Valery Karachik . On solvability of some inverse problems for a nonlocal fourth-order parabolic equation with multiple involution. AIMS Mathematics, 2024, 9(3): 6832-6849. doi: 10.3934/math.2024333
[2]	M. J. Huntul, Kh. Khompysh, M. K. Shazyndayeva, M. K. Iqbal . An inverse source problem for a pseudoparabolic equation with memory. AIMS Mathematics, 2024, 9(6): 14186-14212. doi: 10.3934/math.2024689
[3]	Bauyrzhan Derbissaly, Makhmud Sadybekov . Inverse source problem for multi-term time-fractional diffusion equation with nonlocal boundary conditions. AIMS Mathematics, 2024, 9(4): 9969-9988. doi: 10.3934/math.2024488
[4]	Arivazhagan Anbu, Sakthivel Kumarasamy, Barani Balan Natesan . Lipschitz stability of an inverse problem for the Kawahara equation with damping. AIMS Mathematics, 2020, 5(5): 4529-4545. doi: 10.3934/math.2020291
[5]	Yilihamujiang Yimamu, Zui-Cha Deng, Liu Yang . An inverse volatility problem in a degenerate parabolic equation in a bounded domain. AIMS Mathematics, 2022, 7(10): 19237-19266. doi: 10.3934/math.20221056
[6]	Ahmed M.A. El-Sayed, Eman M.A. Hamdallah, Hameda M. A. Alama . Multiple solutions of a Sturm-Liouville boundary value problem of nonlinear differential inclusion with nonlocal integral conditions. AIMS Mathematics, 2022, 7(6): 11150-11164. doi: 10.3934/math.2022624
[7]	Murat A. Sultanov, Vladimir E. Misilov, Makhmud A. Sadybekov . Numerical method for solving the subdiffusion differential equation with nonlocal boundary conditions. AIMS Mathematics, 2024, 9(12): 36385-36404. doi: 10.3934/math.20241726
[8]	Abdelkader Laiadi, Abdelkrim Merzougui . Free surface flows over a successive obstacles with surface tension and gravity effects. AIMS Mathematics, 2019, 4(2): 316-326. doi: 10.3934/math.2019.2.316
[9]	A. M. A. El-Sayed, W. G. El-Sayed, Somyya S. Amrajaa . On a boundary value problem of arbitrary orders differential inclusion with nonlocal, integral and infinite points boundary conditions. AIMS Mathematics, 2022, 7(3): 3896-3911. doi: 10.3934/math.2022215
[10]	Lin Fan, Shunchu Li, Dongfeng Shao, Xueqian Fu, Pan Liu, Qinmin Gui . Elastic transformation method for solving the initial value problem of variable coefficient nonlinear ordinary differential equations. AIMS Mathematics, 2022, 7(7): 11972-11991. doi: 10.3934/math.2022667

Abstract

1. Introduction

Modern problems of natural science lead to the need to generalize the classical problems of mathematical physics, as well as to the formulation of qualitatively new problems, which include non-local problems for differential equations. Among nonlocal problems, problems with integral conditions are of great interest. Integral conditions are encountered in the study of physical phenomena in the case when the boundary of the process flow region is inaccessible for direct measurements. Inverse problems arise in various fields of human activity, such as seismology, mineral exploration, biology, medical visualization, computed tomography, earth remote sensing, spectral analysis, nondestructive control, etc. Various inverse problems for certain types of partial differential equations have been studied in many works. A more detailed bibliography and a classification of problems are found in ^[1,2,3,4,5]. Inverse problems for one-dimensional pseudo-parabolic equations of third-order were studied in ^[6]. The existence and uniqueness of the solution of the inverse problem for the third order pseudoparabolic equation with integral over-determination condition is studied in ^[7]. Khompysh ^[8] investigated the reconstruction of unknown coefficient in pseudo-parabolic inverse problem with the integral over determination condition and studied the uniqueness and existence of solution by means of method of successive approximations. Studies of wave propagation in cold plasma and magnetohydrodynamics also reduce to the partial differential equations of fourth-order. To the study of nonlocal boundary value problems (including integral conditions) for partial differential equations of the fourth-order are devoted large number of works, see, for example, ^[9,10]. It should be noted that boundary value problems with integral conditions are of particular interest. From physical considerations, the integral conditions are completely natural, and they arise in mathematical modelling in cases where it is impossible to obtain information about the process occurring at the boundary of the region of its flow using direct measurements or when it is possible to measure only some averaged (integral) characteristics of the desired quantity.

In this article, we study the an inverse boundary value problem for a fourth order pseudo parabolic equation with periodic and integral condition to identify the time-dependent coefficients along with the solution function theoretically, i.e. existence and uniqueness.

Statement of the problem and its reduction to an equivalent problem. In the domain $D_{T} = \{ (x, t):0\le x\le 1, \, \, 0\le t\le T\}$ , we consider an inverse boundary value problem of recovering the timewise dependent coefficients $p(t)$ in the pseudo-parabolic equation of the fourth-order

$\begin{equation} u_{t} (x, t)-bu_{txx} (x, t)+a(t)u_{xxxx} (x, t) = p(t)u(x, t)+f(x, t) \end{equation}$

(1.1)

with the initial condition

$\begin{equation} u(x, 0)+\delta u(x, T) = \varphi (x)\, \, \, \, \, \, \, \, (0\le x\le 1), \end{equation}$

(1.2)

boundary conditions

$\begin{equation} u(0, t) = u(1, t), u_{x} (0, t) = u_{x} (1, t), u_{xx} (0, t) = u_{xx} (1, t)\, \, \, \, \, \, \, \, \, (0\le t\le T), \end{equation}$

(1.3)

nonlocal integral condition

$\begin{equation} \int _{0}^{1}u(x, t)dx = 0 (0\le t\le T) \end{equation}$

(1.4)

and with an additional condition

$\begin{equation} u(0, t) = \int _{0}^{t}\gamma (\tau )u(1, \tau )d\tau + h(t)\, \, \, \, \, (0\le t\le T), \end{equation}$

(1.5)

where $b > 0$ , $\delta \ge 0$ -given numbers, $a(t) > 0, \, \, f(x, t), \, \, \varphi (x), \, \, \gamma (\tau), h(t)$ -given functions, $u(x, t)$ and $p(t)$ - required functions.

Denote

$\bar{C}^{4, 1} (D_{T} ) = \left\{u(x, t):u(x, t)\in C^{2, 1} (D_{T} ), u_{txx} , u_{xxxx} \in C(D_{T} )\right\}.$

Definition.By the classical solution of the inverse boundary value problem (1.1)-(1.5)we mean the pair $\left\{u(x, t), p(t)\right\}$ functions $u(x, t)\in \bar{C}^{4, 1} \left(D_{T} \right)$ , $p(t)\in C\left[0, T\right]$ satisfying equation (1.1) in $D_{T}$ , condition (1.2) in [0, 1] and conditions (1.3)-(1.5) in [0, T].

Theorem 1. Let be $b > 0, \, \, \delta \ge 0, \, \, \,$ $\, \, \varphi (x)\in C\, [0, 1], \, \, \,$ $f(x, t)\in C(D_{T})$ , $\int _{0}^{1}f(x, t)dx = 0$ , $\, 0 < a(t)\in C\, [0, T], \,$ $h(t)\in C^{1} [0, T],$ $h(t)\ne 0\,$ $\, (0\le t\le T)$ , $\gamma (t)\in C\, [0, T], \,$ $\delta \gamma (t) = 0$ $\, (0\le t\le T)$ and

$\int _{0}^{1}\varphi (x)dx = 0, \varphi (0) = h(0)+\delta h(T).$

Then the problem of finding a solution to problem (1.1)-(1.5) is equivalent to the problem of determining the functions $u(x, t)\in \bar{C}^{4, 1} (D_{T})$ and $p(t)\in C[0, T]$ , from (1.1)-(1.3) and

$\begin{equation} u_{xxx} (0, t) = u_{xxx} (1, t)\, \, \, \, \, \, \, (0\le t\le T), \end{equation}$

(1.6)

$\gamma (t)u(1, t)+h'(t)-bu_{txx} (0, t)+a(t)u_{xxxx} (0, t) =$

$\begin{equation} = p(t)\left(\int _{0}^{t}\gamma (\tau )u(1, \tau )d\tau + h(t)\, \right)\, +f(0, t)\, \, \, (0\le t\le T). \end{equation}$

(1.7)

Proof. Let be $\{ u(x, t), \, p(t)\}$ is a classical solution to problem (1.1)-(1.5). Integrating equation (1.1) with respect to $x$ from 0 to 1, we get:

$\frac{d}{dt} \int _{0}^{1}u(x, t)dx -b(u_{tx} (1, t)-u_{tx} (0, t))+a(t)(u_{xxx} (1, t)-u_{xxx} (0, t)) =$

$\begin{equation} = p(t)\int _{0}^{1}u(x, t)dx +\int _{0}^{1}f(x, t)dx \, \, \, (0\le t\le T). \end{equation}$

(1.8)

Assuming that $\int _{0}^{1}f(x, t)dx = 0$ , taking into account (1.3) and (1.4), we arrive at the fulfillment of (1.6).

Further, considering $h(t)\in C^{1} [0, T]$ and differentiating with respect to $t$ (1.5), we get:

$\begin{equation} u_{t} (0, t) = \gamma (t)u(1, t)+h'(t)\, \, \, \, \, (0\le t\le T) \end{equation}$

(1.9)

Substituting $x = 0$ into equation (1.1), we have:

$\begin{equation} u_{t} (0, t)-bu_{txx} (0, t)+a(t)u_{xxxx} (0, t) = p(t)u(0, t)+f(0, t)\, \, \, (0\le t\le T). \end{equation}$

(1.10)

Now, suppose that $\{ u(x, t), \, p(t)\}$ is a solution to problem (1.1)-(1.3), (1.6), (1.7). Then from (1.8), taking into account (1.3) and (1.6), we find:

$\begin{equation} \frac{d}{dt} \int _{0}^{1}u(x, t)dx -p(t)\int _{0}^{1}u(x, t)dx = 0\, \, \, (0\le t\le T). \end{equation}$

(1.11)

Due to (1.2) and $\int _{0}^{1}\varphi (x)dx = 0$ , it's obvious that

$\begin{equation} \int _{0}^{1}u(x, 0)dx +\delta \int _{0}^{1}u(x, T)dx = \int _{0}^{1}\varphi (x)dx = 0. \end{equation}$

(1.12)

Obviously, the general solution(1.11) has the form:

$\begin{equation} \int _{0}^{1}u(x, t)dx = ce^{-\int _{0}^{t}p(\tau )d\tau } \, \, \, \, \, (0\le t\le T). \end{equation}$

(1.13)

From here, taking into account (1.12), we obtain:

$\begin{equation} \int _{0}^{1}u(x, 0)dx +\delta \int _{0}^{1}u(x, T)dx = c(1+\delta e^{-\int _{0}^{T}p(\tau )d\tau } ) = 0. \end{equation}$

(1.14)

By virtue of $\delta \ge 0$ , from (1.14) we get that $c = 0$ , and substituting into (1.13) we conclude, that $\int _{0}^{1}u(x, t)dx = 0\, \, \, \, (0\le t\le T)$ . Therefore, condition (1.4) is also satisfied.

Further, from (1.7) and (1.10), we obtain:

$\frac{d}{dt} \left[u(0, t)-\, \left(\int _{0}^{t}\gamma (\tau )u(1, \tau )d\tau + h(t)\, \right)\, \right] =$

$\begin{equation} = p(t)\left[u(0, t)-\, \left(\int _{0}^{t}\gamma (\tau )u(1, \tau )d\tau + h(t)\, \right)\, \right]\, \, (0\le t\le T). \end{equation}$

(1.15)

Let introduce the notation:

$\begin{equation} y(t)\equiv u(0, t)-\, \left(\int _{0}^{t}\gamma (\tau )u(1, \tau )d\tau + h(t)\, \right)\, \, (0\le t\le T) \end{equation}$

(1.16)

and rewrite the last relation in the form:

$\begin{equation} y'(t)+p(t)y(t) = 0\, \, \, \, (0\le t\le T). \end{equation}$

(1.17)

From (1.16), taking into account (1.2), $\delta \gamma (t) = 0$ $\, (0\le t\le T)$ and $\varphi (0) = h(0)+\delta h(T)$ , it is easy to see that

$y(0)+\delta y(T) = \, u(0, 0)-h(0)+\delta \left[u(0, T)-\, \left(\int _{0}^{T}\gamma (\tau )u(1, \tau )d\tau + h(T)\, \right)\, \right] = u(0, 0)+$

$\begin{equation} +\delta u(0, T)-(h(0)+\delta h(T)\, )-\delta \int _{0}^{T}\gamma (\tau )u(1, \tau )d\tau = \varphi (0)-\left(h(0)+\delta h(T)\right) = 0. \end{equation}$

(1.18)

Obviously, the general solution (1.17) has the form:

$\begin{equation} y(t) = ce^{-\int _{0}^{t}p(\tau )d\tau } \, \, \, \, \, (0\le t\le T). \end{equation}$

(1.19)

From here, taking into account (1.18), we obtain:

$\begin{equation} y(0)+\delta y(T) = c(1+\delta e^{-\int _{0}^{T}\frac{a_{0} (\tau )}{a_{1} (\tau )} d\tau } ) = 0. \end{equation}$

(1.20)

By virtue of $\delta \ge 0$ , from (1.20) we get that $c = 0$ , and substituting into (1.19) we conclude that $y(t) = 0\, \, \, \, (0\le t\le T)$ . Therefore, from (1.16) it is clear that the condition (1.5). The theorem has been proven.

2. The existence and uniqueness of the classical solution of the inverse boundary value problem.

It is known ^[5] that the system

$\begin{equation} 1, \cos \lambda _{1} x, \sin \lambda _{1} x, ..., \cos \lambda _{k} x, \sin \lambda _{k} x, ... \end{equation}$

(2.1)

forms the basis of $L_{2} (0, 1)$ , where $\lambda _{k} = 2k\pi \, \, \, \, (k = 0, 1, ...)$ .

Since system (2.1) forms a basis in $L_{2} (0, 1)$ , it is obvious that for each solution $\{ u(x, t), a(t)\} \,$ problem (1.1)–(1.3), (1.6), (1.7):

$\begin{equation} u(x, t) = \sum\limits _{k = 0}^{\infty }u_{1k} (t)\cos \lambda {}_{k} x +\sum\limits _{k = 1}^{\infty }u_{2k} (t)\sin \lambda {}_{k} x \, \, \, (\lambda _{k} = 2\pi k), \end{equation}$

(2.2)

where

$u_{10} (t) = \int _{0}^{1}u(x, t)dx\, \, \, , u_{1k} (t) = 2\int _{0}^{1}u(x, t)\cos \lambda _{k} xdx\, \, \, (k = 1, 2, ...)\, ,$

$u_{2k} (t) = 2\int _{0}^{1}u(x, t)\sin \lambda _{k} xdx\, \, \, (k = 1, 2, ...).$

Applying the formal scheme of the Fourier method, to determine the desired coefficients $u_{1k} (t)\, \, \, \, (k = 0, 1, ...)$ and $u_{2k} (t)\, \, \, \, (k = 1, 2, ...)$ functions $u(x, t)$ from (1.1) and (1.2) we get:

$\begin{equation} u''_{10} (t) = F_{10} (t;u, p)\, \, \, \, (0\le t\le T) , \end{equation}$

(2.3)

$\begin{equation} (1+b\lambda _{k}^{2} )u'_{ik} (t)+a(t)\lambda _{k}^{4} u_{ik} (t) = F_{ik} (t;u, p)\, \, (i = 1, 2;\, 0\le t\le T;\, \, k = 1, 2, ...) , \end{equation}$

(2.4)

$\begin{equation} u_{10} (0)+\delta u_{10} (T) = \varphi _{10} \, \, \, , \end{equation}$

(2.5)

$\begin{equation} u_{ik} (0)+\delta u_{ik} (T) = \varphi _{ik} \, \, \, (\, i = 1, 2;\, k = 1, 2, ...), \end{equation}$

(2.6)

where

$F_{1k} (t;u, a, b) = p(t)u_{1k} (t)+f_{1k} (t)\, \, \, \, (k = 0, 1, ...)\, \, ,$

$f_{10} (t) = \int _{0}^{1}f(x, t)dx \, \, , f_{1k} (t) = 2\int _{0}^{1}f(x, t)\cos \lambda _{k} xdx \, \, (k = 1, 2, ...)\, \, , \, \, \,$

$\varphi _{10} = \int _{0}^{1}\varphi (x)dx, \, \, \, \, \, \, \, \, \, \, \varphi _{1k} = 2\int _{0}^{1}\varphi (x)\cos \lambda _{k} xdx\, \, \, \, \, \, \, \, (k = 1, 2, ...)\, ,$

$F_{2k} (t;u, a, b) = p(t)u_{2k} (t)+f_{2k} (t), \, \,$

$f_{2k} (t) = 2\int _{0}^{1}f(x, t)\sin \lambda _{k} xdx \, \, \, \, \, (k = 1, 2, ...)\, \, , \varphi _{2k} = 2\int _{0}^{1}\varphi (x)\sin \lambda _{k} xdx\, \, \, \, \, \, (k = 1, 2, ...)\, .$

Solving problem (2.3)-(2.6), we find:

$\begin{equation} u_{10} (t) = (1+\delta )^{-1} \left(\varphi _{10} -\delta \int _{0}^{T}F_{0} (\tau ;u, p)d\tau \right)+\int _{0}^{t}F_{10} (\tau ;u, p)d\tau \, \, (0\le t\le T), \end{equation}$

(2.7)

$u_{ik} (t) = \frac{e^{-\int _{0}^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;dss } } \varphi _{ik} +\frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{t}F_{ik} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, -$

$\begin{equation} -\frac{\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } } \frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{T}F_{ik} (\tau ;u, p) e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, \, \, (i = 1, 2;0\le t\le T;k = 1, 2, ...). \end{equation}$

(2.8)

After substituting the expression $u_{1k} (t)\, \, \, \, (k = 0, \, 1, ...),$ $u_{2k} (t)\, \, \, \, (k = 1, 2, ...)$ in (2.2), to define a component $u(x, t)$ solution of problem (1.1)-(1.3), (1.6), (1.7), we obtain:

$u(x, t) = (1+\delta )^{-1} \left(\varphi _{0} -\delta \int _{0}^{T}F_{0} (\tau ;u, p)d\tau \right)+\int _{0}^{t}F_{0} (\tau ;u, p)d\tau \, +$

$+\sum\limits _{k = 1}^{\infty }\left\{\frac{e^{-\int _{0}^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;dss } } \varphi _{1k} \right. _{k} +\frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{t}F_{1k} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau -$

$-\left. \frac{\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } } \frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{T}F_{1k} (\tau ;u, p) e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \right\}\, \, \cos \lambda _{k} x+$

$+\sum\limits _{k = 1}^{\infty }\left\{\frac{e^{-\int _{0}^{t}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;dss } } \varphi _{2k} \right. _{k} +\frac{1}{1+b\lambda\; _{k}^{2} } \int _{0}^{t}F_{2k} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;ds } d\tau -$

$\begin{equation} -\left. \frac{\delta e^{-\int _{0}^{T}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s){\lambda\; _{k}^{4}} }{1+b{\lambda\; _{k}^{2}} } \;ds } } \frac{1}{1+b\lambda\; _{k}^{2} } \int _{0}^{T}F_{2k} (\tau ;u, p) e^{-\int _{\tau }^{t}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;ds } d\tau \right\}\, \, \sin \lambda\; _{k} x. \end{equation}$

(2.9)

Now from (1.7), taking into account (2.2), we have:

$p(t) = \left[h(t)\right]^{-1} \left\{{\mathop{h'(t)-f(0, t)+\gamma (t)u_{10} (t)-}\limits_{}} \right. p(t)\int _{0}^{t}\gamma (\tau )u_{{}_{10} } (\tau )d\tau +$

$\begin{equation} +\sum\limits _{k = 1}^{\infty }\left({\mathop{b\lambda _{k}^{2} u'_{1k} (t)+a(t)\lambda _{k}^{4} u_{1k} (t)+\gamma (t)u_{1k} (t)}\limits_{}} \right. -\left. p(t)\int _{0}^{t}\gamma (\tau )u_{{}_{1k} } (\tau )d\tau \right). \end{equation}$

(2.10)

Further, from (2.4), taking into account (2.8), we obtain:

$b\lambda _{k}^{2} u'_{1k} (t)+a(t)\lambda _{k}^{4} u_{1k} (t)+\gamma (t)u_{1k} (t) = F_{1k} (t;u, p)-u'_{1k} (t)+\gamma (t)u_{1k} (t) =$

$= \, \frac{b\lambda _{k}^{2} }{1+b\lambda _{k}^{2} } F_{1k} (t;u, p)+\left(\frac{a(t)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } +\gamma (t)\right)u_{1k} (t) =$

$\, = \frac{b\lambda\; _{k}^{2} }{1+b\lambda\; _{k}^{2} } F_{k} (t;u, p)+\left(\frac{a(t)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } +\gamma (t)\right)\left[\frac{e^{-\int _{0}^{t}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda\; _{k}^{4} }{1+b\lambda\; _{k}^{2} } \;dss } } \varphi _{{}_{1k} } \right. +$

$+\frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{t}F_{1k} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, -$

$\begin{equation} -\left. \frac{\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } } \frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{T}F_{1k} (\tau ;u, p) e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, \right]\, \, (0\le t\le T;k = 1, 2, ...). \end{equation}$

(2.11)

$p(t) = \left[h(t)\right]^{-1} \left\{{\mathop{h'(t)-f(0, t)+}\limits_{}} \right. \,$

$+\gamma (t))\left[(1+\delta )^{-1} \left(\varphi _{10} -\delta \int _{0}^{T}F_{0} (\tau ;u, p)d\tau \right)+\int _{0}^{t}F_{10} (\tau ;u, p)d\tau \right]-$

$-p(t)\int _{0}^{t}\gamma (\tau )u_{{}_{10} } (\tau )d\tau +\sum\limits _{k = 1}^{\infty }\left[\frac{b\lambda _{k}^{2} }{1+b\lambda _{k}^{2} } F_{1k} (t;u, p)\right. +$

$+\, \left(\frac{a(t)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } +\gamma (t)\right)\left[\frac{e^{-\int _{0}^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } }{1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} }\; dss } } \varphi _{1k} \right. +\frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{t}F_{1k} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, -$

$+\frac{1}{1+b\lambda _{k}^{2} } \int _{0}^{t}F_{1k} (\tau ;u, p)e^{-\int _{\tau }^{t}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } d\tau \, -$

$\begin{equation} \left. \left. +p(t)\int _{0}^{t}\gamma (\tau )u_{1k} (\tau )d\tau \, \right]\, \, \right\}. \end{equation}$

(2.12)

Thus, the solution of problem (1.1)–(1.3), (1.6), (1.7)is reduced to the solution of system (2.9), (2.12) with respect to unknown functions $u(x, t)$ and $p(t)$ .

To study the question of the uniqueness of the solution of problem (1.1)–(1.3), (1.6), (1.7) the following plays an important role.

Lemma 1. If $\{ u(x, t), \, p(t)\}$ -any solution of problem (1.1)–(1.3), (1.6), (1.7), then the functions

$u_{10} (t) = \int _{0}^{1}u(x, t)dx\, \, \, , u_{1k} (t) = 2\int _{0}^{1}u(x, t)\cos \lambda _{k} xdx\, \, \, (k = 1, 2, ...)\, ,$

$u_{2k} (t) = 2\int _{0}^{1}u(x, t)\sin \lambda _{k} xdx\, \, \, (k = 1, 2, ...)$

satisfy the system consisting of equations (27), (28) on $[0, T]$ .

It is obvious that if $u_{10} (t) = \int _{0}^{1}u(x, t)dx\, \, \, ,$ $u_{1k} (t) = 2\int _{0}^{1}u(x, t)\cos \lambda _{k} xdx\, \, \, (k = 1, 2, ...)\, ,$ $u_{2k} (t) = 2\int _{0}^{1}u(x, t)\sin \lambda _{k} xdx\, \, \, (k = 1, 2, ...)$ is a solution to system (2.7), (2.8), then the pair $\left\{u(x, t), p(t)\right\}$ functions $u(x, t) = \sum _{k = 0}^{\infty }u_{1k} (t)\cos \lambda {}_{k} x +\sum _{k = 1}^{\infty }u_{2k} (t)\sin \lambda {}_{k} x \, \, \, (\lambda _{k} = 2\pi k)$ and $p(t)$ is a solution to system (2.9), (2.12).

Consequence. Let system (29), (32) have a unique solution. Then problem (1.1)–(1.3), (1.6), (1.7) cannot have more than one solution, i.e. if problem (1.1)-(1.3), (1.6), (1.7) has a solution, then it is unique.

In order to study the problem (1.1)–(1.3), (1.6), (1.7) consider the following spaces.

Denote by $B_{2, T}^{\alpha }$ ^[6] the set of all functions of the form

$u(x, t) = \sum\limits _{k = 0}^{\infty }u_{1k} (t)\cos \lambda {}_{k} x +\sum\limits _{k = 1}^{\infty }u_{2k} (t)\sin \lambda {}_{k} x \, \, \, (\lambda _{k} = 2\pi k),$

considered in $D_{T}$ , where each of the functions $u_{1k} (t)\, \, \, \, (k = 0, \, 1, ...),$ $u_{2k} (t)\, \, \, \, (k = 1, 2, ...)$ continuous on $[0, T]$ and

$J(u) = \left\| u_{10} (t)\right\| _{C[0, T]} +\left\{\sum\limits _{k = 1}^{\infty }\left(\lambda _{k}^{\alpha } \left\| u_{1k} (t)\right\| _{C[0, T]} \right)^{2} \right\}^{\frac{1}{2} } +\left\{\sum\limits _{k = 1}^{\infty }\left(\lambda _{k}^{\alpha } \left\| u_{2k} (t)\right\| _{C[0, T]} \right)^{2} \right\}^{\frac{1}{2} } < +\infty ,$

$\alpha \ge 0$ . We define the norm in this set as follows:

$\left\| u(x, t)\right\| _{B_{2, T}^{\alpha } } = J(u).$

Through $E_{T}^{\alpha}$ denote the space $B_{2, T}^{\alpha } \times C[0, T]$ vector - functions $z(x, t) = \{ u(x, t), \, p(t)\}$ with norm

$\left\| z(x, t)\right\| _{E_{T}^{\alpha } } = \left\| u(x, t)\right\| _{B_{2, T}^{\alpha } } +\left\| p(t)\right\| _{C[0, T]} .$

It is known that $B_{2, T}^{\alpha }$ and $E_{T}^{\alpha }$ are Banach spaces.

Now consider in space $E_{T}^{5}$ operator

$\Phi (u, p) = \left\{\Phi _{1} (u, p), \Phi _{2} (u, p)\right\}\, ,$

operator

$\Phi _{1} (u, p)) = \tilde{u}(x, t)\equiv \sum\limits _{k = 0}^{\infty }\tilde{u}_{1k} (t)\cos \lambda {}_{k} x +\sum\limits _{k = 1}^{\infty }\tilde{u}_{2k} (t)\sin \lambda {}_{k} x \, \, , \Phi _{2} (u, p) = \tilde{p}(t),$

$\tilde{u}_{10} (t), \tilde{u}_{ik} (t)\, \, \, \, \, (i = 1, 2;\, \, k = 1, 2, ...), \, \, \tilde{p}(t)$ are equal to the right-hand sides of (2.7), (2.8) and (2.12), respectively.

It is easy to see that

$1+b\lambda _{k}^{2} > b\lambda _{k}^{2} , \, \, \, \, 1+\delta \ge 1, . 1+\delta e^{-\int _{0}^{T}\frac{a(s)\lambda _{k}^{4} }{1+b\lambda _{k}^{2} } \;ds } \ge 1.$

Then, we have:

$\begin{equation} \left\| \tilde{u}_{0} (t)\right\| _{C[0, T]} \le \left|\varphi _{10} \right|+(1+\delta )\sqrt{T} \left(\int _{0}^{T}\left|f_{10} (\tau )\right|^{2} d\tau \right)^{\frac{1}{2} } +(1+\delta )T\left\| p(t)\right\| _{C[0, T]} \left\| u_{10} (t)\right\| _{C[0, T]} \, , \end{equation}$

(2.13)

$\left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{5} \left\| \tilde{u}_{ik} (t)\right\| _{C[0, T]} )^{2} \right)^{\frac{1}{2} } \le \sqrt{3} \left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{5} \left|\varphi _{ik} \right|)^{2} \right)^{\frac{1}{2} } \\+\frac{\sqrt{3} (1+\delta )}{b} \sqrt{T} \left(\int _{0}^{T}\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{3} \left|f_{ik} (\tau )\right|)^{2} d\tau \right)^{\frac{1}{2} } +\,$

$\begin{equation} +\frac{\sqrt{3} (1+\delta )}{b} T\left\| p(t)\right\| _{C[0, T]} \left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{5} \left\| u_{ik} (t)\right\| _{C[0, T]} )^{2} \right)^{\frac{1}{2} } \, \, (i = 1, 2) , \end{equation}$

(2.14)

$\left\| \tilde{p}(t)\right\| _{C[0, T]} \le \left\| [h(t)]^{-1} \right\| _{C[0, T]} \left\{\left\| h'(t)-f(0, t)\right\| _{C[0, T]} +\right.$

$+\left\| \gamma (t)\right\| _{C[0, T]} \left[\left|\varphi _{0} \right|+(1+\delta )\sqrt{T} \left(\int _{0}^{T}\left|f_{0} (\tau )\right|^{2} d\tau \right)^{\frac{1}{2} } +(1+\delta )T\left\| p(t)\right\| _{C[0, T]} \left\| u_{0} (t)\right\| _{C[0, T]} \right]+$

$+T\left\| \gamma (t)\right\| _{C[0, T]} \left\| p(t)\right\| _{C[0, T]} \left\| u_{{}_{10} } (t)\right\| _{C[0, T]} +$

$+\left(\sum\limits _{k = 1}^{\infty }\lambda _{k}^{-2} \right)^{\frac{1}{2} } \left[\left(\sum\limits _{k = 1}^{\infty }\left(\lambda _{k} \left\| f_{1k} (t)\right\| \right)_{C[0, T]}^{2} \right)^{\frac{1}{2} } \right. +\left\| p(t)\right\| _{C[0, T]} \left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{3} \left\| u_{1k} (t)\right\| _{C[0, T]} )^{2} \right)^{\frac{1}{2} } +$

$\begin{equation*} +\left(\left\| \gamma (t)\right\| _{C[0, T]} +\frac{1}{b} \left\| a(t)\right\| _{C[0, T]} \right)\left[\left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{3} \left|\varphi _{1k} \right|)^{2} \right)^{\frac{1}{2} } \right. \, +\, \frac{\sqrt{T} (1+\delta )}{b} \left(\int _{0}^{T}\sum\limits _{k = 1}^{\infty }(\lambda _{k} \left|f_{1k} (\tau )\right|)^{2} d\tau \right)^{\frac{1}{2} } + \end{equation*}$

$\begin{equation} +\, \frac{T(1+\delta )}{b} \left\| p(t)\right\| _{C[0, T]} \left. \left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{5} \left\| u_{1k} (t)\right\| _{C[0, T]} )^{2} \right)^{\frac{1}{2} } \right]\, +\left. \\+T\left\| \gamma (t)\right\| _{C[0, T]} \left\| p(t)\right\| _{C[0, T]} \left. \left(\sum\limits _{k = 1}^{\infty }(\lambda _{k}^{5} \left\| u_{1k} (t)\right\| _{C[0, T]} )^{2} \right)^{\frac{1}{2} } \right]\right\} . \end{equation}$

(2.15)

Let us assume that the data of problem (1.1)–(1.3), (1.6), (1.7) satisfy the following conditions:

$1. \varphi (x)\in W_{2} {}^{(5)} (0, 1), \, \, \, \, \varphi (0) = \varphi (1), \varphi '(0) = \varphi '(1), \,$

$\varphi ''(0) = \varphi ''(1), \, \, \varphi '''(0) = \varphi '''(1), \, \, \varphi ^{(4)} (0) = \varphi ^{(4)} (1);$

$2. f(x, t), \, \, \, f_{x} (x, t), \, \, \, f_{xx} (x, t)\in C(D_{T} ), \, \, \, f_{xxx} (x, t)\in L_{2} (D_{T} )\, , \,$

$\, f(0, t) = \, \, \, f(1, t)\, \, , f_{x} (0, t) = \, \, \, f_{x} (1, t)\, , \, \, f_{xx} (0, t) = \, \, \, f_{xx} (1, t)\, \, \, \, \, \, (0\le t\le T);$

$3. b > 0, \, \, \delta \ge 0, \, \, \gamma (t), \, \, \, a(t)\in C[0, T], \, \, \, \, h(t)\in C^{1} [0, T], \, \, \, \, \, h(t)\ne 0\, \, \, \, \, \, (0\le t\le T).$

Then from (2.10)–(2.12), we have:

$\begin{equation} \left\| \tilde{u}(x, t)\right\| _{B_{2, T}^{5} } \le A_{1} (T)+B_{1} (T)\left\| p(t)\right\| _{C[0, T]} \left\| u(x, t)\right\| _{B_{2, T}^{5} } , \end{equation}$

(2.16)

$\begin{equation} \left\| \tilde{p}(t)\right\| _{C[0, T]} \le A_{2} (T)+B_{2} (T)\left\| p(t)\right\| _{C[0, T]} \left\| u(x, t)\right\| _{B_{2, T}^{5} } , \end{equation}$

(2.17)

where

$A_{1} (T) = \left\| \varphi (x)\right\| _{L_{2} (0, 1)} +(1+\delta )\sqrt{T} \left\| f(x, t)\right\| _{L_{2} (D_{T} )} +2\sqrt{3} \left\| \varphi ^{(5)} (x)\right\| _{L_{2} (0, 1)} +$

$\, +\frac{2\sqrt{3} }{b} (1+\delta )\sqrt{T} \left\| f_{xxx} (x, t)\right\| _{L_{2} (D_{T} )} , \, \, \, \, \, \, \, \, \, \, B_{1} (T) = (1+\delta )\left(1+\frac{\sqrt{3} }{b} \right)T,$

$A_{2} (T) = \left\| [h(t)]^{-1} \right\| _{C[0, T]} \left\{\left\| h'(t)-f(0, t)\right\| _{C[0, T]} +\right.$

$+\left\| \gamma (t)\right\| _{C[0, T]} \, \left({\mathop{\left\| \varphi (x)\right\| _{L_{2} (0, 1)} +(1+\delta )\sqrt{T} \left\| f(x, t)\right\| _{L_{2} (D_{T} )} }\limits_{}} \right)+\,$

$+\, \left(\sum\limits _{k = 1}^{\infty }\lambda _{k}^{-2} \right)^{\frac{1}{2} } \left[\left\| \left\| f_{x} (x, t)\right\| _{C[0, T]} \right\| _{L_{2} (0, 1)} \right. \, +$

$\, \, \, \, \, \, \, \, \, \, \, \, \, \, +\, \, \left(\left\| \gamma (t)\right\| _{C[0, T]} +\frac{1}{b} \left\| a(t)\right\| _{C[0, T]} \right)\, \left. \left. \left(\left\| \varphi ^{(3)} (x)\right\| _{L_{2} (0, 1)} +\frac{\sqrt{T} (1+\delta )}{b} \left\| f_{x} (x, t)\right\| _{L_{2} (D_{T} )} \right)\right]\right\},$

$\begin{array}{l} {B_{2} (T) = \left\| [h(t)]^{-1} \right\| _{C[0, T]} \left(\sum _{k = 1}^{\infty }\lambda _{k}^{-2} \right)^{\frac{1}{2} } \left[\left(\left\| \gamma (t)\right\| _{C[0, T]} +\frac{1}{b} \left\| a(t)\right\| _{C[0, T]} \right)\right. \frac{T(2+\delta )}{b} +} \\ {\, \, } \end{array}$

$\left. {\mathop{+T\left\| \gamma (t)\right\| _{C[0, T]} +1}\limits_{}} \right].$

From inequalities (2.16), (2.17) we conclude:

$\begin{equation} \left\| u(x, t)\right\| _{B_{2, T}^{5} } +\left\| \tilde{p}(t)\right\| _{C[0, T]} \le A(T)+B(T)\left\| p(t)\right\| _{C[0, T]} \left\| u(x, t)\right\| _{B_{2, T}^{5} } , \end{equation}$

(2.18)

$A(T) = A_{1} (T)+A_{2} (T), \, \, B(T) = B_{1} (T)+B_{2} (T).$

We can prove the following theorem.

Theorem 2. Let conditions 1-3 be satisfied and

$\begin{equation} (A(T)+2)^{2} B(T) < 1. \end{equation}$

(2.19)

Then problem (1.1)–(1.3), (1.6), (1.7) has in $K = K_{R} (\left\| z\right\| _{E_{T}^{5} } \le R = A(T)+2)$ in the space $E_{T}^{5}$ only one solution.

Proof. In space $E_{T}^{5}$ consider the equation

$\begin{equation} z = \Phi z, \end{equation}$

(2.20)

where $z = \{ u, p\}$ , components $P$ $\Phi _{1} (u, p), \Phi _{2} (u, p)$ of operators $\Phi (u, p)$ are defined by the right-hand sides of equations (2.9) and (2.12).

Consider the operator $\Phi (u, p)$ in a ball $K = K_{R}$ from $E_{T}^{5} \,$ . Similarly to (2.18) we obtain that for any $z = \{ u, p\}$ , $z_{1} = \{ u_{1}, p_{1} \}$ , $z_{2} = \{ u_{2}, p_{2} \} \in K_{R}$ :

$\begin{equation} \left\| \Phi z\right\| _{E_{T}^{5} } \le A(T)+B(T)\left\| p(t)\right\| _{C[0, T]} \left\| u(x, t)\right\| _{B_{2, T}^{5} } , \end{equation}$

(2.21)

$\begin{equation} \left\| \Phi z_{1} -\Phi z_{2} \right\| _{E_{T}^{5} } \le B(T)R(\left\| p_{1} (t)-p_{2} (t)\right\| _{C[0, T]} +\left\| u_{1} (x, t)-u_{2} (x, t)\right\| _{B_{2, T}^{5} } ). \end{equation}$

(2.22)

Then from estimates (2.21), (2.22), taking into account (2.19), it follows that the operator $\Phi$ acts in a ball $K = K_{R}$ and is contractive. Therefore, in the ball $K = K_{R}$ operator $\Phi$ has a single fixed point $\{ u, p\}$ , which is the only one in the ball $K = K_{R}$ solution of equation (2.20), i.e. is the only one solution in the ball $K = K_{R}$ of system (2.9), (2.12) in the ball.

Functions $u(x, t)$ , as an element of space $B_{2, T}^{5}$ is continuous and has continuous derivatives $\, u_{x} (x, t), \, u_{xx} (x, t)$ , $\, u_{xxx} (x, t), \, u_{xxxx} (x, t)$ in $D_{T}$ .

From (2.4), it is easy to see that

$\left(\sum\limits _{k = 1}^{\infty }\left(\lambda _{k}^{} \left\| u'_{ik} (t)\right\| _{C[0, T]} \right)^{2} \right)^{\frac{1}{2} } \le \frac{\sqrt{2} }{b} \left\| a(t)\right\| _{C[0, T]} \left(\sum\limits _{k = 1}^{\infty }\left(\lambda _{k}^{5} \left\| u_{ik} (t)\right\| _{C[0, T]} \right)^{2} \right)^{\frac{1}{2} } +$

$+\frac{\sqrt{2} }{b} \left\| \left\| f_{x} (x, t)+p(t)u_{x} (x, t)\right\| _{C[0, T]} \right\| _{L_{2} (0, 1)} \, \, (i = 1, 2).$

Hence it follows that $u_{t} (x, t)$ and $u_{txx}$ continuous in $D_{T}$ .

It is easy to check that equation (1.1) and conditions (1.2), (1.3), (1.6), (1.7) are satisfied in the usual sense. Consequently, $\{ u(x, t), p(t)\}$ is a solution to problem (1.1)–(1.3), (1.6), (1.7). By the corollary of Lemma 1, it is unique in the ball $K = K_{R}$ . The theorem has been proven.

With the help of Theorem 1, the unique solvability of the original problem (1.1)–(1.5) immediately follows from the last theorem.

Theorem 3. Let all the conditions of Theorem 1 be satisfied, $\int _{0}^{1}f(x, t)dx = 0\, \, (0\le t\le T)$ , $\delta \gamma (t) = 0$ $\, (0\le t\le T)$ and the matching condition is met:

$\int _{0}^{1}\varphi (x)dx = 0, \varphi (0) = h(0)+\delta h(T).$

Then problem (1.1)–(1.5) has in the ball $K = K_{R} (\left\| z\right\| _{E_{T}^{5} } \le R = A(T)+2)$ from $E_{T}^{5} \,$ the only classical solution.

3. Conclusions

The article considered an inverse boundary value problem with a periodic and integral condition, when the unknown coefficient depends on time for a linear pseudoparabolic equation of the fourth order. An existence and uniqueness theorem for the classical solution of the problem is proved.

Conflict of interest

The authors have declared no conflict of interest.

References

[1]	M. Lakestani, B. N. Saray, Numerical solution of Telegraph equation using interpolating scaling functions, Comput. Math. Appl., 60 (2010) 1964–1972. https://doi.org/10.1016/j.camwa.2010.07.030
[2]	X. Yang, H. Zhang, J. Tang, The OSC solver for the fourth-order sub-diffusion equation with weakly singular solutions, Comput. Math. Appl., 82 (2021), 1–12.
[3]	M. Dehghan, M. Lakestani, The use of Chebyshev cardinal functions for solution of the second-order one-dimensional telegraph equation, Numer. Meth. Part. D. E., 25 (2009), 931–938. https://doi.org/10.1002/num.20382 doi: 10.1002/num.20382
[4]	M. Dehghan, A. Shokri, A numerical method for solving the hyperbolic telegraph equation, Numer. Math. Part. D. E., 24 (2008), 1080–1093. https://doi.org/10.1002/num.20306 doi: 10.1002/num.20306
[5]	O. Nikan, Z. Avazzadeh, J. A. T. Machado, M. N. Rasoulizadeh, An accurate localized meshfree collocation technique for the telegraph equation in propagation of electrical signals, Eng. Comput., 39 (2023), 2327–2344. https://doi.org/10.1007/s00366-022-01630-9 doi: 10.1007/s00366-022-01630-9
[6]	S. Sharifi, J. Rashidinia, Numerical solution of hyperbolic telegraph equation by cubic B-spline collocation method, Appl. Math. Comput., 281 (2016), 28–38. https://doi.org/10.1016/j.amc.2016.01.049 doi: 10.1016/j.amc.2016.01.049
[7]	R. M. Hafez, Numerical solution of linear and nonlinear hyperbolic telegraph type equations with variable coefficients using shifted Jacobi collocation method, Comput. Appl. Math., 37 (2018), 5253–5273. https://doi.org/10.1007/s40314-018-0635-1 doi: 10.1007/s40314-018-0635-1
[8]	H. B. Jebreen, Y. C. Cano, I. Dassios, An efficient algorithm based on the multi-wavelet Galerkin method for telegraph equation, AIMS Math., 6 (2020), 1296–1308. https://doi.org/10.3934/math.2021080 doi: 10.3934/math.2021080
[9]	W. M. Abd-Elhameed, E. H. Doha, Y. H. Youssri, M. A. Bassuony, New Tchebyshev-Galerkin operational matrix method for solving linear and nonlinear hyperbolic telegraph type equations, Numer. Meth. Part. D. E., 32 (2016), 1553–1571. https://doi.org/10.1002/num.22074 doi: 10.1002/num.22074
[10]	A. Saadatmandi, M. Dehghan, Numerical solution of hyperbolic telegraph equation using the Chebyshev tau method, Numer. Meth. Part. D. E., 26 (2010), 239–252. https://doi.org/10.1002/num.20442 doi: 10.1002/num.20442
[11]	H. B. Jebreen, I. Dassios, A biorthogonal hermite cubic spline Galerkin method for solving fractional riccati equation, Mathematics, 10 (2022), 1461. https://doi.org/10.3390/math10091461 doi: 10.3390/math10091461
[12]	H. B. Jebreen, C. Cattani, Interpolating scaling functions tau method for solving space-time fractional partial differential equations, Symmetry, 14 (2022), 2463. https://doi.org/10.3390/sym14112463 doi: 10.3390/sym14112463
[13]	M. Asadzadeh, B. N. Saray, On a multiwavelet spectral element method for integral equation of a generalized Cauchy problem, BIT, 62 (2022), 383–1416. https://doi.org/10.1007/s10543-022-00915-1 doi: 10.1007/s10543-022-00915-1
[14]	M. H. Heydari, M. Razzaghi, Highly accurate solutions for space-time fractional Schrödinger equations with non-smooth continuous solution using the hybrid clique functions, Math. Sci., 17 (2023), 31–42. https://doi.org/10.1007/s40096-021-00437-x doi: 10.1007/s40096-021-00437-x
[15]	X. Jiang, J. Wang, W. Wang, H. Zhang, A predictor-corrector compact difference scheme for a nonlinear fractional differential equation, Fractal Fract., 7 (2023), 521. https://doi.org/10.3390/fractalfract7070521 doi: 10.3390/fractalfract7070521
[16]	X. Yang, Q. Zhang, G. Yuan, Z. sheng, On positivity preservation in nonlinear finite volume method for multi-term fractional subdiffusion equation on polygonal meshes, Nonlinear Dynam., 92 (2018), 595–612.
[17]	H. Zhang, X. Yang, D. Xu, An efficient spline collocation method for a nonlinear fourth-order reaction subdiffusion equation, J. Sci. Comput., 85 (2020). https://doi.org/10.1007/s10915-020-01308-8
[18]	A. Iqbal, T. Akram, A numerical study of anomalous electro-diffusion cells in cable sense with a non-singular kernel, Demonstr. Math., 55 (2022), 574–586. https://doi.org/10.1515/dema-2022-0155 doi: 10.1515/dema-2022-0155
[19]	T. Akram, M. Abbas, A. Ali, A. Iqbal, D. Baleanu, A numerical approach of a time fractional reaction-diffusion model with a non-singular kernel, Symmetry, 12 (2020), 1653. https://doi.org/10.3390/sym12101653 doi: 10.3390/sym12101653
[20]	V. R. Hosseini, W. Chen, Z. Avazzadeh, Numerical solution of fractional telegraph equation by using radial basis functions, Eng. Anal. Bound. Elem., 38 (2014), 31–39. https://doi.org/10.1016/j.enganabound.2013.10.009 doi: 10.1016/j.enganabound.2013.10.009
[21]	Y. H. Youssri, W. M. Abd-Elhameed, Numerical spectral Legendre-Galerkin algorithm for solving time fractional telegraph equation, Rom. J. Phys., 63 (2018), 1–16.
[22]	N. Mollahasani, M. M. Mohseni, K. Afrooz, A new treatment based on hybrid functions to the solution of telegraph equations of fractional order, Appl. Math. Model., 40 (2016), 2804–2814. https://doi.org/10.1016/j.apm.2015.08.020 doi: 10.1016/j.apm.2015.08.020
[23]	A. Saadatmandi, M. Mohabbati, Numerical solution of fractional telegraph equation via the tau method, Math. Rep., 17 (2015), 155–166.
[24]	J. Chen, F. Liu, V. Anh, Analytical solution for the time-fractional telegraph equation by the method of separating variables, J. Math. Anal. Appl., 338 (2008), 1364–1377. https://doi.org/10.1016/j.jmaa.2007.06.023 doi: 10.1016/j.jmaa.2007.06.023
[25]	W. Jiang, Y. Lin, Representation of exact solution for the time-fractional telegraph equation in the reproducing kernel space, Commun. Nonlinear Sci., 16 (2011), 3639–3645. https://doi.org/10.1016/j.cnsns.2010.12.019 doi: 10.1016/j.cnsns.2010.12.019
[26]	A. Ali, T. Abdeljawad, A. Iqbal, T. Akram, M. Abbas, On unconditionally stable new modified fractional group iterative scheme for the solution of 2D time-fractional telegraph model, Symmetry, 13 (2021), 2078. https://doi.org/10.3390/sym13112078 doi: 10.3390/sym13112078
[27]	M. Shahriari, B. N. Saray, B. Mohammadalipour, S. Saeidian, Pseudospectral method for solving the fractional one-dimensional Dirac operator using Chebyshev cardinal functions, Phys. Scripta., 98 (2023), 055205. https://doi.org/10.1088/1402-4896/acc7d3 doi: 10.1088/1402-4896/acc7d3
[28]	A. Afarideh, F. D. Saei, M. Lakestani, B. N. Saray, Pseudospectral method for solving fractional Sturm-Liouville problem using Chebyshev cardinal functions, Phys. Scripta, 96 (2021), 125267. https://doi.org/10.1088/1402-4896/ac3c59 doi: 10.1088/1402-4896/ac3c59
[29]	C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral methods fundamentals in single domains, Berlin: Springer-Verlag, 2006.
[30]	A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier, 2006.
[31]	G. Dahlquist, A. Björck, Numerical methods, Englewood Cliffs: Prentice Hall, 1974.
[32]	H. Zhang, Y. Liu, X. Yang, An efficient ADI difference scheme for the nonlocal evolution problem in three-dimensional space, J. Appl. Math. Comput., 69 (2023), 651–674.

This article has been cited by:

1.	Batirkhan Turmetov, Valery Karachik, On solvability of some inverse problems for a nonlocal fourth-order parabolic equation with multiple involution, 2024, 9, 2473-6988, 6832, 10.3934/math.2024333
2.	İrem Bağlan, Ahmet Akdemir, Mustafa Dokuyucu, Inverse coefficient problem for quasilinear pseudo-parabolic equation by fourier method, 2023, 37, 0354-5180, 7217, 10.2298/FIL2321217B

Reader Comments

Your name:*

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