ON THE INVERSE PROBLEM OF THE BITSADZE–SAMARSKII TYPE FOR A FRACTIONAL PARABOLIC EQUATION Текст научной статьи по специальности «Математика»

20

Probl. Anal. Issues Anal. Vol. 12 (30), No3, 2023, pp. 20-40

DOI: 10.15393/j3.art.2023.14370

UDC 517.956.42

R. R. ASHUROV, B. J. KADIRKULOV, B. KH. TURMETOV

ON THE INVERSE PROBLEM OF THE BITSADZE-SAMARSKII TYPE FOR A FRACTIONAL PARABOLIC EQUATION

Abstract. In this paper, the inverse problem of the Bitsadze-Samarsky type is studied for a fractional order equation with a Hadamard-Caputo fractional differentiation operator. The problem is solved using the spectral method. The spectral aspects of the obtained problem are investigated, root functions are found, and their basis property is proved. The conjugate problem is investigated. The uniqueness and existence theorems for a regular solution to this problem are proved.

Key words: Hadamard-Caputo fractional operator, Riesz basis, Le Roy function, inverse problem

2020 Mathematical Subject Classification: 35K65, 34R30,

34K37, 34L10

1. Problem statement. The Hadamard integration operator of order a > 0 is the expression [1]

t a_ 1 h Joty(t) = I (ln ^f y(r)y, t > 0.

T(a)

o

If a = 0, we suppose that h J°ty(t) = y(t).

For a E (m — l,m],m = 1, 2,..., the following expression

i

1 f / + \ m—a— 1 /7

hcDamy(t) = h J£-a(tmy(t)) = -r ln - 5my(r)-L

T(m — a) J \ tJ t

o

is called the Hadamard-Caputo differentiation operator of the order a > 0,

l d_ 'dt

where 6 = tf and 5k = 5 • 5k-\k ^ 1 [2].

© Petrozavodsk State University, 2023

Let 0 < a ^ > 0. Consider the following equation:

t-l3HCD0tu(x, t) = uxx(x, t) + g(x). For the case a = 1, ¡3 = 1, we get

d \ du(t,x)

t-lHC D0tu(t,x) = i-1( t^u(t,x) =

dt

and equation (1) coincides with the classical parabolic equation.

Let Q = {(x,t): 0 < x < 1, 0 <t<T}, where T is a positive real number. For equation (1), consider the following problem in the domain Q:

Problem BS. Find a pair of functions (u(x,t), g(x)) from the class

u(x,t), t-l3HCD0tu(x,t), Uxx(x,t) E C(Q), g(x) E C[0,1], (2)

satisfying in the domain Q equation (1) and the conditions

u(x, 0) = <p(x), 0 ^ x ^ 1, (3)

u(x,T) = ^(x), 0 ^ x ^ 1, (4)

u(0,t) = 0, 0 ^ t ^ T, (5)

u(1,t) = u(x0,t), 0 ^ t ^ T. (6)

Here p(x),ip(x) are the given functions, [3,a,x0 are the given real numbers, such that ¡3 > 0, 0 < x0 < 1.

Among the first works devoted to the issues related to the solvability of nonlocal problems, we note the work of T. Carleman [3] (see also [4]), where the problem with a nonlocal condition, which consists in finding a holomorphic function in a bounded domain that connects the values of this function at different points of the boundary, is studied. This problem was reduced to a singular integral equation with deviation. We should also note the papers [5] and [6], where abstract nonlocal elliptic boundary value problems were studied.

A nonlocal boundary-value problem of a new type for an elliptic differential equation that arises in plasma theory was formulated by A. V. Bit-sadze and A. A. Samarskii [7]. This problem was reduced to an integral Fredholm equation of the second kind. Using the extremum principle for elliptic equations, the uniqueness of the classical solution is proved.

Further, the results on the theory of partial differential equations and functional-differential equations made it possible to study the solvability

problem for a wide class of nonlocal elliptic boundary value problems. In the monograph [8] (see also [9]), a detailed classification of nonlocal elliptic boundary value problems is given, the uniqueness and solvability of such problems in Sobolev spaces and weighted Kondratiev spaces are studied, the properties of the index, the spectral properties of the corresponding operators, the asymptotic behavior of solutions, and the smoothness of generalized solutions are considered. In addition, applications of nonlocal problems to the processes of heat distribution, diffusion, and cooling of aircraft engines are described.

Similar problems with operators of integer or fractional order with Riemann-Liouville, Caputo, and Hadamard-Caputo derivatives were studied in [10], [11], [12], [13], and for parabolic systems in [14], [15].

Note that inverse problems for a parabolic equation of fractional order with the Gerasimov-Caputo operator were also studied in [16], [17], and for degenerate equations, in [18], [19].

In this paper, an inverse problem of the Bitsadze-Samarskii type is studied for a degenerate fractional parabolic equation with the Hadamard-Caputo operator. Using the spectral method, the eigenvalues, as well as the corresponding root functions, are found, and their basis property is proved. The spectral issues of the conjugate problem are also investigated.

2. The Cauchy problem for a one-dimensional fractional differential equation.

In this subsection, we study the Cauchy problem for a one-dimensional fractional differential equation with a Hadamard-Caputo derivative.

Let 0 < a ^ 1, ¡3 > 0. We introduce the following operators:

Bi y(t) = t-Hc D0ty(t),B-3y(t) = J^y]®

and study some properties of these operators.

From the definition of operators h Joty(t) and hcD0ity(t), it follows that

h ) = ^ ,a,p> 0, hc DZt(r) = \0,a ^ a E (0,1]. (7)

\Hat^, ^ > 0,

Lemma 1. Let a E (0,1], p > 0, and f(t) E C[0,4 Then B-3f(t) E C[0, d] and the following estimate holds:

d3

P-13 f (t)L m^ p 11/WUc [o,d],

where ||/(t)^c[0,d] = v}^ \f(t)\-

Proof. Let f (t) E C[0, d]. Then, taking into account the definition of the operator B, we obtain

b-J f (t)

1 f «A t

T

r(a)

ln l)^ f (T) V

<

<

C[0,d]

r(a)

/ f \ a-i (ln-) Tl-1 dr

If (t)lc[Q,d] t1 r(a) f3a

Ta-1e-Tdr =

c[o,d] ti

¡3'

Thus, we get the following estimate:

w'mw^=m*x. \B-/S(f)\« lfmcM

l

0<t<d

n max t1 = —

!3a o^tKd ¡3a

C[0,d].

Lemma 1 is proved. □

Lemma 2. Let 0 < a ^ 1 and f(t) E C 1[0, d], Then hcD^J(t) E C[0, d], HCD0[tf (0) = 0, and the following estimate holds:

II HC'D0tf (i)|| C [0,d] ^11 f'(t)II C[0,d].

df(t). Hence, for f(t) E C1[

get hcD^f(t) E C[0,d]. It is obvious that hcD1 f(t)\t=o = tfj(t)\ =o = 0. d

Let 0 < a < 1. Then

Proof. If a =1, then HcD^tf(t) = tft f(t). Hence, for f(t) E C 1[0, d] we

' d ' dt-

H C

DawJ (t)

1

r(1 -a)

(ln ") Sf(r)

dr

<

<

1

t X — ^ r(1 -a) J Klnr

o

(ln-) \ f(r)\dr^

<

1

C[0,d]

r(1 - a)

(ln -) dr = t\\f

C[0,d].

Thus, we obtain

H hcD,tf(t)\\cM < d\\f'(t)\\cM,HCD0tf(t) |i=o = lino hcD&f(t) = 0. Lemma 2 is proved. □

Lemma 3. Let 0 < a ^ 1, /3 > 0 and f (t) E C 1[0,d]. Then the following equality holds:

B-J3 [-B3[f]} (t) = f(t) — f(0). (8)

Proof. As f (t) E C 1[0,d], Lemma 2 implies that the function B3[f](t) belongs to the class C[0, d]. If a = 1, then

B-3

t t *3w] « = jSB> Wr) 51 = I*(T)'t = f {t) — f {a,

oo For 0 < a < 1, taking into account the formula

HJMHJ3 = HJ0S+3, ^ 0

from [1], we get

t

B,-3 №[f]} (t) = f£j J (ln I)a'1r3B3a[f](r)^ =

0

= hJ3 [.nJ1t-a[rn] (t) = h4 [f] (t) = № — f(0).

Lemma 3 is proved. □

In the domain (0, d), consider the following Cauchy problem:

B3ay(t) = Xy(t), 0 <t<d, (9)

(0) = , (10)

where b is a real positive number.

A solution to this problem is a function from the class y(t) E [0, d], B3y(t) E C(0, d) satisfying equation (9) and condition (10) in the classical sense.

Let a function y(t) be a solution to problem (9), (10). Applying the operator y(t) to both parts of (9), we have

B-3 Hy] (t) = XB-3 [y](t), 0 <t<d.

Hence, taking into account (7) and (8), we obtain y(t) = b + XB-J [y](t), 0 <t<d.

Thus, if y(t) is a solution of problem (9), (10), then it satisfies the Volterra integral equation of the second kind of the type

t

M = K(t, r)y(r)dr + b, (11)

0

where K(t, t) = (in(t) — in(r))a~ V3-1.

To find a solution to the integral equation (11), we apply the method of normalized systems [20]. To do this, denote L1 = E, L2 = XB-3, where E is the unit operator. Then equation (11) can be rewritten as (L1 — L2)y(t) = b. As L1 = E, we get L-1 = E. Let g0 = b. Further, we use the technique used in [20]. According to this method, consider the system

gk (t) = (L-1 •L2)kg o,k = 1, 2,....

For k = 1, we get

a— 1

91(t) = (XB-)go = XB-3 [b] = J^J(]n ^ r3-1dr =

0

X X

— f^-1113-^-(3-1)Hetdi=-X^ [sa-h-Sds = Xa t13.

^ ^ Orv P ( ^ \ I fJrv

T(a) J s paT(a) J ft

00

Let k = 2. Then we get

X2 [ft

a-1

92(t) = (XB-J3)g 1(t) = b—— (in -) r23-1dr

0

OO

b[jaT{a)J ^ 1g ^ 2»[j2a tW. 0

Further, using the method of mathematical induction, for k ^ 2 the function gk (t) can be represented as

gk (t) = ( XB—f )kb

Xkb

( k\)af

ka

tkf.

Lemma 4. Let 0 < a ^ 1 and ft > 0. Then problem (9), (10) has a unique solution, which has the form

y(t)=bY;

Xk

fkf

k=0

( k\)aft

k a

'12)

Proof. Let us first show the convergence of the series (12). To do this, we evaluate the ratio 9k+j. Then we get

9k(t)

9k+i( t) Xk+l tP(k+1) Xktfk

9 k (t) f a(k+1)((k +1)\)a fak (k\)a

->• 0.

XI

Therefore, by the d'Alembert criterion, the series (12) converges uniformly for t E [0,d] (in general, for t ^ 0). Since the functions gk(t), k = 1,2,..., are continuous on [0, the sum of the series is also

continuous in t E [0, d].

Further, applying the operator B^ to the functions gk(t), and taking into account (7), we have

B"a9k (t)=ba

Xk

L( k\)afka

tkf

Bff go(t) = 0,

Xk

f3ka-1 ((k - 1)\)'

7t(k-1)f = Xgk—i(t), k > 1.

Then, formally applying the operator B3 to the series (12), we obtain ^ ^ ^

B^y (t) = Y.Bi9k (t) = x £ gk-! (t) = x £ gk (t) = Xy (t). (13)

k=0

k=l

k=0

This implies a uniform convergence of the series Y1 Bfgk (t) and ful-

k=0 a

fillment of the condition Bfy(t) E C[0, d]. It also follows from (13) that a function (12) satisfies equation (9). It is obvious that y(0) = b, i.e., condition (10) is also satisfied. Lemma 4 is proved. □

Remark. Function (12) can be represented as

where La(z) is written as

00 k

La(z) = ^ ,a> 0, zeR, (14)

k=o ( !)

and is called the Le Roy function [21].

Now, in the domain (0,d), we consider an inhomogeneous equation of the form

B33y (t) = Xy (t) + a, 0 <t<d,a = 0, (15)

and find its solution from the class y(t) E [0, d], B3y(t) E C(0,d), satisfying condition (10). Since the solution to equation (9) is known, then, according to the theory, it is sufficient to find a particular solution to equation (15). It is easy to see that such a solution is a constant function of the type y(t) = — |. Then the solution to problem (15), (10) is written as

/, a \ (\ tß\ a

"{t)=(b+ x)L°{W) -~y (16

XJ a\pa ) X

3. A spectral problem for a second-order differential equation.

To solve the BS problem, we apply the Fourier method, according to which we look for a non-trivial solution of a homogeneous equation ((g(x) = 0)), corresponding to (1) in the form u(x,t) = X(x) • T(t). Substituting it into equation (1) and using conditions (5) and (6) to find the eigenvalues X and eigenfunctions X(x), we obtain a problem for eigenvalues in the form

—X''(x) = XX(x), 0 < x < 1, (17)

X(0) = 0, X(1) = X(x0), 0 <x0 < 1. (18)

Further, we find the eigenfunctions of problem (17), (18). For X ^ 0, this problem has only the trivial solution, so consider the case with X > 0. In this case, we obtain two series of eigenvalues:

((2 n — 1W \2 A ( 2rnr \2

which correspond to eigenfunctions of the form

Xni (x) = sin^/X^x, Xn2 (x) = sin \fX^2x, n E N. (20)

Note that among the two series of eigenvalues from (19), there are coinciding ones. Indeed, comparing Xsi and Xm2 from (19), we obtain the following relation between x0, s and m:

2 s + 2m + 1

xo = '2s + 2m + 1' S'm EN; (21)

in this case, the corresponding values of Xsi and Xm2 coincide, so that the system of eigenfunctions (20) is not complete and the problem of supplementing this system with associated functions [22] arises.

Problem (17), (18) was first studied in [22], where the completeness of the root functions of the differential operator corresponding to this problem was shown.

The spectral problem in a more general formulation than the problem of the form (17), (18), is studied in [23]. A more general basis property criterion is found, which makes it possible to study the basis property of systems of eigenfunctions and associated functions of the above-mentioned problem, as well as of its adjoint problem.

Now, since relation (21) depends on the point x0, let us find out the nature of this point, thereby uttering the results of the aforementioned works.

Clearly, when x0 is an irrational number from the interval (0,1), relation (21) does not hold, and, so, both eigenvalues (19) and the corresponding eigenfunctions are different. The following example shows that (21) does not hold even for some rational values of x0. Indeed, consider the case of x0 = 2. Then, from (21) it follows 2s = 6m + 1. By virtue of the fact that m, s E N, 6m +1 is always odd, and 2s is even, that is, (21) does not occur under any m, E N. So, there are rational fractions in which all eigenfunctions are distinct. Now we give a criterion that shows the relationship between the values of x0, s and m at which the equality (21) takes place, and also gives an algorithm for finding the corresponding values of m, E N.

Let x0 be a rational number from the interval (0,1), such that x0 = |, p < q, and p, q- are coprime natural numbers. Here are some known facts about the number x0:

1) Numbers + and - at the same time are both either even or odd.

2) Let the number q — p (or p + q) be even; then either p + q or q — p is a multiple of 4.

Lemma 5. Let x0 E (0,1) be a rational number, such that x0 = |, p < q, p and q be odd coprime natural numbers, such that q — p is a multiple of 4. Then there exists countable numbers s and m, such that for two series of eigenvalues from (19) the equality XSl = Xm2 is valid.

Proof. Indeed, from the condition on xo, the relation (21) will take the following form:

S = m(P + ^ + 1-. (22)

q — p 2

By condition, q — p is a multiple of 4, then it is represented as q — p = 4r, where r E N. Hence, q = p + 4r. Then (22) takes the following form

s = ^^ + (23)

By condition, p, q are odd numbers. Then, obviously, the number p+2r is also odd, and the numbers p + 2r and 2r have no common divisors. Therefore, the condition E N can be provided only by choosing m, and it is easy to see that it is true if and only if m = k ■ r, where k is any odd natural number. Thus, at these values m from (23) we have:

_ k(p + 2r) 1

S = 2 + 2,

and, since k(p + 2r) is odd, it follows that s E N.

Finally, to find the values s and m, at which (21) is true, we obtained the following formulae:

m(p + q) 1 ,

m = k ■ r, s = ——^ + -, (24)

— 2

where r = , k = 1, 3, 5),.... □

Taking Lemma 1 and formula (21) into account, we divide the rational numbers from the interval (0,1) into two sets Q\ and Q2, where the set Q2 contains all rational numbers satisfying the conditions of Lemma 1, and Qi contains the remaining rational numbers.

Let us consider the case x0 E Q1. In this case, we obtain two series of eigenvalues (19), which correspond to eigenfunctions (20), and all these functions are different.

The spectral problem (17)—(18) is not self-adjoint, and, therefore, the system (20) is not orthogonal. Then, according to the spectral theory of operators, the problem arises of studying the problem adjoint to problem (17)-(18). It is easy to determine that the following problem is adjoint to it:

-Y''(x) = XY(x), x E (0,xo) U (xo, 1), (25)

Y (0) = 0,Y (1) = 0,

(26)

Y (xo + 0)= Y (xo - 0), Y'(1) = Y' (xo + 0) - Y'(xo - 0). (27)

Note that the solution of equation (25) that satisfies conditions (26) and (27) is found uniquely, that is, it has no extra conditions. This problem should be considered as two boundary-value problems with gluing conditions of the form (27).

As in the case of problem (17)-(18), consider the case x0 E Q1. Then, solving the problem (25)-(27), we find that the numbers (19) are also eigenvalues of this problem, and the corresponding eigenfunctions have the form

[Yni (x);Ym (x)},n EN, (28)

where

Yni (x)

1 + xo

sin^/Xnix, 0 ^ x ^ x0,

(xo + 1) cosyXi

sin ^JXni (1 - x) , xo ^ x ^ 1,

Yn2 (x)

0, 0 ^ x ^ xo, 2

(1 - xo) COSyXn~2

sin

^Jxn2 (1 - x) , xo ^ x ^ 1,

Lemma 6. The systems of the functions (20) and (28) are biorthonor-mal, i.e., it takes place

(XSi (x),Ym. (x))

ft

1, s = m,j = 1,

{oA) [0, s = m,j = 1, 2, 1, s = m,j = 2,

(Ys2 (x),XXmj (x))L (o,l) \ n , -TO

2 0, = m, = 1, 2.

4

The proof of Lemma 6 is carried out directly by calculating the corresponding integrals. From Lemma 6 it follows that the system {Yni (x);Yn2 (x)} is an biorthogonal adjoint to the system {Xni (x);Xn2 (x)}, n E N.

Lemma 7. Let x0 E Q1. Then:

1) The system of the root functions of the problems (17), (18) and (25)-(27) consists only of the eigenfunctions (i.e., there are no associated functions);

2) System (20) is complete and minimal;

3) Systems of eigenfunctions (20) and (28) form the Riesz basis in L2(0, 1).

Note that more detailed information about the Riesz bases can be found in [24]. Before proving the lemma, we present some definitions and facts concerning Riesz bases.

Definition 1. [24] A complete and minimal system of the functions {<pn(x)} is called a Bessel system, if for any f E L2(a, b) the series of squared coefficients of its biorthogonal expansion {'pn(x)} converges, i.e.,

f E L2(a, b) implies that E \(f,^n)L2ab)\ < &, where {^n(x)} is the

n=l

system of associated functions.

Definition 2. [24] A complete and minimal system of functions {pn(x)} is called a Hilbert system, if for any sequence of numbers cn, such that

IX

°n < there exists one and only one f E L2(a, b), for which these cn

k=1

are the coefficients of its biorthogonal expansion in {'pn(x)}, i.e.,

Cn = ( f, ^n^ L2(a, b), n =1,

Definition 3. [24] A complete and minimal system is called a Riesz basis, if it is both Bessel and Hilbert system at the same time.

Proof. The proof of part 1) follows from [22] and also from Lemma 1. The completeness of system (20) was proved in [22], and the minimality follows from Lemma 6. Thus, systems (20) and (28) satisfy all three conditions A from [23]. Let us now verify the fulfillment of the conditions for the basis property of the main theorem from [23]. As follows from this theorem, for this it suffices to prove that the following two conditions are satisfied:

a) E 1 <C1,i = 1, 2, for any real A ^ 0,

b) \\Xki ||o ■ \\Yki ||o ^C2,i = 1, 2, for all numbers k.

The fulfillment of the first condition is verified directly. To check the second condition, we calculate the corresponding norms:

where n E N.

Then, since xq is a fixed rational number, the fulfillment of condition b) follows from Lemma 5, as well as from the finiteness of the set of values cos \ni and cos \n2. Thus, systems (20) and (28) satisfy all conditions of the main theorem from [23]. Hence, it follows that these systems form a Riesz basis. Lemma 7 is proved. □

The spectral questions of problems (17), (18) and (25)-(27) are studied in the same way for values x from Q2. Note that in this case, these problems, in addition to eigenfunctions, also have associated functions that correspond to those eigenvalues, whose serial numbers are determined by formulae (24).

4. Existence and uniqueness of a solution to Problem BS.

Now we turn to the study of the existence and uniqueness of the solution of problem BS. According to the theory, we will seek the solution u(x, t), g(x) of problem in the form of an expansion in a specially chosen basis from the system of functions {Xn.}, i= 1, 2 from (20):

<x

u(x,t) (t) • Xm (x) +Un2 (t) • Xn2 (x)), (29)

n=1 <x

9(x) = Y1 (gni

• Xni (x) + 9n2 • Xn2 (x)), (30)

n=1

where uni(t),nn2(t) are unknown functions, gni, gn2 are unknown constants, n E N.

Substituting (29) and (30) into equation (1), we obtain the following equations for finding the functions uni (t) , un2 (t) and constants gni, gn2 :

t-13 HcD^Um (t) + \niUni (t) = gm ,i = 1, 2. (31)

From representation (29), taking into account conditions (3) and the completeness of system (20), we obtain that the unknown functions uni (t), un2 (t) satisfy the conditions

Un,(0) = ^, i=l, 2, (32)

where pni, pn2 are expansion coefficients of the function p(x) in terms of the system of functions (20), which are found by the formulas

= ip(x)Yn. (x)dx, (33)

and Yni (x), i = 1, 2, are functions defined by formulas (28).

Applying the operator B-3 to both parts of equation (31), taking into account formulas (15) and (16), we find that the solution to equation (31) that satisfies the condition (32) has the form

\n,tß

Un, = ^ •L4 - + ^ [1 -L4 - ^jj, ¿ = 1, 2, (34)

where La (z) is the Le-Roy function, which has the form (14).

Now we find the unknown constants gni,gn2. To do this, we expand the function ip(x) into a series in terms of the system of functions (20):

\n,tß

4>{x) = ^2, (^ni ■ Xm (x) + Vpn2 ■ Xn2 (x)),

n=l

where ^ni,vpn2 are expansion coefficients, i.e.,

i

^ni = ^(x)Yni(x)dx, i=l, 2. (35)

vni

0

Further, from representation (29), taking into account condition (4), we obtain that the unknown functions uni(t), i = 1, 2, also satisfy the conditions

u,ni (T) = ^m, i = 1, 2. (36)

i

Then, taking into account the conditions (36) and using (34) to find the constants gni,gn2, we obtain the following equations:

( - ^ ) + £ I1 "M - f )]=*

1, 2.

From here we find

Ki - <pniLa( " ) )

V V ß , x = 1,2. (37)

1- L ( -

Substituting all this into the expression for functions un\(t),un2(t), we obtain

L ( ^¿P \_t ( ^ilP \ i _ L ( ^

Ln\ pa I La\ pa j 1 Ln\

a\ ßT I ßT I ßT .

uni [t) = -7-—JV-+-7-—JY*ni ^ = 11, 2.

1 !ßr-) 1 - -fßT-j

(38)

Thus, a formal solution of the problem is found in the form of series (29) and (30), where the coefficients gni,gn2 and functions uni (t),un2 (t) are determined, respectively, by formulas (37) and (38).

Now we will consider the proof of theorems on the existence and uniqueness of a solution to the B S problem.

5. The uniqueness of a solution to Problem BS. In proving the uniqueness and existence of a solution to problem B S, we use the following estimate for the Le-Roy function La (z) obtained in [21].

Lemma 8. Let 0 < a < ft < 1. Then for any z ^ 0 the estimate is valid:

e-z ^ Lp(-z) < La(-z) ^ .

1 + Z

From here it is easy to obtain the following estimate: Corollary 1. For any z > £ > 0 it holds that

^-^ C < oo.

1 - La(-z)

Theorem 1. If a solution to the BS problem exists, it is unique.

Proof. Let us assume the opposite. Let there be two solutions {u1(x, t), g1(x)} and {u2(x, t),g2(x)} to the BS problem.

Introduce the notation U(x, t) = U1(x, t) — U2(x, t) and g(x) = gi(x) — g2(x). Then the functions U(x, t) and g(x) satisfy the equation

t-l3HcD0;tU = Uxx + g(x) (39)

and conditions

u(x, 0) = 0, U(x, T) = 0, 0 ^x ^ 1, (40)

u(0, t) = 0, u(1, t) = u(x0, t), 0 ^t^T. (41)

Consider the functions

i

Uni (t) = j U(x, t)Yni (x)dx, i = 1, 2, (42)

0

where functions Yni(x),i = 1, 2 are determined by formulas (28).

Applying operator Bto both parts of equality (42) and taking into account equation (39), as well as conditions (40), (41), we conclude that the function Uni (t) and constant gni satisfy the following equation and conditions:

t-1 HcD0>tUni(t) + Xniuni (t) = gni, Urn(0) = 0, Urn(T) = 0, (43)

where

9ni = (9(x),Ym (x))o.

From (31) and (34), we obtain that the solution of this equation, which satisfies the first boundary condition in (43)), has the form

~ gn T1 r / Xmt \

Un ® = XT 1 — L< — ^

, = 1, 2.

Hence, satisfying the second boundary condition in (43), we obtain

OnM Xn

1 — L«( — ^)

P C

0,i = 1, 2.

As La (0) = 1, Xm ^ Xu > 0, we get La( — = 1. Then it follows

that Uni(t) = 0,gni = 0. Consequently, problem (43) has only a trivial solution.

As a result, we obtain that for any fixed t E [0,T] functions u(x, t), g(x) are orthogonal to system (28), which is complete in L2(0,1). Then u(x, t) = 0, g(x) = 0 almost everywhere in Q and [0,1], respectively. Since u E C(Q), f(x) E C[0,1], from here we get that u(x, t) = 0,g(x) = 0. The uniqueness of the solution of problem is proved. □

6. The existence of a solution to Problem BS. Let us prove the existence of a solution to the problem.

Theorem 2. Let the functions <p(x), 1(x) satisfy the conditions

<p(x) E C4[0,1], <p(0) = <p"(0) = 0, V(1) = V(x0), >/(1) = V"(x0),

iKx) E C4[0,1], ^(0) = 4>"(0) = 0, ^(1) = 4>(x0), r(1) = 4>"(x0). Then a solution to problem BS exists.

Proof. Since the system (20) forms the Riesz basis in the space L2(0,1), the functions u(x, t) and g(x) can be represented in the form (29) and (30), where the coefficients gni,gn2 and functions uni (t), un2 (t) are determined, respectively, by formulas (37), (38).

It is easy to show by direct calculation that the functions u(x, t) and g(x), defined by series (29) and (30), satisfy equation (1) and conditions (3)-(6). It remains to show that the functions u(x, t) and g(x) are from class (2).

Let us show that uxx E C(Q). From (29), differentiating twice with respect to the variable x, we obtain

uxx(x,t) = ^2 (-Xni • uni (t)Xni (x) - Xn2 • un2 (t)Xn2 (x)). (44)

n=1

As lX,ni(x) I ^ 1, i = 1, 2, it follows that

<x

| u x x ( x,

t)I ^ ( X ni | uni (t)I +Xn2 I un (t)I). (45)

n=1

Let us estimate the functions uni (t) and un2 (t). Taking Lemma 1 and Corollary 1 into account, from (38) we get

K(t)I < fei I + lini I) , 1=1, 2. (46)

Here and below C is positive (in general, different) constant.

Further, taking into account (46) and the conditions of Theorem 2, integrating by parts the expressions for the coefficients , , % = 1, 2 in (33), (35), we obtain from (45)

k, (x,t)i ^ jr (i ci + ici)+Tl (i ci + id))

n=1 \Ani ^n2 J

where

1 1

C = i fIV (x)YtH (x)dx^ = -L f ijIV (x)Yni (x)dx, i = 1, 2.

Am J Ani J

o 1 0

Thus, the series (44) is majorized by the series

ro

Z 1 (I rf?| + №?| + K>| + K>I).

n=1

whose convergence follows from the Cauchy-Schwarz inequality, as well as from the convergence of the series E Iv^T, E |Ct, ^ = 1, 2.

n=1 n=1

Then, according to the Weierstrass theorem [25], the series (44) converges absolutely and uniformly in the domain Q, and its sum is a continuous function in this domain. In the same way, it is shown that t-fHCD%tu(x, t) e C(Q), and g(x) e C[0,1] follows from the fact that t-fHCD0>tu(x, t) e C(Q), and from equation (1). □

Acknowledgment. The research of the third author is supported by the grant of the Committee of Sciences, Ministry of Education and Science of the Republic of Kazakhstan, project AP09259074.

References

[1] Kilbas A. A, Srivastava H. M., Trujillo J. J. Theory and applications of fractional differential equations. North-Holland Mathematics Studies, 2006.

[2] Jarad F., Abdeljawad T., Baleanu D. Caputo-type modification of the Hadamard fractional derivatives. Adv. Differ. Equat., 2012, vol. 2012, no. 142, pp. 1-8.

DOI: https://doi.org/10.1186/1687-1847-2012-142

[3] Carleman T. La theorie des equations integrales et ses applications. Annales de l'l. H. P., 1930, vol. 1, no. 4, pp. 401-430.

DOI: http://www.numdam.org/item?id=AIHP_1930__1_4_401_0

[4] Garif'yanov F. N., Strezhneva E. V. On a Carleman problem in the case of a doubly periodic group. Probl. Anal. Issues Anal., 2022, vol. 11(29), no 3, pp. 45-55.

DOI: http://10.15393/j3.art.2022.11390

[5] Vyshik M. I. On strongly elliptic system of differential equations., Mat. sb., 1951, vol. 29, no. 3, pp. 615-676 (in Russian).

[6] Browder F. Non-local elliptic boundary value problems. Amer. J. Math. 1964. vol. 86, no. 4, pp. 735-750.

DOI: https://doi.org/10.2307/2373156

[7] Bitsadze A. V., Samarskii A. A. Some elementary generalizations of linear elliptic boundary value problems. Dokl. AN SSSR. 1969,vol. 185, no. 7, pp. 739-740 (Russian).

[8] Skubachevskii A. L. Nonclassical boundary value problems. I. J. Math. Sci., 2008, vol. 155, no. 2, pp. 199-334.

DOI: https://doi.org/10.1007/s10958-008-9218-9

[9] Skubachevskii A. L. Nonclassical boundary value problems. II. J. Math. Sci., 2010, vol. 166, no. 4, pp. 377-561.

DOI: https://doi.org/10.1007/s10958-010-9873-5

[10] Ashyralyev A., F. S. Ozesenli Tetikoglu. it A note on Bitsadze-Samarskii type nonlocal boundary value problems, Well-posedness. Numer. Funct. Anal. Optim. 2013, vol. 34(9), pp. 939-975.

DOI: https://doi.org/10.1080/01630563.2012.738458

[11] Shamsi H. Al., Kadirkulov B. J., Kerbal S. The Bitsadze-Samarskii type problem for mixed type equation in the domain with the deviation from the characteristics. Lobachevskii J. Math., 2020, vol. 41, no. 6, pp. 1021-1030. DOI: https://doi.org/10.1134/S1995080220060025

[12] Khubiev K. U., The Bitsadze-Samarskii problem for some characteristically loaded hyperbolic-parabolic equation. Journal of Samara State Technical University, Ser. physical and mathematical sciences. 2019, vol 23, no. 4, pp. 789-796. DOI: https://doi.org/10.14498/vsgtu1677

[13] Karachik V. V., Turmetov B. Kh. On solvability of some nonlocal boundary value problems for biharmonic equation. Math. Slovaca. 2020, vol. 70, no. 2, pp. 329-342. DOI: https://doi.org/10.1515/ms-2017-0355

[14] Baderko E. A., Cherepova M. Solution of the Bitsadze-Samarskii problem for a parabolic system in a semibounded domain. Differ. Equ., 2017, vol. 53,

no. 10, pp. 1294-1302.

DOI: https://doi.org/10.1134/S0012266117100068

[15] Baderko E. A., Cherepova M. Bitsadze-Samarskii problem for parabolic systems with Dini continuous coefficients. Complex Var. Elliptic Equ., 2018, vol. 64, no. 5, pp. 753-765.

DOI: https://doi.org/10.1080/17476933.2018.1501039

[16] Ashurov R. R., Muhiddinova O. Inverse problem of determining the heat source density for the subdiffusion equation. Differ. Equ., 2020, vol. 56, no. 12, pp. 1550-1563.

DOI: https://doi.org/10.1134/S00122661200120046

[17] Alimov Sh. A., Ashurov R. R. On the backward problems in time for time-fractional subdiffusion equations. Fract. Differ. Calc., 2021, vol. 11, no. 2, pp. 203-217.

DOI: https://dx.doi.org/10.7153/fdc-2021-11-14

[18] Kirane M., Sadybekov M. A., Sarsenbi A. A. On an inverse problem of reconstructing a subdiffusion process from nonlocal data. Math. Methods Appl. Sci., 2019, vol. 42, no. 6, pp. 2043-2052.

DOI: https://doi.org/10.1002/mma.5498

[19] Turmetov B. Kh., Kadirkulov B. J. On a problem for nonlocal mixed-type fractional order equation with degeneration. Chaos Solitons Fractals, 2021, vol. 146, pp. 1-5.

DOI: https://doi .org/10.1016/j.chaos.2021.110835

[20] Turmetov B. Kh. On a method for constructing a solution of integraldifferential equations of fractional order. El. J. Qualit. Theory Differ. Equat.,2018, vol. 25, pp. 1-14.

DOI: https://doi.org/10.14232/ejqtde.2018.1.25

[21] Boudabsa L., Simon T. Some Properties of the Kilbas-Saigo Function. Mathematics. 2021, vol. 9, no. 217, pp. 1-24.

DOI: https://doi .org/10.3390/math9030217

[22] Eroshenkov E.P. The spectrum of the Bitsadze-Samarskii problem. Differ-entsial'nye Uravneniya, 1983, vol. 19, no. 1, pp. 169-171 (in Russian).

[23] Il'in V.A. Necessary and sufficient conditions for being a Riesz basis of root vectors of second-order discontinuous operators. Differentsial'nye Uravneniya, 1986, vol. 22, no. 1, pp. 2059-2071 (in Russian).

[24] Bari N. K. Biorthogonal systems and bases in Hilbert space. Mathemat-ics.Uch. Zap. Mosk. Gos. Univ., 1951, vol 4, 69-107 (in Russian).

[25] Il in V. A., Poznyak E. G. Basics of mathematical analysis, Part II. Fiz-matlit, Moscow, 1973 (in Russian).

Received June 19, 2023. Accepted July 10, 2023. Published online August 6, 2023.

R. R. Ashurov

Institute of Mathematics, Uzbekistan Academy of Science, University street 9 , Tashkent,100174, Uzbekistan School of Engineering, Central Asian University, 264, Milliy Bog St., 111221, Tashkent, Uzbekistan E-mail: ashurovr@gmail.com

B. J. Kadirkulov

Tashkent State University of Oriental Studies, Amir Temur Str. 20, Tashkent,100060, Uzbekistan E-mail: kadirkulovbj@gmail.com

B. Kh.Turmetov

Khoja Akhmet Yassawi International Kazakh-Turkish University, Bekzat Sattarhanov ave., 29, Turkistan, 161200, Kazakhstan E-mail: batirkhan.turmetov@ayu.edu.kz