In this paper, the solvability of some inverse problems for a nonlocal analogue of a fourth-order parabolic equation was studied. For this purpose, a nonlocal analogue of the biharmonic operator was introduced. When defining this operator, transformations of the involution type were used. In a parallelepiped, the eigenfunctions and eigenvalues of the Dirichlet type problem for a nonlocal biharmonic operator were studied. The eigenfunctions and eigenvalues for this problem were constructed explicitly and the completeness of the system of eigenfunctions was proved. Two types of inverse problems on finding a solution to the equation and its righthand side were studied. In the two problems, both of the righthand terms depending on the spatial variable and the temporal variable were obtained by using the Fourier variable separation method or reducing it to an integral equation. The theorems for the existence and uniqueness of the solution were proved.
Citation: Batirkhan Turmetov, Valery Karachik. On solvability of some inverse problems for a nonlocal fourth-order parabolic equation with multiple involution[J]. AIMS Mathematics, 2024, 9(3): 6832-6849. doi: 10.3934/math.2024333
[1] | Yashar Mehraliyev, Seriye Allahverdiyeva, Aysel Ramazanova . On one coefficient inverse boundary value problem for a linear pseudoparabolic equation of the fourth order. AIMS Mathematics, 2023, 8(2): 2622-2633. doi: 10.3934/math.2023136 |
[2] | 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 |
[3] | Fengfei Jin, Baoqiang Yan . Existence and global behavior of the solution to a parabolic equation with nonlocal diffusion. AIMS Mathematics, 2021, 6(5): 5292-5315. doi: 10.3934/math.2021313 |
[4] | 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 |
[5] | 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 |
[6] | Andrey Muravnik . Nonclassical dynamical behavior of solutions of partial differential-difference equations. AIMS Mathematics, 2025, 10(1): 1842-1858. doi: 10.3934/math.2025085 |
[7] | Jeong-Kweon Seo, Byeong-Chun Shin . Reduced-order modeling using the frequency-domain method for parabolic partial differential equations. AIMS Mathematics, 2023, 8(7): 15255-15268. doi: 10.3934/math.2023779 |
[8] | Yaxin Zhao, Xiulan Wu . Asymptotic behavior and blow-up of solutions for a nonlocal parabolic equation with a special diffusion process. AIMS Mathematics, 2024, 9(8): 22883-22909. doi: 10.3934/math.20241113 |
[9] | Anh Tuan Nguyen, Le Dinh Long, Devendra Kumar, Van Thinh Nguyen . Regularization of a final value problem for a linear and nonlinear biharmonic equation with observed data in $ L^{q} $ space. AIMS Mathematics, 2022, 7(12): 20660-20683. doi: 10.3934/math.20221133 |
[10] | Jinhua Feng, Shixi Wang, Hai Bi, Yidu Yang . The a posteriori error estimates of the Ciarlet-Raviart mixed finite element method for the biharmonic eigenvalue problem. AIMS Mathematics, 2024, 9(2): 3332-3348. doi: 10.3934/math.2024163 |
In this paper, the solvability of some inverse problems for a nonlocal analogue of a fourth-order parabolic equation was studied. For this purpose, a nonlocal analogue of the biharmonic operator was introduced. When defining this operator, transformations of the involution type were used. In a parallelepiped, the eigenfunctions and eigenvalues of the Dirichlet type problem for a nonlocal biharmonic operator were studied. The eigenfunctions and eigenvalues for this problem were constructed explicitly and the completeness of the system of eigenfunctions was proved. Two types of inverse problems on finding a solution to the equation and its righthand side were studied. In the two problems, both of the righthand terms depending on the spatial variable and the temporal variable were obtained by using the Fourier variable separation method or reducing it to an integral equation. The theorems for the existence and uniqueness of the solution were proved.
Let Q=Π×(0,T), where Π={x∈Rn:0<xj<pj,j=1,…,n} is a parallelepiped, pj>0, and T>0. Consider the mappings Sj:Rn→Rn,1≤j≤n of the type Sjx=(x1,…,xj−1,pj−xj,xj+1,…,xn). Obviously, the mappings Sj are involutions, i.e., S2j=I, where I is the identity mapping. Let us consider all possible products of mappings Sj, i.e., Sij=SiSi, or Sijk=SiSiSk,.... The total number of such mappings, taking into account the identity mapping S0x=x, is equal to 2n. To number such mappings, we will use the binary number system, namely, if 0≤i<2n in the binary number system, the representation i≡(in…i1)2=i1+2i2+…+2n−1in, where ik=0,1, is valid. Therefore, introducing the vector i=(i1,…,in), we can consider mappings of the type Si≡Si11…Sinn corresponding to the index i. Using these mappings, we introduce the operator
Lxv(x)=2n−1∑i=0aiΔ2v(Six), |
where a0,a1,a2,a3,…,a2n−1 is a set of real numbers, and Δ2 is a biharmonic operator.
Let us consider the following problems in the domain Q.
Problem 1. Find a pair of such functions {u(t,x),f(x)} that are smooth u(t,x)∈C(ˉQ), f(x)∈C(ˉΠ), ut(t,x), Lxu(t,x)∈C(Q), satisfy the equation
ut(t,x)+Lxu(t,x)=f(x)g(t),(t,x)∈Q, | (1.1) |
the boundary conditions
u(t,x)|t∈[0,T],x∈∂Π=0,∂2u(t,x)∂x2i|t∈[0,T],x∈∂Π=0,i=1,2,...,n, | (1.2) |
and the initial and final conditions
u(0,x)=φ(x),x∈ˉΠ, | (1.3) |
u(T,x)=ψ(x),x∈ˉΠ, | (1.4) |
where g(t), φ(x), ψ(x) are the given functions.
Problem 2. Find a pair of such functions {u(t,x),g(t)} that are smooth u(t,x)∈C(ˉQ), g(t)∈C[0,T], ut(t,x),Lxu(t,x)∈C(Q) and satisfy conditions (1.1)–(1.3) and an additional condition
u(t,x0)=h(t),0≤t≤T, | (1.5) |
where x0∈Π is a given fixed point, h(t),φ(x), and f(x) are the given functions, and h(0)=φ(x0).
Differential equations with involution are an important part of the general theory of functional differential equations. As it is known, the first works on differential equations with involution were written by Babbage [1] and continued by Carleman [2]. Various issues in the theory of differential equations with involution were studied in a series of papers written by Przeworska-Rolewicz [3,4,5] The case n=1 is considered in several papers [6,7,8,9,10,11], where the solvability of inverse problems on finding the righthand side of differential equations with involution is studied. Thus, in the work of Nasser Al-Salti et al. [7], in the rectangular domain Ω={(t,x):0<t<T,−π<x<π} for the equation
ut(t,x)−uxx(t,x)−εuxx(t,π−x)=f(x),(t,x)∈Ω, |
inverse problems for determining a pair of functions {u(t,x),f(x)} have been studied. In this work, the authors state that redefined initial and final conditions
u(0,x)=φ(x),u(T,x)=φ(x),−π<x<π, |
as well as one of the boundary conditions, Dirichlet, Neumann, periodic, or antiperiodic conditions, can be considered as additional conditions. The problem studied in the paper is solved using the method of separation of variables. In this case, a spectral problem arises with respect to the spatial variable for the equation with involution. For example, in the case of Dirichlet boundary conditions, we have
−X″(x)−εX″(π−x)=λX(x),−π<x<π,X(−π)=X(π)=0. |
It is proven that the eigenfunctions of this problem are the functions
X1k(x)=cos(k+12)x,k≥0;X2k(x)=sinkx,k≥1, |
and the corresponding eigenvalues are the numbers
λ1k=(1−ε)(k+12)2,k≥0;λ2k=(1+ε)k,k≥1. |
Similar studies in the case of two spatial variables for classical equations were carried out in [12,13,14,15,16,17,18], as well as for equations with involution in [19,20]. The authors of this paper studied inverse problems for equations of parabolic type with an operator of the fourth and higher orders in the spatial variable (see, for example, [21,22,23,24,25,26,27,28]).
The methods used to solve the inverse problems in order to determine the right side of the equation depend on whether the function f(x) or the function g(t) is unknown. In the case when the function f(x) is unknown, the problem under consideration is usually solved using the Fourier method, whereas in the case when the function g(t) is unknown, the problem is reduced to the Volterra integral equation (see, for example, [18]).
As we have already noted, studies of direct and inverse problems for equations with involutive transformed arguments were mainly carried out for equations with one and two spatial variables. For equations with many variables, such problems are insufficiently studied. In this direction, we can only note the work of Kozhanov and Bzheumikhova [25].
When studying Problems 1 and 2, a spectral problem arises for the nonlocal operator Lx. In the one-dimensional case in [7], this problem is solved by finding a general solution to the studied equation and using the boundary conditions to determine the eigenfunctions and eigenvalues. In the two-dimensional case in [19], the corresponding spectral problem was solved by reducing it to four auxiliary problems.
In general, in the n− dimensional case, n≥3, such spectral problems in a parallelepiped have not been considered. When solving this problem, to denote the summation index we used the notation in the binary number system. The transition to this system allowed us to construct eigenfunctions and eigenvalues of this problem explicitly and to prove the completeness of the systems of eigenfunctions in space L2(Π).
It should be also noted that the use of differential equations with involution when modeling a specific physical process of thermal diffusion in the case n=1 is given in [7,9], and in [29,30] in the case n=2, where equations with involutive transformations, which have applications in modeling of optical systems, are considered. In addition, the influence of nonlocality is graphically illustrated in [7] in the one-dimensional case.
In this section, we study the eigenfunctions and eigenvalues of a Dirichlet type problem for a nonlocal biharmonic equation.
Consider the following boundary value problem.
Problem S. Find a function v(x)≠0 from the class v(x)∈C2(ˉΠ)∩C4(Π) that satisfies the equation
Lxv(x)=λv(x),x∈Π, | (2.1) |
and boundary conditions
v(x)|∂Π=∂2v(x)∂x2i|∂Π=0,i=1,2,...,n, | (2.2) |
where λ∈R.
If a0=1, aj=0, j=1,…,2n−1, then the problem coincides with the spectral problem with the Dirichlet condition for the classical biharmonic operator.
Note that in the case of a sphere, a similar spectral problem for the nonlocal Laplace operator was studied in [31]. Various boundary value problems for nonlocal harmonic and biharmonic equations are studied in [32,33,34].
Let us consider a set of functions
vk(x)=vk1k2...kn(x)=C(n,p)n∏j=1sinkjπxjpj, | (2.3) |
where k=(k1k2...kn)∈Nn and C(n,p)=2n/2n∏j=11√pj,p=(p1,…,pn).
The following statement is proved in [35].
Lemma 2.1. The system of functions {vk(x):k∈Nn} is orthonormal and complete in space L2((0,p1)×(0,p2)×...×(0,pn)).
Let us introduce the following numbers εk=2n−1∑i=0(−1)i⋅(k+e)ai, where k=(k1,k2,...,kn), i=(i1,…,in), i⋅k=i1k1+...+inkn, and e=(1,…,1). Note that components of the vector k=(k1,k2,...,kn) can be calculated using mode 2 as in this case the εk value will not change. Taking this into account we can state that the equality k+e=k∗ is valid, where the vector k∗ is conjugate to k. For example, k=(0,1,0,1)⇒k∗=(1,0,1,0). Under these assumptions,
εk=εkmod2=2n−1∑i=0(−1)i⋅k∗ai. |
Theorem 2.1. Let the conditions εk≠0 be satisfied for all k∈Nn, then the system of functions {vk(x):k∈Nn} is a system of eigenfunctions of Problem S. The corresponding eigenvalues are determined by the equalities
λk=εkμ2k,μk=π2n∑j=1k2jp2j,k=(k1,k2,…,kn)∈Nn. | (2.4) |
Proof. It is obvious that the functions vk(x), by their structure, satisfy conditions (2.2) vk(x)|∂Π=0, and as
∂2vk∂x2i=C(n,p)n∏j=1,j≠isinkjπxjpj⋅∂2∂x2i(sinkiπxipi)=−(kiπpi)2C(n,p)n∏j=1,j≠isinkjπxjpj⋅(sinkiπxipi)=−(kiπpi)2vk(x), |
they also satisfy conditions (2.2). Let us check whether Eq (2.1) is fulfilled. Applying the Laplace operator to this function, we obtain
Δvk(x)=−(n∑i=1(kiπpi)2)vk(x), |
which means that
Δ2vk(x)=μ2kvk(x),μk=π2(n∑i=1k2ip2i). |
Hence, for any m∈{1,2,...,n}, we get
Δ2vk(Smx)=μ2kvk(x1,...,xm−1,pm−xm,xm+1,...,xn)=μ2kC(p)sinkmπ(pm−xm)pm⋅n∏j=1,j≠msinkjπxjpj=μ2kC(p)(−1)km+1sinkmπxmpm⋅n∏j=1,j≠msinkjπxjpj=μ2kC(p)(−1)km+1n∏j=1sinkjπxjpj=(−1)km+1μ2kvk(x). |
Let 0≤i≤2n−1, which corresponds to vector i=(i1,…,in). If im=1, then
Δ2vk(Simmx)=μ2kC(p)sinkmπ(pm−xm)pm⋅n∏j=1,j≠msinkjπxjpj=(−1)km+1μ2kvk(x). |
Obviously, this equality is also valid for the case im=0. Thus, we get
Δ2vk(Simmx)=(−1)im(km+1)μ2kvk(x). | (2.5) |
In the general case, the following equality holds:
Δ2vk(Six)=Δ2vk(Si11…Sinnx)=(−1)i1(k1+1)+i2(k2+1)+...+in(kn+1)μ2kvk(x) |
=(−1)i⋅(k+1)μ2kvk(x)=(−1)i⋅k∗μ2kvk(x), |
where k=(k1,k2,...,kn), i=(i1,…,in), i⋅k=i1k1+...+inkn and e=(1,…,1).
Now, if we apply the operator Lx to the function vk, then from the previous equalities it follows that
Lxvk(x)=2n−1∑i=0aiΔ2vk(Six)=2n−1∑i=0ai(−1)i⋅k∗μ2kvk(x)=μ2kvk(x)(2n−1∑i=0(−1)i⋅k∗ai)=εkmod2μ2kvk(x). |
Thus, in addition to conditions (2.2), the function vk(x) also satisfies the equality Lxvk(x)=λkvk(x), i.e., Eq (2.1). The theorem is proved.
Remark 2.1. If the condition εkmod2>0 is satisfied for any k∈Nn, then all eigenvalues of Problem S are positive.
Example 2.1. Let n=2, then Eq (2.1) takes the form
a0Δ2v(x1,x2)+a1Δ2v(p1−x1,x2)+a2Δ2v(x1,p2−x2)+a2Δ2v(p1−x1,p2−x2)=λv(x1,x2), |
and the boundary conditions are written as
v(0,x2)=v(p1,x2)=0,0≤x2≤p2;v(x1,0)=v(x1,p2)=0,0≤x1≤p1,vx1x1(0,x2)=vx1x1(p1,x2)=0,0≤x2≤p2;vx2x2(x1,0)=vx2x2(x1,p2)=0,0≤x1≤p1. |
Eigenfunctions are specified in accordance with (2.3) as
v(m,k)(x)=2√p1p2sinmπx1p1sinkπx2p2,m,k=1,2,..., |
and the corresponding eigenvalues are
λ(m,k)=ε(m,k)π4(m2p21+k2p22)2,m,k=1,2,…, |
where ε(m,k) is
ε(m,k)=ε(m,k)mod2=a0+(−1)m+1a1+(−1)k+1a2+(−1)k+ma3=a0+(−1)m∗a1+(−1)k∗a2+(−1)k∗+m∗a3. |
More precisely, ε(m,k) can be written as
ε(2m−1,2k−1)=a0+a1+a2+a3;ε(2m−1,2k)=a0+a1−a2−a3;ε(2m,2k−1)=a0−a1+a2−a3;ε(2m,2k)=a0−a1−a2+a3, |
where m,k=1,2,....
Let
h(x)=∞∑k1=1...∞∑kn=1hk1...knvk1...kn(x)=∑k∈Nnhkvk(x) | (3.1) |
be the Fourier series expansion of the function h(x) by the system {vk(x):k∈Nn}, where
hk=(h,vk)≡∫p10...∫pn0h(x1,x2,...,xn)vk1...kn(x1,x2,...,xn)dx1...dxn. | (3.2) |
Further, we will use the symbol C to denote an arbitrary constant whose value does not affect our conclusions.
Lemma 3.1. Let the function h(x) be continuous in a closed domain ˉΠ and have continuous partial derivatives ∂jh(x)∂x1...∂xj, 1≤j≤n in ˉΠ. If the conditions
h(0,x2,...,xn)=h(p1,x2,...,xn)=0,0≤xj≤pj,j=2,...,n,h(x1,0,...,xn)=h(x1,p2,...,xn)=0,0≤xj≤pj,j=1,...,n,j≠2,…h(x1,...,xn−1,0)=h(x1,...,xn−1,pn)=0,0≤xj≤pj,j=1,...,n−1 | (3.3) |
are satisfied, then the number series ∑k∈Nn|hk| converges.
Proof. If h(x)∈C(ˉΠ) and the function ∂h(x)∂x1 is continuous, then integrating the integral in (3.2) by parts over the variable x1 and taking into account equality (3.3), we obtain
hk=hk1...kn=1k1h1,0,..,0k1...kn, |
where
h1,0,..,0k1...kn=C∫p10...∫pn0∂h(x1,x2,...,xn)∂x1v1,0,..,0k1...kn(x1,x2,...,xn)dx1...dxn,v1,0,..,0k1...kn(x1,x2,...,xn)=cosk1πx1p1⋅n∏j=2sinkjπxjpj. |
Applying this process to all j∈{2,3,...,n}, we get
hk=hk1...kn=1k1k2...knh1,1,..,1k1...kn, |
where
h1,1,..,1k1...kn=C∫p10...∫pn0∂nh(x1,x2,...,xn)∂x1...∂xnv1,1,..,1k1...kn(x1,x2,...,xn)dx1...dxn,v1,1,..,1k1...kn(x1,x2,...,xn)=n∏j=1coskjπxjpj. |
Using the Cauchy-Bunyakovsky inequality, we obtain
∑k∈Nn|hk|≤∞∑k1=1...∞∑kn=1|1k1...knh1,1,..,1k1...kn|≤√∞∑k1=1...∞∑kn=11k21...k2n√∞∑k1=1...∞∑kn=1|h1,1,..,1k1...kn|2. |
As the system {v1,1,..,1k1...kn(x1,x2,...,xn)} is orthogonal in space L2(Π) and ∂nh(x1,x2,...,xn)∂x1...∂xn∈L2(Π), then due to Bessel's inequality, the series ∞∑k1=1...∞∑kn=1|h1,1,..,1k1...kn|2 converges. Moreover, the series
∞∑k1=1...∞∑kn=11k21...k2n=∞∑k1=11k21...∞∑kn=11k2n<∞ |
also converges. This implies the assertion of the lemma.
The following assertion is proved in a similar way.
Lemma 3.2. Let a function h(x) belong to a class C4(ˉΠ) and have continuous partial derivatives of the form ∂n+4h(x)∂x1...∂x5j...∂xn for j∈{1,2,...,n} in ˉΠ. If the functions h(x),∂2h(x)∂x2j, and ∂4h(x)∂x4j satisfy conditions (3.3), the number series ∞∑k1=1...∞∑kj=1...∞∑kn=1k4j|hk1...kj...kn| converges.
If the functions h(x),∂2h(x)∂x2j, and ∂4h(x)∂x4j satisfy conditions (3.3), the number series
∞∑k1=1...∞∑kj=1...∞∑kn=1k4j|hk1...kj...kn| |
converges.
Proof. Let the functions h(x) and ∂2h(x)∂x2j,∂4h(x)∂x4j satisfy conditions (3.3), then integrating the integral four times over the variable xj in the equality
hk1...kn=Ck1...kn∫p10...∫pn0∂nh(x1,x2,...,xn)∂x1...∂xnv1,1,..,1k1...kn(x1,x2,...,xn)dx1...dxn, |
we get
hk1...kn=1k1k2...k5j...knh1,1,..,5,...,1k1k2...kj...kn, |
where
h1,1,..,5,...,1k1k2...kj...kn=C∫p10...∫pn0∂n+4h(x1,x2,...,xn)∂x1...∂x5j...∂xnv1,1,..,5,...,1k1k2...kj...kn(x1,x2,...,xn)dx1...dxn,v1,1,..,5,...,1k1k2...kj...kn(x1,x2,...,xn)=n∏j=1coskjπxjpj. |
Using the Cauchy-Bunyakovsky inequality for ∞∑k1=1...∞∑kn=1|hk1...kn|, we obtain
∞∑k1=1...∞∑kj=1...∞∑kn=1k4j|hk1...kn|=∞∑k1=1...∞∑kn=1k4jk1k2...k5j...kn|h1,1,..,3,...,1k1k2...kj...kn|≤√∞∑k1=1...∞∑kn=11k21...k2j...k2n√∞∑k1=1...∞∑kn=1|h1,1,..,3,...,1k1k2...kj...kn|2. |
As the system {v1,1,..,1k1...kn(x1,x2,...,xn)} is orthogonal in space L2(Π) and ∂nh(x1,x2,...,xn)∂x1...∂xn∈L2(Π), then due to Bessel's inequality, the series ∞∑k1=1...∞∑kn=1|h1,1,..,1k1...kn|2 converges. This implies the statement of the lemma.
Corollary 3.1. Let the conditions of Lemma 3.2 be satisfied, then the number series
∑k∈Nnλk|hk|=∞∑k1=1...∞∑kj=1...∞∑kn=1λk1...kj...kn|hk1...kj...kn| |
converges.
The proof of this assertion follows from the representation of eigenvalues in the form λk1k2...kn=εk1k2...knπ4(n∑j=1k2jp2j)2 and the assertion of Lemma 3.2.
By definition, the solution to Problem 1 is expressed as functions u(t,x),f(x). If such functions exist, they must be the elements of space L2(Π) and, therefore, they can be represented as series
u(t,x)=∞∑k∈Nnuk(t)vk(x)=∞∑k1=1∞∑kn=1uk1...kn(t)vk1...kn(x), | (4.1) |
f(x)=∑k∈Nnfkvk(x)=∞∑k1=1…∞∑kn=1fkvk(x), | (4.2) |
where uk(t)=uk1...kn(t) and fk=fk1...kn are the coefficients to be determined.
Substituting (4.1) and (4.2) into Eq (1.1) and equating the coefficients for vk(x), we obtain the following problem for the coefficients uk(t):
∂uk(t)∂t=−λkuk(t)+fkg(t), | (4.3) |
uk(0)=φk,uk(T)=ψk. | (4.4) |
The general solution to Eq (4.3) is the function
uk(t)=Ck⋅e−λkt+fk∫t0e−λk(t−τ)g(τ)dτ,k=(k1,...,kn)∈Nn, | (4.5) |
where Ck are arbitrary constants. Let us introduce the notation
gk(t)=∫t0e−λk(t−τ)g(τ)dτ, |
then uk(t) from (4.5) is represented in the form
uk(t)=Ck⋅e−λkt+fkgk(t),k=(k1,...,kn)∈Nn. |
Now, we will consider two cases separately: g(t)=1 and g(t)≠1.
(Ⅰ) Suppose g(t)=1, then
gk(t)=∫t0e−λk(t−τ)dτ=1λk(1−e−λkt) |
and uk(t) are represented as
uk(t)=Ck⋅e−λkt+fkλk(1−e−λkt). |
If we use condition (4.4), we obtain
Ck=uk(0)=φk,ψk=uk(T)=φke−λkT+fkλk(1−e−λkT), |
then we find
fk=λkψk−φke−λkT1−e−λkT. | (4.6) |
The solution to problem (4.3), (4.4) is represented as
uk(t)=e−λkt−e−λkT1−e−λkTφk+1−e−λkt1−e−λkTψk. | (4.7) |
Hence, the solution to Problem 1 can be written as
u(t,x)=∞∑k1=1...∞∑kn=1[e−λkt−e−λkT1−e−λkTφk+1−e−λkt1−e−λkTψk]vk(x), | (4.8) |
f(x)=∞∑k1=1...∞∑kn=1λkψk−φke−λkT1−e−λkTvk(x). | (4.9) |
By construction and also due to the properties of functions vk(x), the function u(t,x) from (4.8) formally satisfies conditions (1.1)–(1.4). Let us examine the smoothness of the functions u(t,x) and f(x).
Let the functions φ(x) and ψ(x) satisfy the conditions of Corollary 3.1, then the number series
∞∑k1=1...∞∑kn=1λk(|φk|+|ψk|) |
converges. As λk=εkμ2k and εk>0 does not tend to 0, the functions 11−e−λkT and e−λkT1−e−λkT are bounded. For the coefficients fk from equality (4.6), we obtain
|fk|≤λk11−e−λkT|ψk|+λke−λkT1−e−λkT|φk|≤C(λk|φk|+λk|ψk|). |
Based on this estimate and the uniform boundedness of the moduli of the eigenfunctions |vk1...kn(x)| from (2.3), we obtain absolute and uniform convergence of the functional series (4.9) in the closed domain ˉΠ. Hence, the sum of this series, i.e., the function f(x), is continuous in the domain ˉΠ. For all 0≤t≤T, there are estimates
|e−λkt−e−λkT1−e−λkT|=e−λkT|eλk(T−t)−11−e−λkT|≤C,|1−e−λkt1−e−λkT|≤C. |
The series (4.8) satisfies the estimate
|∞∑k1=1...∞∑kn=1[e−λkt−e−λkT1−e−λkTφk+1−e−λkt1−e−λkTψk]vk(x)|≤C∞∑k1=1...∞∑kn=1[|φk|+|ψk|]<∞. |
This means that this series converges absolutely and uniformly in a closed domain ˉQ and, therefore the function u(t,x), the sum of this series, belongs to the class C(ˉQ).
Differentiating series (4.8) with respect to the variable t, we obtain
ut(t,x)=∞∑k1=1...∞∑kn=1[−λke−λkt1−e−λkTφk+λke−λkt1−e−λkTψk]vk(x). | (4.10) |
If f(λ)=λe−λt, then maxλ≥0f(λ)=f(1/t)=1te−1, which means that λe−λt≤1δe for t≥δ. Therefore, the series (4.10) satisfies the estimate
|ut(t,x)|≤C∞∑k1=1...∞∑kn=1[|φk|+|ψk|]. |
If the condition of Lemma 3.1 is satisfied, the last series converges, then for an arbitrary δ>0, series (4.10) converges absolutely and uniformly in the closed domain ˉQδ=[δ≤t≤T]×ˉΠ⊂Q. Hence, ut(t,x)∈C(ˉQδ). Thus, by virtue of arbitrariness of δ>0, it follows that ut(t,x)∈C(Q). The inclusion Lxu(t,x)∈C(Q) is proved in a similar way. Thus, functions (4.8) and (4.9) for g(t)=1 satisfy all the conditions of Problem 1.
If in Problem 1 homogeneous conditions (4.2) are specified, then for the coefficients fk1...kn in equality (4.6), we obtain that fk1...kn≡(f,vk1k2...k1)=0. Hence, the function f(x) is orthogonal to all elements of the system {vk(x)}k∈Nn. Due to the completeness of this system and continuity of the function f(x), we obtain f(x)≡0,x∈ˉΠ.
Similarly, for the coefficients uk(t) from equality (4.5) we obtain uk(t)≡(u,vk)=0. Hence, the equality u(t,x)=0, x∈ˉΠ is valid for almost all t∈[0,T]. Due to the continuity of the function u(t,x) in ˉQ, we obtain that u(t,x)≡0, (t,x)∈ˉQ. This implies the uniqueness of the solution to Problem 1.
Thus, we proved the following assertion.
Theorem 4.1. Let the functions φ(x) and ψ(x) in Problem 1 satisfy the conditions of Corollary 3.1, g(t)=1 and let the coefficients ai,i=0,1,...,2n−1 be such that the conditions εk>0 are satisfied. The solution to the problem exists, is unique, and is represented in the form of series (4.8) and (4.9).
(Ⅱ) Let us study Problem 1 for the case g(t)≠1. If we search for a solution to the problem in the form of (4.1) and (4.2), then for unknown coefficients uk(t), we obtain problem (4.3) and (4.4). Moreover, the general solution to Eq (4.3) is determined by equality (4.5). Substituting this function into condition (4.4), we have
φk=uk(0)=Ck,ψk=uk(T)=φk⋅e−λkT+fkgk(T). |
Thus, if for all k∈Nn, the condition gk(T)≠0, then
fk=1gk(T)[ψk−φk⋅e−λkT], | (4.11) |
and
uk1...kn(t)=[1−gk1...kn(t)gk1...kn(T)e−λk1...kn(T−t)]e−λk1...kntφk1...kn+gk1...kn(t)gk1...kn(T)ψk1...kn. | (4.12) |
If in Problem 1, as in the case of g(t)=1, homogeneous conditions (4.2) are given, we obtain f(x)≡0, x∈ˉΠ and u(t,x)≡0,(t,x)∈ˉQ. Therefore, if the conditions gk(T)≠0, k∈Nn are satisfied, the solution to Problem 1 is unique. If for some m=(m1,m2,...,mn) the equality gm(T)=0 holds, then the homogeneous Problem 1 has a nonzero solution. Let us show that if this condition is satisfied, a pair of functions
u(t,x)=fmgm(t)vm(x),f(x)=fmvm(x) |
will be a solution to homogeneous Problem 1, where fm is an arbitrary constant. Indeed, applying the operators ∂∂t and Lx to the function u(t,x), we get
ut(t,x)=fm∂gm(t)∂tvm(x)=fmg(t)vm(x)−λmfmgm(t)vm(x),Lxu(t,x)=fmgm(t)Lvm(x)=−λmfmgm(t)vm(x). |
Hence,
ut(t,x)−Lu(t,x)=fmg(t)vm(x)=f(x)g(t). |
It is obvious that this function satisfies homogeneous conditions (1.2) and (1.3). Thus, we proved the following assertion.
Theorem 4.2. If a solution to Problem 1 exists, then it is unique if, and only if, the conditions gk(T)≠0 are satisfied for all k∈Nn.
Regarding the existence of a solution to Problem 1 in the case g(t)≠1, the following assertion is valid.
Theorem 4.3. Let in Problem 1 the functions φ(x) and ψ(x) satisfy the conditions of Corollary 3.1, and the coefficients ai, i=0,1,...,2n−1 are such that the conditions εk>0 are satisfied. A solution to the problem exists and is represented in the form of series (4.8) and (4.9), where the coefficients fk and uk are respectively determined by equalities (4.11) and (4.12).
Proof. If g(t)∈C[0,T],|g(t)|≥g0≡C, then according to the mean value theorem there is a point ξ∈[0,T] such that
gk(T)=∫T0e−λk(T−τ)g(τ)dτ=g(ξ)∫T0e−λk(T−τ)dτ=g(ξ)1−e−λkTλk. |
This gives us the following lower estimate:
|gk1...kn(T)|=|g(ξ)|1−e−λk1...knTλk1...kn≥Cλk1...kn. | (4.13) |
Using inequality (4.13) for fk from equality (4.11), we get
|fk|≤1|gk(T)||ψk−φk⋅e−λkT|≤Cλk(|φk|+|ψk|). |
Similarly from equality (4.12), we obtain for uk(t),
|uk(t)|≤|e−λkt−gk(t)gk(T)e−λkT||φk|+|gk(t)gk(T)||ψk|. |
As the function g(t) is continuous on the interval [0,T], it follows that
|gk(t)|≤∫t0e−λk(t−τ)|g(τ)|dτ≤max0≤τ≤T|g(τ)|∫T0e−λk(t−τ)dτ. |
Hence, |gk(t)gk(T)|≤C. Therefore, the estimate
|uk(t)|≤C(|φk|+|ψk|) |
is valid for uk(t). The estimate
|∂uk(t)∂t|≤Cλk(|φk|+|ψk|) |
is proved in a similar way. From these estimates and the convergence of the series
∞∑k1=1...∞∑kn=1(|φk|+|ψk|),∞∑k1=1...∞∑kn=1λk(|φk|+|ψk|) |
absolute and uniform convergence of series (4.8) and (4.9) follows. For the sums of these series we get f(x)∈C(ˉΠ) and u(t,x)∈C(ˉQ).
The absolute and uniform convergence of the series
∂u(t,x)∂t=∞∑k1=1...∞∑kn=1∂uk(t)∂tvk(x),Lxu(t,x)=∞∑k1=1...∞∑kn=1uk(t)Lvk(x) |
in an arbitrary closed domain ˉQδ⊂Q and δ>0 is proved similar to the case g(t)=1. Therefore, the inclusions ∂u(t,x)∂t, Lu(t,x)∈C(Q) are valid. The theorem is proved.
The main assertion regarding Problem 2 is the following theorem.
Theorem 5.1. Let εk>0, the function φ(x), satisfy the conditions of Corollary 3.1, then if f(x0)≠0, h(t)∈C1[0,T], h(0)=φ(x0), the solution to Problem 2 exists and is unique.
Proof. If we assume that the function g(t) is known, then the solution to Problem 2 can be represented as
u(t,x)=∞∑k1=1...∞∑kn=1[φk1...kn⋅e−λkt+fk∫t0e−λk(t−τ)g(τ)dτ]vk(x), | (5.1) |
where
fk=∫p10...∫pn0f(x)vk(x)dx1...dxn. |
Let us suppose that x=x0 in (5.1), then
h(t)=∞∑k1=1...∞∑kn=1[φk⋅e−λkt+fk∫t0e−λk(t−τ)g(τ)dτ]vk(x0)=φ0(t)+∞∑k1=1...∞∑kn=1[fk∫t0e−λk(t−τ)g(τ)dτ]vk(x0)=∞∑k1=1...∞∑kn=1φk1...kn⋅e−λktvk(x0)+∞∑k1=1...∞∑kn=1[fk∫t0e−λk(t−τ)g(τ)dτ]vk(x0), |
where
φ0(t)=∞∑k1=1...∞∑kn=1φk⋅e−λktvk(x0). | (5.2) |
Denote r(t)=h(t)−φ0(t) and
K(t,τ)=∞∑k1=1...∞∑kn=1fke−λk(t−τ)vk(x0). | (5.3) |
For the function g(t), we obtain the following Volterra integral equation of the first kind
∫t0K(t,τ)g(τ)dτ=r(t). | (5.4) |
Lemma 5.1. If the function φ(x) satisfies the conditions of Corollary 3.1, then the function φ0(t) from (5.2) is continuous and has a continuous derivative on the interval [0,T].
Proof. Let us differentiate the series (5.2)
φ′0(t)=−∞∑k1=1...∞∑kn=1λkφk⋅e−λktvk(x0). |
As for the points of the domain ˉQ, we have the estimate |e−λktvk(x0)|≤C, then
|φ′0(t)|≤∞∑k1=1∞∑kn=1λk|φk|. |
If the condition of Corollary 3.1 is satisfied, the last number series converges, then series (5.2) converges absolutely and uniformly. Therefore, the sum of this series represents a continuous function, i.e., φ0(t)∈C1[0,T]. The lemma is proved.
Let us study the properties of the kernel K(t,τ). Differentiating series (5.3) with respect to t, we get
Kt(t,τ)=−∞∑k1=1...∞∑kn=1λkfke−λk(t−τ)vk(x0). | (5.5) |
For series (5.3) and (5.5), we obtain the estimates
|K(t,τ)|≤∞∑k1=1...∞∑kn=1|fk|,|Kt(t,τ)|≤∞∑k1=1...∞∑kn=1λk|fk|. |
If the function f(x) satisfies the conditions of Corollary 3.1, then the number series
∞∑k1=1...∞∑kn=1|fk|,∞∑k1=1...∞∑kn=1λk|fk| |
converges. The series (5.3) and (5.5) converge absolutely and uniformly in the closed domain [0,T]×[0,T]. Therefore, the functions K(t,τ) and Kt(t,τ) are continuous in this domain. Further, differentiating equality (5.4), we obtain
K(t,t)g(t)+∫t0Kt(t,τ)g(τ)dτ=r′(t). | (5.6) |
As
K(t,t)=∞∑k1=1...∞∑kn=1fkvk(x0)=f(x0)≠0, |
the equality (5.6) is an integral Volterra equation of the second kind with a continuous kernel and a continuous righthand side. According to the general theory, such an equation has a unique solution g(t) from the class C[0,T].
If we substitute this function into equality (5.1), then the pair of functions {u(t,x),g(t)} satisfies all the conditions of Problem 2. The smoothness of the function and the uniqueness of the solution are proved as in the case of Problem 1. The theorem is proved.
In this paper, the solvability of some inverse problems for a nonlocal analogue of the fourth-order parabolic equation is studied. The nonlocal operator is introduced using involutive mappings. Unlike previous works of the authors, in this paper the problems are studied in the n-dimensional case. The considered problems are solved by applying the Fourier method and reducing them to the Volterra integral equation. In this case, a spectral problem arises for the nonlocal analogue of the biharmonic operator. The eigenfunctions and eigenvalues for this problem are found explicitly and the completeness of the system of eigenfunctions is proved. Solutions to the main problems are constructed in the form of series using a system of eigenfunctions. Further, it is planned to continue the study of inverse problems for high order differential equations with involutions.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research has been funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan (grant No. AP19677926).
The authors declare no conflicts of interest in this paper.
[1] | C. Babbage, An essay towards the calculus of functions, Philos. Trans. R. Soc. Lond., 105 (1815), 389–423. |
[2] | T. Carleman, La théorie des équations intégrales singulières et ses applications, Ann. Inst. H. Poincaré, 1 (1930), 401–430. |
[3] |
D. Przeworska-Rolewicz, On equations with several involutions of different orders and its applications to partial differential-difference equations, Stud. Math., 32 (1969), 101–113. https://doi.org/10.4064/sm-32-2-101-113 doi: 10.4064/sm-32-2-101-113
![]() |
[4] |
D. Przeworska-Rolewicz, On equations with rotations, Stud. Math., 35 (1970), 51–68. https://doi.org/10.4064/sm-35-1-51-68 doi: 10.4064/sm-35-1-51-68
![]() |
[5] | D. Przeworska-Rolewicz, Equations with transformed argument: an algebraic approach, Polish Scientific Publishers and Elsevier Scientific Publishing Company, 1973. |
[6] |
A. Ahmad, M. Ali, S. A. Malik, Inverse problems for diffusion equation with fractional Dzherbashian-Nersesian operator, Fract. Calc. Appl. Anal., 24 (2021), 1899–1918. https://doi.org/10.1515/fca-2021-0082 doi: 10.1515/fca-2021-0082
![]() |
[7] | N. Al-Salti, M. Kirane, B. T. Torebek, On a class of inverse problems for a heat equation with involution perturbation, Hacet. J. Math. Stat., 48 (2019), 669–681. |
[8] |
M. Ali, S. Aziz, S. A. Malik, Inverse source problems for a space-time fractional differential equation, Inverse Probl. Sci. Eng., 28 (2020), 47–68. https://doi.org/10.1080/17415977.2019.1597079 doi: 10.1080/17415977.2019.1597079
![]() |
[9] |
M. Kirane, M. A. Sadybekov, A. A. Sarsenbi, On an inverse problem of reconstructing a subdiffusion process from nonlocal data, Math. Methods Appl. Sci., 42 (2019), 2043–2052. https://doi.org/10.1002/mma.5498 doi: 10.1002/mma.5498
![]() |
[10] |
E. Mussirepova, A. Sarsenbi, A. Sarsenbi, The inverse problem for the heat equation with reflection of the argument and with a complex coefficient, Bound. Value Probl., 2022 (2022), 1–13. https://doi.org/10.1186/s13661-022-01675-1 doi: 10.1186/s13661-022-01675-1
![]() |
[11] |
A. Sarsenbi, A. Sarsenbi, Boundary value problems for a second-order differential equation with involution in the second derivative and their solvability, AIMS Math., 8 (2023), 26275–26289. https://doi.org/10.3934/math.20231340 doi: 10.3934/math.20231340
![]() |
[12] |
R. Brociek, A. Wajda, D. Slota, Inverse problem for a two-dimensional anomalous diffusion equation with a fractional derivative of the Riemann-Liouville type, Energies, 14 (2021), 1–17. https://doi.org/10.3390/en14113082 doi: 10.3390/en14113082
![]() |
[13] |
Z. C. Deng, F. L. Liu, L. Yang, Numerical simulations for initial value inversion problem in a two-dimensional degenerate parabolic equation, AIMS Math., 6 (2021), 3080–3104. https://doi.org/10.3934/math.2021187 doi: 10.3934/math.2021187
![]() |
[14] |
D. K. Durdiev, A. A. Rahmonov, Z. R. Bozorov, A two-dimensional diffusion coefficient determination problem for the time-fractional equation, Math. Methods Appl. Sci., 44 (2021), 10753–10761. https://doi.org/10.1002/mma.7442 doi: 10.1002/mma.7442
![]() |
[15] |
S. Kerbal, B. J. Kadirkulov, M. Kirane, Direct and inverse problems for a Samarskii-Ionkin type problem for a two dimensional fractional parabolic equation, Prog. Fract. Differ. Appl., 4 (2018), 147–160. https://doi.org/10.18576/pfda/040301 doi: 10.18576/pfda/040301
![]() |
[16] |
M. Kirane, S. A. Malik, M. A. Al-Gwaiz, An inverse source problem for a two dimensional time fractional diffusion equation with nonlocal boundary conditions, Math. Methods Appl. Sci., 36 (2013), 1056–1069. https://doi.org/10.1002/mma.2661 doi: 10.1002/mma.2661
![]() |
[17] |
S. A. Malik, S. Aziz, An inverse source problem for a two parameter anomalous diffusion equation with nonlocal boundary conditions, Comput. Math. Appl., 73 (2017), 2548–2560. https://doi.org/10.1016/j.camwa.2017.03.019 doi: 10.1016/j.camwa.2017.03.019
![]() |
[18] |
K. B. Sabitov, A. R. Zainullov, Inverse problems for a two-dimensional heat equation with unknown right-hand side, Russ. Math., 65 (2021), 75–88. https://doi.org/10.3103/S1066369X21030087 doi: 10.3103/S1066369X21030087
![]() |
[19] |
B. K. Turmetov, B. J. Kadirkulov, An inverse problem for a parabolic equation with involution, Lobachevskii J. Math., 42 (2021), 3006–3015. https://doi.org/10.1134/S1995080221120350 doi: 10.1134/S1995080221120350
![]() |
[20] |
B. K. Turmetov, B. J. Kadirkulov, On the solvability of an initial-boundary value problem for a fractional heat equation with involution, Lobachevskii J. Math., 43 (2022), 249–262. https://doi.org/10.1134/S1995080222040217 doi: 10.1134/S1995080222040217
![]() |
[21] |
A. Ahmad, D. Baleanu, On two backward problems with Dzherbashian-Nersesian operator, AIMS Math., 8 (2023), 887–904. https://doi.org/10.3934/math.2023043 doi: 10.3934/math.2023043
![]() |
[22] |
R. R. Ashurov, A. T. Muhiddinova, Inverse problem of determining the heat source density for the subdiffusion equation, Differ. Equ., 56 (2020), 1550–1563. https://doi.org/10.1134/S00122661200120046 doi: 10.1134/S00122661200120046
![]() |
[23] |
A. S. Berdyshev, B. J. Kadirkulov, On a nonlocal problem for a fourth-order parabolic equation with the fractional Dzhrbashyan-Nersesyan operator, Differ. Equ., 52 (2016), 122–127. https://doi.org/10.1134/S0012266116010109 doi: 10.1134/S0012266116010109
![]() |
[24] |
M. Kirane, A. A. Sarsenbi, Solvability of mixed problems for a fourth-order equation with involution and fractional derivative, Fractal Fract., 7 (2023), 1–12. https://doi.org/10.3390/fractalfract7020131 doi: 10.3390/fractalfract7020131
![]() |
[25] |
A. I. Kozhanov, O. I. Bzheumikhova, Elliptic and parabolic equations with involution and degeneration at higher derivatives, Mathematics, 10 (2022), 1–10. https://doi.org/10.3390/math10183325 doi: 10.3390/math10183325
![]() |
[26] |
K. Li, P. Wang, Properties for fourth order discontinuous differential operators with eigenparameter dependent boundary conditions, AIMS Math., 7 (2022), 11487–11508. https://doi.org/10.3934/math.2022640 doi: 10.3934/math.2022640
![]() |
[27] |
Y. Mehraliyev, S. Allahverdiyeva, A. Ramazanova, On one coefficient inverse boundary value problem for a linear pseudoparabolic equation of the fourth order, AIMS Math., 8 (2023), 2622–2633. https://doi.org/10.3934/math.2023136 doi: 10.3934/math.2023136
![]() |
[28] |
M. Muratbekova, B. Kadirkulov, M. Koshanova, B. Turmetov, On solvability of some inverse problems for a fractional parabolic equation with a nonlocal biharmonic Operator, Fractal Fract., 7 (2023), 1–18. https://doi.org/10.3390/fractalfract7050404 doi: 10.3390/fractalfract7050404
![]() |
[29] |
A. A. Kornuta, V. A. Lukianenko, Dynamics of solutions of nonlinear functional differential equation of parabolic type, Iz. VUZ Appl. Nonlinear Dyn., 30 (2022), 132–151. https://doi.org/10.18500/0869-6632-2022-30-2-132-151 doi: 10.18500/0869-6632-2022-30-2-132-151
![]() |
[30] |
A.V. Razgulin, A class of parabolic functional-differential equations of nonlinear optics, Differ. Equ., 36 (2000), 449–456. https://doi.org/10.1007/BF02754466 doi: 10.1007/BF02754466
![]() |
[31] |
B. Turmetov, V. Karachik, On eigenfunctions and eigenvalues of a nonlocal Laplace operator with multiple involution, Symmetry, 13 (2021), 1–20. https://doi.org/10.3390/sym13101781 doi: 10.3390/sym13101781
![]() |
[32] |
V. V. Karachik, A. M. Sarsenbi, B. K. Turmetov, On the solvability of the main boundary value problems for a nonlocal Poisson equation, Turkish J. Math., 43 (2019), 1604–1625. https://doi.org/10.3906/mat-1901-71 doi: 10.3906/mat-1901-71
![]() |
[33] |
B. Turmetov, V. Karachik, M. Muratbekova, On a boundary value problem for the biharmonic equation with multiple involutions, Mathematics, 9 (2021), 1–23. https://doi.org/10.3390/math9172020 doi: 10.3390/math9172020
![]() |
[34] |
V. Karachik, B. Turmetov, On solvability of some nonlocal boundary value problems for biharmonic equation, Math. Slovaca, 70 (2020), 329–342. https://doi.org/10.1515/ms-2017-0355 doi: 10.1515/ms-2017-0355
![]() |
[35] | K. Rektorys, Variational methods in mathematics, science and engineering, Dordrecht: Springer, 1977. https://doi.org/10.1007/978-94-011-6450-4 |
1. | Mohammed Elamine Beroudj, Abdelaziz Mennouni, Carlo Cattani, Hermite solution for a new fractional inverse differential problem, 2024, 0170-4214, 10.1002/mma.10516 | |
2. | Ying-Qing Liang, Fan Yang, Xiao-Xiao Li, A hybrid regularization method for identifying the source term and the initial value simultaneously for fractional pseudo-parabolic equation with involution, 2024, 1017-1398, 10.1007/s11075-024-01944-3 |