Об обобщенной разрешимости смешанной задачи для нелинейного уравнения с псевдопараболическим оператором высокой степени
Автор: Юлдашев Т.К.
Журнал: Сибирский аэрокосмический журнал @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 ^м