Next Article in Journal
Degenerate Derangement Polynomials and Numbers
Previous Article in Journal
Solutions of Bernoulli Equations in the Fractional Setting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inverse Problem for a Partial Differential Equation with Gerasimov–Caputo-Type Operator and Degeneration

by
Tursun K. Yuldashev
1,*,† and
Bakhtiyar J. Kadirkulov
2,†
1
Uzbek-Israel Joint Faculty of High Technology and Engineering Mathematics, National University of Uzbekistan, Tashkent 100174, Uzbekistan
2
Department of Mathematics and Information Technologies, Tashkent State University of Oriental Studies, Tashkent 100060, Uzbekistan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 5 May 2021 / Revised: 10 June 2021 / Accepted: 17 June 2021 / Published: 21 June 2021
(This article belongs to the Section Mathematical Physics)

Abstract

:
In the three-dimensional open rectangular domain, the problem of the identification of the redefinition function for a partial differential equation with Gerasimov–Caputo-type fractional operator, degeneration, and integral form condition is considered in the case of the 0 < α 1 order. A positive parameter is present in the mixed derivatives. The solution of this fractional differential equation is studied in the class of regular functions. The Fourier series method is used, and a countable system of ordinary fractional differential equations with degeneration is obtained. The presentation for the redefinition function is obtained using a given additional condition. Using the Cauchy–Schwarz inequality and the Bessel inequality, the absolute and uniform convergence of the obtained Fourier series is proven.

1. Introduction

When the boundary of a domain of a physical process is impossible to study, a nonlocal condition of integral form can be obtained as additional information sufficient for the unique solvability of the problem. Therefore, in recent years, research has intensified on the nonlocal direct and inverse boundary value problems for differential and integro-differential equations with integral conditions (see, for example, [1,2,3,4,5,6,7,8,9,10,11,12,13,14]). The many problems of gas dynamics, the theory of elasticity, and the theory of plates and shells have been described by high-order partial differential equations.
Fractional calculus plays an important role in mathematical modeling in many scientific and engineering disciplines [15,16,17,18]. In [19] where problems of continuum and statistical mechanics are considered. The construction of various models of theoretical physics using fractional calculus is described in ([20], Vol. 4, 5), [21,22]. A detailed review of the application of fractional calculus to solving problems in applied sciences is provided in ([23], Vol. 6–8), [24]. In [25], an inverse problem to determine the right-hand side for a mixed type integro-differential equation with fractional order Gerasimov–Caputo operators is considered. The problem of determining the source function for a degenerate parabolic equation with the Gerasimov–Caputo operator was investigated [26]. In [27], the solvability of the nonlocal boundary problem for a mixed-type differential equation with a fractional-order operator and degeneration is studied. In the applications of fractional derivatives to solving partial differential equations, interesting results have been obtained [28,29,30,31,32].
We recall some materials from the theory of fractional order integro-differentiation. Let 0 ; T be an interval on the set of non-negative real numbers, 0 < T < . The Riemann–Liouville 0 < α -order fractional integral for the function η ( t ) has the form
I 0 + α η ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 η ( s ) d s , α > 0 , t 0 ; T ,
where Γ ( α ) is the Gamma function.
For the case n 1 < α n , n N , the Riemann–Liouville α -order fractional derivative for the function η ( t ) is defined as follows:
D 0 + α η ( t ) = d n d t n I 0 + n α η ( t ) , t 0 ; T .
The Gerasimov–Caputo α -order fractional derivative for the function η ( t ) is defined by the following formula
* D 0 + α η ( t ) = I t 0 + n α η ( n ) ( t ) = 1 Γ ( n α ) 0 t η ( n ) ( s ) d s ( t s ) α n + 1 , t 0 ; T .
These derivatives are reduced to the nth order derivatives for α = n N :
D 0 + n η ( t ) = * D 0 + n η ( t ) = d n d t n η ( t ) , t 0 ; T .
In this paper, for the case of the 0 < α 1 order, we study the regular one value solvability of the inverse boundary value problem for the Gerasimov–Caputo-type fractional partial differential equation with degeneration. This partial differential equation is a fractional-order ordinary differential equation with respect to the first argument and is a higher even-order partial differential equation with respect to spatial arguments. The stability of the solution on the given functions is proved.
So, in the three-dimensional open domain Ω = { ( t , x , y ) | 0 < t < T , 0 < x , y < l } , a partial differential equation of the following form is considered
D t , x , y β , α , 4 , 4 [ U ( t , x , y ) ] = a ( t ) b ( x , y )
with a nonlocal condition on the integral form
U ( T , x , y ) + I 0 + ρ U ( t , x , y ) t = T = φ ( x , y ) , 0 x , y l ,
where ρ , T, and l are given positive real numbers,
D t , x , y β , α , 4 , 4 [ U ] = * D 0 + α + ε D 0 + α 4 k x 4 k + 4 k y 4 k +
+ ω t β 4 k x 4 k + 4 k y 4 k U ( t , x , y ) ,
where ω and β are non-negative parameters, ε is a positive parameter, ε > δ > 0 , δ = const , 0 < α 1 , k is a given positive integer, a ( t ) C ( Ω T ) , Ω T [ 0 ; T ] , Ω l [ 0 ; l ] , b ( x , y ) C Ω l 2 is a known function, and φ ( x , y ) is a redefinition function, Ω l 2 Ω l × Ω l . We assume that for the given functions, the following boundary conditions are true
φ ( 0 , y ) = φ ( l , y ) = φ ( x , 0 ) = φ ( x , l ) = 0 ,
b ( 0 , y ) = b ( l , y ) = b ( x , 0 ) = b ( x , l ) = 0 .
Statement of Problem. We find the pair of functions { U ( t , x , y ) ; φ ( x , y ) } , the first of which satisfies the partial differential Equation (1), nonlocal integral condition (2), and boundary value conditions
U ( t , 0 , y ) = U ( t , l , y ) = U ( t , x , 0 ) = U ( t , x , l )
= 2 x 2 U ( t , 0 , y ) = 2 x 2 U ( t , l , y ) = 2 x 2 U ( t , x , 0 ) = 2 x 2 U ( t , x , l )
= 2 y 2 U ( t , 0 , y ) = 2 y 2 U ( t , l , y ) = 2 y 2 U ( t , x , 0 ) = 2 y 2 U ( t , x , l ) = . . .
= 4 k 2 x 4 k 2 U ( t , 0 , y ) = 4 k 2 x 4 k 2 U ( t , l , y ) = 4 k 2 x 4 k 2 U ( t , x , 0 ) = 4 k 2 x 4 k 2 U ( t , x , l )
= 4 k 2 y 4 k 2 U ( t , 0 , y ) = 4 k 2 y 4 k 2 U ( t , l , y ) = 4 k 2 y 4 k 2 U ( t , x , 0 ) = 4 k 2 y 4 k 2 U ( t , x , l ) = 0 ,
class of functions
U ( t , x , y ) C ( Ω ¯ ) , * D 0 + α U ( t , x , y ) C x , y 4 k , 4 k ( Ω ) C x , y 4 k + 0 ( Ω ) C x , y 0 + 4 k ( Ω )
and additional condition
U ( t 1 , x , y ) = ψ ( x , y ) , 0 < t 1 < T , 0 x , y l ,
φ ( x , y ) C [ 0 ; l ] 2 , where ψ ( x , y ) are a given smooth function and
ψ ( 0 , y ) = ψ ( l , y ) = ψ ( x , 0 ) = ψ ( x , l ) = 0 ,
C x , y 4 k + 0 ( Ω ) is the class of continuous functions 4 k U ( t , x , y ) x 4 k on Ω , whereas C x , y 0 + 4 k ( Ω ) is the class of continuous functions 4 k U ( t , x , y ) y 4 k on Ω , 4 k 2 y 4 k 2 U ( t , x , l ) we understand as 4 k 2 y 4 k 2 U ( t , x , y ) y = l , Ω ¯ = { ( t , x , y ) | 0 t T , 0 x , y l } .

2. Cauchy Problem for a Fractional Ordinary Differential Equation with Degeneration

It is well-known that the two-parametric Mittag–Leffler function is defined as (see, for example, [33])
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , α > 0 , z , β C .
The generalized Mittag–Lefler (Kilbas–Saigo)-type function was defined for real α , m R and complex l C by Kilbas and Saigo in the following form [33]
E α , m , l ( z ) = k = 0 c k z k , c 0 = 1 , c k = j = 0 k 1 Γ ( α [ j m + l ] + 1 ) Γ ( α [ j m + l + 1 ] + 1 ) , k = 1 , 2 ,
These functions belong to the class of entire functions on the complex plane.
Let us consider the Cauchy problem for an ordinary differential equation of fractional order with degeneration
* D 0 + α u ( t ) = λ t β u ( t ) + f ( t ) , t ( 0 ; T ) , u ( 0 ) = u 0 ,
where β , λ , u 0 R , f ( t ) is a given continuous function.
Let γ [ 0 ; 1 ) . Then, we consider the class of following functions ([34], p. 4, 205):
C γ Ω T = g ( t ) : t γ g ( t ) C Ω T ,
C γ α Ω T = g ( t ) C Ω T : * D 0 + α g ( t ) C γ Ω T .
Lemma 1.
Let γ [ 0 ; α ] , β 0 . Then, for all f ( t ) C γ Ω T , there exists a unique solution u ( t ) C γ α Ω T of the Cauchy problem (8). This solution has the following form
u ( t ) = u 0 E α , 1 + β α , β α λ t α + β + 0 t K ( t , τ ) f ( τ ) d τ ,
where
K ( t , τ ) = i = 1 K i ( t , τ ) ,
K 0 ( t , τ ) = 1 Γ ( α ) ( t τ ) α 1 , K i ( t , τ ) = λ Γ ( α ) τ t s β ( t s ) α 1 K i 1 ( s , τ ) d s ,
where E α , 1 + β α , β α λ t α + β is the Kilbas–Saigo function, defined by (7).
Proof. 
The uniqueness of the solution u ( t ) C γ α Ω T of the problem (8) was proven in [34], p. 205. In this paper, the existence of solution for the case f ( t ) = 0 is also proved. So, we consider the inhomogeneous problem (8) and the solution to this problem as the sum of two functions
u ( t ) = v ( t ) + w ( t ) ,
where the functions v ( t ) and w ( t ) , respectively, are solutions to the following two problems:
* D 0 + α v ( t ) = λ t β v ( t ) , t ( 0 ; T ) , v ( 0 ) = u 0 ,
* D 0 + α w ( t ) = λ t β w ( t ) + f ( t ) , t ( 0 ; T ) , w ( 0 ) = 0 .
As implied by [34], p. 233, in their Theorem 4.4, the problem (13) has a unique solution v ( t ) C γ α Ω T of the form
v ( t ) = u 0 E α , 1 + β α , β α λ t α + β .
We consider the problem (14). According to [34] in their Corollary 3.24, p. 202, the problem is equivalent to the one value solvability of theVolterra-type integral equation of the second kind
w ( t ) = λ I 0 + α t β w ( t ) + I 0 + α f ( t ) , t ( 0 ; T ) .
We apply the method of successive approximations to solve the integral Equation (16). In this order, we place
w 0 ( t ) = I 0 + α f ( t ) , w m ( t ) = w 0 ( t ) + λ I 0 + α t β w m 1 ( t ) , m = 1 , 2 ,
With the convergence of the iterative process for the integral Equation (16), i.e., the existence of a solution to Equation (16), we can similarly provide the proof of the corresponding part of Theorem 3.25 in [34], p. 202. Hence, it is not difficult to determine that the solution to (16) has the form
w ( t ) = 0 t K ( t , τ ) f ( τ ) d τ ,
where the kernel K ( t , τ ) is defines by the Formulas (10) and (11). According to (12), from the representations (15) and (17), we have (9). Lemma 1 is proved. □
Lemma 2.
For the case of γ [ 0 ; α ] , β 0 , there holds the following estimate
| ( t τ ) 1 α K ( t , τ ) | E α , α | λ | t β ( t τ ) α .
Proof. 
By virtue of (10) and (11), we obtain that there the following estimates hold:
| K ( t , τ ) | = i = 1 K i ( t , τ ) | K 0 ( t , τ ) | + | K 1 ( t , τ ) | + | K 2 ( t , τ ) | + + | K i ( t , τ ) | + ,
K 0 ( t , τ ) = 1 Γ ( α ) ( t τ ) α 1 ,
K 1 ( t , τ ) = λ Γ ( α ) τ t s β ( t s ) α 1 K 0 ( s , τ ) d s | λ | Γ ( α ) τ t s β ( t s ) α 1 K 0 ( s , τ ) d s
| λ | Γ 2 ( α ) τ t s β ( t s ) α 1 ( s τ ) α 1 d s | λ | t β Γ 2 ( α ) τ t ( t s ) α 1 ( s τ ) α 1 d s .
Through the change in variables s = τ + ( t τ ) z , we obtain
τ t ( t s ) ρ 1 ( s τ ) σ 1 d s = Γ ( ρ ) Γ ( σ ) Γ ( ρ + σ ) ( t τ ) ρ + σ 1 .
Hence, considering ρ = σ = α , for (21), we derive
K 1 ( t , τ ) | λ | t β Γ 2 ( α ) Γ 2 ( α ) Γ ( 2 α ) ( t τ ) 2 α 1 = | λ | Γ ( 2 α ) t β ( t τ ) 2 α 1 .
Now, we estimate the next kernel K 2 ( t , τ ) :
K 2 ( t , τ ) = λ Γ ( α ) τ t s β ( t s ) α 1 K 1 ( s , τ ) d s | λ | Γ ( α ) τ t s β ( t s ) α 1 K 1 ( s , τ ) d s
| λ | Γ 2 ( α ) | λ | Γ ( 2 α ) τ t s 2 β ( t s ) α 1 ( s τ ) 2 α 1 d s | λ | 2 t 2 β Γ ( α ) Γ ( 2 α ) τ t ( t s ) α 1 ( s τ ) 2 α 1 d s .
Hence, considering ρ = α , σ = 2 α , we derive
K 2 ( t , τ ) | λ | 2 Γ ( 3 α ) t 2 β ( t τ ) 3 α 1 .
Through the induction method, we obtain
K m ( t , τ ) | λ | m Γ ( ( m + 1 ) α ) t m β ( t τ ) ( m + 1 ) α 1 .
Considering (6) and the estimates (20), (22)–(24), for (19), we have
K ( t , τ ) 1 Γ ( α ) ( t τ ) α 1 + | λ | Γ ( 2 α ) t β ( t τ ) 2 α 1 + | λ | 2 Γ ( 3 α ) t 2 β ( t τ ) 3 α 1
+ + | λ | m Γ ( m α + α ) t m β ( t τ ) ( m + 1 ) α 1 +
= ( t τ ) α 1 1 Γ ( α ) + | λ | Γ ( 2 α ) t β ( t τ ) α + | λ | 2 Γ ( 3 α ) t 2 β ( t τ ) 2 α
+ + | λ | m Γ ( m α + α ) t m β ( t τ ) m α + = ( t τ ) α 1 E α , α | λ | t β ( t τ ) α .
Hence, we obtain the estimate (18). Lemma 2 is proved. □

3. Expansion of the Solution into Fourier Series

Nontrivial solutions of the inverse problem are sought as a Fourier series
U ( t , x , y ) = n , m = 1 u n , m ( t ) ϑ n , m ( x , y ) ,
where
u n , m ( t ) = 0 l 0 l U ( t , x , y ) ϑ n , m ( x , y ) d x d y ,
ϑ n , m ( x , y ) = 2 l sin π n l x sin π m l y , n , m = 1 , 2 , . . .
We also suppose that the following function is expand to Fourier series
b ( x , y ) = n , m = 1 b n , m ϑ n , m ( x , y ) ,
where
b n , m = 0 l 0 l b ( x , y ) ϑ n , m ( x , y ) d x d y .
Substituting Fourier series (25) and (27) into partial differential Equation (1), we obtain the countable system of ordinary fractional differential equations of 0 < α < 1 -order with degeneration
* D 0 + α u n , m ( t ) + ω λ n , m 2 k t β u n , m ( t ) = a ( t ) b n , m 1 + ε μ n , m 4 k ,
where
λ n , m 2 k = μ n , m 4 k 1 + ε μ n , m 4 k , μ n , m k = π l k n 2 k + m 2 k .
According to Lemma 1, the general solution of the countable system of differential Equation (29) has the form
u n , m ( t ) = C n , m E α , 1 + β α , β α ω λ n , m 2 k t α + β + b n , m h n , m ( t ) ,
where
h n , m ( t ) = 1 1 + ε μ n , m 4 k 0 t K ( t , τ ) a ( τ ) d τ ,
C n , m is arbitrary constant, function K ( t , τ ) is defined by the formula (11).
By Fourier coefficients (26), the integral condition (2) is rewritten in the form
u n , m ( T ) + I 0 + ρ u n , m ( t ) t = T
= 0 l 0 l U ( T , x , y ) + I 0 + ρ U ( t , x , y ) t = T ϑ n , m ( x , y ) d x d y
= 0 l 0 l φ ( x , y ) ϑ n , m ( x , y ) d x d y = φ n , m .
To find the unknown coefficients C n , m in (30), we use condition (31), and from (30), we have
C n , m = 1 σ 0 n , m φ n , m b n , m σ 1 n , m ,
where
σ 0 n , m = E α , 1 + β / α , β / α ω λ n , m 2 k T α + β + I 0 + ρ E α , 1 + β / α , β / α λ n , m 2 k ω t α + β t = T ,
σ 1 n , m = h n , m ( T ) + I 0 + ρ h n , m ( t ) t = T .
Note that σ 0 n , m 0 , which is separate from zero and finite. This easily follows from the complete monotonicity of the function E α , m , l ( z ) , z > 0 at α 1 and from the finiteness of λ n , m 2 k [35]. Substituting the defined coefficients (32) into presentation (30), we derive that
u n , m ( t ) = φ n , m A n , m ( t ) + b n , m B n , m ( t ) ,
where
A n , m ( t ) = 1 σ 0 n , m E α , 1 + β / α , β / α ω λ n , m 2 k t α + β ,
B n , m ( t ) = h n , m ( t ) σ 1 n , m σ 0 n , m E α , 1 + β / α , β / α ω λ n , m 2 k t α + β .
Substituting the presentation of Fourier coefficients (33) of the main unknown function into Fourier series (25), we obtain
U ( t , x , y ) = n , m = 1 ϑ n , m ( x , y ) φ n , m A n , m ( t ) + b n , m B n , m ( t ) .
We consider the Fourier series (34) as a formal solution of the direct problem (1)–(4).

4. Determination of the Redefinition Function

Using the additional condition (5) and considering presentation (33), we obtain from the Fourier series (34) the following countable system for the Fourier coefficients of the redefinition function
φ n , m A n , m ( t 1 ) + b n , m B n , m ( t 1 ) = ψ n , m ,
where
ψ n , m = 0 l 0 l ψ ( x , y ) ϑ n , m ( x , y ) d x d y .
From the relation (35), we find the redefinition function as
φ n , m = ψ n , m χ 1 n , m + b n , m χ 2 n , m ,
where
χ 1 n , m = 1 A n , m ( t 1 ) , χ 2 n , m = B n , m ( t 1 ) A n , m ( t 1 ) ,
A n , m ( t 1 ) = 1 σ 0 n , m E α , 1 + β / α , β / α λ n , m 2 k ω t 1 α + β 0 , 0 < t 1 < T .
Since φ n , m are Fourier coefficients defined by (31), we substitute presentation (35) into the Fourier series
φ ( x , y ) = n , m = 1 ϑ n , m ( x , y ) ψ n , m χ 1 n , m + b n , m χ 2 n , m .
We prove the absolute and uniform convergence of the Fourier series (38) of the redefinition function. We need this to use the concepts of the following well-known Banach spaces and the Hilbert coordinate space 2 of number sequences φ n , m n , m = 1 with norm
φ 2 = n , m = 1 φ n , m 2 < .
The space L 2 ( Ω l 2 ) of square-summable functions on the domain Ω l 2 = Ω l × Ω l has norm
ϑ ( x , y ) L 2 Ω l 2 = 0 l 0 l ϑ ( x , y ) 2 d x d y < .
Smoothness conditions. For functions
ψ ( x , y ) , b ( x , y ) C 4 k ( Ω l 2 ) ,
let there exist piecewise continuous 4 k + 1 -order derivatives. Then, by integrating the functions (28) and (36) 4 k + 1 times over every variable x , y in parts, we obtain the following relations
ψ n , m = l π 8 k + 2 ψ n , m ( 8 k + 2 ) n 4 k + 1 m 4 k + 1 , b n , m = l π 8 k + 2 b n , m ( 8 k + 2 ) n 4 k + 1 m 4 k + 1 ,
ψ n , m ( 8 k + 2 ) 2 2 l 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 ,
b n , m ( 8 k + 2 ) 2 2 l 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 ,
where
ψ n , m ( 8 k + 2 ) = 0 l 0 l 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 ϑ n , m ( x , y ) d x d y ,
b n , m ( 8 k + 2 ) = 0 l 0 l 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 ϑ n , m ( x , y ) d x d y .
By obtaining estimates for the solution, we use the above properties of the Kilbas–Saigo function and Lemma 2. Then, it is easy to see that
σ 2 = max n , m χ 1 n , m ; χ 2 n , m   < ,
where
χ 1 n , m = 1 A n , m ( t 1 ) , χ 2 n , m = B n , m ( t 1 ) A n , m ( t 1 ) , 0 < t 1 < T .
Theorem 1.
Suppose that the conditions of smoothness are fulfilled. Then, Fourier series (38) is absolute and uniform convergent.
Proof. 
We use Formulas (39)–(41) and Estimate (42). Using the Cauchy–Schwartz inequality for series (38), we obtain the estimate
φ ( x , y ) n , m = 1 ϑ n , m ( x , y ) · ψ n , m χ 1 n , m + b n , m χ 2 n , m
2 l σ 2 n , m = 1 ψ n , m + n , m = 1 b n , m
2 l l π 8 k + 2 σ 2 n , m = 1 ψ n , m ( 8 k + 2 ) n 4 k + 1 m 4 k + 1 + n , m = 1 b n , m ( 8 k + 2 ) n 4 k + 1 m 4 k + 1
2 l l π 8 k + 2 σ 2 C 01 ψ n , m ( 8 k + 2 ) 2 + b n , m ( 8 k + 2 ) 2
γ 1 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 + 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 < ,
where
γ 1 = σ 2 C 0 1 2 l 2 l π 8 k + 2 , C 0 1 = n , m = 1 1 n 8 k + 2 m 8 k + 2 < .
Estimate (43) implies the absolute and uniform convergence of Fourier series (38). Theorem 1 is proved. □

5. Determination of the Main Unknown Function

So, the redefinition function is determined as a Fourier series (38). Now, the redefinition function is known. Substituting representation (37) into the Fourier series (34), the main unknown function can be presented as
U ( t , x , y ) = n , m = 1 ϑ n , m ( x , y ) ψ n , m P n , m ( t ) + b n , m Q n , m ( t ) ,
where
P n , m ( t ) = χ 1 n , m A n , m ( t ) , Q n , m ( t ) = χ 2 n , m A n , m ( t ) + B n , m ( t ) .
To establish the uniqueness of the function U ( t , x , y ) , we suppose that there are two functions, U 1 and U 2 , satisfying the given conditions (1)–(5). Then, their difference U = U 1 U 2 is a solution of the differential Equation (1), satisfying conditions (2)–(5) with function ψ ( x , y ) 0 . By virtue of relations (36), we have ψ n , m = 0 . Hence, we obtain from formulas (26), (28), and (44) in the domain Ω that the zero identity follows
0 l 0 l U ( t , x , y ) ϑ n , m ( x , y ) d x d y 0 .
Hence, by virtue of the completeness of the systems of eigenfunctions 2 l sin π n l x , 2 l sin π m l y in L 2 Ω l 2 , we deduce that U ( t , x , y ) 0 for all x Ω l 2 [ 0 , l ] 2 and t Ω T [ 0 ; T ] .
Since U ( t , x , y ) C Ω ¯ , then U t , x , y 0 in the domain Ω ¯ . Therefore, the solution of the problem (1)–(5) is unique in the domain Ω ¯ .
Theorem 2.
Let the conditions of the Theorem 1 be fulfilled. Then, the series (44) converges. At the same time, their term-by-term differentiation is possible.
Proof. 
By virtue of conditions of Theorem 1 and the properties of the Mittag–Leffler and Kilbas–Saigo functions, the functions P n , m ( t ) , Q n , m ( t ) are uniformly bounded on the segment [ 0 ; T ] . So, for any positive integers n , m , there exists a finite constant σ 3 , so the following estimate is true:
max n , m max 0 t T P n , m ( t ) ; max 0 t T Q n , m ( t ) σ 3 .
Using estimation Formulas (39)–(41) and (45), analogously to the case of estimate (43), for series (44), we obtain
U ( t , x , y ) n , m = 1 ϑ n , m ( x , y ) · ψ n , m P n , m ( t ) + b n , m Q n , m ( t )
γ 2 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 + 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 < ,
where γ 2 = C 01 σ 3 2 l 2 l π 8 k + 2 .
The estimate (46) implies the absolute and uniform convergence of the Fourier series (44). We differentiate this function (44) the required number of times
4 k x 4 k U ( t , x , y ) = n , m = 1 π n l 4 k ϑ n , m ( x , y ) ψ n , m P n , m ( t ) + b n , m Q n , m ( t ) ,
4 k y 4 k U ( t , x , y ) = n , m = 1 π m l 4 k ϑ n , m ( x , y ) ψ n , m P n , m ( t ) + b n , m Q n , m ( t ) .
The expansions of the following functions into Fourier series are defined similarly
* D 0 + α U ( t , x , y ) , 4 k x 4 k * D 0 + α U ( t , x , y ) , 4 k y 4 k * D 0 + α U ( t , x , y ) .
We show the convergence of series (47) and (48). As in the case of estimate (43), applying the Cauchy–Schwarz inequality, we obtain
4 k x 4 k U ( t , x , y ) n , m = 1 π n l 4 k u n , m ( t ) · ϑ n , m ( x , y )
2 l π l 4 k σ 3 n , m = 1 n 4 k ψ n , m + n , m = 1 n 4 k b n , m
2 l l π 4 k + 2 σ 3 n , m = 1 ψ n , m ( 8 k + 2 ) n m 4 k + 1 + n , m = 1 b n , m ( 8 k + 2 ) n m 4 k + 1
2 l l π 4 k + 2 σ 3 C 02 ψ n , m ( 8 k + 2 ) 2 + b n , m ( 8 k + 2 ) 2
γ 3 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 + 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 < ,
where γ 3 = 2 l 2 l π 4 k + 2 σ 3 C 02 , C 0 2 = n , m = 1 1 n m 8 k + 2 < ;
4 k y 4 k U ( t , x , y ) n , m = 1 π m l 4 k u n , m ( t ) · ϑ n , m ( x , y )
2 l π l 4 k σ 3 n , m = 1 m 4 k ψ n , m + n , m = 1 m 4 k b n , m
2 l l π 4 k + 2 σ 3 n , m = 1 ψ n , m ( 8 k + 2 ) n 4 k + 1 m + n , m = 1 b n , m ( 8 k + 2 ) n 4 k + 1 m
2 l l π 4 k + 2 σ 3 C 03 ψ n , m ( 8 k + 2 ) 2 + b n , m ( 8 k + 2 ) 2
γ 4 8 k + 2 ψ ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 + 8 k + 2 b ( x , y ) x 4 k + 1 y 4 k + 1 L 2 Ω l 2 < ,
where γ 4 = 2 l 2 l π 4 k + 2 σ 3 C 03 , C 0 3 = n , m = 1 1 n 8 k + 2 m < .
The convergence of the Fourier series for functions
* D 0 + α U ( t , x , y ) , 4 k x 4 k * D 0 + α U ( t , x , y ) , 4 k y 4 k * D 0 + α U ( t , x , y )
is easy to prove because the necessary estimates are obtained similarly for the cases of estimates (43), (49) and (50). Therefore, the function U ( t , x , y ) belongs to the (4) class of functions. Theorem 2 is proved. □

6. Stability of the Solution U (t, x, y) with Respect to the Given Functions

Theorem 3.
Suppose that all the conditions of Theorem 2 are fulfilled. Then, the function U ( t , x , y ) as a solution of the problem (1)–(5) is stable with respect to given function ψ ( x , y ) .
Proof. 
We show that the solution of the differential Equation (1) U ( t , x , y ) is stable with respect to a given function ψ ( x , y ) . Let U 1 ( t , x , y ) and U 2 ( t , x , y ) be two different solutions of the inverse boundary value problem (1)–(5), corresponding to two different values of the function ψ 1 ( x , y ) and ψ 2 ( x , y ) , respectively.
We state that ψ 1 n , m ψ 2 n , m < δ n , m , where 0 < δ n , m is a sufficiently small positive quantity and the series n , m = 1 δ n , m is convergent. Then, as such, given the conditions of the theorem, from the Fourier series (44), we obtain
U 1 ( t , x , y ) U 2 ( t , x , y ) C ( Ω ¯ ) 2 l σ 3 n , m = 1 ψ 1 n , m ψ 2 n , m < 2 l σ 3 n , m = 1 δ n , m < .
If we set ε = 2 l σ 3 n , m = 1 δ n , m < , then, from the last estimate, we derive assertions about the stability of the solution of the differential Equation (1) with respect to a given function ψ ( x , y ) in (5). Theorem 3 is proved. □
Similarly, we proved that the following two theorems hold.
Theorem 4.
Suppose that all the conditions of Theorem 2 are fulfilled. Then, the function U ( t , x , y ) as a solution of the problem (1)–(5) is stable with respect to a given function b ( x , y ) in the right-hand side of Equation (1).
Theorem 5.
Suppose that all the conditions of Theorem 2 are fulfilled. Then, the function U ( t , x , y ) , as a solution of the problem (1)–(5), is stable with respect to the source function φ ( x , y ) .

7. Conclusions

In three-dimensional domain, an inverse boundary value problem (1)–(5) of the identification of the redefinition function φ ( x , y ) for a partial differential equation with degeneration and integral form condition was studied in this paper. The case of a 0 < α 1 -order Gerasimov–Caputo-type fractional operator was considered. The solution of the partial differential equation was studied in the class of regular functions. Equation (1) depends on three independent arguments: t , x , y . The first argument is the time argument, and with respect to this argument, Equation (1) is a fractional Gerasimov–Caputo-type ordinary differential equation with degeneration. The other two variables x , y are spatial, and Equation (1) with respect to them is a partial differential equation of higher even order. The Fourier series method was used and a countable system of ordinary differential equations was obtained. Using the given additional condition (5), the presentation for the redefinition function was obtained. When conditions of smoothness are fulfilled, then using the Cauchy–Schwarz inequality and the Bessel inequality, the absolute and uniform convergence of the obtained Fourier series was proved. The stability of main unknown function u ( t , x , y ) of the problem (1)–(5) with respect to given functions b ( x , y ) , ψ ( x , y ) and redefinition function φ ( x , y ) was studied.
This work is theoretical in nature and develops the theory of differential equations with fractional operators. We studied the unique solvability in the classical sense of a nonlocal inverse problem for a partial differential equation. In the future, we intend to continue our research in the direction of superposition of several fractional-order operators.
Remark 1.
Equation (1) is a generalization of the Barenblatt–Zheltov–Kochina differential equation. The Barenblatt–Zheltov–Kochina equation simulates the filtration of a viscoelastic fluid in fractured porous media [36].

Author Contributions

Conceptualization, T.K.Y. and B.J.K. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ashurov, R.R.; Mukhiddinova, A.T. Inverse problem of determining the heat source density for the subdiffusion equation. Differ. Equ. 2020, 56, 1550–1563. [Google Scholar] [CrossRef]
  2. Ashyralyyev, C.; Cay, A. Numerical solution to elliptic inverse problem with Neumann-type integral condition and overdetermination. Bull. Karaganda Univ. Math. Ser. 2020, 99, 5–17. [Google Scholar] [CrossRef]
  3. Kostin, A.B. The inverse problem of recovering the source in a parabolic equation under a condition of nonlocal observation. Sbornik. Math. 2013, 204, 1391–1434. [Google Scholar] [CrossRef]
  4. Mamedov, K.R. Uniqueness of the solution to the inverse problem of scattering theory for the Sturm-Liouville operator with a spectral parameter in the boundary condition. Math. Notes 2003, 74, 136–140. [Google Scholar] [CrossRef]
  5. Yuldashev, T.K. Nonlocal inverse problem for a pseudohyperbolic-pseudoelliptic type integro-differential equations. Axioms 2020, 9, 45. [Google Scholar] [CrossRef]
  6. Assanova, A.T. An integral-boundary value problem for a partial differential equation of second order. Turk. J. Math. 2019, 43, 1967–1978. [Google Scholar] [CrossRef]
  7. Assanova, A.T.; Imanchiyev, A.E.; Kadirbayeva, Z.M. A nonlocal problem for loaded partial differential equations of fourth order. Bull. Karaganda Univ. Math. 2020, 97, 6–16. [Google Scholar] [CrossRef]
  8. Il’in, V.A. Uniqueness of generalized solutions of mixed problems for the wave equation with nonlocal boundary conditions. Differ. Equ. 2008, 44, 692–700. [Google Scholar] [CrossRef]
  9. Yuldashev, T.K. Nonlocal mixed-value problem for a Boussinesq-type integro-differential equation with degenerate kernel. Ukr. Math. J. 2016, 68, 1278–1296. [Google Scholar] [CrossRef]
  10. Yuldashev, T.K. Mixed problem for pseudoparabolic integr-odifferential equation with degenerate kernel. Differ. Equ. 2017, 53, 99–108. [Google Scholar] [CrossRef]
  11. Yuldashev, T.K. Solvability of a boundary value problem for a differential equation of the Boussinesq type. Differ. Equ. 2018, 54, 1384–1393. [Google Scholar] [CrossRef]
  12. Yuldashev, T.K. On a boundary-value problem for Boussinesq type nonlinear integro-differential equation with reflecting argument. Lobachevskii J. Math. 2020, 41, 111–123. [Google Scholar] [CrossRef]
  13. Yuldashev, T.K.; Kadirkulov, B.J. Boundary value problem for weak nonlinear partial differential equations of mixed type with fractional Hilfer operator. Axioms 2020, 9, 68. [Google Scholar] [CrossRef]
  14. Yuldashev, T.K.; Kadirkulov, B.J. Nonlocal problem for a mixed type fourth-order differential equation with Hilfer fractional operator. Ural Math. J. 2020, 6, 153–167. [Google Scholar] [CrossRef]
  15. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  16. Gerasimov, A.N. Generalization of laws of the linear deformation and their application to problems of the internal friction. Prikl. Mat. Mekh. 1948, 12, 251–260. [Google Scholar]
  17. Novozhenova, O.G. Life and science of Alexey Gerasimov, one of the pioneers of fractional calculus in Soviet Union. Fract. Calc. Appl. Anal. 2017, 20, 790–809. [Google Scholar] [CrossRef]
  18. Rossikhin, Y.A. Reflections on two parallel ways in the progress of fractional calculus in mechanics of solids. Appl. Mech. Rev. 2010, 63, 010701. [Google Scholar] [CrossRef]
  19. Mainardi, F. Fractional calculus: Some basic problems in continuum and statistical mechanics. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: Wien, Austria, 1997. [Google Scholar]
  20. Handbook of Fractional Calculus with Applications; Tenreiro Machado, J.A. (Ed.) Walter de Gruyter GmbH: Berlin, Germany, 2019; Volume 8, pp. 47–85. [Google Scholar]
  21. Kumar, D.; Baleanu, D. Fractional Calculus and Its Applications in Physics. Front. Phys. 2019, 7. [Google Scholar] [CrossRef]
  22. Sun, H.; Chang, A.; Zhang, Y.; Chen, W. A review on variable-order fractional differential equations: Mathematical foundations, physical models, numerical methods and applications. Fract. Calc. Appl. Anal. 2019, 22, 27–59. [Google Scholar] [CrossRef] [Green Version]
  23. Hilfer, R. Application of Fractional Calculus in Physics; World Scientific Publishing Company: Singapore, 2000. [Google Scholar]
  24. Patnaik, S.; Hollkamp, J.P.; Semperlotti, F. Applications of variable-order fractional operators: A review. Proc. R. Soc. 2020, 476, 1–32. [Google Scholar] [CrossRef] [Green Version]
  25. Yuldashev, T.K.; Karimov, E.T. Inverse problem for a mixed type integro-differential equation with fractional order Caputo operators and spectral parameters. Axioms 2020, 9, 121. [Google Scholar] [CrossRef]
  26. 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, 42, 2043–2052. [Google Scholar] [CrossRef]
  27. Turmetov, B.K.; Kadirkulov, B.J. On a problem for nonlocal mixed-type fractional order equation with degeneration. Chaos Solitons Fractals 2021, 146, 1–5. [Google Scholar] [CrossRef]
  28. Berdyshev, A.S.; Kadirkulov, J.B. On a nonlocal problem for a fourth-order parabolic equation with the fractional Dzhrbashyan–Nersesyan operator. Differ. Equ. 2016, 52, 122–127. [Google Scholar] [CrossRef]
  29. Malik, S.A.; Aziz, S. An inverse source problem for a two parameter anomalous diffusion equation with nonlocal boundary conditions. Comput. Math. Appl. 2017, 73, 2548–2560. [Google Scholar] [CrossRef]
  30. Serikbaev, D.; Tokmagambetov, N. A source inverse problem for the pseudo-parabolic equation with the fractional Sturm-Liouville operator. Bull. Karaganda Univ. Math. Ser. 2020, 99, 143–151. [Google Scholar] [CrossRef]
  31. Yuldashev, T.K.; Abdullaev, O.K. Unique solvability of a boundary value problem for a loaded fractional parabolic-hyperbolic equation with nonlinear terms. Lobachevskii J. Math. 2021, 42, 1113–1123. [Google Scholar] [CrossRef]
  32. Yuldashev, T.K.; Kadirkulov, B.J. Inverse boundary value problem for a fractional differential equations of mixed type with integral redefinition conditions. Lobachevskii J. Math. 2021, 42, 649–662. [Google Scholar] [CrossRef]
  33. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag–Leffler Functions, Related Topics and Applications; Springer Monographs in Mathematics; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar] [CrossRef]
  34. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and applications of fractional differential equations. In North-Holland Mathematical Studies; Elsevier Science B. V.: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  35. Boudabsa, L.; Simon, T.; Vallois, P. Fractional extreme distributions. Math. Probab. 2021, 1–46. [Google Scholar]
  36. Barenblatt, G.I.; Zheltov, Y.P.; Kochina, I.N. On the basic concepts of the theory of filtration in fractured media. Appl. Math. Mech. 1960, 24, 58–73. (In Russian) [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuldashev, T.K.; Kadirkulov, B.J. Inverse Problem for a Partial Differential Equation with Gerasimov–Caputo-Type Operator and Degeneration. Fractal Fract. 2021, 5, 58. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020058

AMA Style

Yuldashev TK, Kadirkulov BJ. Inverse Problem for a Partial Differential Equation with Gerasimov–Caputo-Type Operator and Degeneration. Fractal and Fractional. 2021; 5(2):58. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020058

Chicago/Turabian Style

Yuldashev, Tursun K., and Bakhtiyar J. Kadirkulov. 2021. "Inverse Problem for a Partial Differential Equation with Gerasimov–Caputo-Type Operator and Degeneration" Fractal and Fractional 5, no. 2: 58. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020058

Article Metrics

Back to TopTop