Next Article in Journal
Frac-Vector: Better Category Representation
Next Article in Special Issue
Studying and Simulating the Fractional COVID-19 Model Using an Efficient Spectral Collocation Approach
Previous Article in Journal
Numerical Simulation of Nonlinear Stochastic Analysis for Measles Transmission: A Case Study of a Measles Epidemic in Pakistan
Previous Article in Special Issue
The Right Equivalent Integral Equation of Impulsive Caputo Fractional-Order System of Order ϵ∈(1,2)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvability of Mixed Problems for a Fourth-Order Equation with Involution and Fractional Derivative

by
Mokhtar Kirane
1,† and
Abdissalam A. Sarsenbi
2,*,†
1
Département de Mathématiques, Université de La Rochelle, 17042 La Rochelle, France
2
Department of Mathematics, M. Auezov South Kazakhstan University, Shymkent 160000, Kazakhstan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2023, 7(2), 131; https://doi.org/10.3390/fractalfract7020131
Submission received: 13 December 2022 / Revised: 19 January 2023 / Accepted: 26 January 2023 / Published: 30 January 2023
(This article belongs to the Special Issue Recent Advances in Time/Space-Fractional Evolution Equations)

Abstract

:
In the present work, two-dimensional mixed problems with the Caputo fractional order differential operator are studied using the Fourier method of separation of variables. The equation contains a linear transformation of involution in the second derivative. The considered problem generalizes some previous problems formulated for some fourth-order parabolic-type equations. The basic properties of the eigenfunctions of the corresponding spectral problems, when they are defined as the products of two systems of eigenfunctions, are studied. The existence and uniqueness of the solution to the formulated problem is proved.

1. Introduction and Problem Statement

The transformation S of the function f x , x 1 , 1 , which satisfies the condition S 2 f x = f x , is called an involution. We consider equations containing an involution of the form S f x = f x .
Differential equations with involution are of theoretical and practical interest, as they occupy a special place among differential equations with deviating arguments and functional differential equations. More information on the theoretical developments and applications of such equations can be found in [1,2,3,4,5].
The solvability of direct, inverse problems for the second-order differential equations with involution is considered in [4,5,6,7,8,9,10,11,12,13] (see also references therein). Most of these studies are based on the method of separation of variables, which leads to spectral problems with involution. Regarding spectral problems for differential equations with involution, good references are [14,15,16,17,18,19,20,21,22,23].
The boundary value problems for the fourth-order partial differential equations are studied in (see [24], and references therein). However, in the scientific literature, there are no works devoted to the study of the initial-boundary value problems for the fourth-order partial differential equations with involution.
This paper is devoted to the study of the existence and uniqueness of a solution to mixed problems for a two-dimensional fourth-order equation with involution
C D 0 t β u + 4 x 4 u x , y , t + α 4 x 4 u x , y , t + 4 y 4 u x , y , t = f x , y , t ,
where 1 < α < 1 . Here C D 0 t β u = I 0 t 1 β u t , β 0 , 1 is the Caputo derivative of order β [25] where for any integrable function g the left-hand fractional Riemann–Liouville integral of order β > 0 is defined as
I 0 β g t = 1 Γ β 0 t g τ t τ 1 β d τ ,
where Γ β is the Euler gamma function.
Let Q = x , y , t : 1 < x , y < 1 and Q ¯ be its closure. The space C x , y k , l Q ¯ consists of all functions u x , y , t having continuous derivatives of the order k with respect to x and continuous derivatives of the order l with respect to y in the domain Q ¯ , t > 0 .
The purpose of the work is to find, for given functions f x , y , t and φ x , y , a function that satisfies the following conditions:
(1) u x , y , t C x , y 3 , 3 Q ¯ C x , y 4 , 4 Q , t 0 ;
(2) u x , y , t satisfies Equation (1) in the domain Q, t > 0 ;
(3) u x , y , t satisfies the boundary conditions
u x , y , 0 = φ x , y , 1 x , y 1 ,
u 1 , y , t = 0 , u 1 , y , t = 0 , 2 x 2 u x = 1 = 2 x 2 u x = 1 = 0 , 1 y 1 , t 0 ,
l j u a j y u x , 1 , t + b j y u x , 1 , t + a j 1 u x , 1 , t + b j 1 u x , 1 , t = 0 , l ˜ j u a ˜ j 3 y 3 u x , 1 , t + b ˜ j 3 y 3 u x , 1 , t + a ˜ j 1 2 y 2 u x , 1 , t + b ˜ j 1 2 y 2 u x , 1 , t = 0 ,
where for j = 1 , 2 , coefficients a j , a j 1 , b j , b j 1 are some given numbers, x , t 1 , 1 × 0 , .
As can be seen from conditions (4), we consider boundary conditions of the most general type.
The paper consists of three sections and a list of references. Section 2 is devoted to the study of basic properties of the fourth-order differential operators. In Section 3, we give a solution to the stated problem.

2. Properties of Eigenfunctions of Spectral Problems

Consider a homogeneous equation
C D 0 t β u + 4 x 4 u x , y , t + α 4 x 4 u x , y , t + 4 y 4 u x , y , t = 0
with boundary conditions (2)–(4). For simplicity and without loss of generality, we seek for a nonzero solution of Equation (5) in the separated form u x , y , t = w x , y c t . Substituting this into Equation (5), we obtain the spectral problem
4 x 4 w x , y + α 4 x 4 w x , y + 4 y 4 w x , y = σ w x , y , 1 < x , y < 1 ,
w 1 , y = w 1 , y = 0 , 2 x 2 w x = 1 = 2 x 2 w x = 1 = 0 , 1 y 1 ,
l j w = 0 , j = 1 , 2 , l ˜ j w = 0 , j = 1 , 2 .
The solution to problem (6)–(8) is sought in the form w x , y = X x Y y . Then, the spectral problem (6)–(8) generates the following two one-dimensional spectral problems.
Spectral problem 1.
X x + α X x = λ X x , 1 < x < 1 ,
X 1 = X 1 = X 1 = X 1 = 0 .
Spectral problem 2.
Y y = μ Y y , 1 < y < 1 ,
l j Y a j Y 1 + b j Y 1 + a j 1 Y 1 + b j 1 Y 1 = 0 , j = 1 , 2 , l ˜ j θ Y a ˜ j Y 1 + b ˜ j Y 1 + a ˜ j 1 Y 1 + b ˜ j 1 Y 1 = 0 , j = 1 , 2 ,
where μ = σ λ .
Spectral problem 1. To study problems (9), (10), we use the following notation: α 0 = 1 1 + α 4 , α 1 = 1 1 α 4 , ρ = λ 4 , (here, r e i φ 1 4 = r 1 4 e i φ 4 , π 2 < φ 3 π 2 ).
The solution to Equation (9) can be written as
X x = c 1 e α 0 ρ x + e α 0 ρ x + c 2 e i α 0 ρ x + e i α 0 ρ x + c 3 e α 1 ρ x e α 1 ρ x + c 4 e i α 1 ρ x e i α 1 ρ x .
Using boundary conditions (10), we find the eigenvalues of the first series
λ k 1 = 1 + α k 1 2 4 π 4 , k = 1 , 2 , ,
whose corresponding eigenfunctions are
X k 1 = cos k 1 2 π x .
The second series of eigenvalues
λ k 2 = 1 α k π 4 , k = 1 , 2 , 3 , ,
whose corresponding eigenfunctions are X k 2 = sin k π x .
The system of eigenfunctions of spectral problems (9) and (10) can be written as
X k x = X 2 k x = cos k 1 2 π x , X 2 k 1 x = sin k π x .
Lemma 1. 
System (13) forms a complete orthonormal system in L 2 1 , 1 .
The lemma is proved by the direct calculation of the Fourier coefficients of an arbitrary function f x L 2 1 , 1 .
Spectral problem 2. To study the fourth-order spectral problem (11), (12), we introduce the operator L Y y = Y θ y , where θ = ± 1 , with the domain D ( L ) = { Y y W 2 2 1 , 1 : l j Y = 0 } , where W 2 2 1 , 1 is the Sobolev space [26]. Hereafter, the operator L acts on functions depending on the variable y. The operator L commutes with the operator S, i.e., L S = S L , and the equality L 2 Y y = Y y is valid if L Y y D L . The domain of the operator L 2 is
D L 2 = Y y W 4 2 1 , 1 : l j Y = 0 , l ˜ j θ Y = 0 .
In this case, the coefficients of the expression l ˜ j Y in (12) depend on the number θ . Consider the spectral problem for the second-order equation
Y θ y = ν Y y , 1 < y < 1 , θ = ± 1 ,
l j Y = 0 , j = 1 , 2 ,
where (15) is the first equality in (12). Equation (14) is self-conjugate, whereas boundary conditions (15) may be non-self-conjugate. The conjugate boundary conditions are denoted by l j * Y . For θ = 1 , we have a classical spectral problem for an ordinary second-order differential operator. For θ = 1 , problem (14), (15) was studied in [27,28]. In the case θ = 1 , the result is well known [29].
Theorem 1 
([29]). Let θ = 1 . Then, the system of eigenfunctions of the spectral problem (14), (15) forms a Riesz basis in L 2 1 , 1 if one of the following three conditions A:
(1) 
a 1 b 2 a 2 b 1 0 ;
(2) 
a 1 = a 2 = b 1 = b 2 = 0 , a 11 b 21 a 21 b 11 0 ;
(3) 
a 1 b 2 a 2 b 1 = 0 , a 1 + b 1 > 0 , a 1 2 b 1 2 a 21 2 b 21 2 0 ;
 is satisfied.
In the case θ = 1 , properties of eigenfunctions in problems (14), (15) were studied in [27,28]. Based on the results of [28], the following theorem is formulated.
Theorem 2 
([28]). Let θ = 1 . Then, the system of eigenfunctions of spectral problems (14) and (15) forms a Riesz basis in L 2 1 , 1 if one of the three conditions A is satisfied. If a 1 b 2 a 1 b 2 = 0 , a 1 + b 1 > 0 , and one of the following four conditions:
( 1 ) a 1 = b 1 , a 21 + b 21 0 ; ( 2 ) a 1 = b 1 , a 21 b 21 0 ; ( 3 ) a 21 = b 21 , a 1 + b 1 0 ; ( 4 ) a 21 = b 21 , a 1 b 1 0 ,
is satisfied, then the system of eigenfunctions of the spectral problem (14), (15) also forms a Riesz basis in L 2 1 , 1 .
It follows from the results of [28] that the system of eigenfunctions of problem (14) does not form a basis in L 2 1 , 1 for boundary conditions of the type Y 1 = ± Y 1 ,   Y 1 = Y 1 .
It was shown in [27] that if L Y y D L , then the system of eigenfunctions of the second-order spectral problems (14) and (15) coincides with the system of eigenfunctions of the fourth-order spectral problems (11) and (12). Therefore, Theorems 1 and 2 are also valid for the fourth-order spectral problems (11) and (12).
Remark 1. 
If θ = 1 , then the expressions l ˜ j θ Y in (12) can be represented in the form
l ˜ j θ Y a j Y 1 + b j Y 1 + a j 1 Y 1 + b j 1 Y 1 = 0 , j = 1 , 2 .
According to this, the second boundary condition in (4) takes the form
l ˜ j u a j 3 y 3 u x , 1 , t + b j 3 y 3 u x , 1 , t + a j 1 2 y 2 u x , 1 , t + b j 1 2 y 2 u x , 1 , t = 0 .
Remark 2. 
If θ = 1 , then the expressions l ˜ j θ Y in (12) can be represented in the form
l ˜ j θ Y b j Y 1 + a j Y 1 b j 1 Y 1 a j 1 Y 1 = 0 , j = 1 , 2 .
According to this, the second boundary condition in (4) takes the form
l ˜ j u b j 3 y 3 u x , 1 , t + a j 3 y 3 u x , 1 , t b j 1 2 y 2 u x , 1 , t a j 1 2 y 2 u x , 1 , t = 0 .
Let us formulate theorems on the basis property of the eigenfunctions of the spectral problems (11) and (12) for a fourth-order equation.
Theorem 3. 
Let θ = 1 and L Y y D L . Then, the system of eigenfunctions of the spectral problems (11) and (12) forms the Riesz basis in L 2 1 , 1 if the coefficients of the boundary conditions (12) satisfy one of the three conditions A.
Theorem 4. 
Let θ = 1 and L Y y D L . Then, the system of eigenfunctions of the spectral problems (11) and (12) forms the Riesz basis in the following cases:
(1) if the coefficients of the boundary conditions (12) satisfy one of the three conditions A;
or
(2) if a 1 b 2 a 2 b 1 = 0 , a 1 + b 1 > 0 , and for the coefficients of the boundary conditions (12), one of the four conditions in (16) is satisfied.
Note that when one of the three conditions A is satisfied, all eigenvalues of the spectral problem for the second-order Equations (14) and (15) are single. If one of the conditions (16) is satisfied, the boundary conditions (12) for the fourth-order Equation (11) are regular in the sense of Birkhoff [30] (but are not strongly regular), the eigenvalues of the fourth order spectral problem (11), (12) may be multiple. For example [27], the boundary value problem (14) with boundary conditions
Y 1 = Y 1 = 0 , Y 1 = Y 1 , Y 1 = Y 1
( a 1 = b 1 , a 21 = 1 , all other coefficients are equal to zero) is not self-conjugate, has an infinite number of double eigenvalues, has no associated functions, and the system of eigenfunctions forms a Riesz basis in L 2 1 , 1 . Hence, it appears that the system of eigenfunctions of the conjugate problem
Y * 1 = 0 , Y * 1 = 0 , Y * 1 = Y * 1 , Y * 1 = Y * 1 ,
also forms the Riesz basis in L 2 1 , 1 .
Thus, in the case θ = 1 , we cover a larger range of boundary conditions for the boundary value problem (11), (12).
All the systems we consider are Riesz bases, and they remain Riesz bases for any enumeration of system elements (unconditional bases). By virtue of [31], the systems we consider are uniformly bounded. Therefore, the coefficients of bi-orthonormalization, in the case of the non-orthogonality of the studied systems, are also uniformly bounded. This means that the bi-orthnormalization coefficients do not affect either the convergence of expansions in terms of eigenfunctions or the convergence of series formed by the expansion coefficients.
Note that boundary conditions (10) are self-conjugate, whereas boundary conditions (12) may be non-self-conjugate. Therefore, problem (6)–(8) in the general case is a non-self-conjugate problem. For specific boundary conditions of the form (8), the conjugate boundary conditions are written in the standard way.
Let us write a problem conjugate to problem (6)–(8):
4 x 4 w * x , y + α 4 x 4 w * x , y + 4 y 4 w * x , y = σ ¯ w * x , y , 1 < x , y < 1 ,
w * 1 , y = w * 1 , y = 0 , 2 w * x 2 x = 1 = 2 w * x 2 x = 1 = 0 , 1 y 1 ,
a j * y w * x , 1 + b j * y w * x , 1 + a j 1 * w * x , 1 + b j 1 * w * x , 1 = 0 , j = 1 , 2 , a ˜ j * 3 y 3 w * x , 1 + b ˜ j * 3 y 3 w * x , 1 + a ˜ j 1 * 2 y 2 w * x , 1 + b ˜ j 1 * 2 y 2 w * x , 1 = 0 , j = 1 , 2 .
Taking the above into account, we assume that all eigenvalues of problems (9)–(12) are numbered as follows: λ 2 n = 1 + α n 1 2 4 π 4 , λ 2 n 1 = 1 α n π 4 , n = 1 , 2 , , μ k μ k + 1 , for any number k, σ n k = λ n + μ k . If an eigenvalue μ k is double, then μ k + 1 = μ k . The eigenfunctions wiare numbered respectively.
We denote the system of eigenfunctions of problem (6)–(8) as
w k n x , y = X k x Y n y ,
where X k x are eigenfunctions of problem (9), (10), and Y n y are eigenfunctions of problem (11), (12). We denote the system of eigenfunctions of the corresponding conjugate problem (19)–(21) as
w * k n x , y = X k x Y n * y ,
where Y n * y are eigenfunctions of the problem conjugate to problems (11) and (12).
Systems (22) and (23) are bi-orthogonal as systems of eigenfunctions of mutually conjugate spectral problems (6)–(8) and (19)–(21). Since the system X k x is an ortho-normal basis and the system Y n y is a Riesz basis, then each of the systems (22) and (23) is a Riesz basis [32] in L 2 1 , 1 × 1 , 1 .

3. Existence and Uniqueness of the Problem Solution

In this section, we assume that all the conditions of Theorem 4 are satisfied, and the second boundary condition in (4) is of the form (17) or (18).
Then, each of the systems (22) and (23) forms a Riesz basis in L 2 1 , 1 × 1 , 1 . To solve the inhomogeneous problem (1)–(4), we first consider the homogeneous Equation (5). The solution to problems (5), (2)–(4) is sought in the form of a Fourier series
u x , y , t = k , n = 0 C k n t w k n x , y
using the Riesz basis w k n x , y , where the unknown functions C k n t are defined by the relations
C k n t = Ω x y u x , y , t w * k n x , y d x d y = 1 1 1 1 u x , y , t X k x Y n * y d x d y .
From Equations (5) and (25) and condition (2), we obtain the problems
C D 0 t β C k n t + σ n k C k n t = 0 , C k n 0 = φ k n ,
whose solutions are written as [25]
C k n t = φ k n E β σ k n t β ,
where φ k n = φ x , y , w * k n x , y in the sense of the scalar product in L 2 1 , 1 × 1 , 1 . Here, E β z is the Mittag–Leffler function [25] satisfying the estimate
E β z M 1 + z .
Using relation (26), expression (24) can be written as
u x , y , t = k , n = 1 φ k n E β σ k n t β X k x Y n y .
It is necessary to show that the formal solution (27) is a classical solution to problems (5), (2)–(4).
Lemma 2. 
Let φ x , y L 2 Q ¯ . Then, series (27) converges absolutely and uniformly in the domain Q ¯ , t > 0 .
Proof. 
For any t t 0 > 0 , the Mittag–Leffler function satisfies the estimate
E β σ k n t β M 1 + σ k n t β .
Since σ n k = λ n + μ k , λ k 1 = 1 + α k 1 2 4 π 4 , λ k 2 = 1 α k π 4 , μ k = O k 4 , then
M 1 + σ k n t β M σ k n t β M σ k n t 0 β M t 0 β λ k + μ k M t 0 β 2 λ k μ k .
Then
u x , y , t = k , n = 1 φ k n E β σ k n t β X k x Y n y k , n = 1 1 1 1 1 φ x , y X k x Y n * y d x d y E β σ k n t β X k x Y n y k , n = 1 M 0 M t 0 β 2 λ k μ k 1 1 1 1 φ x , y X k x Y n * y d x d y .
Applying the Cauchy–Schwarz inequality and the Bessel inequality, we obtain
k , n = 0 φ k n M 1 φ L 2 Q .
The lemma is proved. □
Note that the method of separation of variables imposes stricter requirements on the initial function. However, in this work, the problem of reducing the smoothness of the initial function was not posed.
The main result of the work is the following theorem.
Theorem 5. 
Let the function φ x , y be such that φ x , y C x , y 3 , 0 Q ¯ C x , y 0 , 3 Q ¯ C x , y 4 , 0 Q C x , y 0 , 4 Q C x , y 2 , 2 Q ¯ , L φ x , y D L , L 2 φ x , y D L , t > 0 and satisfy self-conjugate conditions (3) with respect to the variable x. Then, problem (5), (2)–(4) has a unique solution, which can be represented in the form (27).
Proof. 
Let us first prove the uniqueness. If the problem has two solutions u 1 , u 2 , then the function u = u 1 u 2 satisfies u x , y , 0 = 0 , Equation (5), and boundary conditions (3), (4). Then, from (26), we obtain C k n t = 0 , and from (25) it follows that u x , y , t = 0 in Q ¯ due to the completeness of the system w * k n x , y in L 2 1 , 1 × 1 , 1 . By virtue of Lemma 2, from the continuity of the solution in the entire closed domain, we obtain u x , y , t = 0 in Q ¯ . The uniqueness of the solution is proved.
Let us prove the existence of a solution. Formula (27) represents a formal solution to problems (5), (2)–(4). In order to justify that the obtained formal solution is indeed a true solution, we need to show that the series C D 0 t β u , p u x p , p u y p ,   p = 1 , 2 , 3 and 4 u x 4 , 4 u y 4 converge uniformly in Q ¯ and Q, t > 0 , respectively.
It is obvious that for the elements of system (13), the estimates
d p X 2 k d x p k 1 2 π p ,
d p X 2 k 1 d x p k π p , p = 1 , 2 , 3 , 4
are valid. Now, we must prove the uniform convergence of the series 4 u x , y , t x 4 in Q, t > 0 . Formally differentiating series (27) four times with respect to x and integrating by parts four times with respect to the variable x, taking into account (13), we get
4 x 4 u x , y , t = k , n = 1 φ k n E β σ k n t β X k 4 x Y n y
= k , n = 1 1 1 1 1 φ x , y X k x Y n * y d x d y E β σ k n t β X k 4 x Y n y
= k , n = 1 1 1 1 1 φ x , y X k x Y n * y d x d y E β σ k n t β λ k X k x Y n y
α 2 k , n = 1 1 1 1 1 φ x 4 x , y X k x λ k Y n * y d x d y E β σ k n t β λ k X k x Y n y ,
where α 2 = m a x α 0 1 , α 1 1 . Further, applying the Lemma 2 we get
4 u x , y , t x 4 M 0 φ x 4 L 2 Q .
where φ x 4 means the fourth derivative with respect to x. Thus, the uniform convergence of the series 4 u x , y , t x 4 in Q is proved.
The convergence of the series 4 u x , y , t y 4 is proved similarly by replacing the fourth-order derivative Y 4 y with the help of Equation (11). In this case, the function inside the integral Y n * y is represented as its fourth derivative using Equation (11). Here are the corresponding calculations:
4 y 4 u x , y , t = k , n = 1 φ k n E β σ k n t β X k x Y n 4 y
= k , n = 1 1 1 1 1 φ x , y X k x Y n * y d x d y E β σ k n t β X k x μ n Y n y
= k , n = 1 1 1 1 1 φ x , y X k x Y n * 4 y μ n d x d y E β σ k n t β X k x μ n Y n y
= k , n = 1 1 1 1 1 φ y 4 x , y X k x Y n * y μ n d x d y E β σ k n t β X k x μ n Y n y .
Hence, as in the previous case, we obtain the estimate
4 u x , y , t y 4 M 0 φ y 4 L 2 Q .
Thus, the uniform convergence of the series 4 u x , y , t y 4 in Q , t > 0 is proved. The uniform convergence of the series 2 u x , y , t y 2 in Q ¯ , t > 0 is proved in a similar way.
Using estimates [33]
Y n s y c o n s t 1 + μ n 4 s max Y n y , s = 1 , 2 , 3 ,
for the eigenfunctions of the fourth-order differential operator, as above, we prove the uniform convergence of the series p u y p , p = 1 , 3 , in Q ¯ .
Each term of series (27) is a solution of problems (5), (2)–(4). Therefore, boundary conditions (3) and (4) are satisfied. The function (27) satisfies the relations of the type (2). To do this, it sufficient to show that the series (27) converges uniformly in Q ¯ , t 0 . For t 0 , the inequality
u x , y , t M k , n = 1 1 1 1 1 φ x , y X k x Y n * y d x d y X k x Y n y
is satisfied. The function inside the integral Y n * y is represented as its second derivative using Equation (14). Integrating by parts two times over the variable x and two times over the variable y, applying the Cauchy–Schwarz inequality and the Bessel inequality, we conclude that the series (27) converges absolutely and uniformly in Q ¯ , t 0 . The theorem is proved. □
Remark 3. 
The solution to problem (2)–(4) for the inhomogeneous Equation (1) is sought in the form of a series (24). We represent the function f x , y , t as an expansion in the basis (22)
f x , y , t = k , n = 0 f k n t w k n x , y , f k n t = f x , y , t , w * k n x , y .
The functions C k n t have the form
C k n t = φ k n E β σ k n t β + 0 t t τ β 1 E β , β σ k n t τ β f k n τ d τ ,
where φ k n = φ x , y , w * k n x , y in the sense of the scalar product in L 2 1 , 1 × 1 , 1 , E β , γ z is the Mittag–Leffler function [25] defined by formula E β , γ z = k = 0 z k Γ β k + γ and satisfying the estimate
E β , γ z M 1 + z , 0 < β < 2 , γ R , δ arg z π , δ β π 2 , min π ; β π .
Substituting (28) into (24), we obtain the solution of problem (1)–(4).
Example 1. 
Let the boundary conditions (4) have the form (case θ = 1 )
l 1 u y u x , 1 , t + y u x , 1 , t = 0 , l 2 u u x , 1 , t = 0 , l ˜ 1 u 3 y 3 u x , 1 , t + 3 y 3 u x , 1 , t = 0 , l ˜ 2 u 2 y 2 u x , 1 , t = 0 .
Then, the spectral problem (11), (12) takes the form
Y y = μ Y y , 1 < y < 1 ,
Y 1 + Y 1 = 0 , Y 1 = 0 , Y 1 + Y 1 = 0 , Y 1 = 0 .
where the eigenvalues μ k = k 1 2 4 π 4 correspond to two series of eigenfunctions
Y k 1 = sin k 1 2 π y + 1 k + 1 e k 1 2 π y + e k 1 2 π y e k 1 2 π + e k 1 2 π , Y k 2 = cos k 1 2 π y ,
k = 1 , 2 , The conjugate boundary conditions are
Y * 1 = 0 , Y * 1 + Y * 1 = 0 , Y * 1 = 0 , Y * 1 + Y * 1 = 0 .
The eigenfunctions of the conjugate problem are
Y k 1 * = sin k 1 2 π y , Y k 2 * = cos k 1 2 π y + 1 k + 1 e k 1 2 π y e k 1 2 π y e k 1 2 π + e k 1 2 π .
Each of systems (30) and (31) forms a Riesz basis in L 2 1 , 1 [27].
Theorem 5 is formulated in the following form.
Theorem 6. 
Let the function φ x , y be such that φ x , y C x , y 3 , 0 Q ¯ C x , y 0 , 3 Q ¯ C x , y 4 , 0 Q C x , y 0 , 4 Q C x , y 2 , 2 Q ¯ , satisfy self-conjugate conditions (3) with respect to the variable x and conditions (29) with respect to the variable y. Then problem (5), (2), (3), (29) has a unique solution, which can be represented in the form (27).
Example 2. 
Considering the problem (5), (2), (3), (29). Let in (2) φ x , y = sin π x cos π 2 y . Then in (27) coefficient φ 11 = 1 , all other coefficients φ k n , k , n 1 , equal to zero. Since λ 1 = 1 α π 4 , μ 1 = 1 16 π 4 , σ 11 = 17 16 α π 4 , from Equation (26), we get C 11 = E β 17 16 α π 4 t β . The solution in the form (27) of problem (5), (2), (3), (29) takes the form
u x , y , t = E β 17 16 α π 4 t β sin π x cos π 2 y .

4. Conclusions

The existence and uniqueness of regular solutions of a certain class of problems for two-dimensional parabolic equations of fractional order is established. With respect to one variable, the equation contains an involution transformation. The completeness of the eigenfunctions of the boundary value problem for the fourth-order differential equation with involution is proved. The Riesz basis property of the eigenfunctions of boundary value problems for the fourth-order ordinary differential equation is shown for a special class of boundary conditions. In this case, the properties of eigenfunctions of boundary value problems for the second-order differential equation with involution are used.

Author Contributions

M.K. and A.A.S. contributed equally to the conceptualization and the formal analysis of the problem discussed in the paper. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by a grant from the Ministry of Science and Education of the Republic of Kazakhstan (grant No. AP13068539).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data is present within the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wiener, J. Generalized Solutions of Functional Differential Equations, 1st ed.; World Scientific: Singapore; River Edge, NJ, USA; London, UK; Hong Kong, China, 1993; ISBN 981-02-1207-0. [Google Scholar]
  2. Przeworska-Rolewicz, D. Equations with Transformed Argument. An Algebraic Approach, 1st ed.; Elsevier-PWN: Amsterdam, The Netherlands; Warszawa, Poland, 1973. [Google Scholar]
  3. Cabada, A.; Tojo, F.A.F. Differential Equations with Involutions, 1st ed.; Atlantis Press: Paris, France, 2015; ISBN 978-94-6239-120-8. [Google Scholar]
  4. 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]
  5. Mussirepova, E.; Sarsenbi, A.A.; Sarsenbi, A.M. Solvability of mixed problems for the wave equation with reflection of the argument. Math. Methods Appl. Sci. 2022, 45, 11262–11271. [Google Scholar] [CrossRef]
  6. Al-Salti, N.; Kerbal, S.; Kirane, M. Initial-boundary value problems for a time-fractional differential equation with involution perturbation. Math. Model. Nat. Phenom. 2019, 14, 315. [Google Scholar] [CrossRef] [Green Version]
  7. Ashyralyev, A.; Sarsenbi, A.M. Well-posedness of a parabolic equation with involution. Num. Funct. Anal. Optim. 2017, 38, 1295–1304. [Google Scholar] [CrossRef]
  8. Karachik, V.V.; Sarsenbi, A.M.; Turmetov, B.K. On the solvability of the main boundary value problems for a nonlocal Poisson equation. Turk. J. Math. 2019, 43, 1604–1625. [Google Scholar] [CrossRef]
  9. Ilyas, A.; Malik, S.A.; Saif, S. Inverse problems for a multi-term time fractional evolution equation with an involution. Inverse Probl. Sci. Eng. 2021, 29, 3377–3405. [Google Scholar] [CrossRef]
  10. Yarka, U.; Fedushko, S.; Vesely, P. The Dirichlet problem for the perturbed elliptic equation. Mathematics 2020, 8, 2108. [Google Scholar] [CrossRef]
  11. Sarsenbi, A.A. A solvability conditions of mixed problems for equations of parabolic type with involution. Bull. Karaganda Univ.-Math. 2018, 4, 87–93. [Google Scholar] [CrossRef] [Green Version]
  12. Usmanov, K.; Turmetov, B.; Nazarova, K. On the Solvability of Some Boundary Value Problems for the Nonlocal Poisson Equation with Boundary Operators of Fractional Order. Fractal Fract. 2022, 6, 308. [Google Scholar] [CrossRef]
  13. Kozhanov, A.I.; Bzheumikhova, O.I. Elliptic and Parabolic Equations with Involution and Degeneration at Higher Derivatives. Mathematics 2022, 10, 3325. [Google Scholar] [CrossRef]
  14. Burlutskaya, M.S.; Khromov, A.P. Fourier method in an initial-boundary value problem for a first-order partial differential equation with involution. Comput. Math. Math. Phys. 2011, 51, 2102–2114. [Google Scholar] [CrossRef]
  15. Baskakov, A.G.; Krishtal, I.A.; Uskova, N.B. Spectral Properties of Classical Dirac Operators and Operators with Involution in Homogeneous Function Spaces. Diff. Equ. 2021, 57, 1273–1278. [Google Scholar] [CrossRef]
  16. Kritskov, L.V.; Sarsenbi, A.M. Basicity in Lp of root functions for differential equations with involution. Electron. J. Differ. Eq. 2015, 278, 1–9. [Google Scholar]
  17. Kritskov, L.V.; Sadybekov, M.A.; Sarsenbi, A.M. Properties in Lp of root functions for a nonlocal problem with involution. Turk. J. Math. 2019, 43, 393–401. [Google Scholar] [CrossRef]
  18. Kritskov, L.V.; Sadybekov, M.A.; Sarsenbi, A.M. Nonlocal spectral problem for a second-order differential equation with an involution. Bull. Karaganda Univ.-Math. 2018, 3, 53–60. [Google Scholar] [CrossRef]
  19. Sarsenbi, A.A.; Sarsenbi, A.M. On Eigenfunctions of the boundary value problems for second order differential equations with involution. Symmetry 2021, 13, 1972. [Google Scholar] [CrossRef]
  20. Kalenyuk, P.I.; Baranetskij, Y.O.; Kolyasa, L.I. A nonlocal problem for a differential operator of even order with involution. J. Appl. Anal. 2020, 26, 297–307. [Google Scholar] [CrossRef]
  21. Bondarenko, N.P. Inverse spectral problems for functional-differential operators with involution. J. Differ. Equ. 2022, 318, 169–186. [Google Scholar] [CrossRef]
  22. Polyakov, D.M. Spectral Asymptotics of Two-Term Even Order Operators with Involution. J. Math. Sci. 2022, 260, 810–819. [Google Scholar] [CrossRef]
  23. Sarsenbi, A.A. Unconditional basicity of eigenfunctions’ system of Sturm-Liouville operator with an involutional perturbation. Bull. Karaganda Univ.-Math. 2018, 3, 117–127. [Google Scholar] [CrossRef]
  24. Kerbal, S.; Kadirkulov, B.J.; Kirane, M. Direct and Inverse Problems for a Samarskii-Ionkin Type Problem for a Two Dimensional Fractional Parabolic Equation. Progr. Fract. Differ. Appl. 2018, 4, 147–160. [Google Scholar] [CrossRef]
  25. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Application of Fractional Differential Equations, 1st ed.; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006; ISBN 0444518320/9780444518323. [Google Scholar]
  26. Adams, R.A. Sobolev Spaces; Academic Press: New York, NY, USA; San Francisco, CA, USA; London, UK, 1975. [Google Scholar]
  27. Sarsenbi, A.M. Unconditional bases related to a nonclassical second-order differential operator. Differ. Equ. 2010, 46, 509–514. [Google Scholar] [CrossRef]
  28. Sadybekov, M.A.; Sarsenbi, A.M. Criterion for the basis property of the eigenfunction system of a multiple differentiation operator with an involution. Differ. Equ. 2012, 48, 1112–1118. [Google Scholar] [CrossRef]
  29. Keselman, G.M. On non-conditional convergence of expansions by eigenfunctions of some differential operators. Izv. Vuzov. Matem. 1964, 39, 82–93. [Google Scholar]
  30. Naimark, M.A. Linear Differential Operators; Ungar: New York, NY, USA, 1968. [Google Scholar]
  31. Sarsenbi, A.M. Necessary and sufficient conditions for the Riesz basis property of the eigen- and associated functions of high-order differential operators on an interval. Dokl. Math. 2008, 77, 290–292. [Google Scholar] [CrossRef]
  32. Il’in, V.A.; Barnovska, M. Riesz basis of a spectral problem with an infinite-to-one eigenvalues. Math. Slovaca 1985, 35, 165–167. [Google Scholar]
  33. Kerimov, N.B. On the basis property and uniform minimality of systems of root functions of differential operators, I. Differ. Uravn. 1996, 32, 321–327. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kirane, M.; Sarsenbi, A.A. Solvability of Mixed Problems for a Fourth-Order Equation with Involution and Fractional Derivative. Fractal Fract. 2023, 7, 131. https://doi.org/10.3390/fractalfract7020131

AMA Style

Kirane M, Sarsenbi AA. Solvability of Mixed Problems for a Fourth-Order Equation with Involution and Fractional Derivative. Fractal and Fractional. 2023; 7(2):131. https://doi.org/10.3390/fractalfract7020131

Chicago/Turabian Style

Kirane, Mokhtar, and Abdissalam A. Sarsenbi. 2023. "Solvability of Mixed Problems for a Fourth-Order Equation with Involution and Fractional Derivative" Fractal and Fractional 7, no. 2: 131. https://doi.org/10.3390/fractalfract7020131

APA Style

Kirane, M., & Sarsenbi, A. A. (2023). Solvability of Mixed Problems for a Fourth-Order Equation with Involution and Fractional Derivative. Fractal and Fractional, 7(2), 131. https://doi.org/10.3390/fractalfract7020131

Article Metrics

Back to TopTop