Next Article in Journal / Special Issue
A Rosenzweig–MacArthur Model with Continuous Threshold Harvesting in Predator Involving Fractional Derivatives with Power Law and Mittag–Leffler Kernel
Previous Article in Journal
Local Spectral Theory for R and S Satisfying RnSRn = Rj
Previous Article in Special Issue
Finite Series of Distributional Solutions for Certain Linear Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inverse Problem for a Mixed Type Integro-Differential Equation with Fractional Order Caputo Operators and Spectral Parameters

by
Tursun K. Yuldashev
1,* and
Erkinjon T. Karimov
2
1
Uzbek-Israel Joint Faculty of High Technology and Engineering Mathematics, National University of Uzbekistan, Tashkent 100174, Uzbekistan
2
V. I. Romanovskiy Institute of Mathematics, Uzbekistan Academy of Sciences, Tashkent 100174, Uzbekistan
*
Author to whom correspondence should be addressed.
Submission received: 25 September 2020 / Revised: 15 October 2020 / Accepted: 16 October 2020 / Published: 20 October 2020

Abstract

:
The questions of the one-value solvability of an inverse boundary value problem for a mixed type integro-differential equation with Caputo operators of different fractional orders and spectral parameters are considered. The mixed type integro-differential equation with respect to the main unknown function is an inhomogeneous partial integro-differential equation of fractional order in both positive and negative parts of the multidimensional rectangular domain under consideration. This mixed type of equation, with respect to redefinition functions, is a nonlinear Fredholm type integral equation. The fractional Caputo operators’ orders are smaller in the positive part of the domain than the orders of Caputo operators in the negative part of the domain under consideration. Using the method of Fourier series, two systems of countable systems of ordinary fractional integro-differential equations with degenerate kernels and different orders of integro-differentation are obtained. Furthermore, a method of degenerate kernels is used. In order to determine arbitrary integration constants, a linear system of functional algebraic equations is obtained. From the solvability condition of this system are calculated the regular and irregular values of the spectral parameters. The solution of the inverse problem under consideration is obtained in the form of Fourier series. The unique solvability of the problem for regular values of spectral parameters is proved. During the proof of the convergence of the Fourier series, certain properties of the Mittag–Leffler function of two variables, the Cauchy–Schwarz inequality and Bessel inequality, are used. We also studied the continuous dependence of the solution of the problem on small parameters for regular values of spectral parameters. The existence and uniqueness of redefined functions have been justified by solving the systems of two countable systems of nonlinear integral equations. The results are formulated as a theorem.

1. Introduction

Fractional calculus plays an important role in the mathematical modeling of many natural and engineering processes (see [1]). We can gladly refer to many examples of applied research works, where fractional integro-differential operators are successfully and widely used. For example, in [2] some applications of basic problems in continuum and statistical mechanics are considered. In [3], the mathematical problems of an Ebola epidemic model by fractional order equations are studied. In [4,5], fractional models of the dynamics of tuberculosis infection and novel coronavirus (nCoV-2019) are studied, respectively. The construction of various models for studying problems of theoretical physics by the aid of fractional calculus is described in [6] (vol. 4, 5), [7,8]. A specific physical interpretation of the fractional derivatives, describing the random motion of a particle moving on the real line at Poisson-paced times with finite velocity, is given in [9]. A detailed review of the applications of fractional calculus in solving practical problems is given in [6] (vol. 6–8), [10]. More detailed information, as well as a bibliography related to the theory of fractional integro-differentiation and fractional derivatives, can also be found in [11,12,13,14,15,16,17,18].
We also note the special role of generalized special functions, such as polynomials, in solving fractional differential equations. In [19], using Hermite polynomials of higher and fractional order, some operational techniques to find general solutions of extended forms to d’Alembert and Fourier equations. In [20], the solutions of various generalized forms of the Heat Equation, by means of different tools ranging from the use of Hermite–Kampé de Fériet polynomials of higher and fractional order to operational techniques, are discussed. In [21], the combined use of integral transforms and special polynomials provides a powerful tool to deal with fractional derivatives and integrals. The real need to know the properties of such special functions in solving direct and inverse problems for fractional partial differential equations has been shown in [22].
Applications for equations of mixed type are studied in the works of many researchers. For example, in [23], an example of gas motion in a channel surrounded by a porous medium was studied, with the gas motion in a channel being described by a wave equation, while—outside the channel—a diffusion equation was posed. In [24], a problem related to the propagation of electric oscillations in compound lines, when the losses on a semi-infinite line were neglected and the rest of the line was treated as a cable with no leaks, was investigated. This reduced the problem under consideration to a mixed parabolic–hyperbolic type equation. In [25], a hyperbolic–parabolic system, in relation to pulse combustion, is investigated. Mixed type fractional differential equations are studied in many works by scientists—particularly in [26,27,28,29,30,31,32,33,34,35].
The theories of integral and integro-differential equations are important in studying the large directions of the general theory of equations of mathematical physics. The presence of an integral term in differential equations of the first and second order has an important role in the theory of dynamical systems of automatic control [36,37]. Boundary value problems for integro-differential equations with spectral parameters have singularities in studying the questions of one-value solvability [38,39]. Mixed type integer order integro-differential equations with degenerate kernels and spectral parameters are studied in [40,41].
To find the solutions of direct mixed and boundary value problems of mathematical physics, it is required to set the coefficients of the equation, the boundary of the domain under consideration, and the initial and boundary data. It usually happens that, in solving the applied problems experimentally, the quantitative characteristics of the object under study are not available for direct observation, or it is impossible to carry out the experiment itself for one reason or another. Then, in practice, the researchers can obtain some indirect information and draw a conclusion about the properties of the studied object. This information is determined by the nature of the object under study and here requires mathematical processing and the interpretation of research results. Nonlocal integral conditions often arise when the experiment gives averaged information about this object. When the structure of the mathematical model of the studying process is known, the problem of redefining the mathematical model is posed. Such problems belong to the class of inverse problems. By inverse problems we mean problems whose solution consists of determining the parameters of a model based on the available observation results and other experimental information. Inverse problems for equations of mixed type are studied relatively rarely due to the complexity of the studying process.
In the present paper, we study the questions of the one-value solvability of an inverse boundary value problem for a mixed type integro-differential equation with Caputo operators of different fractional orders and spectral parameters in a multidimensional rectangular domain.
The rest of this paper is organized as follows. In Section 2, we state the problem, which we will investigate in this work. Section 3 is devoted to formally expanding the solution of the direct problem into Fourier series. In Section 4, we formally determine the redefinition functions. Section 5 contains the proof of existence and uniqueness of Fourier coefficients of redefinition functions from a countable system of nonlinear integral equations. Section 6 is devoted to the justification of convergence and the possibility of the term by term differentiation of the obtained Fourier series. Section 7 contains the proof of the continuous dependence of the solution on the small parameter. In the last Section 8, as a conclusion, we formulate the theorem, which we have proved in this paper.

2. Statement of the Problem

We recall that the Caputo differential operator of fractional order m 1 < α < m has the form
C D a t α f ( t ) = 1 Γ ( m α ) a t ( t s ) m α 1 f ( m ) ( s ) d s ,
where Γ ( z ) is Euler gamma function.
In the multidimensional domain Ω = { T < t < T , 0 < x 1 , , x m < l } , a mixed type integro-differential equation of the following form is considered:
A ε ( U ) B ω ( U ) = ν 0 T K 1 ( t , s ) U ( s , x ) d s + F 1 ( t , x ) , t > 0 , ν T 0 K 2 ( t , s ) U ( s , x ) d s + F 2 ( t , x ) , t < 0 ,
where
F i ( t , x ) = k i ( t ) g i ( x ) + f i x , Ω l m Θ i ( y ) g i ( y ) d y , i = 1 , 2 ,
A ε ( U ) = 1 + sgn ( t ) 2 C D 0 t α 1 ε i = 1 m 2 x i x i C D 0 t β 1 U ( t , x )
+ 1 sgn ( t ) 2 C D 0 t α 2 ε i = 1 m 2 x i x i C D 0 t β 2 U ( t , x ) ,
B ω ( U ) = i = 1 m U x i x i , t > 0 , ω 2 i = 1 m U x i x i , t < 0 ,
T and l are given as positive real numbers, ω is a positive spectral parameter, ε is a positive small parameter, ν is a real non-zero spectral parameter, 0 K j ( t , s ) = a j ( t ) b j ( s ) , a j ( t ) C 2 [ T ; T ] , b j ( s ) C [ T ; T ] , f i C x 2 Ω l m × R , Ω l m Θ i ( y ) d y < , Ω l m Θ i ( y ) d y = 0 l 0 l Θ i ( y ) d y 1 · · d y m , i , j = 1 , 2 , k 1 ( t ) C 2 [ 0 ; T ] , k 2 ( t ) C 2 [ T ; 0 ] , while g 1 ( x ) and g 2 ( x ) are redefinition functions, R ( ; ) , x Ω l m [ 0 ; l ] m , 0 < β 1 < α 1 1 , 1 < β 2 < α 2 2 .
Problem. Find in the domain Ω a triple of unknown function
U ( t , x ) C ( Ω ¯ ) C 0 , 1 Ω C α 1 , 2 ( Ω + ) C α 2 , 2 ( Ω ) C t , x α 1 + 2 ( Ω + ) C t , x α 2 + 2 ( Ω ) C t , x 1 , x 2 , , x m α 1 + 2 + 0 + + 0 ( Ω + ) C t , x 1 , x 2 , , x m α 2 + 2 + 0 + + 0 ( Ω ) C t , x 1 , x 2 , x 3 , , x m α 1 + 0 + 2 + 0 + + 0 ( Ω + ) C t , x 1 , x 2 , x 3 , , x m α 2 + 0 + 2 + 0 + + 0 ( Ω ) C t , x 1 , , x m 1 , x m α 1 + 0 + + 0 + 2 ( Ω + ) C t , x 1 , , x m 1 , x m α 2 + 0 + + 0 + 2 ( Ω )
and redefinition functions g i ( x ) C Ω l m , i = 1 , 2 , satisfying the mixed integro-differential Equation (1) and the following boundary conditions
U ( T , x ) = φ 1 ( x ) , C D 0 t θ U ( T , x ) = φ 2 ( x ) , x Ω l m ,
U ( t , 0 ) = U ( t , l ) = 0 , T < t < T
and additional conditions
0 T Φ 1 ( t ) U ( t , x ) = ψ 1 ( x ) , x Ω l m ,
T 0 Φ 2 ( t ) U ( t , x ) = ψ 2 ( x ) , x Ω l m ,
where 0 < θ < 1 , φ i ( x ) , ψ i ( x ) are given smooth functions, φ i ( 0 ) = φ i ( l ) = 0 , ψ i ( 0 ) = ψ i ( l ) = 0 , i = 1 , 2 , C r ( Ω ) is a class of functions U ( t , x 1 , , x m ) with continuous derivatives r U t r , r U x 1 r , , r U x m r in Ω , C t , x r , s ( Ω ) is a class of functions U ( t , x 1 , , x m ) with continuous derivatives r U t r , s U x 1 s , , s U x m s in Ω , C t , x 1 , x 2 , , x m r + r + 0 + + 0 ( Ω ) is a class of functions U ( t , x 1 , , x m ) with continuous derivative 2 r U t r x 1 r in Ω , , C t , x 1 , , x m 1 , x m r + 0 + + 0 + r ( Ω ) is a class of functions U ( t , x 1 , , x m ) with continuous derivative 2 r U t r x m r in Ω , r , s are positive real numbers, Ω ¯ = T t T , x Ω l m , Ω = Ω { x 1 , , x m = 0 } { x 1 , , x m = l } , Ω = { T < t < 0 , 0 < x 1 , , x m < l } , Ω + = { 0 < t < T , 0 < x 1 , , x m < l } .

3. Expansion of the Solution of the Direct Problem (1)–(4) into Fourier Series

Our investigation is based on the application of sine Fourier series to the mixed type integro-differential Equation (1) of the complicated form. Hence, the solution of the mixed integro-differential Equation (1) in domain Ω is sought in the form of the following Fourier series
U ( t , x ) = n 1 , , n m = 1 u n 1 , , n m ± ( t ) ϑ n 1 , , n m ( x ) ,
where
u n 1 , , n m ± ( t ) = u n 1 , , n m + ( t ) = Ω l m U ( t , x ) ϑ n 1 , , n m ( x ) d x , t > 0 , u n 1 , , n m ( t ) = Ω l m U ( t , x ) ϑ n 1 , , n m ( x ) d x , t < 0 ,
Ω l m U ( t , x ) ϑ n 1 , , n m ( x ) d x = 0 l 0 l U ( t , x ) ϑ n 1 , , n m ( x ) d x 1 · · d x m ,
ϑ n 1 , , n m ( x ) = 2 l m sin π n 1 l x 1 · · sin π n m l x m , n 1 , , n m = 1 , 2 ,
In this order, we also suppose that the redefinition functions and nonlinear functions on the right-hand side of the integro-differential Equation (1) are representable as the following Fourier series
g i ( x ) = n 1 , , n m = 1 g i n 1 , , n m ϑ n 1 , , n m ( x ) , f i ( x , V i ) = n 1 , , n m = 1 f i n 1 , , n m ( V i ) ϑ n 1 , , n m ( x ) ,
where
g i n 1 , , n m = Ω l m g i ( x ) ϑ n 1 , , n m ( x ) d x , f i n 1 , , n m ( V i ) = Ω l m f i ( y , V i ) ϑ n 1 , , n m ( y ) d y ,
f i ( y , V i ) = f i y , Ω l m Θ i ( z ) g i ( z ) d z , i = 1 , 2 .
Substituting series (7) and (9) into mixed Equation (1), we obtain two fractional countable systems of ordinary integro-differential equations
C D 0 t α 1 u n 1 , , n m + ( t ) + ε μ n 1 , , n m 2 C D 0 t β 1 u n 1 , , n m + ( t ) + μ n 1 , , n m 2 u n 1 , , n m + ( t ) = ν 0 T a 1 ( t ) b 1 ( s ) u n 1 , , n m + ( s ) d s + F 1 n 1 , , n m ( t ) , t > 0 ,
C D 0 t α 2 u n 1 , , n m ( t ) + ε μ n 1 , , n m 2 C D 0 t β 2 u n 1 , , n m ( t ) + μ n 1 , , n m 2 ω 2 u n 1 , , n m ( t ) = ν T 0 a 2 ( t ) b 2 ( s ) u n 1 , , n m ( s ) d s + F 2 n 1 , , n m ( t ) , t < 0 ,
where μ n 1 , , n m = π l n 1 2 + + n m 2 ,
F i n 1 , , n m ( t ) = k i ( t ) g i n 1 , , n m + f i n 1 , , n m ( V i ) , i = 1 , 2 .
We use the method of degenerate kernels. In this order, by the aid of designations
τ n 1 , , n m + = 0 T b 1 ( s ) u n 1 , , n m + ( s ) d s ,
τ n 1 , , n m = T 0 b 2 ( s ) u n 1 , , n m ( s ) d s
we present the countable systems of ordinary integro-differential Equations (10) and (11) as follows
C D 0 t α 1 u n 1 , , n m + ( t ) + ε μ n 1 , , n m 2 C D 0 t β 1 u n 1 , , n m + ( t ) + μ n 1 , , n m 2 u n 1 , , n m + ( t ) = ν a 1 ( t ) τ n 1 , , n m + + F 1 n 1 , , n m ( t ) , t > 0 ,
C D 0 t α 2 u n 1 , , n m ( t ) + ε μ n 1 , , n m 2 C D 0 t β 2 u n 1 , , n m ( t ) + μ n 1 , , n m 2 ω 2 u n 1 , , n m ( t ) = ν a 2 ( t ) τ n 1 , , n m + F 2 n 1 , , n m ( t ) , t < 0 .
The solutions of the countable systems of differential Equations (15) and (16), satisfying conditions
u n 1 , , n m + ( 0 ) = C 1 n 1 , , n m + , u n 1 , , n m ( 0 ) = C 1 n 1 , , n m , d d t u n 1 , , n m ( 0 ) = C 2 n 1 , , n m
have the following form:
u n 1 , , n m + ( t ) = ν τ n 1 , , n m + Ψ 11 n 1 , , n m ( t , ε ) + Ψ 12 n 1 , , n m ( t , ε ) + C 1 n 1 , , n m + Ψ 13 n 1 , , n m ( t , ε ) , t > 0 ,
u n 1 , , n m ( t ) = ν τ n 1 , , n m Ψ 21 n 1 , , n m ( t , ε , ω ) + Ψ 22 n 1 , , n m ( t , ε , ω ) + C 1 n 1 , , n m Ψ 23 n 1 , , n m ( t , ε , ω ) C 2 n 1 , , n m Ψ 24 n 1 , , n m ( t , ε , ω ) , t < 0 ,
where C 1 n 1 , , n m + , C i n 1 , , n m , ( i = 1 , 2 ) are for unknown constants to be uniquely determined,
Ψ 11 n 1 , , n m ( t , ε ) = 0 t a 1 ( t s ) s α 1 1 E ( α 1 β 1 , α 1 ) , α 1 ε μ n 1 , , n m 2 s α 1 β 1 , μ n 1 , , n m 2 s α 1 d s ,
Ψ 12 n 1 , , n m ( t , ε ) = 0 t F 1 n 1 , , n m ( t s ) s α 1 1 E ( α 1 β 1 , α 1 ) , α 1 ε μ n 1 , , n m 2 s α 1 β 1 , μ n 1 , , n m 2 s α 1 d s ,
Ψ 13 n 1 , , n m ( t , ε ) = E ( α 1 β 1 , α 1 ) , 1 ε μ n 1 , , n m 2 t α 1 β 1 , μ n 1 , , n m 2 t α 1 ,
Ψ 21 n 1 , , n m ( t , ε , ω ) = t 0 a 2 ( s t ) ( s ) α 2 1 Ψ 25 n 1 , , n m ( t , ε , ω ) d s ,
Ψ 22 n 1 , , n m ( t , ε , ω ) = t 0 F 2 n 1 , , n m ( s t ) ( s ) α 2 1 Ψ 25 n 1 , , n m ( t , ε , ω ) d s ,
Ψ 23 n 1 , , n m ( t , ε , ω ) = E ( α 2 β 2 , α 2 ) , 1 ε μ n 1 , , n m 2 ( t ) α 2 β 2 , μ n 1 , , n m 2 ω 2 ( t ) α 2 ,
Ψ 24 n 1 , , n m ( t , ε , ω ) = t E ( α 2 β 2 , α 2 ) , 2 ε μ n 1 , , n m 2 ( t ) α 2 β 2 , μ n 1 , , n m 2 ω 2 ( t ) α 2 ,
Ψ 25 n 1 , , n m ( t , ε , ω ) = E ( α 2 β 2 , α 2 ) , α 2 ε μ n 1 , , n m 2 ( s ) α 2 β 2 , μ n 1 , , n m 2 ω 2 ( s ) α 2 ,
The function E ( α , β ) , γ ( z 1 , z 2 ) is a Mittag–Leffler function of two variables:
E ( α , β ) , γ ( z 1 , z 2 ) = m 1 , m 2 = 0 z 1 m 1 z 2 m 2 Γ ( γ + α m 1 + β m 2 ) ,
where z i , α , β , γ C , Re ( α ) > 0 , Re ( β ) > 0 .
From the statement of the problem (properties in (2)), it follows that the continuous conjugation condition is fulfilled for the main unknown function: U ( 0 + 0 , x ) = U ( 0 0 , x ) . Therefore, by taking Formula (6) into account, we have the conditions for Fourier coefficients of the main unknown function
u n 1 , , n m + ( 0 + 0 ) = Ω l m U ( 0 + 0 , x ) ϑ n 1 , , n m ( x ) d x = Ω l m U ( 0 0 , x ) ϑ n 1 , , n m ( x ) d x = u n 1 , , n m ( 0 0 ) .
We put
φ i n 1 , , n m = Ω l m φ i ( x ) ϑ n 1 , , n m ( x ) d x , i = 1 , 2 .
Then, taking (8) into account, from the conditions in (3), we obtain
u n 1 , , n m ( T ) = Ω l m U ( T , x ) ϑ n 1 , , n m ( x ) d x = Ω l m φ 1 ( x ) ϑ n 1 , , n m ( x ) d x = φ 1 n 1 , , n m ,
C D 0 t θ u n 1 , , n m ( T ) = Ω l m C D 0 t θ U ( T , x ) ϑ n 1 , , n m ( x ) d x = Ω l m φ 2 ( x ) ϑ n 1 , , n m ( x ) d x = φ 2 n 1 , , n m .
By the aid of the continuous conjugation condition (19) from (17) and (18), we have the relation that C 1 n 1 , , n m + = C 1 n 1 , , n m . To find the unknown coefficients of the integration C 1 n 1 , , n m and C 2 n 1 , , n m in (18), we use the conditions (20) and (21) and deduce the following system of linear algebraic equations:
ν τ n 1 , , n m Ψ 21 n 1 , , n m ( T , ε , ω ) + Ψ 22 n 1 , , n m ( T , ε , ω ) + + C 1 n 1 , , n m Ψ 23 n 1 , , n m ( T , ε , ω ) C 2 n 1 , , n m Ψ 24 n 1 , , n m ( T , ε , ω ) = φ 1 n 1 , , n m , ν τ n 1 , , n m D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω ) + D 0 t θ Ψ 22 n 1 , , n m ( T , ε , ω ) + + C 1 n 1 , , n m D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) C 2 n 1 , , n m D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) = φ 2 n 1 , , n m ,
where by D 0 t θ Ψ ( T ) is denoted D 0 t θ Ψ ( t ) t = T . We assume that
σ n 1 , , n m ( ω ) = Ψ 24 n 1 , , n m ( T , ε , ω ) · D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) Ψ 23 n 1 , , n m ( T , ε , ω ) · D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) 0 .
If the condition (23) is fulfilled, then the system (22) with respect to C 1 n 1 , , n m and C 2 n 1 , , n m is uniquely solvable. By solving this system (22), we arrive at the following presentations for these unknown coefficients
C 1 n 1 , , n m = 1 σ n 1 , , n m ( ω )
× φ 1 n 1 , , n m D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) + φ 2 n 1 , , n m Ψ 24 n 1 , , n m ( T , ε , ω ) ν τ n 1 , , n m
× Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω ) Ψ 21 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω )
+ Ψ 22 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 22 n 1 , , n m ( T , ε , ω ) ,
C 2 n 1 , , n m = 1 σ n 1 , , n m ( ω )
× φ 1 n 1 , , n m D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) + φ 2 n 1 , , n m Ψ 23 n 1 , , n m ( T , ε , ω ) ν τ n 1 , , n m
× Ψ 23 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω ) Ψ 21 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω )
+ Ψ 22 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) Ψ 23 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 22 n 1 , , n m ( T , ε , ω ) .
By substituting these results into (18) and taking into account C 1 n 1 , , n m + = C 1 n 1 , , n m in (17) and designation (12), we obtain the following representations for the Fourier coefficients of the main unknown functions in the positive and negative parts of the domain:
u n 1 , , n m + ( t , ε , ω , ν ) = φ 1 n 1 , , n m + φ 2 n 1 , , n m N 11 n 1 , , n m ( t , ε , ω ) + ν τ n 1 , , n m + N 12 n 1 , , n m ( t , ε ) ν τ n 1 , , n m N 13 n 1 , , n m ( t , ε , ω ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) × N 14 n 1 , , n m ( t , ε ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) N 15 n 1 , , n m ( t , ε , ω ) , t > 0 ,
u n 1 , , n m ( t , ε , ω , ν ) = φ 1 n 1 , , n m N 21 n 1 , , n m ( t , ε , ω ) + φ 2 n 1 , , n m N 22 n 1 , , n m ( t , ε , ω ) + ν τ n 1 , , n m N 23 n 1 , , n m ( t , ε , ω ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) N 24 n 1 , , n m ( t , ε , ω ) , t < 0 ,
where
N 11 n 1 , , n m ( t , ε , ω ) = 1 σ n 1 , , n m ( ω ) Ψ 13 n 1 , , n m ( t , ε ) Ψ 24 n 1 , , n m ( T , ε , ω ) ,
N 12 n 1 , , n m ( t , ε ) = Ψ 11 n 1 , , n m ( t , ε ) ,
N 13 n 1 , , n m ( t , ε , ω ) = 1 σ n 1 , , n m ( ω ) Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω )
Ψ 21 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) Ψ 13 n 1 , , n m ( t , ε ) ,
N 14 n 1 , , n m ( t , ε ) = Ψ ¯ 12 n 1 , , n m ( t , ε ) ,
N 15 n 1 , , n m ( t , ε , ω ) = 1 σ n 1 , , n m ( ω ) Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω )
Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) Ψ 13 n 1 , , n m ( t , ε ) ,
N 21 n 1 , , n m ( t , ε , ω ) = 1 σ n 1 , , n m ( ω ) Ψ 23 n 1 , , n m ( t , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω )
Ψ 24 n 1 , , n m ( t , ε , ω ) D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) ,
N 22 n 1 , , n m ( t , ε , ω ) = 1 σ n 1 , , n m ( ω ) Ψ 23 n 1 , , n m ( t , ε , ω ) Ψ 24 n 1 , , n m ( T , ε , ω )
+ Ψ 24 n 1 , , n m ( t , ε , ω ) Ψ 23 n 1 , , n m ( T , ε , ω ) , N 23 n 1 , , n m ( t , ε , ω ) = Ψ 21 n 1 , , n m ( t , ε , ω )
1 σ n 1 , , n m ( ω ) Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω )
Ψ 21 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) Ψ 23 n 1 , , n m ( t , ε , ω )
+ 1 σ n 1 , , n m ( ω ) Ψ 23 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 21 n 1 , , n m ( T , ε , ω )
Ψ 21 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) Ψ 24 n 1 , , n m ( t , ε , ω ) ,
N 24 n 1 , , n m ( t , ε , ω ) = Ψ ¯ 22 n 1 , , n m ( t , ε , ω )
+ 1 σ n 1 , , n m ( ω ) Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω )
Ψ 24 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) Ψ 23 n 1 , , n m ( t , ε , ω )
+ 1 σ n 1 , , n m ( ω ) Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω )
Ψ 23 n 1 , , n m ( T , ε , ω ) D 0 t θ Ψ ¯ 22 n 1 , , n m ( T , ε , ω ) Ψ 24 n 1 , , n m ( t , ε , ω ) ,
Ψ ¯ 12 n 1 , , n m ( t , ε ) = 0 t k 1 ( t s ) s α 1 1 E ( α 1 β 1 , α 1 ) , α 1 ε μ n 1 , , n m 2 s α 1 β 1 , μ n 1 , , n m 2 s α 1 d s ,
Ψ ¯ 22 n 1 , , n m ( t , ε , ω ) = t 0 k 2 ( s t ) ( s ) α 2 1 Ψ 25 n 1 , , n m ( t , ε , ω ) d s .
According to the degenerate kernels method, we substitute these presentations, (24) and (25), into designations (13) and (14):
τ n 1 , , n m + 1 ν χ 12 n 1 , , n m ( ε , ω ) + ν τ n 1 , , n m χ 13 n 1 , , n m ( ε , ω ) = φ 1 n 1 , , n m + φ 2 n 1 , , n m χ 11 n 1 , , n m ( ε , ω ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) χ 14 n 1 , , n m ( ε , ω ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) χ 15 n 1 , , n m ( ε , ω ) ,
τ n 1 , , n m 1 ν χ 23 n 1 , , n m ( ε , ω ) = φ 1 n 1 , , n m χ 21 n 1 , , n m ( ε , ω ) + φ 2 n 1 , , n m χ 22 n 1 , , n m ( ε , ω ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) χ 24 n 1 , , n m ( ε , ω ) ,
where
χ 1 i n 1 , , n m ( ε , ω ) = 0 T b 1 ( s ) N 1 i n 1 , , n m ( s , ε , ω ) d s , i = 1 , 5 ¯ ,
χ 2 i n 1 , , n m ( ε , ω ) = T 0 b 2 ( s ) N 2 i n 1 , , n m ( s , ε , ω ) d s , i = 1 , 4 ¯ .
We solve the linear algebraic Equations (26) and (27) as a system of algebraic equations with respect to quantities τ n 1 , , n m + and τ n 1 , , n m . If the following conditions are fulfilled
ν χ 12 n 1 , , n m ( ε , ω ) 1 , ν χ 23 n 1 , , n m ( ε , ω ) 1 ,
then, from (26) and (27), we derive
τ n 1 , , n m + = φ 1 n 1 , , n m M 11 n 1 , , n m ( ε , ω ) + φ 2 n 1 , , n m M 12 n 1 , , n m ( ε , ω ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) M 13 n 1 , , n m ( ε , ω ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) M 14 n 1 , , n m ( ε , ω ) ,
τ n 1 , , n m = φ 1 n 1 , , n m M 21 n 1 , , n m ( ε , ω ) + φ 2 n 1 , , n m M 22 n 1 , , n m ( ε , ω ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) M 23 n 1 , , n m ( ε , ω ) ,
where
M 1 i n 1 , , n m ( ε , ω ) = 1 1 ν χ 12 n 1 , , n m ( ε , ω ) χ 11 n 1 , , n m ( ε , ω ) ν χ 13 n 1 , , n m ( ε , ω ) χ 2 i n 1 , , n m ( ε , ω ) 1 ν χ 23 n 1 , , n m ( ε , ω ) ,
i = 1 , 2 , M 13 n 1 , , n m ( ε , ω ) = χ 14 n 1 , , n m ( ε , ω ) 1 ν χ 12 n 1 , , n m ( ε , ω ) ,
M 14 n 1 , , n m ( ε , ω ) = 1 1 ν χ 12 n 1 , , n m ( ε , ω ) χ 15 n 1 , , n m ( ε , ω ) ν χ 13 n 1 , , n m ( ε , ω ) χ 24 n 1 , , n m ( ε , ω ) 1 ν χ 23 n 1 , , n m ( ε , ω ) ,
M 2 i n 1 , , n m ( ε , ω ) = χ 2 i n 1 , , n m ( ε , ω ) 1 ν χ 23 n 1 , , n m ( ε , ω ) , i = 1 , 2 , M 23 n 1 , , n m ( ε , ω ) = χ 24 n 1 , , n m ( ε , ω ) 1 ν χ 23 n 1 , , n m ( ε , ω ) .
Substituting presentations (29) and (30) of τ n 1 , , n m ± into (24) and (25), we derive
u n 1 , , n m + ( t , ε , ω , ν ) = φ 1 n 1 , , n m Q 11 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m Q 12 n 1 , , n m ( t , ε , ω , ν ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) Q 13 n 1 , , n m ( t , ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
u n 1 , , n m ( t , ε , ω , ν ) = φ 1 n 1 , , n m Q 21 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m Q 22 n 1 , , n m ( t , ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 ,
where
Q 1 i n 1 , , n m ( t , ε , ω , ν ) = N 11 n 1 , , n m ( t , ε , ω ) + ν N 12 n 1 , , n m ( t , ε , ω ) M 1 i n 1 , , n m ( ε , ω )
ν N 13 n 1 , , n m ( t , ε , ω ) M 2 i n 1 , , n m ( ε , ω ) , i = 1 , 2 ,
Q 13 n 1 , , n m ( t , ε , ω , ν ) = N 14 n 1 , , n m ( t , ε , ω ) + ν N 12 n 1 , , n m ( t , ε , ω ) M 13 n 1 , , n m ( ε , ω ) ,
Q 14 n 1 , , n m ( t , ε , ω , ν ) = N 15 n 1 , , n m ( t , ε , ω ) + ν N 12 n 1 , , n m ( t , ε , ω ) M 14 n 1 , , n m ( ε , ω )
ν N 13 n 1 , , n m ( t , ε , ω ) M 23 n 1 , , n m ( ε , ω ) ,
Q 2 i n 1 , , n m ( t , ε , ω , ν ) = N 2 i n 1 , , n m ( t , ε , ω ) + ν N 23 n 1 , , n m ( t , ε , ω ) M 2 i n 1 , , n m ( ε , ω ) , i = 1 , 2 ,
Q 23 n 1 , , n m ( t , ε , ω , ν ) = N 24 n 1 , , n m ( t , ε , ω ) + ν N 23 n 1 , , n m ( t , ε , ω ) M 23 n 1 , , n m ( ε , ω ) .
Now, we substitute presentations (31) and (32) into the Fourier series (7) and obtain the following formal solution of the direct problem (1)–(4)
U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) φ 1 n 1 , , n m Q 11 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m Q 12 n 1 , , n m ( t , ε , ω , ν ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) Q 13 n 1 , , n m ( t , ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) φ 1 n 1 , , n m Q 21 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m × Q 22 n 1 , , n m ( t , ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 .
We suppose that the conditions of (23) were violated for some values of spectral parameter ω . So, we have to consider the algebraic equation with respect to spectral parameter ω
σ n 1 , , n m ( ω ) = Ψ 24 n 1 , , n m ( T , ε , ω ) · D 0 t θ Ψ 23 n 1 , , n m ( T , ε , ω ) Ψ 23 n 1 , , n m ( T , ε , ω ) · D 0 t θ Ψ 24 n 1 , , n m ( T , ε , ω ) = 0 .
The set of positive solutions of this algebraic Equation (35) with respect to the spectral parameter ω , we denote by 1 . We call these values ω 1 as irregular values and, for these values, the condition (23) is violated. Another set Λ 1 = ( 0 ; ) 1 is called the set of regular values of the spectral parameter ω and, for these regular values, the condition (23) is fulfilled.
Now, we assume that the conditions in (28) are violated ν χ 12 n 1 , , n m ( ε , ω ) = 1 , ν χ 23 n 1 , , n m ( ε , ω ) = 1 . Hence, we have
ν 1 = 1 χ 12 n 1 , , n m ( ε , ω ) , ν 2 = 1 χ 23 n 1 , , n m ( ε , ω ) .
For regular values ω Λ 1 there hold χ 12 n 1 , , n m ( ε , ω ) 0 , χ 23 n 1 , , n m ( ε , ω ) 0 . So, we denote the set { ν 1 , ν 2 } by 2 . Then a set Λ 2 = ( ; 0 ) ( 0 ; ) 2 is called the set of regular values of the spectral parameter ν . Therefore, for all values of ν Λ 2 , condition (28) is satisfied. We use the following notation = n 1 , , n m N ; ω Λ 1 ; ν Λ 2 , where N is the set of natural numbers. This is the set on which all values of the spectral parameters ω and ν are regular. Therefore, in this case, we study the solution of the direct problem (1)–(4) in the domain Ω as Fourier series (33) and (34).

4. Redefinition Functions

Suppose that the functions ψ i ( x ) expand in the Fourier series
ψ i ( x ) = n 1 , , n m = 1 ψ i n 1 , , n m ϑ n 1 , , n m ( x ) ,
where
ψ i n 1 , , n m = Ω l m ψ i ( x ) ϑ n 1 , , n m ( x ) d x , i = 1 , 2 , n 1 , , n m = 1 , 2 ,
By virtue of series (33), (34) and (36), we apply conditions (5) and (6):
n 1 , , n m = 1 ψ 1 n 1 , , n m ϑ n 1 , , n m ( x ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x )
× 0 T Φ 1 ( t ) φ 1 n 1 , , n m Q 11 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m Q 12 n 1 , , n m ( t , ε , ω , ν )
+ g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) Q 13 n 1 , , n m ( t , ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 14 n 1 , , n m ( t , ε , ω , ν ) d t ,
n 1 , , n m = 1 ψ 2 n 1 , , n m ϑ n 1 , , n m ( x ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x )
× T 0 Φ 2 ( t ) φ 1 n 1 , , n m Q 21 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m Q 22 n 1 , , n m ( t , ε , ω , ν )
+ g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Q 23 n 1 , , n m ( t , ε , ω , ν ) d t .
Hence, we obtain
ψ 1 n 1 , , n m = φ 1 n 1 , , n m Υ 11 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Υ 12 n 1 , , n m ( ε , ω , ν ) + g 1 n 1 , , n m + f 1 n 1 , , n m ( V 1 ) Υ 13 n 1 , , n m ( ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Υ 14 n 1 , , n m ( ε , ω , ν ) ,
ψ 2 n 1 , , n m = φ 1 n 1 , , n m Υ 21 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Υ 22 n 1 , , n m ( ε , ω , ν ) + g 2 n 1 , , n m + f 2 n 1 , , n m ( V 2 ) Υ 23 n 1 , , n m ( ε , ω , ν ) ,
where
Υ 1 i n 1 , , n m ( ε , ω , ν ) = 0 T Φ 1 ( t ) Q 1 i n 1 , , n m ( t , ε , ω , ν ) d t , i = 1 , 4 ¯ ,
Υ 2 i n 1 , , n m ( ε , ω , ν ) = T 0 Φ 2 ( t ) Q 2 i n 1 , , n m ( t , ε , ω , ν ) d t , i = 1 , 3 ¯ .
The relations of (37) and (38) we consider as a system of functional algebraic equations with respect to coefficients of redefinition functions. By solving this system, we obtain the following representations
g 1 n 1 , , n m ( ε , ω , ν ) + f 1 n 1 , , n m ( V 1 ) = ψ 1 n 1 , , n m Δ 11 n 1 , , n m ( ε , ω , ν ) + ψ 2 n 1 , , n m Δ 12 n 1 , , n m ( ε , ω , ν ) + φ 1 n 1 , , n m Δ 13 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Δ 14 n 1 , , n m ( ε , ω , ν ) ,
g 2 n 1 , , n m ( ε , ω , ν ) + f 2 n 1 , , n m ( V 2 ) = ψ 2 n 1 , , n m Δ 21 n 1 , , n m ( ε , ω , ν ) + φ 1 n 1 , , n m Δ 22 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Δ 23 n 1 , , n m ( ε , ω , ν ) ,
where
Δ 11 n 1 , , n m ( ε , ω , ν ) = Υ 13 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 12 n 1 , , n m ( ε , ω , ν ) = Υ 23 n 1 , , n m ( ε , ω , ν ) Υ 13 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 13 n 1 , , n m ( ε , ω , ν ) = Υ 11 n 1 , , n m ( ε , ω , ν ) + Υ 21 n 1 , , n m ( ε , ω , ν ) Υ 23 n 1 , , n m ( ε , ω , ν ) Υ 13 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 14 n 1 , , n m ( ε , ω , ν ) = Υ 12 n 1 , , n m ( ε , ω , ν ) + Υ 22 n 1 , , n m ( ε , ω , ν ) Υ 23 n 1 , , n m ( ε , ω , ν ) Υ 13 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 21 n 1 , , n m ( ε , ω , ν ) = Υ 23 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 22 n 1 , , n m ( ε , ω , ν ) = Υ 21 n 1 , , n m ( ε , ω , ν ) Υ 23 n 1 , , n m ( ε , ω , ν ) 1 ,
Δ 23 n 1 , , n m ( ε , ω , ν ) = Υ 22 n 1 , , n m ( ε , ω , ν ) Υ 23 n 1 , , n m ( ε , ω , ν ) 1 .
We rewrite Formulas (39) and (40) in the form of countable systems of nonlinear integral equations (CSNIE)
g i n 1 , , n m ( ε , ω , ν ) = I i g i n 1 , , n m c i n 1 , , n m ( ε , ω , ν ) Ω l m f i y , Ω l m Θ i ( z ) n 1 , , n m = 1 g i n 1 , , n m ( ε , ω , ν ) ϑ n 1 , , n m ( z ) d z ϑ n 1 , , n m ( y ) d y , i = 1 , 2 ,
where
c 1 n 1 , , n m ( ε , ω , ν ) = ψ 1 n 1 , , n m Δ 11 n 1 , , n m ( ε , ω , ν ) + ψ 2 n 1 , , n m Δ 12 n 1 , , n m ( ε , ω , ν )
+ φ 1 n 1 , , n m Δ 13 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Δ 14 n 1 , , n m ( ε , ω , ν ) ,
c 2 n 1 , , n m ( ε , ω , ν ) = ψ 2 n 1 , , n m Δ 21 n 1 , , n m ( ε , ω , ν )
+ φ 1 n 1 , , n m Δ 22 n 1 , , n m ( ε , ω , ν ) + φ 2 n 1 , , n m Δ 23 n 1 , , n m ( ε , ω , ν ) .

5. Unique Solvability of CSNIE (41)

We use the concepts of the following well-known Banach spaces, including a Hilbert coordinate space 2 of number sequences b n 1 , , n m n 1 , , n m = 1 with the norm
b 2 = n 1 , , n m = 1 b n 1 , , n m 2 < .
We also use the space L 2 ( Ω l m ) of square-summable functions on the domain Ω l m with the norm
ϑ ( x ) L 2 ( Ω l m ) = Ω l m ϑ ( x ) 2 d x < .
In the process of proofing the unique solvability of CSNIE (41), we need the following conditions.
Smoothness conditions. Let functions
φ i ( x ) , ψ i ( x ) C 2 ( Ω l m ) , f i x , Ω l m Θ i ( y ) g i ( y ) d y C x 2 ( Ω l m × R ) , i = 1 , 2
in the domain Ω l m have piecewise continuous third order derivatives.
Then, by integrating them in parts three times over all variables x 1 , x 2 , , x m , we obtain the following formulas [41]
n 1 , , n m = 1 φ i n 1 , , n m ( 3 m ) 2 2 l m Ω l m 3 m φ i ( x ) x 1 3 x 2 3 x m 3 2 d x ,
n 1 , , n m = 1 ψ i n 1 , , n m ( 3 m ) 2 2 l m Ω l m 3 m ψ i ( x ) x 1 3 x 2 3 x m 3 2 d x ,
φ i n 1 , , n m = l π 3 m φ i n 1 , , n m ( 3 m ) n 1 3 n m 3 , ψ i n 1 , , n m = l π 3 m ψ i n 1 , , n m ( 3 m ) n 1 3 n m 3 ,
where
φ i n 1 , , n m ( 3 m ) = Ω l m 3 m φ i ( x ) x 1 3 x 2 3 x m 3 ϑ n 1 , , n m ( x ) d x ,
ψ i n 1 , , n m ( 3 m ) = Ω l m 5 m ψ i ( x ) x 1 3 x 2 3 x m 3 ϑ n 1 , , n m ( x ) d x , i = 1 , 2 .
We obtain also that
f i n 1 , , n m ( V i ) = l π 3 m f i n 1 , , n m ( 3 m ) ( x , V i ) n 1 3 n m 3 ,
n 1 , , n m = 1 f i n 1 , , n m ( 3 m ) ( V i ) 2 2 l m Ω l m 3 m f i ( x , V i ) x 1 3 x 2 3 x m 3 2 d x ,
where
f i n 1 , , n m ( 3 m ) ( V i ) = Ω l m 3 m f i ( x , V i ) x 1 3 x 2 3 x m 3 ϑ n 1 , , n m ( x ) d x , i = 1 , 2 .
We use also the following well known properties of the Mittag–Leffler function:
( 1 ) For all k > 0 , α 0 , β 0 , γ 0 ( 0 ; 2 ] , α 0 β 0 γ 0 , t 0 the function t β 0 1 E α 0 , β 0 , γ 0 k t α , k t β is complete and monotonous and there holds
( 1 ) s t β 0 1 E ( α 0 , β 0 ) , γ 0 k t α 0 , k t β 0 ( s ) 0 , s = 0 , 1 , 2 ,
( 2 ) For all α 0 , β 0 ( 0 , 2 ) , γ R and arg z 1 = π , there hold the following estimates
E ( α 0 , β 0 ) , γ 0 z 1 , z 2 C 1 1 + | z 1 | ,
E ( α 0 , β 0 ) , γ 0 ε 1 z 1 , z 2 E ( α 0 , β 0 ) , γ 0 ε 2 z 1 , z 2 ε 1 ε 2 C 2 1 + | z 1 | ,
where 0 < C i = const does not depend on z, ε i 0 ; ε 0 , 0 < ε 0 = const , i = 1 , 2 .
According to the properties of the Mittag–Leffler function (Formulas (47) and (48)) the quantities Δ 1 , j n 1 , , n m ( ε , ω , ν ) ( j = 1 , 4 ¯ ) and Δ 2 , j n 1 , , n m ( ε , ω , ν ) ( j = 1 , 3 ¯ ) are uniformly bounded. So, for any positive integers n 1 , , n m , there exist finite constant numbers C 0 i ( i = 1 , 2 ) , by which the following estimates take place
max n 1 , , n m N max j Δ 1 , j n 1 , , n m ( ε , ω , ν ) C 01 , j = 1 , 4 ¯ ,
max n 1 , , n m / m a t h b b N max j Δ 2 , j n 1 , , n m ( ε , ω , ν ) C 02 , j = 1 , 3 ¯ ,
where 0 < C 0 i = const , i = 1 , 2 .
Lemma 1.
Suppose that the smoothness conditions are fulfilled and
f i ( x , V 1 i ) f i ( x , V 2 i ) K 1 i ( x ) V 1 i V 2 i , ρ < 1 ,
where ρ = C 04 γ 3 Θ i ( x ) L 2 ( Ω l m ) , γ 3 = C 03 l π 3 m 2 l m ,
C 03 = n 1 , , n m = 1 1 n 1 6 n m 6 < ; max i 3 m K 1 i ( x ) x 1 3 x 2 3 x m 3 L 2 ( Ω l m ) C 04 < , i = 1 , 2 .
Then, for regular values of spectral parameters ω and ν, CSNIE (41) is uniquely solvable in the space 2 . In this case, successive approximations are defined as follows:
g i n 1 , , n m 0 ( ε , ω , ν ) = c i n 1 , , n m , g i n 1 , , n m k + 1 ( ε , ω , ν ) = I i n 1 , , n m ( g i k ) , i = 1 , 2 .
Proof. 
We apply the method of successive approximations and the method of compressive mappings. We use Formulas (42)–(44) and estimates (50) and (51). By the aid of the Cauchy–Schwartz inequality and the Bessel inequality for the zeroth approximation of the coefficients of the redefinition functions from successive approximations (52), we obtain
g 1 0 ( ε , ω , ν ) 2 n 1 , , n m = 1 c 1 n 1 , , n m ( ε , ω , ν ) n 1 , , n m = 1 ψ 1 n 1 , , n m Δ 1 1 n 1 , , n m ( ε , ω , ν ) + n 1 , , n m = 1 ψ 2 n 1 , , n m Δ 12 n 1 , , n m ( ε , ω , ν ) + n 1 , , n m = 1 φ 1 n 1 , , n m Δ 13 n 1 , , n m ( ε , ω , ν ) + n 1 , , n m = 1 φ 2 n 1 , , n m Δ 14 n 1 , , n m ( ε , ω , ν ) C 01 n 1 , , n m = 1 ψ 1 n 1 , , n m + n 1 , , n m = 1 ψ 2 n 1 , , n m + n 1 , , n m = 1 φ 1 n 1 , , n m + n 1 , , n m = 1 φ 2 n 1 , , n m C 01 l π 3 m n 1 , , n m = 1 ψ 1 n 1 , , n m ( 3 m ) n 1 3 n m 3 + n 1 , , n m = 1 ψ 2 n 1 , , n m ( 3 m ) n 1 3 n m 3 + n 1 , , n m = 1 φ 1 n 1 , , n m ( 3 m ) n 1 3 n m 3 + n 1 , , n m = 1 φ 2 n 1 , , n m ( 3 m ) n 1 3 n m 3 C 01 l π 3 m n 1 , , n m = 1 1 n 1 6 n m 6 ψ 1 ( 3 m ) 2 + ψ 2 ( 3 m ) 2 + φ 1 ( 3 m ) 2 + φ 2 ( 3 m ) 2 γ 1 3 m ψ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m < ,
where γ 1 = C 01 C 03 2 l m l π 3 m , C 03 = n 1 , , n m = 1 1 n 1 6 n m 6 < ;
g 2 0 ( ε , ω , ν ) 2 n 1 , , n m = 1 c 2 n 1 , , n m ( ε , ω , ν ) n 1 , , n m = 1 ψ 2 n 1 , , n m Δ 21 n 1 , , n m ( ε , ω , ν ) + n 1 , , n m = 1 φ 1 n 1 , , n m Δ 22 n 1 , , n m ( ε , ω , ν ) + n 1 , , n m = 1 φ 2 n 1 , , n m Δ 23 n 1 , , n m ( ε , ω , ν ) C 02 l π 3 m n 1 , , n m = 1 ψ 2 n 1 , , n m ( 3 m ) n 1 3 n m 3 + n 1 , , n m = 1 φ 1 n 1 , , n m ( 3 m ) n 1 3 n m 3 + n 1 , , n m = 1 φ 2 n 1 , , n m ( 3 m ) n 1 3 n m 3 γ 2 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m < ,
where γ 2 = C 02 C 03 2 l m l π 3 m .
By Formulas (45) and (46), using the Cauchy–Schwartz inequality and Bessel inequality for the first difference of approximation (52), we obtain
g i 1 ( ε , ω , ν ) g i 0 ( ε , ω , ν ) 2 n 1 , , n m = 1 Ω l m f i x , Ω l m Θ i ( y ) g i 0 ( y , ε , ω , ν ) d y ϑ n 1 , , n m ( x ) d x l π 3 m n 1 , , n m = 1 f i n 1 , , n m ( 3 m ) x , V i 0 n 1 3 n m 3 γ 3 3 m f i x , V i 0 x 1 3 x 2 3 x m 3 L 2 Ω l m < ,
where γ 3 = C 03 l π 3 m 2 l m , V i 0 = Ω l m Θ i ( x ) g i 0 ( x , ε , ω , ν ) d x , i = 1 , 2 .
Analogously, by the condition of the lemma and expansion (9), using the Cauchy–Schwartz inequality and Bessel inequality for an arbitrary difference of approximation (52), we obtain
g i k + 1 ( ε , ω , ν ) g i k ( ε , ω , ν ) 2 γ 3 3 m x 1 3 x 2 3 x m 3 f i ( x , V i k ) f i ( x , V i k 1 ) L 2 ( Ω l m ) γ 3 Ω l m Θ i ( y ) · g i k ( y , ε , ω , ν ) g i k 1 ( y , ε , ω , ν ) d y 3 m K 1 i ( x ) x 1 3 x 2 3 x m 3 L 2 ( Ω l m ) C 04 γ 3 Ω l m Θ i ( y ) n 1 , , n m = 1 g i n 1 , , n m k ( ε , ω , ν ) g i n 1 , , n m k 1 ( ε , ω , ν ) · ϑ n 1 , , n m ( y ) d y ρ g i k ( ε , ω , ν ) g i k 1 ( ε , ω , ν ) 2 , i = 1 , 2 ,
where
ρ = C 04 γ 3 Θ i ( x ) L 2 Ω l m , max i 3 m K 1 i ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m C 04 < , i = 1 , 2 .
By the condition of the lemma, ρ < 1 . Therefore, it follows from estimate (56) that the operators on the right-hand side of (41) are contracting. From the estimates (53)–(56), it is implied that there exists a unique pair of fixed points g 1 n 1 , , n m ( ε , ω , ν ) ; g 2 n 1 , , n m ( ε , ω , ν ) , which is a solution of CSNIE (41) in the space 2 . The Lemma 1 is proved.  □

6. Convergence of Fourier Series (57)

Now, we determine the redefinition functions. In this order, we substitute representations (41) into the Fourier series (9) and obtain
g i ( x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) c i n 1 , , n m ( ε , ω , ν ) Ω l m f i y , Ω l m Θ i ( z ) n 1 , , n m = 1 g i n 1 , , n m ( ε , ω , ν ) ϑ n 1 , , n m ( z ) d z ϑ n 1 , , n m ( y ) d y , i = 1 , 2 .
We prove that the following lemma holds.
Lemma 2.
Assume that the conditions of Lemma 1 are satisfied. Then for regular values of spectral parameters ω and ν, the series (57) converge absolutely.
Proof. 
We use estimates (53)–(55). Using the Cauchy–Schwartz inequality and Bessel inequality for series (57), we obtain the following estimates
g 1 ( x , ε , ω , ν ) n 1 , , n m = 1 ϑ n 1 , , n m ( x ) [ c 1 n 1 , , n m ( ε , ω , ν ) + Ω l m f 1 y , Ω l m Θ 1 ( z ) n 1 , , n m = 1 g 1 n 1 , , n m ( ε , ω , ν ) ϑ n 1 , , n m ( z ) d z ϑ n 1 , , n m ( y ) d y C 03 2 l 3 m 2 l π 3 m C 01 ψ 1 ( 3 m ) 2 + C 01 ψ 2 ( 3 m ) 2 + C 01 φ 1 ( 3 m ) 2 + C 01 φ 2 ( 3 m ) 2 + f 1 ( 3 m ) ( V 1 ) 2 γ 4 3 m ψ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m f 1 ( x , V 1 ) x 1 3 x 2 3 x m 3 L 2 ( Ω l m ) < ,
g 2 ( x , ε , ω , ν ) n 1 , , n m = 1 ϑ n 1 , , n m ( x ) [ c 2 n 1 , , n m ( ε , ω , ν ) + Ω l m f 2 y , Ω l m Θ 2 ( z ) n 1 , , n m = 1 g 2 n 1 , , n m ( ε , ω , ν ) ϑ n 1 , , n m ( z ) d z ϑ n 1 , , n m ( y ) d y γ 4 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 L 2 Ω l m + 3 m f 2 ( x , V 2 ) x 1 3 x 2 3 x m 3 L 2 Ω l m < ,
where
3 m f i x , V 2 x 1 3 x 2 3 x m 3 L 2 Ω l m = Ω l m 3 m f i x , V i x 1 3 x 2 3 x m 3 2 d x ,
V i = Ω l m Θ i ( y ) n 1 , , n m = 1 g i n 1 , , n m ( ε , ω , ν ) ϑ n 1 , , n m ( z ) d z , i = 1 , 2 ,
γ 4 = C 03 C 05 2 l 3 m 2 l π 3 m , C 05 = max C 01 ; C 02 ; 1 .
From (58) and (59) the convergence of series (57) is implied. Lemma 2 is proved.  □
So, we determined the redefinition functions as a Fourier series (57). Using representations (39) and (40),we can present Fourier series (33) and (34) of the main unknown functions as
U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) ψ 1 n 1 , , n m W 11 n 1 , , n m ( t , ε , ω , ν ) + ψ 2 n 1 , , n m W 12 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 13 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) ψ 2 n 1 , , n m W 21 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 22 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 ,
where
W i 1 n 1 , , n m ( t , ε , ω , ν ) = Δ i 1 n 1 , , n m ( ε , ω , ν ) Q i 3 n 1 , , n m ( t , ε , ω , ν ) , i = 1 , 2 ,
W 12 n 1 , , n m ( t , ε , ω , ν ) = Δ 12 n 1 , , n m ( ε , ω , ν ) Q 13 n 1 , , n m ( t , ε , ω , ν ) + Δ 21 n 1 , , n m ( ε , ω , ν ) Q 14 n 1 , , n m ( t , ε , ω , ν ) ,
W 1 j n 1 , , n m ( t , ε , ω , ν ) = Q 1 j 2 n 1 , , n m ( t , ε , ω , ν ) + Δ 1 j n 1 , , n m ( ε , ω , ν ) Q 13 n 1 , , n m ( t , ε , ω , ν )
+ Δ 2 j 1 n 1 , , n m ( ε , ω , ν ) Q 14 n 1 , , n m ( t , ε , ω , ν ) , j = 3 , 4 ,
W 2 k n 1 , , n m ( t , ε , ω , ν ) = Q 2 k 1 n 1 , , n m ( t , ε , ω , ν ) + Δ 2 k n 1 , , n m ( ε , ω , ν ) Q 23 n 1 , , n m ( t , ε , ω , ν ) , k = 2 , 3 .
To establish the uniqueness of the function U ( t , x , ε , ω , ν ) , we suppose that there are two solutions U 1 and U 2 to this problem. Then, their difference U = U 1 U 2 is a solution of Equation (1), satisfying conditions (2)–(6) with functions φ i ( x ) 0 , ψ i ( x ) 0 ( i = 1 , 2 ) . Then, for φ i n 1 , , n m = ψ i n 1 , , n m = 0 ( i = 1 , 2 ) , it follows from Formulas (60) and (61) in the domain Ω that
Ω l m U ( t , x , ε , ω , ν ) ϑ n 1 , , n m ( x ) d x = 0 .
Hence, by virtue of the completeness of the systems of eigenfunctions 2 l sin π n 1 l x 1 , 2 l sin π n 2 l x 2 , , 2 l sin π n m l x m in L 2 Ω l m , we deduce that U ( t , x , ε , ω , ν ) 0 for all x Ω l m [ 0 ; l ] m and t [ T ; T ] .
Therefore, for regular values of spectral parameters ω and ν , the function U ( t , x , ε , ω , ν ) is a unique solution tp the mixed type integro-differential Equation (1) with conditions (2)–(6), if this function exists in the domain Ω .
Lemma 3.
Let smoothness conditions hold. Then, for regular values of spectral parameters ω and ν, series (60) and (61) converge. At the same time, their term by term differentiation is possible.
Proof. 
According to the properties of the Mittag–Leffler function (Formulas (47) and (48)), the functions W 1 i n 1 , , n m ( t , ε , ω , ν ) ( i = 1 , 4 ¯ ) and W 2 j n 1 , , n m ( t , ε , ω , ν ) ( j = 1 , 3 ¯ ) are uniformly bounded on the segment [ T ; T ] . So, for any positive integers n 1 , , n m , there exist finite constant numbers C 1 k ( k = 1 , 2 ) ; then, the following estimates take place
max n 1 , , n m N max i = 1 , 4 ¯ W 1 i n 1 , , n m ( t , ε , ω , ν ) C 11 , max n 1 , , n m N max j = 1 , 3 ¯ W 1 j n 1 , , n m ( t , ε , ω , ν ) C 12 ,
where C 1 k = const , k = 1 , 2 .
Analogously to the estimates (58) and (59), by applying estimates (62), the Cauchy–Schwartz inequality and Bessel inequality for series (60) and (61), we
U ( t , x , ε , ω , ν ) n 1 , , n m = 1 u n 1 , , n m + ( t , ε , ω , ν ) · ϑ n 1 , , n m ( x ) 2 l m C 11 n 1 , , n m = 1 ψ 1 n 1 , , n m + ψ 2 n 1 , , n m + φ 1 n 1 , , n m + φ 2 n 1 , , n m γ 5 Ω l m 3 m ψ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x < ,
where γ 5 = 2 l 3 m 2 C 11 C 03 l π 3 m ,
U ( t , x , ε , ω , ν ) n 1 , , n m = 1 u n 1 , , n m ( t , ε , ω , ν ) · ϑ n 1 , , n m ( x ) γ 6 Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x < ,
where γ 6 = 2 l 3 m 2 C 12 C 03 l π 3 m .
It follows from estimates (63) and (64) that the series (60) and (61) are convergent absolutely and uniformly in the domain Ω ¯ for the
( n 1 , , n m , ω , ν ) = n 1 , , n m N ; ω Λ 1 ; ν Λ 2 .
Therefore, for the ( n 1 , , n m , ω , ν ) functions, (63) and (64) formally differentiate in Ω ¯ the required number of times
C D 0 t α 1 U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) × ψ 1 n 1 , , n m C D 0 t α 1 W 11 n 1 , , n m ( t , ε , ω , ν ) + ψ 2 n 1 , , n m C D 0 t α 1 W 12 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m C D 0 t α 1 W 13 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m C D 0 t α 1 W 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
C D 0 t α 2 U ( t , x , ε , ω , ν ) = n 1 , , n m = 1 ϑ n 1 , , n m ( x ) ψ 2 n 1 , , n m C D 0 t α 2 W 21 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m C D 0 t α 2 W 22 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m C D 0 t α 2 W 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 ,
U x 1 x 1 ( t , x , ε , ω , ν ) = n 1 , , n m = 1 π n 1 l 2 ϑ n 1 , , n m ( x ) ψ 1 n 1 , , n m W 11 n 1 , , n m ( t , ε , ω , ν ) + ψ 2 n 1 , , n m W 12 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 13 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
U x 1 x 1 ( t , x , ε , ω , ν ) = n 1 , , n m = 1 π n 1 l 2 ϑ n 1 , , n m ( x ) ψ 2 n 1 , , n m W 21 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 22 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 ,
U x 2 x 2 ( t , x , ε , ω , ν ) = n 1 , , n m = 1 π n 2 l 2 ϑ n 1 , , n m ( x ) ψ 1 n 1 , , n m W 11 n 1 , , n m ( t , ε , ω , ν ) + ψ 2 n 1 , , n m W 12 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 13 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 14 n 1 , , n m ( t , ε , ω , ν ) , t > 0 ,
U x 2 x 2 ( t , x , ε , ω , ν ) = n 1 , , n m = 1 π n 2 l 2 ϑ n 1 , , n m ( x ) ψ 2 n 1 , , n m W 21 n 1 , , n m ( t , ε , ω , ν ) + φ 1 n 1 , , n m W 22 n 1 , , n m ( t , ε , ω , ν ) + φ 2 n 1 , , n m W 23 n 1 , , n m ( t , ε , ω , ν ) , t < 0 .
The expansions of the following functions into Fourier series are defined in the domain Ω in a similar way
U x 3 x 3 ( t , x , ε , ω , ν ) , , U x m x m ( t , x , ε , ω , ν ) , C D 0 t α 1 U x 1 x 1 ( t , x , ε , ω , ν ) ,
C D 0 t α 2 U x 1 x 1 ( t , x , ε , ω , ν ) , C D 0 t α 1 U x 2 x 2 ( t , x , ε , ω , ν ) , , C D 0 t α 2 U x 2 x 2 ( t , x , ε , ω , ν ) , ,
C D 0 t α 1 U x m x m ( t , x , ε , ω , ν ) , C D 0 t α 2 U x m x m ( t , x , ε , ω , ν ) .
The convergence of series (65) and (66) is proved similarly to the proof of the convergence of series (60) and (61). So, it is enough to show the convergence of series (67) and (70). Taking into account Formulas (42)–(44) and estimates (62) and applying the Cauchy–Schwartz inequality and Bessel inequality, we obtain
U x 1 x 1 ( t , x , ε , ω , ν ) n 1 , , n m = 1 π n 1 l 2 u n 1 , , n m + ( t , ε , ω , ν ) · ϑ n 1 , , n m ( x )
2 l m π l 2 C 11 n 1 , , n m = 1 n 1 2 ψ 1 n 1 , , n m + ψ 2 n 1 , , n m + φ 1 n 1 , , n m + φ 2 n 1 , , n m
2 l m C 11 l π 3 m 2 n 1 , , n m = 1 1 n 1 n 2 3 n m 3 ψ 1 n 1 , , n m ( 3 m ) + n 1 , , n m = 1 1 n 1 n 2 3 n m 3 ψ 2 n 1 , , n m ( 3 m )
+ n 1 , , n m = 1 1 n 1 n 2 3 n m 3 φ 1 n 1 , , n m ( 3 m ) + n 1 , , n m = 1 1 n 1 n 2 3 n m 3 φ 2 n 1 , , n m ( 3 m )
γ 7 Ω l m 3 m ψ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x
+ Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x < ,
where γ 7 = 2 l 3 m 2 C 11 C 06 l π 3 m 2 , C 06 = n 1 , , n m = 1 1 n 1 2 n 2 6 n m 6 ;
U x 2 x 2 ( t , x , ε , ω , ν ) n 1 , , n m = 1 π n 2 l 2 u n 1 , , n m ( t , ε , ω , ν ) · ϑ n 1 , , n m ( x )
2 l m π l 2 C 12 n 1 , , n m = 1 n 2 2 ψ 2 n 1 , , n m + φ 1 n 1 , , n m + φ 2 n 1 , , n m
2 l m C 12 l π 3 m 2 n 1 , , n m = 1 1 n 1 3 n 2 n 3 3 n m 3 ψ 2 n 1 , , n m ( 3 m )
+ n 1 , , n m = 1 1 n 1 3 n 2 n 3 3 n m 3 φ 1 n 1 , , n m ( 3 m ) + n 1 , , n m = 1 1 n 1 3 n 2 n 3 3 n m 3 φ 2 n 1 , , n m ( 3 m )
γ 8 Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x
+ Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x < ,
where γ 8 = 2 l 3 m 2 C 12 C 07 l π 3 m 2 , C 07 = n 1 , , n m = 1 1 n 1 6 n 2 2 n m 6 .
The convergence of series (68) and (69) is similar to the convergence of series (67) and (70). The convergence of Fourier series for functions
U x 3 x 3 ( t , x , ε , ω , ν ) , , U x m x m ( t , x , ε , ω , ν ) , C D 0 t α 1 U x 1 x 1 ( t , x , ε , ω , ν ) ,
C D 0 t α 2 U x 1 x 1 ( t , x , ε , ω , ν ) , C D 0 t α 1 U x 2 x 2 ( t , x , ε , ω , ν ) , , C D 0 t α 2 U x 2 x 2 ( t , x , ε , ω , ν ) , ,
C D 0 t α 1 U x m x m ( t , x , ε , ω , ν ) , C D 0 t α 2 U x m x m ( t , x , ε , ω , ν )
is proved in a similar way in the domain Ω . It follows from these last estimates that functions (60) and (61) possess the properties of (2) for the regular values of spectral parameters ω and ν .  □

7. Continuous Dependence of Solution on the Small Parameter

We consider the continuous dependence of the solution to the problem (1)–(4) on small-parameter ε > 0 for regular values of spectral parameters ω and ν . Let ε 1 and ε 2 be two different values of small positive parameter ε . It is easy to check from (47)–(49) that the following estimates hold
max n 1 , , n m N max t [ 0 ; T ] W 1 i n 1 , , n m ( t , ε 1 , ω , ν ) W 1 i n 1 , , n m ( t , ε 2 , ω , ν ) C 21 ε 1 ε 2 , i = 1 , 4 ¯ ,
max n 1 , , n m N max t [ T ; 0 ] W 2 i n 1 , , n m ( t , ε 1 , ω , ν ) W 2 i n 1 , , n m ( t , ε 2 , ω , ν ) C 22 ε 1 ε 2 , i = 1 , 3 ¯ ,
where 0 < C 2 i = const , ε i 0 ; ε 0 , 0 < ε 0 = const , i = 1 , 2 .
Then, taking estimates (63), (64), (71) and (72) into account and applying the Cauchy–Schwartz inequality and Bessel inequality, from series (60) and (61), we obtain
U ( t , x , ε 1 , ω , ν ) U ( t , x , ε 2 , ω , ν n 1 , , n m = 1 u n 1 , , n m + ( t , ε 1 , ω , ν ) u n 1 , , n m + ( t , ε 2 , ω , ν ) · ϑ n 1 , , n m ( x ) 2 l m C 21 ε 1 ε 2 n 1 , , n m = 1 ψ 1 n 1 , , n m + ψ 2 n 1 , , n m + φ 1 n 1 , , n m + φ 2 n 1 , , n m γ 9 ε 1 ε 2 Ω l m 3 m ψ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x = ε 1 ε 2 · C 31 ,
where γ 9 = 2 l 3 m 2 C 21 C 03 l π 3 m , 0 < C 31 = const < ;
U ( t , x , ε 1 , ω , ν ) U ( t , x , ε 2 , ω , ν n 1 , , n m = 1 u n 1 , , n m ( t , ε 1 , ω , ν ) u n 1 , , n m ( t , ε 2 , ω , ν ) · ϑ n 1 , , n m ( x ) γ 10 ε 1 ε 2 Ω l m 3 m ψ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x + + Ω l m 3 m φ 1 ( x ) x 1 3 x 2 3 x m 3 2 d x + Ω l m 3 m φ 2 ( x ) x 1 3 x 2 3 x m 3 2 d x = ε 1 ε 2 · C 32 ,
where γ 10 = 2 l 3 m 2 C 22 C 03 l π 3 m , 0 < C 32 = const < .
It follows from estimates (73) and (74) that U ( t , x , ε 1 , ω , ν ) U ( t , x , ε 2 , ω , ν is small if ε 1 ε 2 is small in the domain Ω ¯ for the ( n 1 , , n m , ω , ν ) .

8. Conclusions and Statement of the Theorem

In the present paper, we study the questions of the one-value solvability of an inverse boundary value problem (1)–(6) for a mixed type integro-differential equation with Caputo operators of different fractional orders and spectral parameters in a multidimensional rectangular domain. For ( n 1 , , n m , ω , ν ) , we proved four lemmas under the following conditions A: Let functions
φ i ( x ) , ψ i ( x ) C 2 ( Ω l m ) , f i x , Ω l m Θ i ( y ) g i ( y ) d y C x 2 ( Ω l m × R ) , i = 1 , 2
in the domain Ω l m have piecewise continuous third order derivatives.
We will formulate a theorem as generalizing the above four proved lemmas. Thus, the following theorem is true.
Theorem 1.
Let the conditions ofAbe fulfilled. Then, for the possible numbers n 1 , , n m and regular values of spectral parameters ω and ν from the set ℵ, the inverse boundary value problem (1)–(6) is uniquely solvable in the domain Ω and the triple of solutions is represented in the form of series (57), (60) and (61). Moreover, it is true that
lim ε 0 U ( t , x , ε , ω , ν ) = U ( t , x , 0 , ω , ν ) ,
where U ( t , x , 0 , ω , ν ) is the solution of the mixed type fractional integro-differential equation of the form
A 0 ( U ) B ω ( U ) = ν 0 T K 1 ( t , s ) U ( s , x ) d s + F 1 ( t , x ) , t > 0 , ν T 0 K 2 ( t , s ) U ( s , x ) d s + F 2 ( t , x ) , t < 0 ,
A 0 ( U ) = 1 + sgn ( t ) 2 C D 0 t α 1 + 1 sgn ( t ) 2 C D 0 t α 2 U ( t , x ) , B ω ( U ) = i = 1 m U x i x i , t > 0 , ω 2 i = 1 m U x i x i , t < 0
with boundary value conditions (3)–(6) under consideration,
F i ( t , x ) = k i ( t ) g i ( x ) + f i x , Ω l m Θ i ( y ) g i ( y ) d y , i = 1 , 2 .
As a conclusion, we say that the numerical methods for solving fractional differential equations are important in the implementation of applied problems. In the future, we will also try to consider the applications of the numerical solution to the problems that we are solving. There are many methods for the numerical implementation of fractional differential equations. In this regard, we note the work done in [42]. In this paper, a new class of ( C , G f ) -invex functions is introduced and given nontrivial numerical examples, which justify the existence of such functions. Moreover, we construct generalized convexity definitions (such as, ( F , G f ) -invexity, C-convex etc.).

Author Contributions

Conceptualization, T.K.Y. and E.T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  2. 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]
  3. Area, I.; Batarfi, H.; Losada, J.; Nieto, J.J.; Shammakh, W.; Torres, A. On a fractional order Ebola epidemic model. Adv. Differ. Equ. 2015, 278. [Google Scholar] [CrossRef] [Green Version]
  4. Hussain, A.; Baleanu, D.; Adeel, M. Existence of solution and stability for the fractional order novel coronavirus (nCoV-2019) model. Adv. Differ. Equ. 2020, 384. [Google Scholar] [CrossRef] [PubMed]
  5. Ullah, S.; Khan, M.A.; Farooq, M.; Hammouch, Z.; Baleanu, D. A fractional model for the dynamics of tuberculosis infection using Caputo-Fabrizio derivative. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 975–993. [Google Scholar] [CrossRef] [Green Version]
  6. Tenreiro Machado, J.A. Handbook of Fractional Calculus with Applications in 8 Volumes; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2019. [Google Scholar]
  7. Kumar, D.; Baleanu, D. Fractional Calculus and Its Applications in Physics. Front. Phys. 2019, 7, 81. [Google Scholar] [CrossRef]
  8. 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]
  9. Saxena Ram, K.; Garra, R.; Orsingher, E. Analytical solution of space-time fractional telegraph-type equations involving Hilfer and Hadamard derivatives. Integral Transform. Spec. Funct. 2015, 27, 30–42. [Google Scholar] [CrossRef]
  10. Patnaik, S.; Hollkamp, J.P.; Semperlotti, F. Applications of variable-order fractional operators: A review. Proc. R. Soc. 2020, 476, 2234. [Google Scholar] [CrossRef] [Green Version]
  11. Garra, R.; Gorenflo, R.; Polito, F.; Tomovski, Ž. Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef] [Green Version]
  12. Tenreiro Machado, J.A. (Ed.) Handbook of Fractional Calculus with Applications; Walter de Gruyter GmbH: Berlin, Germany, 2019; Volume 8. [Google Scholar]
  13. Hilfer, R. Application of Fractional Calculus in Physics; World Scientific Publishing Company: Singapore, 2000. [Google Scholar]
  14. Hilfer, R. On fractional relaxation. Fractals 2003, 11, 251–257. [Google Scholar] [CrossRef]
  15. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  16. Klafter, J.; Lim, S.C.; Metzler, R. Fractional Dynamics, Recent Advances; World Scientific: Singapore, 2011; Chapter 9. [Google Scholar]
  17. Sandev, T.; Tomovski, Z. Fractional Equations and Models: Theory and Applications; Springer Nature Switzerland AG: Cham, Switzerland, 2019. [Google Scholar]
  18. Xu, C.; Yu, Y.; Chen, Y.Q.; Lu, Z. Forecast analysis of the epidemic trend of COVID-19 in the United States by a generalized fractional-order SEIR model. medRxiv 2020. [Google Scholar] [CrossRef]
  19. Cesarano, C. Generalized special functions in the description of fractional diffusive equations. Commun. Appl. Ind. Math. 2019, 10, 31–40. [Google Scholar] [CrossRef] [Green Version]
  20. Assante, D.; Cesarano, C.; Fornaro, C.; Vazquez, L. Higher Order and Fractional Diffusive Equations. J. Eng. Sci. Technol. Rev. 2015, 8, 202–204. [Google Scholar] [CrossRef]
  21. Dattoli, G.; Cesarano, C.; Ricci, P.; Vazquez, L. Special Polynomials and Fractional Calculus. Math. Comput. Model. 2003, 37, 729–733. [Google Scholar] [CrossRef]
  22. Restrepo, J.; Ruzhansky, M.; Suragan, D. Explicit representations of solutions for linear fractional differential equations with variable coefficients. arXiv 2020, arXiv:2006.1535v1. [Google Scholar]
  23. Gel’fand, I.M. Some questions of analysis and differential equations. Uspekhi Mat. Nauk. 1959, 14, 3–19. (In Russian) [Google Scholar]
  24. Uflyand, Y.S. On oscillation propagation in compound electric lines. Inzhenerno-Phizicheskiy Zhurnal 1964, 7, 89–92. (In Russian) [Google Scholar]
  25. Terlyga, O.; Bellout, H.; Bloom, F. A hyperbolic-parabolic system arising in pulse combustion: Existence of solutions for the linearized problem. Electron. J. Differ. Equ. 2013, 2013, 1–42. [Google Scholar]
  26. Abdullaev, O.K.; Sadarangani, K. Nonlocal problems with integral gluing condition for loaded mixed type equations involving the Caputo fractional derivative. Electron. J. Differ. Equ. 2016, 2016, 1–10. Available online: http://ejde.math.txstate.edu (accessed on 25 September 2020).
  27. Agarwal, P.; Berdyshev, A.S.; Karimov, E.T. Solvability of a nonlocal problem with integral transmitting condition for mixed type equation with Caputo fractional derivative. Results Math. 2017, 71, 1235–1257. [Google Scholar] [CrossRef]
  28. Zarubin, A.N. Boundary value problem for a differential-difference mixed-compound equation with fractional derivative and with functional delay and advance. Differ. Equ. 2019, 55, 220–230. [Google Scholar] [CrossRef]
  29. Karimov, E.T.; Al-Salti, N.; Kerbal, S. An inverse source non-local problem for a mixed type equation with a Caputo fractional differential operator. East Asian J. Appl. Math. 2017, 7, 417–438. [Google Scholar] [CrossRef]
  30. Karimov, E.T.; Kerbal, S.; Al-Salti, N. Inverse Source Problem for Multi-Term Fractional Mixed Type Equation. In Advanes in Real and Complex Analysis with Applications; Springer Nature Singapore Pte Ltd.: Singapore, 2017; pp. 289–301. [Google Scholar] [CrossRef]
  31. Repin, O.A. Nonlocal problem with Saigo operators for mixed type equation of the third order. Russ. Math. 2019, 63, 55–60. [Google Scholar] [CrossRef]
  32. Repin, O.A. On a problem for a mixed-type equation with fractional derivative. Russ. Math. 2018, 62, 38–42. [Google Scholar] [CrossRef]
  33. Salakhitdinov, M.S.; Karimov, E.T. Uniqueness of an inverse source non-local problem for fractional order mixed type equations. Eurasian Math. J. 2016, 7, 74–83. Available online: http://mi.mathnet.ru/rus/emj/v7/i1/p74 (accessed on 25 September 2020).
  34. 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]
  35. 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]
  36. Yuldashev, T.K. Nonlocal boundary value problem for a nonlinear Fredholm integro-differential equation with degenerate kernel. Differ. Equ. 2018, 54, 1646–1653. [Google Scholar] [CrossRef]
  37. Yuldashev, T.K. On the solvability of a boundary value problem for the ordinary Fredholm integrodifferential equation with a degenerate kernel. Comput. Math. Math. Phys. 2019, 59, 241–252. [Google Scholar] [CrossRef]
  38. Yuldashev, T.K. Spectral features of the solving of a Fredholm homogeneous integro-differential equation with integral conditions and reflecting deviation. Lobachevskii J. Math. 2019, 40, 2116–2123. [Google Scholar] [CrossRef]
  39. 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]
  40. Yuldashev, T.K. On an integro-differential equation of pseudoparabolic-pseudohyperbolic type with degenerate kernels. Proc. YSU Phys. Math. Sci. 2018, 52, 19–26. Available online: http://mi.mathnet.ru/rus/uzeru/v52/i1/p1914 (accessed on 25 September 2020). [CrossRef]
  41. Yuldashev, T.K. Nonlocal inverse problem for a pseudohyperbolic-pseudoelliptic type integro-differential equations. Axioms 2020, 9, 45. [Google Scholar] [CrossRef]
  42. Dubey, R.; Mishra, L.N.; Cesarano, C. Multiobjective fractional symmetric duality in mathematical programming with (C,Gf)-invexity assumptions. Axioms 2019, 8, 97. [Google Scholar] [CrossRef] [Green Version]
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.; Karimov, E.T. Inverse Problem for a Mixed Type Integro-Differential Equation with Fractional Order Caputo Operators and Spectral Parameters. Axioms 2020, 9, 121. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040121

AMA Style

Yuldashev TK, Karimov ET. Inverse Problem for a Mixed Type Integro-Differential Equation with Fractional Order Caputo Operators and Spectral Parameters. Axioms. 2020; 9(4):121. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040121

Chicago/Turabian Style

Yuldashev, Tursun K., and Erkinjon T. Karimov. 2020. "Inverse Problem for a Mixed Type Integro-Differential Equation with Fractional Order Caputo Operators and Spectral Parameters" Axioms 9, no. 4: 121. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop