CC BY
32
MSC 35-02, 34G20, 35C09, 35Q35
DOI: 10.14529/mmp220103
SEMILINEAR SOBOLEV TYPE MATHEMATICAL MODELS
A.A. Zamyshlyaeva1, E.V. Bychkov1
1 South Ural State University, Chelyabinsk, Russian Federation E-mail: zamyshliaevaaa@susu.ru, bychkovev@susu.ru
The article contains a review of the results obtained in the scientific school of Georgy Sviridyuk in the field of semilinear Sobolev type mathematical models. The paper presents results on solvability of the Cauchy and Showalter-Sidorov problems for semilinear Sobolev type equations of the first, the second and higher orders, as well as examples of non-classical models of mathematical physics, such as the generalized Oskolkov model of nonlinear filtering, propagation of ion-acoustic waves in plasma, propagation waves in shallow water, which are studied by reduction to one of the above abstract problems. Methods for studying the semilinear Sobolev type equations are based on the theory of relatively p-bounded operators for equations of the first order and the theory of relatively polynomially bounded operator pencils for equations of the second and higher orders in the variable t. The paper uses the phase space method, which consists in reducing a singular equation to a regular one defined on some subspace of the original space, to prove existence and uniqueness theorems, and the Galerkin method to construct an approximate solution.
Keywords: Oskolkov equation; equation of ion-acoustic waves in plasma; modified Boussinesq equation; semilinear Sobolev type equation; relatively p-bounded operators; relatively polynomially bounded operator pencils; Galerkin method; *-weak convergence.
Dedicated to the 70th anniversary of the Teacher Professor Georgy Anatolyevich Sviridyuk
Introduction
Mathematical models based on semilinear Sobolev type equations are called semilinear Sobolev type mathematical models.
The necessity to study semilinear Sobolev type mathematical models caused by the need to study applied problems related to the dynamics of viscoelastic fluids [16], vibrations in the DNA molecule [6], the theory of metal creep [5], wave propagation in shallow water [25], the propagation of ion-acoustic waves in plasma [1], the theory of electric circuits [14], the theory of heat conduction with two temperatures [4], filtration in a fractured-porous medium [2,28] and others [19,39].
The study of semilinear Sobolev type equations was initiated in [29,37] and the concept of a quasistationary trajectory was introduced. Later, on the basis of abstract results, the Oskolkov [33,35] and Hoff [32] mathematical models were investigated, and the structure of the phase space of the Hoff and Oskolkov equations was studied. In parallel with the solution of the problem of existence and uniqueness, the theory of optimal control of solutions to the semilinear Sobolev type equations arose and was developed [23,36,44], the theory of stability of solutions to the semilinear Sobolev type equations [22,34] was studied and the phenomenon of non-uniqueness of solutions was explained [18, 31].
Probably, the first work devoted to the study of equations unsolvable with respect to the highest time derivative belongs to A. Poincare (1885) [20]. However, the regular study of initial-boundary value problems for such equations began with the works of
S.L. Sobolev [27]. At present, the theory of Sobolev type equations is actively developing both in breadth and in depth, as evidenced by many scientific directions around which the scientific schools have developed [1,7,10,14,15,21,24,26]. In this paper, we trace the stages in the development of the theory of semilinear Sobolev type equations and consider its applications to the study of three semilinear Sobolev type mathematical models.
The first one is the Oskolkov mathematical model (a mathematical model of nonlinear filtration in a fractured-porous medium). Let Q C Rn be a bounded domain with boundary dQ of class CIn a cylinder Q x R, consider the equation
ut — ^Aut = v A« — K (u) (1)
with the Cauchy-Dirichlet conditions
u(x, 0) = u0(x), x G Q, (2)
u(x,t) = 0, (x,t) G dQ x R. (3)
Equation (1) describes many processes and phenomena occurring with the participation of a viscoelastic fluid, including filtration. It was obtained by A.P. Oskolkov [19]. The nonlinear term in (1) is such that K(0) = 0, (K(u),u) > 0 ((-, •) is an inner product in L2(Q)). In particular, it can take the form K(u) = u2m+1 or K(u) = shu. Generally speaking, the nonlinearity can be represented by the series
<x
K(u) = a™u2m+\ am G R+.
m=0
The parameters v G R+ characterize the elastic and viscous properties of the fluid, respectively.
The second one is a mathematical model of ion-acoustic waves in plasma. Let Q = (0,a) x (0,b) x (0, c) C R3. In a cylinder Q x R consider the equation
d2u
(A - A)utttt + (A - A')utt + Q!^ = A(m3) (4)
5x3
with the Cauchy-Dirichlet conditions
u(x, 0) = u0(x), ut(x, 0) = u1(x), uti(x, 0) = u2(x), uttt(x, 0) = u3(x), x G Q, (5)
u(x,t) = 0, (x,t) G dQ x R.
Equation (4) describes the ion-acoustic waves in a plasma in an external magnetic field. Here function u is a generalized potential of the electric field, the constants A, A', a characterize ion gyrofrequency, Langmuir frequency and Debye radius.
The third one is a mathematical model of wave propagation in shallow water. It is based on a modified Boussinesq equation. Let Q C Rn be a domain with boundary dQ of class CT G R+. In a cylinder C = Q x (0, T), consider the modified Boussinesq equation
(A — A)utt — a2Au + u3 = 0, (x,t) G Q x (0,T) (6)
with the Cauchy-Dirichlet conditions
u(x, t) = 0, (x, t) G dQ x (0, T) (7)
u(x, 0) = u0(x), ut(x, 0) = u1(x), x G Q, (8)
where A, a G R. Other modifications of the Boussinesq equation with nonlinearities of the form A|u|pu have also become widespread [4,5,17]. The equation has many applications in various fields of natural science. For example, it models the propagation of waves in shallow water, taking into account capillary effects. In this case, the function u = u(x,t) determines the height of the wave.
The article, in addition to the Introduction, Conclusion and References, includes seven paragraphs. The first section contains the main results of the theory of p-bounded operators obtained by G.A. Sviridyuk, necessary for further presentation. The second section presents results on the solvability of the Cauchy problem for semilinear first-order equations obtained by V.O. Kazak. The third paragraph contains some results of the study of Oskolkov generalized mathematical model. The fourth section presents the main statements of the theory of relatively polynomially bounded operator pencils, obtained by A.A. Zamyshlyaeva. The fifth section contains results on the solvability of the Cauchy problem for high-order semilinear Sobolev-type equations obtained in the authors' papers. The sixth paragraph is devoted to the study of the semilinear model of ion-acoustic waves. The seventh paragraph contains the results of the study of the mathematical model of shallow water wave propagation.
1. Relatively p-Bounded Operators
A detailed exposition of the theory of relatively p-bounded operators can be found in [30]. Let U, F be Banach spaces and operators L, M G L(U;F)-
Definition 1. The set
pL(M) = (ß G C : (ßL - M)-1 G L(F;U)}
is called a resolvent set of the operator M with respect to operator L (in short, L-resolvent set of the operator M). The set C\pL(M) = aL(M) is called a spectrum of the operator M with respect to the operator L (in short, the L-spectrum of the operator M).
Definition 2. Operator-functions
(ßL - M)-1, RL = № - M)-1L, LL = L(ßL - M)-1
with domain pL(M) are called respectively resolvent, right resolvent, left resolvent of the operator M with respect to the operator L (in short, the L-resolvent, right L-resolvent, left L-resolvent of the operator M).
Theorem 1. [30] Let the operator L G L(U,F), and the operator M : dom M C U ^ F be linear and closed. Then the L-resolvent, the right and the left L-resolvent are analytic in the set pL(M).
Definition 3. The operator M is said to be spectrally bounded with respect to the operator L (in short, (L, a)-bounded), if
3a> 0 Vß G C : (|ß| > a) ^ (ß G pL(M)).
Lemma 1. [30] Let the operator M be (L, a)-bounded. Then the operators
P = — [ Rx(M)dX andQ = — [ L{(M)d\ 2n J 2ni J
r r
are projectors. Here r = (A G C : |A| = r > a}.
Set U0 = ker P, F0 = ker Q, U1 = im P, F1 = imQ. Denote by Lk(Mk) the restriction of the operator L (M) to the subspace Uk, k = 0,1.
Theorem 2. [30] Let the operator M be (L, a)-bounded. Then
(i) Lk,Mk : Uk ^ Fk, k = 0,1;
(ii) Mo-1 G L(F0,U0);
(iii) operator L-1 G ¿(F1, U1) exists;
(iv) operator M1 G ¿(U1, F1) exists.
Let G ker L\(0} be an eigenvector of the operator L.
Definition 4. An ordered set • • •} C imL is called a chain of M-adjoined vectors
of an eigenvector if
L^q+1 = M^q, q = 0,1, 2 .q G ker L for q = 1, 2, • • •
The chain is said to be finite if there exists an M-adjoined vector such that either <fp G domM or G imL. The power of the final chain is called its length. The linear span of all eigenvectors and M-adjoined vectors of the operator L is called the M-root lineal of the operator L. Under the conditions of the theorem 2, we construct the operators H = M0-1 L0 G L(U0) and S = L-1M1 G ¿(U1). Since the operator-function (^L0 - M0)-1 is an entire function, it can therefore be expanded into a Taylor series
OuLo - Mo)-1 = (^H - I)-1M0-1 = ( — J] ^kHk ) M-1
k=0
absolutely and uniformly convergent on any compact set in C. Let's do the same with the
-1
operator-function (^L1 — M1 )"
(^¿1 — M1)-1 = (^I — S )-1L-1 = ^-1(I — )-1L-
£ S k L-1,
-k k
k=o
where ^ G p(S) or what is the same as ^ G pL(M). Hence, for the (L, a)-bounded operator, by virtue of the last two expansions, we have
(œ \ œ
— J] H M Mo-1 (I — Q) + £ S k-1L-1Q. (9)
k=0 / k=1
Let the operator M be (L, a)-bounded.
œ
œ
Definition 5. The infinity point of the L-resolvent of the operator M is called
• a removable singular point if H = O;
• a pole of order p if Hp = O, Hp+1 = O, p G N;
• an essentially singular point if Hq = O, Vq G N.
Further the removable singular point will be called a pole of order zero.
Remark 1. In what follows the (L, a)-bounded operator M will be called (L,p)-bounded, if the point to is a pole of order p G {0} U N of its L-resolvent.
Theorem 3. [30] Let L be a Fredholm operator (that is, ind L = 0). Then the following statements are equivalent:
1. The operator M is (L, 0)-bounded;
2. Any eigenvector of the operator L does not have M-adjoined vectors.
2. Semilinear Sobolev Type Equations of the First Order
Let the operators L, M G L(U; F), N G Ck(U; F), k G N U {to}, and the operator M be (L,p)-bounded, p G {0} U N. Consider the Cauchy problem
u(0) = u0 (10)
for a semilinear Sobolev type equation
Lu = Mu + N (u). (11)
The vector-function u G Ck((-T,T);U) is called a solution to equation (11) if for some T G R+ it satisfies this equation. The solution u = u(t) of equation (11) is called a solution to problem (10), (11) if it satisfies the initial condition (10).
Example 1. Let U = F = Rsn), the operators L, M, and N be defined by formulas
L =(0J) , M = I, N : u ), u =(£,n)-
Then the Cauchy problem with u0 = (0, 0), for equation (11) will have two solutions (0, 0) and (t/2, t2/4). If instead of the operator N in this case we take the operator N : u ^
2 ^, then the same problem will not have a solution at all. This simple example
shows the need to narrow down the definition of a solution to equations (11).
By virtue of theorem 2 equation (11) can be reduced to an equivalent system
HU0 = u0 + M0-1(I _ Q)N(u), (12)
U1 = Su1 + L-1QN (u), (13)
where u1 = Pu, u0 = u — u1.
Definition 6. Solution u = u(t) of problem (10), (11) is called a quasi-stationary trajectory for equation (11) passing through the point u0, if Hu(t) = 0 for all t G (—T, T).
Obviously, any stationary solution of problem (10), (11) is a quasi-stationary trajectory, however the opposite is not true. In example 1 a stationary solution is the only quasi-stationary trajectory passing through the point (0, 0), and in this sense such solution is unique. Further only quasi-stationary trajectories are considered.
To this end, introduce a set which is called the phase space of the equation (11):
By virtue of theorem 2 and (12) any quasi-stationary trajectory u = u(t) lies in M, that is u(t) G M for all t G (—T, T).
Let the point u0 G M, put u^ = Pu0 G U1. The set M at the point u0 is a Banach Ck-manifold if there are neighborhoods OM C M and O0 C U1 of the points u0 and u0, respectively, and a Ck-diffeomorphism 8 : O1 ^ OM such that 8-1 is equal to the contraction of P onto OM. The set M is called a Banach Ck-manifold modelled by the space U1, if it is a Banach Ck-manifold at each point of U1.
Theorem 4. [29] Let the set M be a Banach Ck-manifold at the point u0. Then there is a unique quasi-stationary trajectory of equation (11) passing through the point u0.
3. Mathematical Model of Oskolkov
Reduce problem (1) - (3) to the Cauchy problem (10) for the semilinear Sobolev type
equation (11). To do this, put U = W2m+2n W2, F = W2m, m G N. All functional spaces are defined on the domain Q. The operators L and M are defined by the formulas L =1 — ^A, M = vA. Obviously, L, M G L(U; F), and L is a Fredholm operator for all ^ G R \ (0}.
Lemma 2. [33] For all v G R \ (0} the operator M is (L, 0)-bounded.
Denote by (Ak} the set of eigenvalues of the homogeneous Dirichlet problem for the Laplace operator A in the domain Q, numbered in non-increasing order, taking into account multiplicity, and by (<^k} denote the set of orthonormal (in the sense of L2) of the corresponding eigenvectors. Then
In case ker L = (0} there is an operator L 1 G L(F; U), and so the operator M G L(U; F) (L, 0)-bounded. If ker L = (0} take the vector ^ G ker L \ (0}, i.e.
M = {u G U : (I — Q)(Mu + N(u)) = 0}.
œ 1 =A;
œ 1=A;
Since
œ 1 =Xi
then the vector ^ has no M-adjoint vectors, and the assertion follows from theorem 3. In what follows, we need regularity theorem [12]:
Lemma 3. Let f G C~(R) and m > n/2. Then F G Cwhere the operator F : u S f (u).
Lemma 4. Let the function K G C^(R) and let m + 2 > n/2. Then the operator N : u S K(u) belongs to the class C^ (U; F).
By virtue of lemma 3 the operator N G C^ (U), and in view of the continuity of embedding U S F , N G C~(U; F).
So, reduction of problem (1) - (3) to (10), (11) is finished. By lemma 1 construct a projector
I, if G {Ai},
P =1 i- E ^^
where (•, •) is an inner product in L2. The projector Q has the same form, but is defined on the space F. Fix m > n/2 — 2 and construct the set
M
and the space
U1
U, if ж-1 G {Ai},
{u G U : (Mu _ N(u), ^i) = 0, ж-1 = Ai}
U, if ж-1 G {Ai},
{u G U : (Mu _ N(u),^i) = 0, ж-1 = Ai}.
In the case of ^ 1 G {A^} the set M is obviously a smooth Banach C^-manifold. In the case of G {A^} this is yet to be proven.
Theorem 5. [33] (i) For any / {Ak}, v G R \ {0}, m > n/2 — 2, u0 G U and some T G R+ there exists a unique solution u G C^((—T,T); U) of problem (1) - (3).
(ii) Let for ^ G {Ak}, v G R \ {0}, m > n/2 — 2 set M at the point u0 be a Banach C^-manifold. Then for some T G R+ there exists a unique solution u G C^((—T, T); M) of problem (1) - (3).
Theorem 5 follows directly from theorem 4. We only note that in our case the operator H = O, and therefore, any solution to problem (1) - (3) necessarily turns out to be quasi-stationary trajectory.
4. Relatively Polynomially Bounded Operator Pencils
The statements presented in this paragraph were obtained in the works of A.A. Zamyshlyaeva [38,45]. Let U, F be Banach spaces and operators A, B0, B1,..., Bn-1 G L(U; F). By B denote the pencil formed by operators Bn-1,..., B1, B0. The sets pA(S) = {¡1 G C : (nnA-nn~lBn_l-...-nBl-B0)-1 G and aA(S) = C\pA(B) are
called an A-resolvent set and an A-spectrum of the pencil Sb respectively. The operator-function of a complex variable RA^(S) = (^nA — ^n-1Bn-1 — ... — ^B1 — B0)-1 with the domain pA(S) is called an A-rcsolvent of the pencil S.
Definition 7. The operator pencil B is called polynomially bounded with respect to an operator A (or polynomially A-bounded) if
3a e R+ fy e C (|^| > a) ^ (RA(B) e L(F;U)).
Remark 2. If there exists an operator A-1 e L(F; U) then the pencil B is A-bounded. Lemma 5. [38] Let the operator pencil B be polynomially A-bounded and condition
[ /rA= O, k = 0,1,..., n - 2, (14)
where the circuit y = G C : = r > a}, be fulfilled. Then the operators
P = 1hJ Q = ¿Ï /
Y Y
are projectors in spaces U and F respectively.
Denote U0 = ker P, F0 = ker Q, U1 = im P, F1 = im Q. According to lemma 5 U = U0 © U1, F = F0 © F1- By Ak (Bf ) denote restriction of operators A (B) onto Uk, k = 0,1; l = 0,1,...,n - 1.
Theorem 6. [38] Let the operator pencil B be polynomially A-bounded and condition (14) be fulfilled. Then
(i) Ak G L(Uk;Fk), k = 0,1;
(ii) Bf G L(Uk ; Fk ), k = 0,1, l = 0,1,...,n - 1;
(iii) operator (A1)-1 G ^(F1; U1) exists;
(iv) operator (B°)-1 G L(F0; U0) exists.
Using theorem 6 construct operators H0 = (B°)-1A0 G L(U0), H1 = (B°)-1B0 G L(U0),..., Hra-1 = (B00)-1B£-1 G L(U0) and S0 = (A1 )-1B0 G ¿(U1), S = (A1)-1Bj G L(U ), ..Sra-1 = (A1)-1B^-1 G ¿(U1).
Definition 8. Define the family of operators (Kj, K?,..., K^1} as follows:
K0s = O, s = n, K£ = I, K1 = H0, k? = -H1,..., Ks = -tfs_1,..., Kn = Hra-1, Kq1 = Kqn_1H0, K? = K1_1 - K_1#1, . . . , Kqs = K^1 - Kn_1#-_1, . . . ,
K = K-11 - Kq-iHn-1, q =1, 2,.... The A-resolvent can be represented by a Laurent series
ro
(^A - ^n-1B„-i - ... - ^Bi - B0)-1 = -J] ^^^(Bo0)-1(I - Q)+
q=0
ro
+ J] (^n-1S„-1 + ■ ■ ■ + + So)q A1-1Q. 9=1
Using this representation we classify the character of the infinity point of the A-resolvent of the operator pencil B.
Definition 9. The point to is called
• a removable singular point of an A-resolvent of the pencil B, if Kf = O, s = 1, 2,...,n;
• a pole of order p G N of an A-resolvent of the pencil B, if 3p such that Kpf = O, s = 1, 2,...,n, but Kp+1 = O, s = 1, 2,...,n;
• an essential singular point of an A-resolvent of the pencil B, if Kq1 = O for all q G N.
Further a removable singular point of an A-resolvent of the pencil B will be called a pole of order 0, for brevity. If the operator pencil B is polynomially A-bounded and the point to is a pole of order p G {0} U N of an A-resolvent of the pencil B then the operator pencil B is called polynomially (A,p)-bounded.
Theorem 7. [42] Let A, Bn-1, ..., B1; B0 G L(U,F) and A be a Fredholm operator. Then the following statements are equivalent:
(i) The lengths of all chains of the B-adjoined vectors of the operator A are bounded by number (p + n — 1) G {0} U N and the chain of length (p + n — 1) exists.
(ii) The operator pencil B is polynomially (A,p)-bounded.
5. Semilinear Sobolev Type Equations of Higher Order
Consider the Cauchy problem
u(k)(0) = uk, k = 0,1,..., n — 1, (15)
for a semilinear Sobolev type equation of higher order
Au(n) = Bn-1u(n-1) + Bn-2u(n-2) + ... + B0u + N (u), (16)
where operators A, Bn-1, Bn-2,..., B0 G L(U; F),N G C^(U; F), and U, F are Banach spaces.
Definition 10. If a vector-function u G C^((—t, t); U), t G R+ satisfies equation (16) then it is called a solution of this equation. If the vector-function satisfies in addition condition (15) then it is called a solution of problem (15), (16).
Definition 11. The set P is called a phase space of (16), if
(i) for all (u0,u1,... ,un-1) G Tn-1p there exists a unique solution of (15), (16);
(ii) a solution u = u(t) of (16) lies in P as a trajectory, i.e. u(t) G P for all t G (—t, t).
If ker A = {0} then equation (16) can be reduced to an equivalent equation
u(n) = F (u,u ,...,u(n-1)),
where F(u,u,... ,u(n-1)) = A-1(Bn-1u(n-1) + Bn-2u(n-2) +... + B0u + N(u)) is a mapping of class C^ by construction. The existence of a unique solution u of (15), (16) for all (u0,u1,... ,un-1) follows from the classical Cauchy theorem.
Let ker A = {0} and operator pencil B be (A, 0)-bounded, then by theorem 6 equation (16) can be reduced to an equivalent system of equations
0 = (I — Q)(B0 + N )(u0 + u1),
dn dn-1 dn-2 (17)
d „.1 _ A-1^(T3 " l D " l l D l AT\t„.0 i „,1\ ^ '
-u
Ar1^-!—-T + ? + ... + B0 + N){u° + ul)
where u1 = Pu, u0 = (I — P)u.
Now consider a set M = {u G U : (I — Q)(B0u + N(u)) = 0}. Let the set M be not empty, i.e. there is a point u0 G M. Denote u01 = Pu G U1. The set M is called a Banach Ck-manifold at point u0 if there exist neighborhoods O C M and O1 C U1 of points u0 and u0 respectively and a Ck-diffeomorphism 5 : O1 ^ O such that 5-1 is a restriction of projector P on O. The set M is called a Banach Ck-manifold modelled by the space U1 if it is a Banach Ck-manifold at any point.
Let the following condition be fulfilled
(I — Q)(B + NU0 ) : U0 ^ F0 is a toplinear isomorfism. (18)
According to the implicit function theorem there exist neighborhoods O0 C U0 and O1 C U1 of points u° = (I — P)u0, u0 = Pu0 respectively and the operator B G Cœ(O1; O0) such that u0 = B^). Lets construct an operator 5 = I + B : O1 ^ M, 5^) = u0. Then the operator 5-1 together with the set O1 makes a map of M and is a restriction of P on 5[O1] = O C M. Thus, we proved
Lemma 6. [42] The set M = {u G U : (I — Q)(B0u + N(u)) = 0} under condition (18) is a Cœ-manifold at point u0.
Lets act with the Frechet derivative 5(n) i i , of order n on the second equation of system (17). Since 5(u1) = u and
n
tl i i = i- (i(U*))
we obtain equation u(n) = F (u, u,..., u(n 1)), where
F(u,u,... ,u(n-1)) = ¿^i,...^^)A-1Q(Bn-1 u(n-1) + Bn-2u(n-2) + ...
+B0u + N(u)) G C~(U).
Therefore, we get
Theorem 8. [42] Let the operator pencil B be (A, 0)-bounded, N G C^(U; F) and condition (18) be fulfilled. Then for any (u0,u1,... ,un-1) G Tn-1M there exists a unique solution of problem (15), (16) lying in M as trajectory.
6. Mathematical Model of Ion-Acoustic Waves in Plasma
As a model example, consider problem (4), (5). In order to reduce mathematical model (4), (5) to problem (15), (16) set
U = {u G W2+2(Q) : u(x) = 0,x G dQ}, F = Wl(Q).
Define operators A = A — A, B2 = (A' — A), B0 = -a-j^, B3 = Bx = O. Operators
A, B3, B2, B0 are G L(U; F) for all l G {0} U N. Denote the eigenvectors of the Dirichlet
f . nkx1 . nmx^ nnx31 ,
problem (5) for the Lapiace operator by U)kmn = < sin-sin—-—sin->, where
[a b c J
/ i 2 / n'm\2 / nn\2
k, m, n G N and denote the eigenvalues by Akmn = — \ ( — j + + \—J • The
spectrum a(A) is negative, discrete, finite and tends only to —to. Since } C C^(Q)
we obtain
/A — ^3B3 — ^ — ^Bi — Bo
0
X—T 4 ' 2 /n^ 2
/ y [(Afcmn — A)/i + (Akmn — A )fJ> — a ( -J ] < ipkmni ' > <-Pkmn,
fc,m,ra=1
where < •, • > is an inner product in L2(Q).
Remark 3. In the case when (i) A G a (A) the A-spectrum of pencil B aA(B) = {^mn : r, m, n G N, j = 1,..., 4}, where ^rmn are the roots of equation
(nn \ 2
—J = 0. (19)
In the case when (ii) (A G a (A)) A (A = A') the A-spectrum of pencil B aA(B) = : k G N}, where are the roots of equation (19) with A = Aj. In the case when (iii) (A G a(A)) A (A = A') the A-spectrum of pencil B (B) = {j : k G N, k = l}.
Check condition (14). In case (i) there exists A-1 G ^(F1;U1) therefore condition (14) is fulfilled. In case (ii)
1 f ^ r J
1 / V^ ^ < ^fcmra, • > ^fcmra"^
E
7 , (At™. - A)/«4 + (A*™, - AOm2 - Q*2?)2
Y k , m , n—1 c
_ 1 f \ ^ ¡^ < ^Pkrnm ' > ^kmndu .
when r =1, therefore condition (14) is not fulfilled and this case is excluded from further considerations. In case (iii) (A G 0"(Д)) Л (A = A') condition (14) is fulfilled.
Lemma 7. [42] Let (i) A G ^(Д)) or (ii) (A G а(Д)) Л (A = A'). Then the pencil B is polynomially (A, 0)-bounded.
In case (i) ker A = {0} that is, the operator A has no eigenvectors and, by remark 2 the pencil B is (A, 0)-bounded.
In case (ii) A G а(Д) and A = A' construct the chain of B-adjoined vectors of an eigenvector = akmn^kmn G ker A \ {0}. Since B3 = B1 = O the first three
A—^fcmn
B-adjoined vectors can be taken equal to zero. On the fourth we obtain
Bo<~Po = Bq( ^^ (Ikmn^kmn) = _Oi ^^ 0>kmn<Pkmn G IHlA,
A Akmn A Akmn
since £ |afcmra | > 0.
A—^fcmn
Therefore the eigenvector doesn't have a B-adjoined vector of order four, the length of the chains of B-adjoined vectors of operator A is bounded by three, and the chain of length three exists.
Construct projectors. In case (i) P = I and Q = I. In case (ii)
P 1 ^ ^ < ^fcmra j ' > ^fcmraj
and the projector Q has the same form but it is defined on space F. Construct the set
(t)!
wi={ue 11: < a (—J u +A(u3),tpkmn > tpkmn = 0}.
^ — ^fcmn
By theorem 8 we have
Theorem 9. [42] (i) Let A £ 0"(A), (u0, ui,..., un-1) G Un. Then for some t = t(u0, «i,... , un-1) > 0 there exists a unique solution u G Cn ((—t, t), U) of problem (4), (5). (ii) Let (A G 0"(A)) A (A = A'), (u0,u1,... ,un-1) G Tn-1M and condition (18) be fulfilled. Then for some t = t(u0, u1,..., un-1) > 0 there exists a unique solution u G Cn((—t, t), M) of problem (4), (5).
7. Mathematical Model of Waves Propagation in Shallow Water
In some particular cases of a non-linear term in equation, one can not only answer the question of the existence and uniqueness of a solution, such as [41], but also find this solution. A detailed algorithm is described in [3], in this section we present only the main steps in finding a solution to problem (6) - (8).
For the solution, we need several function spaces. Let Q C be a domain with the boundary dQ of class C^, denote Q = Q x (0,T). Define spaces L4(Q), H¿(Q) and denote B = L4 (Q) n H¿(Q), D = H 1(Q) n coim L (where coim L = H 1(Q) © ker l). The operator A : H 1(Q) ^ H-1(Q) is given by formula
(Au,v) = — y (VuVv)dx. n
Introduce the notation of the operators
(Lu,v) = J(VuVv + Auv)dx, (Mu,v) = a2 J(VuVv)dx, (N(u),v) = J u3vdx. n n n
In addition, define distribution spaces (functions with values in a Banach space) (0,T; B) and (0,T; L2(H)). Construct dual spaces using the Dunford-Pettis theorem: (P°°(0, T; B))* ~ P^O, T; pt)(Q) U ii_1(n)) and (P°°(0, T; D))* ~ P^O, T; P>*).
Let Ak be the eigenvalues of the homogeneous Dirichlet problem (7)) for the operator A, numbered in nonincreasing order with multiplicity taken into account, and be the corresponding eigenvectors. Moreover, the linear span spanj^i, <^2,..., is dense in B for m —y to and is orthonormal (in the sense of an inner product in L2 (H)).
Theorem 10. [3] Let A G [Ax, +to), (u0,ui) G TP, where u0 G B = H(1(Q) n L4(Q) and U G L2(Q) n coimL. Then there exists a solution u = u(x,t) of problem (6) - (8) such that u G L~(0,T; H1(Q) n L4(Q)) and U G L~(0,T; L2(Q) n coimL).
Below we give a sketch of the proof. The solution of problem (6) - (8) will be sought in the form of the Galerkin approximation
um(t) = £ «Г(t)^fc. (20)
fc=i
Substituting the approximate solution (20) into equation (6) and multipliyng scalarly by the basis functions }m=1, we get
(Lum^fc)- a2(Aum^fc) + ((um)3,^fc> = 0, 1 < k < m. (21)
Using the series expansions of the initial functions in terms of basis functions, we obtain the initial conditions for the system of algebraic differential equations (21)
am(0) = em, am(0) = Yfcm 1 < k < m, (22)
m m
where «m = emVfc ^ u0 in B for m ^ то, and «m = S YfcVfc ^ U in L2(Q) as fc=i fc=i m ^ то. The existence of a unique local solution um = um(t,x), t G [0, tm] was proved.
After that, a priori estimates were obtained as follows. Multiplying equation (21) by am(t) (1 < k < m) and summing over k from 1 to m, we get
(Lum,Um) - a2(Aum, иm) + ((um)3, Um) = 0.
Introduce a norm in the space D (L2(Q) = coimL ф ker L) |U|Hi = (Lu, и). By the Courant principle, this norm is equivalent to the norm induced by the space H 1(П). Using the self-adjointness of L, A and integrating it on the segment [0, t], t < tm we obtain
•m 2 | „2ц„ m\\2 . -^II mi4 ^ ^ ino\
\u \ + a ||и \\Hi + -||и \\L4 < 6. (23)
The constant C does not depend on tm and hence tm = T.
Remark 4. Due to inequality (23) for m ^ то, the sequence of functions um is bounded in the space (0,T; L2(Q)) and и is bounded in L~(0,T; B).
Since the sequence {(umi)3} is bounded in the space L^(0,T; L4/3(Q)), we have
(umi)3 ^ z *-weakly in L~(0, T; L4/3(Q)). (24)
Moreover, it can be shown that z = и3.
Now we can pass term by term to the limit in equality (21), setting m^ = l. Let k be fixed and l > k, we get
(i) due to the density of the system of functions }m=1 in the space B for m ^ то, and the arbitrariness of the choice of , we have equality for arbitrary v G B
d2
— + Q!2(VM, vv) + (w3,w) = 0. (25) _dt/__
Вестник !Ю"УрГ"У. Серия «Математическое моделирование 55
и программирование» (Вестник ЮУрГУ ММП). 2022. Т. 15, № 1. С. 43-59
(ii) u(0) = u0;
(iii) (U(0), <£fc) = (ui, ), V k.
Remark 5. Under the conditions of the theorem 10 and the Rellich-Kondrashov theorem, the solution to problem (6) - (8) is unique.
Conclusion
Instead of the Cauchy condition in all mathematical models, one can consider the Showalter-Sidorov condition
P (u(k)(0) — uk ) = 0, k = 0,1,...,n — 1, (26)
where P is a projector along the kernel of the operator at the highest derivative with respect to t. Condition (26) is a natural generalization of the Cauchy condition for Sobolev type equations.
Further directions of development are seen in the study of semilinear Sobolev type equations with additive "white noise" [8,9,13], nonlinear inverse problems [43], as well as the study of multipoint initial-final problems [11,23,40].
Acknowledgments. The authors express their sincere gratitude to Professor G.A. Sviridyuk for posing the problem, interest in the work, and provided opportunities.
This work was supported by a grant from Ministry of Science and Higher Education of the Russian Federation No. FENU-2020-0022 (2020072GZ).
References
1. Al'shin A.B., Korpusov M.O., Sveshnikov A.G. Blow-up in Nonlinear Sobolev Type Equations. Berlin, Walter De Gruyter, 2011. DOI: 10.1515/9783110255294
2. Barenblatt G.I., Zheltov Yu.P., Kochina I.N. Basic Concepts in the Theory of Seepage of Homogeneous Fluids in Fissurized Rocks. Journal of Applied Mathematics and Mechanics, 1960, vol. 24, no. 5, pp. 1286-1303.
3. Bychkov E.V. Analytical Study of the Mathematical Model of Wave Propagation in Shallow Water by the Galerkin Method. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2021, vol. 14, no. 1, pp. 26-38. DOI: 10.14529/mmp210102
4. Chen P.J., Gurtin M.E. On a Theory of Heat Conduction Involving Two Temperatures. Zeitschrift fiir angewandte mathematik und physik, 1968, vol. 19, pp. 614-627. DOI: 10.1007/BF01594969
5. Clarkson P.A., Leveque R.J., Saxton R., Solitary Wave Interactions in Elastic Rods. Studies in Applied Mathematics, 1986, vol. 75, no. 1, pp. 95-122. DOI: 10.1002/sapm198675295
6. Cristiansen P.L., Muto V., Lomdahl P.S. On a Toda Lattice Model with a Transversal Degree of Freedom. Nonlinearity, 1991, vol. 4, no. 2, pp. 477-501. DOI: 10.1088/0951-7715/4/2/012
7. Demidenko G.V., Uspenskii S.V. Partial Differential Equations and Systems Not Solvable with Respect to the Highest Order Derivative. N.Y., Basel, Hong Kong, Marcel Dekker, 2003.
8. Favini A., Sviridyuk G., Manakova N. Linear Sobolev Type Equations with Relatively p-Sectorial Operators in Space of "Noises". Abstract and Applied Analysis, 2015, vol. 2015, article ID: 697410.
9. Favini A., Sviridyuk G., Sagadeeva M. Linear Sobolev Type Equations with Relatively p-Radial Operators in Space of "Noises". Mediterranean Journal of Mathematics, 2016, vol. 13, no. 6, pp. 4607-4621. DOI: 10.1007/s00009-016-0765-x
10. Favini A., Yagi A. Degenerate Differential Equations in Banach Spaces. N.Y., Basel, Hong Kong, Marcel Dekker, 1999.
11. Favini A., Zagrebina S.A., Sviridyuk G.A. Multipoint Initial-Final Value Problems for Dynamical Sobolev-type Equations in the Space of "Noises". Electronic Journal of Differential Equations, 2018, vol. 2018, article ID: 128.
12. Hassard B.D. Theory and Application of Hopf Bifurcation. Cambridge University Press, Camdribge, 1981.
13. Keller A.V., Shestakov A.L., Sviridyuk G.A., Khudyakov Yu.V. The Numerical Algorithms for the Measurement of the Deterministic and Stochastic Signals. Semigroups of Operators - Theory and Applications. Heidelberg, N.Y., Dordrecht, London, Springer Int. Publ. Switzerland, 2015, pp. 183-195. DOI: 10.1007/978-3-319-12145-1_11
14. Korpusov M.O., Lukyanenko D.V. Instantaneous Blow-up Versus Local Solvability for one Problem of Propagation of Nonlinear Waves in Semiconductors. Journal of Mathematical Analysis and Applications, 2018, vol. 459, no. 1, pp. 159-181. DOI: 10.1016/j.jmaa.2017.10.062
15. Kozhanov A.I. On a Nonlocal Boundary Value Problem with Variable Coefficients for the Heat Equation and the Aller Equation. Differential Equations, 2004, vol. 40, no. 6, pp. 815-826. DOI: 10.1023/B:DIEQ.0000046860.84156.f0
16. Ladyzhenskaya O.A. Matematicheskie voprosy dinamiki vyazkouprugoy neszhimaemoy zhidkosti [Mathematical Problems in the Dynamics of a Viscoelastic Incompressible Fluid]. Moscow, Fizmatgiz, 1961. (in Russian)
17. Lions J.L. Sur quelques methodes de resolution des problemes aux limites non linears. Paris, Dunod, Gauthier Villars, 1969. (in French)
18. Manakova N.A., Gavrilova O.V. About Nonuniqueness of Solutions of the Showalter-Sidorov Problem for One Mathematical Model of Nerve Impulse Spread in Membrane. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2018, vol. 11, no. 4, pp. 161-168. DOI: 10.14529/mmp180413
19. Oskolkov A.P. Nonlocal Problems for Some Class Nonlinear Operator Equations Arising in the Theory Sobolev Type Equations. Journal of Soviet Mathematics, 1993, vol. 64, no. 1, pp. 724-735. DOI: 10.1007/BF02988478
20. Poincare H. Sur l'equilibre d'une massuide animee d'un mouvement de rotation. Acta Mathematica, 1885, vol. 7, pp. 259-380. DOI: 10.1007/BF02402204 (in French)
21. Pyatkov S.G. Operator Theory. Nonclassical Problems. Utrecht, Boston, Köln, Tokyo, VSP, 2002. DOI: 10.1515_9783110900163
22. Sagadeeva M.A., Sviridyuk G.A. The Nonautonomous Linear Oskolkov Model on a Geometrical Graph: The Stability of Solutions and the Optimal Control. Semigroups of Operators - Theory and Applications. Heidelberg, N.Y., Dordrecht, London, Springer Int. Publ. Switzerland, 2015, pp. 257-271. DOI: 10.1007/978-3-319-12145-1_16
23. Sagadeeva M.A., Zagrebina S.A., Manakova N.A. Optimal Control of Solutions of a Multipoint Initial-final Problem for Non-autonomous Evolutionary Sobolev Type Equation. Evolution Equations and Control Theory, 2019, vol. 8, no. 3, pp. 473-488. DOI: 10.3934/eect.2019023
24. Showalter R.E. Hilbert Space Methods for Partial Differential Equations. Pitman, London, San Francisco, Melbourne, 1977.
25. Shubin Wang, Guowang Chen. Small Amplitude Solutions of the Generalized IMBq Equation. Journal of Mathematical Analysis and Applications, 2002, vol. 274, no. 2, pp. 846-866. DOI: 10.1016/S0022-247X(02)00401-8
26. Sidorov N., Loginov B., Sinithyn A., Falaleev M. Lyapunov-Shmidt Methods in Nonlinear Analysis and Applications. Dordrecht, Boston, London, Kluwer Academic Publishers, 2002.
27. Sobolev S.L. On a New Problem of Mathematical Physics. Izv. Akad. Nauk SSSR. Ser. Mat., 1954, vol. 18, no. 1, pp. 3-50. (in Russian)
28. Sviridyuk G.A. A Problem of Generalized Boussinesq Filtration Equation. Soviet Math., 1989, vol. 33, no. 2, pp. 62-73.
29. Sviridyuk G.A. Quasistationary Trajectories of Semilinear Dynamical Equations of Sobolev Type. Russian Academy of Sciences. Izvestiya Mathematics, 1994, vol. 42, no. 3, pp. 601-614. DOI: 10.1070/IM1994v042n03ABEH001547
30. Sviridyuk G.A., Fedorov V.E. Linear Sobolev Type Equations and Degenerate Semigroups of Operators. Utrecht, Boston, Köln, Tokyo, VSP, 2003. DOI: 10.1515/9783110915501
31. Sviridyuk G.A., Karamova A.F. On the Phase Space Fold of a Nonclassical Equation. Differential Equation, 2005, vol. 41, no. 10, pp. 1476-1481. DOI: 10.1007/s10625-005-0300-5
32. Sviridyuk G.A., Kazak V.O. The Phase Space of an Initial-Boundary Value Problem for the Hoff Equation. Mathematical Notes, 2002, vol. 71, no. 2, pp. 262-266. DOI: 10.1023/A:1013919500605
33. Sviridyuk G.A., Kazak V.O. The Phase Space of a Generalized Model of Oskolkov. Siberian Mathematical Journal, 2003, vol. 44, no. 5, pp. 877-882. DOI: 10.1023/A:1026080506657
34. Sviridyuk G.A., Kitaeva O.G. Invariant Manifolds of the Hoff Equation. Mathematical Notes, 2006, vol. 79, no. 3, pp. 408-412. DOI: 10.4213/mzm2713
35. Sviridyuk G.A., Manakova N.A. Phase Space of the Cauchy-Dirichlet Problem for the Oskolkov Equation of Nonlinear Filtration. Russian Mathematics, 2003, vol. 47, no. 9, pp. 33-38.
36. Sviridyuk G.A., Manakova, N.A. An Optimal Control Problem for the Hoff Equation. Journal of Applied and Industrial Mathematics, 2007, vol. 1, no. 2, pp. 247-253. DOI: 10.1134/S1990478907020147
37. Sviridyuk G.A., Sukacheva T.G. The Cauchy Problem for a Class of Semilinear Equations of Sobolev Type. Siberian Mathematical Journal, 1990, vol. 31, no. 5, pp. 794-802. DOI: 10.1007/BF00974493
38. Sviridyuk G.A., Zamyshlyaeva A.A. The Phase Spaces of a Class of Linear HigherOrder Sobolev Type Equations. Differential Equations, 2006, vol. 42, no. 2, pp. 269-278. DOI: 10.1134/S0012266106020145
39. Temam R. Navier-Stokes Equations. Theory and Numerical Analysis. Amsterdam, N.Y., Oxford, North Holland Publ. Co., 1979.
40. Zagrebina S.A., Konkina A.S. The Multipoint Initial-Final Value Condition for the Navier-Stokes Linear Model. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2015, vol. 8, no. 1, pp. 132-136. DOI: 10.14529/mmp150111
41. Zamyshlyaeva A.A., Bychkov E.V. The Cauchy Problem for the Second Order Semilinear Sobolev Type Equation. Global and Stochastic Analysis, 2015, vol. 2, no. 2, pp. 159-166.
42. Zamyshlyaeva A.A., Bychkov E.V. The Cauchy Problem for the Sobolev Type Equetion of Higher Order. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2018, vol. 11, no. 1, pp. 5-14. DOI: 10.14529/mmp180101
43. Zamyshlyaeva A., Lut A. Inverse Problem for the Sobolev Type Equation of Higher Order. Mathematics, 2021, vol. 9, no. 14, article ID: 1647. DOI: 10.3390/math9141647
44. Zamyshlyaeva A.A., Manakova N.A., Tsyplenkova O.N. Optimal Control in Linear Sobolev Type Mathematical Models. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2020, vol. 13, no. 1, pp. 5-27. DOI: 10.14529/mmp200101
45. Zamyshlyaeva A.A., Sviridyuk G.A. Nonclassical Equations of Mathematical Physics. Linear Sobolev Type Equations of Higher Order. Bulletin of the South Ural State University. Series: Mathematics. Mechanics. Physics, 2016, vol. 8, no. 4, pp. 5-16. DOI: 10.14529/mmph160401
Received October 8, 2021
УДК 517.9 DOI: 10.14529/mmp220103
ПОЛУЛИНЕЙНЫЕ МАТЕМАТИЧЕСКИЕ МОДЕЛИ СОБОЛЕВСКОГО ТИПА
А.А. Замышляева1, Е.В. Бычков1
1Южно-Уральский государственный университет, г. Челябинск, Российская Федерация
Статья содержит обзор результатов, полученных в научной школе Георгия Анатольевича Свиридюка, в области полулинейных математических моделей соболевского типа. В работе приведены результаты о разрешимости задачи Коши и Шоуолтера -Сидорова для полулинейных уравнений соболевского типа первого, второго и высокого порядков, а также примеры неклассических моделей математической физики, такие, как обобщенная модель нелинейной фильтрации Осколкова, распространения ионно-акустических волн в плазме, распространения волн на мелкой воде, которые исследуются путем редукции к одной из вышеперечисленных абстрактных задач. Методы исследования полулинейных уравнений соболевского типа базируется на теории относительно р-ограниченных операторов для уравнений первого порядка по переменной t и теории относительно полиномиально ограниченных пучков операторов для уравнений второго и высокого порядка по переменной t. В работе применяется метод фазового пространства, заключающийся в редукции сингулярного уравнения к регулярному, определенному на некотором подпространстве исходного пространства, для доказательства теорем существования и единственности и метод Галеркина для построения приближенного решения.
Ключевые слова: уравнение Осколкова; модифицированное уравнение Буссинеска; уравнение ионно-звуковых волн в плазме; полулинейные уравнения соболевского типа; относительно р-ограниченные операторы; относительно полиномиально ограниченные пучки операторов; метод Галеркина; *-слабая сходимость.
Алена Александровна Замышляева, доктор физико-математических наук, профессор, кафедра «Прикладная математика и программирование:», ЮжноУральский государственный университет (г. Челябинск, Российская Федерация), zamyshliaevaaa@susu.ru.
Евгений Викторович Бычков, кандидат физико-математических наук, кафедра «Уравнения математической физики>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), bychkovev@susu.ru.
Поступила в редакцию 8 октября 2021 г.
