Processing math: 100%
Research article Special Issues

On solvability of some inverse problems for a nonlocal fourth-order parabolic equation with multiple involution

  • 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

    Related Papers:

    [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 Π={xRn:0<xj<pj,j=1,,n} is a parallelepiped, pj>0, and T>0. Consider the mappings Sj:RnRn,1jn of the type Sjx=(x1,,xj1,pjxj,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 0i<2n in the binary number system, the representation i(ini1)2=i1+2i2++2n1in, where ik=0,1, is valid. Therefore, introducing the vector i=(i1,,in), we can consider mappings of the type SiSi11Sinn corresponding to the index i. Using these mappings, we introduce the operator

    Lxv(x)=2n1i=0aiΔ2v(Six),

    where a0,a1,a2,a3,,a2n1 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),0tT, (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,k0;X2k(x)=sinkx,k1,

    and the corresponding eigenvalues are the numbers

    λ1k=(1ε)(k+12)2,k0;λ2k=(1+ε)k,k1.

    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, n3, 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,,2n1, 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)nj=1sinkjπxjpj, (2.3)

    where k=(k1k2...kn)Nn and C(n,p)=2n/2nj=11pj,p=(p1,,pn).

    The following statement is proved in [35].

    Lemma 2.1. The system of functions {vk(x):kNn} is orthonormal and complete in space L2((0,p1)×(0,p2)×...×(0,pn)).

    Let us introduce the following numbers εk=2n1i=0(1)i(k+e)ai, where k=(k1,k2,...,kn), i=(i1,,in), ik=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=2n1i=0(1)ikai.

    Theorem 2.1. Let the conditions εk0 be satisfied for all kNn, then the system of functions {vk(x):kNn} is a system of eigenfunctions of Problem S. The corresponding eigenvalues are determined by the equalities

    λk=εkμ2k,μk=π2nj=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

    2vkx2i=C(n,p)nj=1,jisinkjπxjpj2x2i(sinkiπxipi)=(kiπpi)2C(n,p)nj=1,jisinkjπ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)=(ni=1(kiπpi)2)vk(x),

    which means that

    Δ2vk(x)=μ2kvk(x),μk=π2(ni=1k2ip2i).

    Hence, for any m{1,2,...,n}, we get

    Δ2vk(Smx)=μ2kvk(x1,...,xm1,pmxm,xm+1,...,xn)=μ2kC(p)sinkmπ(pmxm)pmnj=1,jmsinkjπxjpj=μ2kC(p)(1)km+1sinkmπxmpmnj=1,jmsinkjπxjpj=μ2kC(p)(1)km+1nj=1sinkjπxjpj=(1)km+1μ2kvk(x).

    Let 0i2n1, which corresponds to vector i=(i1,,in). If im=1, then

    Δ2vk(Simmx)=μ2kC(p)sinkmπ(pmxm)pmnj=1,jmsinkjπ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(Si11Sinnx)=(1)i1(k1+1)+i2(k2+1)+...+in(kn+1)μ2kvk(x)
    =(1)i(k+1)μ2kvk(x)=(1)ikμ2kvk(x),

    where k=(k1,k2,...,kn), i=(i1,,in), ik=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)=2n1i=0aiΔ2vk(Six)=2n1i=0ai(1)ikμ2kvk(x)=μ2kvk(x)(2n1i=0(1)ikai)=ε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 kNn, 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(p1x1,x2)+a2Δ2v(x1,p2x2)+a2Δ2v(p1x1,p2x2)=λv(x1,x2),

    and the boundary conditions are written as

    v(0,x2)=v(p1,x2)=0,0x2p2;v(x1,0)=v(x1,p2)=0,0x1p1,vx1x1(0,x2)=vx1x1(p1,x2)=0,0x2p2;vx2x2(x1,0)=vx2x2(x1,p2)=0,0x1p1.

    Eigenfunctions are specified in accordance with (2.3) as

    v(m,k)(x)=2p1p2sinmπ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)ma1+(1)ka2+(1)k+ma3.

    More precisely, ε(m,k) can be written as

    ε(2m1,2k1)=a0+a1+a2+a3;ε(2m1,2k)=a0+a1a2a3;ε(2m,2k1)=a0a1+a2a3;ε(2m,2k)=a0a1a2+a3,

    where m,k=1,2,....

    Let

    h(x)=k1=1...kn=1hk1...knvk1...kn(x)=kNnhkvk(x) (3.1)

    be the Fourier series expansion of the function h(x) by the system {vk(x):kNn}, 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, 1jn in ˉΠ. If the conditions

    h(0,x2,...,xn)=h(p1,x2,...,xn)=0,0xjpj,j=2,...,n,h(x1,0,...,xn)=h(x1,p2,...,xn)=0,0xjpj,j=1,...,n,j2,h(x1,...,xn1,0)=h(x1,...,xn1,pn)=0,0xjpj,j=1,...,n1 (3.3)

    are satisfied, then the number series kNn|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=Cp10...pn0h(x1,x2,...,xn)x1v1,0,..,0k1...kn(x1,x2,...,xn)dx1...dxn,v1,0,..,0k1...kn(x1,x2,...,xn)=cosk1πx1p1nj=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=Cp10...pn0nh(x1,x2,...,xn)x1...xnv1,1,..,1k1...kn(x1,x2,...,xn)dx1...dxn,v1,1,..,1k1...kn(x1,x2,...,xn)=nj=1coskjπxjpj.

    Using the Cauchy-Bunyakovsky inequality, we obtain

    kNn|hk|k1=1...kn=1|1k1...knh1,1,..,1k1...kn|k1=1...kn=11k21...k2nk1=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...xnL2(Π), 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...knp10...pn0nh(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=Cp10...pn0n+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)=nj=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...k2nk1=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...xnL2(Π), 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

    kNnλ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(nj=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)=kNnuk(t)vk(x)=k1=1kn=1uk1...kn(t)vk1...kn(x), (4.1)
    f(x)=kNnfkvk(x)=k1=1kn=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)=Ckeλkt+fkt0eλ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)=Ckeλ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(1eλkt)

    and uk(t) are represented as

    uk(t)=Ckeλkt+fkλk(1eλkt).

    If we use condition (4.4), we obtain

    Ck=uk(0)=φk,ψk=uk(T)=φkeλkT+fkλk(1eλkT),

    then we find

    fk=λkψkφkeλkT1eλkT. (4.6)

    The solution to problem (4.3), (4.4) is represented as

    uk(t)=eλkteλkT1eλkTφk+1eλkt1eλkTψk. (4.7)

    Hence, the solution to Problem 1 can be written as

    u(t,x)=k1=1...kn=1[eλkteλkT1eλkTφk+1eλkt1eλkTψk]vk(x), (4.8)
    f(x)=k1=1...kn=1λkψkφkeλkT1eλ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 11eλkT and eλkT1eλkT are bounded. For the coefficients fk from equality (4.6), we obtain

    |fk|λk11eλkT|ψk|+λkeλkT1eλ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 0tT, there are estimates

    |eλkteλkT1eλkT|=eλkT|eλk(Tt)11eλkT|C,|1eλkt1eλkT|C.

    The series (4.8) satisfies the estimate

    |k1=1...kn=1[eλkteλkT1eλkTφk+1eλkt1eλkTψk]vk(x)|Ck1=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λkt1eλkTφk+λkeλkt1eλkTψk]vk(x). (4.10)

    If f(λ)=λeλt, then maxλ0f(λ)=f(1/t)=1te1, which means that λeλt1δe for tδ. Therefore, the series (4.10) satisfies the estimate

    |ut(t,x)|Ck1=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δ=[δtT]×ˉΠ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)}kNn. 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,...,2n1 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)=φkeλkT+fkgk(T).

    Thus, if for all kNn, the condition gk(T)0, then

    fk=1gk(T)[ψkφkeλkT], (4.11)

    and

    uk1...kn(t)=[1gk1...kn(t)gk1...kn(T)eλk1...kn(Tt)]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, kNn 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)=fmgm(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 kNn.

    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,...,2n1 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)|g0C, 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(ξ)1eλkTλk.

    This gives us the following lower estimate:

    |gk1...kn(T)|=|g(ξ)|1eλk1...knTλk1...knCλk1...kn. (4.13)

    Using inequality (4.13) for fk from equality (4.11), we get

    |fk|1|gk(T)||ψkφkeλkT|Cλk(|φk|+|ψk|).

    Similarly from equality (4.12), we obtain for uk(t),

    |uk(t)||eλktgk(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=1uk(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...kneλkt+fkt0eλ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[φkeλkt+fkt0eλk(tτ)g(τ)dτ]vk(x0)=φ0(t)+k1=1...kn=1[fkt0eλk(tτ)g(τ)dτ]vk(x0)=k1=1...kn=1φk1...kneλktvk(x0)+k1=1...kn=1[fkt0eλk(tτ)g(τ)dτ]vk(x0),

    where

    φ0(t)=k1=1...kn=1φkeλ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φkeλktvk(x0).

    As for the points of the domain ˉQ, we have the estimate |eλktvk(x0)|C, then

    |φ0(t)|k1=1kn=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
  • This article has been cited by:

    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
  • Reader Comments
  • © 2024 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

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1287) PDF downloads(69) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog