Об обобщенной разрешимости смешанной задачи для нелинейного уравнения с псевдопараболическим оператором высокой степени
Автор: Юлдашев Т.К.
Журнал: Сибирский аэрокосмический журнал @vestnik-sibsau
Рубрика: Математика, механика, информатика
Статья в выпуске: 2 (48), 2013 года.
Бесплатный доступ
Доказываются теоремы об обобщенной разрешимости смешанной задачи для нелинейного дифференциального уравнения с псевдопараболическим оператором произвольной натуральной степени.
Обобщенная разрешимость, интегральное тождество, система нелинейных интегральных уравнений, метод последовательных приближений
Короткий адрес: https://sciup.org/148177042
IDR: 148177042
Текст научной статьи Об обобщенной разрешимости смешанной задачи для нелинейного уравнения с псевдопараболическим оператором высокой степени
В области D рассматривается уравнение
^ 2 m +1 *\ 4 m +1 -л 4 m Л
— + (-1) m V—--—+ УЦ—--—+ -д— X d t д t д X2 m д t д X4 m д X4 m J (1)
X u ( t , X ) = f ( t , X , u ( t , X ) )
= J ф <4 nm - 2)( y ) dy = 0, J = 1, n , 0
D = D t x D i , D t = [ 0, T ] , D i = [ 0, l ] ,
0 < l
u (t , X) |t=0 = Ф1( x), d j-1
—-1 u(t,X) t=0 = Фj (x), j = 2,n d tJ 1
и граничными условиями
u (t , X)|X =0 = u xx (t , X)|X =0 =. . . = d 2(2nm-1)
= . 2(2 nm-1) u (t , X ) x=0 = J u (t , У) dy =
О x0
малые параметры, n, m - натуральные числа.
Следует отметить, что изучению разного типа линейных и нелинейных дифференциальных уравнений в частных производных и их систем посвящено много работ и при этом применены разные методы (см., например [1–4]). В данной работе, как и в работах [5; 6], используется метод разделения переменных, основанный на поиске решения смешанной задачи (1)–(3) в виде
ll
= J uyy (t, У) dy =... = J
5 2(2nm-1)
---—---- u (t, y) dy = 0, d y 2(2nm-1) v 7 17
го
u(t, x) = E ai(t) ■ bi(x), (4)
i =1
где bi (x) = . - sin X i x , X i =
где f (t, X, u) e С(D X R), VJ (X) e C 2nm+1(Dz),
ФД x) x=0 =Ф j "(x )| x=0 = ... = Ф J4 nm-2)(x) x=0 =
J ФJ(y)dy =J ФJ ”(y)dy = ... = J Ф(4nm-2)(y)dy = 0 ,
00 0
2 i n l
Множество {a(t) = (ai(t))| ai(t) e C(DT), i = 1, 2, ...}
введением нормы
a(t) B2(T)
го
I i=1
, . Л2
max ai(t) teDt 'J
становится банаховым пространством и обозначается через B2(T).
Наряду с пространством B2(T) рассмотрим и про-
l
= J Ф1(y)
д n "1
---- Ф + n д tn-1
д n+4 m-2
д tn-2д y4 m
Ф +
странство B2N (T) с нормой
+ n(n - 1) _
2 d t
d n + 4 m-1
n-3 d y 4 m+2
Ф +
N (
B2N(T)
= I
max t е Dt
. 12
I ai(t )|
.
+
.
Для каждого элемента
a (t) е B 2 (T) определяется
/ ix 4 nm-3 ^4 nm-2
n (n -1) d _ д _
.. +---л---т Ф + n -------— Ф +
2 дt дy4nm-4 дy4nm-2
| д n+2 m-1 д n+6 m-2
+v Ф + n Ф +
(д tn-1д y2 m д tn-2д y6 m
w оператор: Qa (t) = Iai (t)■ bt (x).
i =1
n (n -1) дn+6 m-1
+ 2 д t"-3дy6 m+2
_ n (n -1)
Ф+ . . . +—----X
Обозначено через E2 (D) множество значений
оператора Q . Здесь очевидно, что Q :B2 (T) ^ E2(D)
d 4 nm+2 m-3
X dt dy 4nm+2m-4
Ф + n
4 nm+2 m-2A дI
--;--z—- Ф + дy 4nm+2m-2I
и E 2( D) с L2( D).
Обозначается через W2(k)(D) множество функций
( n +4 m-1
I d
+vp
( dtn-1dy4m
Ф + n
д n+8 m-2
a tn-2a y8 m
Ф +
а2
Ф(t, x) таких, что Ф(t, x), —- Ф (t, x), дxx
..
d 2(2 nm-1)
., n Ф(t,x) dx 2(2 nm-1)
n (n -1) дn+8 m-1 n (n -1)
+---7----„ э Ф + . . . +--X
2 д tn-3дy8 m+2 2
при фиксированном t е DT принадлежат области оп-
d 4nm-2
ределения оператора ^ 4nm-2 , имеют производные
d 4 nm+4 m-3
X------".----л---T Ф + n dt dy 4nm +4m-4
d 4 nm+4 m-2 Ф дy 4nm+4m-2
dy - t=0
.
порядка k по t , принадлежащие L2 (D), и обраща-
ются в нуль при t > T -8 (0 < 8 - зависит от Ф (t, x)).
Определение. Если функция u (t, x) е E2(D )
l
. . +Jф n-1(y)
| d 2 m+1
v
(d t dy2 m
d d4 m
— Ф + n--— Ф + a t a y4 m
а6 m ।
Ф + n ——Ф +
6m dy )
удовлетворяет интегральному тождеству
t 1 Г
JJiu(t,y) 0 0 I
д n
---Ф + n a tn
д n + 4 m-1
д tn-1дy4 m
Ф +
n (n-1) dn+4 m
--Ф + . . . +
2 d tn-2ay4 m+2
n (n -1) д4 nm-2
2 д t2дy 4 nm-4
Ф + n
^4 nm-1 ^4 nm
----: гФ + :— dt ay4nm-2 ay4nm
Ф +
д n+2 m
+Ч atnay2m
Ф + n
д n+6 m-1
a tn-1dy6 m
Ф +
( 4 m+1^8
I dd
+ vp ---— Ф + n—— Ф
(dt dy4 md
1 Г
Xdy-Jv n(y) Ф + v 0 _
X
Jt=0
d 2 m d 4 m
-— Ф + vp ——Ф a y2 m a y4 m
dy
-1=0
для любого Ф (t, x) е W2(n) (D), то она называется обобщенным решением смешанной задачи (1)–(3).
Покажем, что коэффициенты разложения ai (t)
решения смешанной задачи (1)–(3) удовлетворяют следующей счетной системе нелинейных интегральных уравнений (ССНИУ):
n (n -1) an+6 m
2 d tn-2a y6 m+2
n (n -1) Ф+ ... +—----X
ai(t ) = wi(t) +
д 4nm+2m-2
д t 2 д y 4nm+2m-4
Ф + n
4 nm+2 m-1 A d I
Ф + d t d y4nm+2m-2 J
tl (5)
+ JJ f ( 5,y,Q a(5)) bi(y)Pi(t,5)dy ds,
где
д n+4 m
+w( atnay4m n (n -1) дn+8 m
Ф + n
д n+8 m-1
д tn-1дy8 m
Ф +
n w (t) = I Ф ki
tk-1
(k -1)
n
! I6 хk (v, m)X
д tn-2дy8 m+2
n (n-1) Ф+ . . . +—----X
t j - k
X ^7-1)7 ■exp {-61i(v,p)t},
d 4 nm+4 m-2 dt 2 dy 4nm+4m-4
d 4 nm+4 m-1
Ф + n ----------Ф dt dy4nm+4m-2
- f Ф^ dydt =
Pi(t,») =
(n-1)!( t - 5)n-1
6 o,:(v, Ц)
■ exP {—6 1 i(v, Ц)(t-») },
X 4 nm 0; (v, ц=—-—, 1 0 n0i(v, ц)
0 0,(v, ц) = ( 1 + vx 2 m +
n
+ n (n 1) X 4 nm+2 m—4a ,,(t) + n X 4n m+2 m—2a '(t)) +
2 i i i i I
+v ц(X4ma(n)(t) + nX8ma(n—1)(t) +
Ф ki = Jф —(У) bi(У) dy-о
Нас интересует укороченная система нелинейных интегральных уравнений (УСНИУ):
aN (t) = wN (t) +
+ n(n—1)
+ n(n—1)
—
8 m+2a(n—2)(t) + - - - +
4nm + 4m —4 . ил 4nm + 4m—2 f/,\ i ai (t)+ nXi ai (t)
■ Jf (t,y, Qa(t) )■ b(y)dy dt = 0
0 J
—
t1 . . (6)
+ JJf( s,У,QN a(s)) bi(У)Pi(t,s)dy ds,
0 0
где n k——1 n wN(t) = ЕфN7T"ni S0Juk(v,ц) x
— =1 (k 1)! j=— t j - k x jk)7 ■exp {—0 i(v ’ц) t} ’
P(t, s) = (n Jn^t \'— exP {- 01,(v, Ц)(t - s)}, 00i(v, Ц)
4 nm
X •
01i(v,P) = „ , ,
00 i(v, Ц)
0 n0i(v, ц) = ( 1 + vx 2 m +
n
,
QN a (t) = uN (t,.) = £>;'' (tyb, (x)- i =1
Так как h (t) - любая функция, удовлетворяющая указанным выше условиям, то at (t) имеет обобщенные производные порядка n по t в смысле Соболева на отрезке DT -
Из (7) следует, что л (n) А \ 4m (n—1)/v\ . n (n 1) л 4m+2 (n—2)
ai(t)+nXi ai (t)+—2— Xi ai (t)+ --- +
+ n (n2—1) x4nm—4ai"(t) + nX4nm—2ai'(t) + X4nm ai (t) + v ( X2ma(n)(t) + nX6m ain—0(t) + n^n.—1)X6m+2a(n—2)(t) + , n(n— 1)(n— 2)^6 m + 4 „( n—3)/^ , ,
+ 3! X i ai (t)+ - - - + n (n — 1)(n — 2) 4nm+2m—6 „ m/A, + -------7;-------X i ai (t) +
+ n (n2 1) X 4nm+2m—4a i „( t) + n X 4n m+2 m—2a i (t) | + +vц (X,дma(n)(t) + nXima(n—1)(t) +
С учетом (4), в силу того что Ф = Ф j (t, x) = = h (t) bj (x) e W2(n) (D) в определении и bj (x) полны и ортонормированы в L2(Dl), где h (t) e Сn (DT), j = 1,2,--., следует
T
J h (t) [ a(n)(t) + n X 4 ma(n—1)( t) + 0
n(n— 1). 4 m+2 Г n—2)(a nn(n— 1)(n— 2)
+ 2 X i ai (t)+ xX4m+4 a(n—3)(t) + --- + n(n — ^n — 2) X4nm—6 ai"'(t) +
+ n^ 11X4nm—4a,"(t) + nx4nm—2a,'(t) + X4nm a, (t) +
2 m (n )/,\ „ л 6 m Г n-1). n (n — 1) л 6 m+2 „(n-2).
v(Xi ai(t)+nXi ai(t)+—2— X, a(t)+
, n (n — 1)(n — 2)^6m+4 „(n —3)^з,,
+ 3! X i ai (t)+ - - - + n (n — 1)(n — 2) 4nm+2m—6।
+-------3!-------Xi ai(t)+
+ n(n2 1) xfm+2a(n—2)(t) + - - - +
+ n(n-1)x4nm+4m—4ai.(t) + nX4nm+4m—2ai'(t)| =
= Jf (t,y,Qa(t) )■ b(y)dy -
Решая систему (8) методом вариации произвольных постоянных, получаем ai (t , v , ц) = ( С1 i + C2 it + C3 it 2 + C4 it 3 + --- + Cnit" 1 )x tl x exp{—01i(v, ц) t}+ J J f (s, y, Qa (s) )x (9)
0 0
xЬ,(y)P(t,s)dyds, t e DT-
Используя условия ai (0) =Ф1 i , ai '(0) = Ф2i, ai "(0) = ф3i, ---, a(n—1)(0) =фп, для определения коэффициентов Cji (j = 1, n), из (9) получаем (6)-
Теорема 1. Пусть выполняются следующие условия:
T
решения: a(t) е B2(T) и 9(t) е B2(T). Тогда для их разности справедлива оценка
| dt< A < м;
L 2( Dl)
f(t, x, и) е Lip{ L(t,x) |u } ,
где 0 <j||L(s,x)| 0
3 II w(t)!IBN(T)
। L 2( Dl) d s ';
< м,
где II w(t)|IBN(T)
N
= z
f .)2
max ai(t)
t е DT
.
Тогда УСНИУ (6)
имеет единственное решение
в пространстве B2 (T).
Доказательство. Используется метод последовательных приближений:
a0 (t) = w (t), t е Dt , a k+1 (t) = w, (t) + +jjf (s,У, QNak(s))b(У)Pt(t,s)dy ds,к = 0 0
= 0,1,2,3
,...,
t е DT.
Учет условий теоремы в (10) дает оценки
II a 1(t) - a0(t)|| < M1M2 lq A ,
II II B 2 ( T )
I a2(t)-a1(t )l IBn ( t )
2 t
M12nM23,N lq j||L(s,x)|Il2(Dl) ds , 0
где
M1, N =11 P(t , s ) ||B N ( T) , M 2, N = ||b(x) ||B N (l) .
Продолжая этот процесс для произвольного натурального числа к , подобно (12) получаем
IIak+1(t)-ak(t)|IBN(t) <
M
1 Л к+1 l2
M2 k+1
M 2,N
jlIL(s , x )| L2( Dl) ds
A—---------------:
k
к!
.
Существование решения из оценки (13), так как при
УСНИУ (6) следует к ^ м последователь-
к м ность функций { a (t)} сходится равномерно по t к функции a (t) е B2 (T). Для проведения доказательства единственности этого решения в пространстве B2 (T) предполагается, что УСНИУ (6) имеет два
II a(t)-9(t)||Bn(t)
XjllL(s,x )| Il2( Dl )lla (s)-9 (s) I IB 2n ( t ) ds .
Применение к (14) неравенства Гронуолла-Белл-мана дает, что || a(t)-9(t)||bn(T) ^ 0 для всех t е DT. Отсюда следует единственность решения УСНИУ (6) в пространстве B 2 (T).
Рассмотрим формулу (4) как следующее предельное соотношение:
N u (t,x) = lim u N(t,x) = lim V ai(t)■ bi(x). (4‘)
N^м N^м , i =1
Подставляя ССНИУ (5) в предел (4’), получим формальное решение смешанной задачи (1)-(3):
N
tl
u(t, x) = lim Vr wi(t)+f [f (s, y, QN a (s) )x
N ^”i=1 L0o ' (15)
x bi (y)Pi(t,s)dy ds]■ bi (x).
Теорема 2. Пусть выполняются условия теоремы 1 и II w(t)|IB (T) <м. Если a(t)еB2(T) является решением УСНИУ (6), то (15) будет обобщенным решением смешанной задачи (1)-(3).
Доказательство. Так как a (t) е B 2 (T), то из ра-
N венства lim uN (t,x) = lim Vai (t)■ bi (x) = u (t,x) N^м N^м i =1
в силу условий теоремы следует, что lim f ( t, x, uN (t, x) ) = f (t, x, u (t, x)) (16) N ^м
в смысле метрики L2(D).
Рассмотрим последовательность функционалов:
Vn =jj1 uN(t, У)
0 0 [
д n
---Ф + n a tn
d n+4 m-1
----;---Т" Ф + д tn-1dy4 m
+
+n (n -1) 2
n (n -1) а4 nm-2
д n+4 m a tn-2a y4 m+2
2 dt2дy4 nm-4
( д n+2m
Ф + n
+v ^ дtnay2m
Ф +
... +
д 4 nm-1
дt ay4nm-2
4 nm
Ф + 1—Ф + дy4nm
д n+6 m-1
Ф + n ---:---—Ф + д tn-1 д y6 m
n (n -1) dn+6 m
+ 2 d tn-2dy6 m+2
. n (n -1) Ф+ ... + —----X
d 4 nm+2 m-2 ^ 4 nm+2 m-1
X--:---:---:--гФ + n -----;---^Ф + dt 2 dy 4nm+2m-4 dt dy 4nm+2m-2I
+VЦ
д n + 4 m
д tn д y4 m
д n+8 m-1
Ф + n -----:----;
д t
д y8 m
Ф +
+n (n-1)
X
д
д t
n (n -1)
д n+8 m
д tn-2д y86m + 2
д 4nm+4m-2
д t 2 д y 4nm+4m-4
Ф + n
д n+4 m-1
д tn-3
----л---T Ф + . . . + д y4 m +2
_ n (n -1) Ф+ ... +—----X
n (n -1) д4 nm-3
+ * дt д y4nm-4
Ф + n
д4nm-2
д y4 nm-2
Ф +
д 4 nm+4 m-1
д t д y 4nm+4m-2
Ф
l
+f (t, y, uN (t, y)) ф} dydt-/фГ (У) 0
Ф + n
д n +4 m-2
д tn-2дy4 m
Ф +
n (n -1) дn+4 m-1
д tn-3
+
+v l д t
д n+2 m-1
д y2 m
Ф + n
д n+6 m-2
д tn-2д y6 m
Ф +
+n (n -1)
д n+6 m-1
д tn-3
---2—ГТ Ф + д y6 m + 2
Ф + дy4 m + 2
+
.
. . +
n (n -1) д4nm+2m-3
дt дy4nm+2m-4
Ф + n
д 4 nm+2 m-2
дy 4nm+2m-2
Ф +
n (n-1) д4 nm-3
+... + дt дУ4nm-4
Ф + n
д 4 nm-2
д y4nm-2
Ф +
+VЦ
д n +4 m-1
д tn-1д y
—Ф + n
4m
д n+8 m-2
д tn-2д y8 m
Ф +
д n+2 m-1
+vl дtn-1дy
Ф + n
2m
д n+6 m-2
д tn-2д y6 m
Ф +
+n (n -1)
д n+8 m-1
д tn-3
-----ф + д y 8 m + 2
n (n -1) д4nm+2m-3
+... + 2 дt дy4nm+2m-4
Ф + n
д 4 nm+2 m-2
дy 4nm+2m-2
Ф +
n (n -1) д4 nm+4 m-3
. . . * д t д y 4nm+4m-4
Ф +
+ и(и-1)
д n+6 m-1
д tn-3
--—7Ф + vц д y6 m + 2
д n+4 m-1
д tn-1дУ4 m
Ф +
д 4 nm+4 m-2
+ n ----;----;---Ф д y 4nm+4m-2
d y -
J t=0
д n+8 m-2
n д tn-2 д y 8
n (n -1) Ф + —---
д n+8 m-1
д tn-3
-----Ф + д y 8 m + 2
.
l, ( N )
. . +JI фn-1(У)-Хф(n-1).Ь,(У) lX
i=1
n (n-1) д4nm+4m-3
... * дt дy 4nm + 4m-4
Ф +
X
д
— Ф + n д t
д4 m
д y4 m
Ф + v
д 2 m+1
д t д y2 m
Ф + n
д6 m
д y6 m
Ф +
д 4nm+4m-2
n ---;----;---Ф дy 4nm+4m-2
Jt=0
-
l
/фN-1(У)
д 2 m +1
+ v -----5—
( д t д y2 m
Ф + n
X
dy + ...
-
+ v^
д 4 m +1
д t д y4 m
Ф + n
д 8m
д У8 m
Ф
d y -
Jt=0
д
—Ф + n д t
д4 m
д y4m
Ф +
lr ( N )
-JI Фn(У)-Ефnibi(У) lX
д6 m
д y6 m
Ф l + v^
д 4 m+1
д t д y4 m
Ф +
X
Ф
i =1
д2 m
д y2 m
Ф + vц
д4 m
д y4 m
Ф
dy +
д8 m
+n--Ф дy8 m
Ф + v
д2 m
дy2 m
Jt =0
dy +J Ф N (У) X
Ф + vц
д 4 m
Tl
+JJ Ф(t, У)[f(t,У,Qa (t)) 00
Jt=0
дy4 m
Ф
dy .
Jt=0
■ ]T Jf (t,z,QNa(t))• bi(z)dz bi(y)dydt. i=1 0 _
Интегрируя по частям отдельные слагаемые в (17) и учитывая условия теоремы и начальные условия
Очевидно, что первые n интегралы в (18) стре-
получаем
а,:(0)=Ф1/, а,: '(0)= Ф2i,
ai "(0)= Ф 3 i ,
a (n-1)(0) = Ф ni,
l- ( N )
VN =JI Ф 1(У)-^Ф 1 -b-^У) lX
i=1
X
д
д t
Ф + n
д n+4 m-2
д tn-2д y4 m
Ф +
мятся к нулю при N ^м , так как фj (x )еL2(Dl),
j = 1, n. Сходимость последней разности в (18)
при N ^м следует из (16). Отсюда ясно, что lim VN = 0 . Это и доказывает теорему.
N ^м