Асимптотика решения первой краевой задачи для сингулярно возмущенного дифференциального уравнения в частных производных второго порядка параболического типа

Автор: Кожобеков Кудайберди Гапаралиевич, Шооруков Асылбек Абдибахапович, Турсунов Дилмурат Абдиллажанович

Журнал: Вестник Южно-Уральского государственного университета. Серия: Математика. Механика. Физика @vestnik-susu-mmph

Рубрика: Математика

Статья в выпуске: 1 т.14, 2022 года.

Бесплатный доступ

Строится полное равномерное асимптотическое разложение по малому параметру решения первой краевой задачи. Первая краевая задача ставится для сингулярно возмущенного линейного неоднородного дифференциального уравнения в частных производных второго порядка с двумя независимыми переменными параболического типа. Задача исследуется на прямоугольнике. Особенности задачи - присутствие малого параметра перед оператором теплопроводности, существование угловых пограничных слоев на нижних углах прямоугольника. Требуется построить равномерное асимптотическое разложение решения первой краевой задачи на прямоугольнике, с любой степенью точности, при стремлении малого параметра к нулю. Асимптотическое разложение решения по малому параметру строится методом Вишика-Люстерника. При решении поставленной задачи нами используются: методы интегрирования обыкновенных дифференциальных уравнений, классический метод малого параметра, метод пограничных функций Вишика-Люстерника и принцип максимума. Как обычно, задача решается в двух этапах: в первом этапе строится формальное разложение решения первой краевой задачи, а во втором этапе оценивается остаточный член полученного разложения и этим доказывается, что полученное разложение действительно является асимптотическим на всем прямоугольнике. В первом этапе формальное асимптотическое решение ищется в виде суммы шести функций (решений): внешнее решение, определенное на всем прямоугольнике, погранслойное решение в малой окрестности нижней стороны прямоугольника, два боковых погранслойных решения в малой окрестности боковых сторон прямоугольника и два угловых погранслойных решения в окрестностях нижних вершин прямоугольника. Все эти погранслойные решения экспоненциально убывают вне пограничных слоев.

Еще

Асимптотическое решение, малый параметр, сингулярновозмущенная задача, первая краевая задача, уравнение теплопроводности, погранслойное решение

Короткий адрес: https://sciup.org/147235833

IDR: 147235833   |   DOI: 10.14529/mmph220103

Текст научной статьи Асимптотика решения первой краевой задачи для сингулярно возмущенного дифференциального уравнения в частных производных второго порядка параболического типа

Постановка задачи . В прямоугольнике рассмотрим первую краевую задачу [1, 2]:

d z ( t , x )     2 d z ( t , x )      2  , x 6z ( t , x )     z X Z X Z-Z X Z X

s-- a 2-----7.— + s 2 p ( t , x )------ - + q ( t , x ) z ( t , x ) = f ( t , x ),   ( t , x ) gQ ,

( dt         dx2  J             dx

z (0, x ) = ф (x ),   x g [0,1],

z(t,0) = Д1(t), z(t,1) = д,(t),  t g[0, T], где s - малый положительный параметр, 0< a - const, Q = {(t,x)|0 < t < T,0 < x < 1}, q, p, f g C " (Q), Фе C "[0,1], ^, д, g C "[0, T ], q (t, x) > 0: (t, x) gQ, ^(0) = ^1(0), ^(1) = ^,(0).

Требуется построить равномерное асимптотическое приближение первой краевой задачи (1)– (3), при s ^ 0.

Решение задачи (1)–(3) будем искать в виде рядов по степеням малого параметра [3–5]:

z (t , x ) = U ( t , x ) + V ( t , x ) + П 1 ( t , П 1 ) + П 2 ( t , П 2 ) + W 1 ( Т , П 1 ) + W 2 ( т , П 2 )            (4)

где t = 11 s , n 1 = x I X , n 2 = (1 - x )/ X , X = V J .

Подставляя соотношение (4) в задачу (1)–(3), получаем следующие задачи:

дU(t,x)    ?d2U(t,x))    2 / хдU(t,x)     ....

s--a1------— 1 + s2p(t,x)-------- + q(t,x)U(t,x) = f (t,x), (t,x) gQ;(5)

( д t           дx2   Jд дV(t,x)     2 д2V(t,x)    2 / хдV(т,x)

--sa2------—- + s2p(t, x)-------- + q(t, x)V(t, x) = 0, (t, x) g Q1;          (6) дт           дx2д

7 д 2 n j ( t, n j )                             9дп j ( t, n j )     3          дп j ( t, n j )

a2----. 2     - q(t, nj )пj (t, nj ) = A----5-----+ AP(t, nj )---а------, (t, nj ) g Q2j;

дп2                                 д t              Jдn д Wj (T,n j)    2д 2Wj (T,n j)                            4         д Wj (т,п j)

——~-a2---a 2    + q(т,nj)Wj(т,nj) = A 'P(т,nj) — ---—, (т,nj) gQ3j, дт           дn2                                    дnj где           Q1 = {(t, x )|0 < т < ц-T, 0 < x < 1},           Q2 j = {Му )|0 < t < T, 0 < n, < A"1},

Q 3 j = {( т , П j )l 0 т T s -1, 0 < n j A"1 }, j = 1,2.

Подставляя соотношение (4) в начальное условие (2), получаем:

^ ( x ) = U (0, x ) + V (0, x ) + П 1 (0, n 1 ) + П 2 (0, n 2 ) + W (0, n 1 ) + W 2 (0, П з ) ^

V (0, x) = ^( x) - U (0, x),(9)

Wj(0,nj) =-Пj(0,nj), j = 1,2.(10)

Теперь подставляя соотношение (4) в граничные условия (3), имеем:

M 1 ( t ) = U ( t , 0) + V ( t s -1 ,0) + П 1 ( t , 0) + П 2 ( t , a 1 ) + W 1 ( t s -1 ,0) + W 2 ( t s -1 , A ^1 ) , учитывая, что П 2( t , A -1) = 0, W 2( t s -1 , A -1) = 0 - условие для погранслойных функций, получаем:

П1( t ,0)=м t) - u (t ,0),(11)

W1(t,0) = -V (t,0),(12)

аналогично

^ 2 ( t ) = U ( t , 1) + V ( t s -1 ,1) + П 2 ( t , 0) + W 2 ( t s -1 ,0) ^

П1( t ,0) = ^2( t) - U (t ,0),(13)

Wi(t,0) = -V (t,1) .(14)

В результате мы получили шесть задач:

– из уравнения (5) методом малого параметра однозначно определяем U ( t , x );

  • -    из (6) и (9) определяем V ( т , x );

  • -    из (7), (11) и (13) определяем П j ( t , n j ), j = 1, 2;

  • - из (8), (10), (12) и (14) определяем W j ( т , п , ), j = 1, 2.

Начнем с уравнения (5).

Лемма 1. Для решения U(t, x) уравнения (5) справедливо формальное асимптотическое раз- ложение

U ( t , x ) = s k u k ( t , x ),                                        (15)

k =0

где U k g C ( Q ), k = 0,1,2,... - конкретизируются при доказательстве леммы 1.

Доказательство. Формально подставляя (15) в (5) и приравнивая коэффициенты при одина- ковых степенях е, имеем:

u k ( t , x ) =

1 q ( t , x )

д u k -1 ( t , x )

д t

f ( t , x ) u 0 ( t , x ) =   2 A ,

q ( t , x )

- a

2 д 2 u k -1 ( t , x )

д x 2

p t z . x ) a u k -i tM), k g n , u 1( t , x ) , 0, q ( t , x )      д x                       ’

заметим, что u k g C ( Q ), k = 0,1,2,.... Лемма 1 доказана.

Кожобеков К.Г., Шооруков А.А.,                   Асимптотика решения первой краевой задачи

Турсунов Д.А.                  для сингулярно возмущенного дифференциального уравнения…

Перейдем к задаче (6), (9). Пусть

V ( т , X ) = ]^ skV k ( т , X ),                                      (18)

k =0

где v k ( т , x) - пока неизвестные функций.

Подставляя (18) в (6) и (9), имеем:

k дvk ( т , X )      2 д V, ( т , X )    2V* j г х^- - j , X )   X" j г х         х п

^g  —-sa --—--+ g ^гP(x)-------^Д/р(x)vk-j(т,x) = 0, k=0   (   дт             dx2        j=o             dx      j=o

E sk Vk (0, x) = ф( x) - £ gkuk (0, x), k=0

где q j ( x ) = д j q (0, x )/ дт , q0(x ) = q (0, x ), p j ( x ) = д j p (0, x )/ дт .

Отсюда следует, что

^(Px) + q0(x)V0(т,x) = 0 V0(0,x) = ф(x)-u0(0,x);(19)

дт дvk (т, x)      , ..

—k-----+ q 0 (x) vk (т,x) = дт

= a 2 д vk -H ^ x ) + ^ т j P j ( x )   k - j -2 ( , ) + ^ т j q j ( x ) v k j ( т , x ), v k (0, x ) = _ u k (0, x^ k e N .

Решение задачи (19) существует, единственно и представимо в виде: v 0 ( т , x ) = ( ф ( x ) - u 0 (0, x )) e - q ( x ) т .

Для v 1( т , x ) имеем:

д v l ( т , x ) . z x z_ x 2 д2 v o( т , x ) . _ Z X        X ZA X       ZA X

i+ q 0 ( x ) V i ( т , x ) = a ----+ т q 1 ( x ) v 0 ( т , x ), v 1 (0, x ) = - u 1 (0, x ).

дтд

Справедлива

Лемма 2 . Решение задачи

^v^Ii^) + q0 (x)V(т, x) = e"q0 (x)т (P0 (x) +... + Pn (x)тп ), (т, x) e Qi, V(0, x) = V0 (x), x e [0,1 ] дт                               x1

существует, единственно и представимо в виде

V ( т , x ) = e q 0 ( x ) r V 0 ( x ) + e

f               T            n+ +1 )

--q0 (x '7  P0 (x)т + P1(x) —- + ... + Pn (x)----- где q0(x) > 0 x e[0,1], q0,pj,V0 e C”[0,1].

Доказательство . Уравнение

9 V ( т , x ) дт запишем в виде

q 0 ( x ) V( т , x ) = e q 0 ( x ) т ( P 0 ( x ) + P 1 ( x ) т + ... + P n ( x)T ) , ( т , x ) eQ 1 ,

(V(т,x)eq0(x)т ) ^ = (P0(x) + P1 (x)r +... + Pn(x)T ), полученное выражение интегрируем по т, учитывая начальное условие:

т

V(т,x)eq0(x)т -V0(x) = j(P0(x) +... + Pn(x)sn)ds  ^ V(т,x) = e"q0(x)тр0(x) + Pn+1,т(т,x)e"q0(x)т, где Pn+1т(т,x) = P0(x)т + P1(x)— +... + Pn(x)тп+1 /n +1. Лемма 2 доказана.

С помощью леммы 2 доказывается существование и единственность решений задач (20). Кроме этого, из леммы 2 следует, что эти решения экспоненциально стремятся к нулю при стремлении т к бесконечности, т. е.:

v k ( т ,x ) = O ( e q 0( x ) т ) , т ^ro , q 0( x ) 0: x e [0,1], k = 0,1,2,...

Перейдем к задачам (7), (11) и (13). Здесь две задачи относительно функций П1^ , 7 1 ) и П 2(t , 7 2) — аналогичные, поэтому достаточно рассмотреть одну. Мы рассмотрим задачу относительно функций П 1 ( t , 7 1 ):

2 д 2 П1( t 7j)         All z. .   .2 5П] (t ,7)   ,)4 Z, х^П (t ,7)    z. z о a ---— q(t,71)П1(t,71) = 2---V—- + ^p(t,71)---1г—-,  (t,71) е^21, дП2                           ^t

П1 (t, 0) = ^1 (t) - U(t, 0), t e [0, T].(22)

Пусть

П1( t ,71)  L 2П k (t ,71),(23)

к =0 где п 1 k (t , 7 1 ) — пока неизвестные функций, причем lim п 1 k ( t , 7 1 ) = 0, t e [0, T ]. ,                                                                                                      7 1 >z ,

Подставляя (23) в (21) и (22), имеем:

лк 2д П 1, k ( t , 7 1 )         j zA            x        ^+2 d n 1, k ( t , 7 1 )            j ( A^k-j k - j ( t , 7 1 )

L2 a —r^l/(tПk-j(t>71) = L2  —^—+2 ^7^(t)—i, k=0   (       д71       j=0                    J  k=0     (    dt          j=0              dxJ

L ^Xк (t,0) = A(t) - L 22kuk(t,0), t e [0, T], к=0

. .    дjq(t,0)       , .    дjp(t,0)      , ., где qj(t)=———, p} (t) =      j , q0 (t) = q(t,0) • j       дп1      j       дп/

Отсюда для n 1, 0( t, п 1 ) имеем:

д П 10

a ---- к" - q 0 ( t ) п ю = 0, ( t, П 1 ) е^ 21 , П ю( t ,0) = р х (t ) - u 0 ( t ,0), lim п ^( t, П 1 ) = 0, t e [0, T ].

дп                ,                                ,                                   П 1 >z

Решение этой задачи можно представить в виде

4q 0 ( t )

n l

П1,0( t ,П1) = ( ^1( t) - u 0( t ,0)) e   a

Для n 1 k (t , n 1 ), k e N имеем:

? д 2 n 1 k                                                         д п 1 k - j - 2( t, п ? )  д п 1 k - 2( t , n 1 )

a -^г-q0(t)п1,k =Ln1jqy(t)п1,k-j(t,П1)+Ln\Pj(t)—^------+—^—L,(t,П1) e^21, дП1              j=1                      j=0               дx

П 12 k (t , 0) = - U k ( t , 0), П 12 k -1 ( t , 0) = 0, lim П 1 k ( t , П 1 ) = 0, t e [0, T ].

,                                 ,                      П1 ^да

Справедлива

Лемма 3. Решение задачи

2                        ^q 0 ( t )

д П a -—q-q0(t)>n = e a (P1(t)71 +... + Pn(t)П1 ), (t,71)eQ21, дП1

n ( t ,0) = n 0 ( t ), lim n ( t , 7 1 ) = 0, t e [0, T ],

71 ^ro существует, единственно и представимо в виде

- ( t ) 7 1           - ( t ) 7 1

n( t ,71) = e   a n( t) + e   a    71 ( p1( t 7 +... + pn (t ))T]n ), где q0(t) > 0 t e[0,T], q0,pj,n° e Cro[0,T].

Лемма 3 доказывается прямым интегрированием, как и лемма 2.

Кожобеков К.Г., Шооруков А.А.,                   Асимптотика решения первой краевой задачи

Турсунов Д.А.                 для сингулярно возмущенного дифференциального уравнения…

С помощью леммы 3 доказывается существование, единственность решений задач (24). Для решений задач (24) справедливы оценки:

n (t , n 1 ) = O e

' q 0 ( t)

------П 1

, n 1 ^x , t g [0, T ], 0 a , 0 q 0( t ).

Перейдем к задачам (8), (10), (12) и (14) для определения угловых погранслойных функций W j ( т , п j ) , j = 1, 2.

Здесь тоже достаточно рассмотреть одну из них, либо задачу для W1( т , n 1) , либо для W2( т , n 2) , второе исследуется аналогично. Рассмотрим задачу для W1( т , n 1) :

д W ( т , п )    2 d2 W i( т , п )    , XTT,z х       / x^ W( т , П )             х z-x

—я—L - a 2---я 2 1 + q (т,П1Ж(т,П1) + £ p (т,П1)—г---L = 0,  (т,П1) 6 Он, дт          5^1

W1(0,rK) = -П1(0,П1),  W1(t,0) = -V (т,0).(26)

Пусть

X ^Ст,П1) = ^ Zkw1.к(Т’П1), к=0

где w 1, к ( т , п 1 ) - пока неизвестные функций.

Подставляя (27) в (25) и (26), получаем: д^    . д2 w_

—a       + q 0 w j = 0, У ,^) е^ з1 , j = 0,1;

дтдп

5w1,к    2 52 wк .         _ Z                         X   Z    X ~.

—--a — + q 0 w 1 к =ф к ( w 1 к -1 , w 1 к-2 ,..., W 10 , T , П 1 ),   ( т , П 1 ) 6^ 31 , к = 2,3,...

дт5п w1,к (0, П1) = -п1,к (0, П1), w1,2к (т,0) = -vk (t,0) w1,2к+1(т,0) = 0 к = 0Л> —, дn+mq(0,0)

где q =----------, 0 < q0 = q(0,0), Ф* - линейно зависят от предыдущих wXj (j

производных от т/1, полиномиально зависят от т, и щ•

Если ввести обозначение w1,к (т,п1) = e-q0тук (т,п1), то ew1, к (т,^1) _ _ а - q 0т            - q0т дУк (т,П1) . д2 w1, к (т,^1) _ - q 0т д" Ук (т,П1)

=  q0e    ук (т,Ч1) + e         _      ;      _ 2     = e         _ 2

дт                                 дт         дп^2

и рассматриваемая задача примет вид:

дУк     2 д2Ук                                        х Z     X                                     ZZXZXX

— -a —^ = Фк(У1 к-1,У1 к-2,...,У\0-т,П1)-   (тП1) е^з1, к = 0,1,...

дт       дп2           ,       ,,

Ук(0,П1) = -п1,к (0,П1),У2к(т,0)= -eq^'vk(т,0), У2к+1(т,0)= 0, к = 0,1,... где Фj = 0, j = 0,1.

Решения задач (28), (29) существуют, единственны и представимы в виде [6]:

X                         тт

Ук (т, П1) = - j П1, к (0, ^) G (т, ^, П1) d^ + j Ук (t) H (П1, т -1) dt + j j Ф к (t, x) G (П1, x, т -1) dxdt, 0                          000

где

H (Л\-т) =

П1     e~4a?

2 a 4пт3/2

G (т, x ,П1) =

2 a 4пт

(    ( exp

(П- x)21

4a2т J

- exp

_ (П1 + x)  1

, 2          I

V    4 a2т  J J

У2 к (т) = -eq°тук(т,0), У2 к-1(т) = 0,  Фк -Пх) = Фк ( У1, к -1, У1, к-2,..., У1,0,т, П1).

Отсюда следует, что функции w1,к (т,п1) экспоненциально убывают при т + п1 ^да

Таким образом, нами построены все функции, входящие в правую часть равенства (4).

Обоснование. Оценим остаточный член разложения (4).

Пусть z(t, x) = zn (t, x) + R(t, x), где R(t, x) - остаточный член разложения, zn (t,x) = Un (t,x) + Vn (t, x) + П1,2n+1(t ,П1) +П2,2 n+1(t, n2 ) + W1,2 n+1(^,П1) + W2,2n+1UnU n                             n2

Un (t,x) = E ^uk(t,x)’ Vn (r’x) = E ekvk(т,x)’ Пj,2n+1(t,Пj) = E ^nj,k(t, Пj), k=0                       k=0

2 n+1

Wj,2n+1( T nj ) = E ^wj,k(r,Пj ) . k=0

Тогда для остаточного члена получим следующую задачу: 2nz, dR(t,x)     2 d R(t,x)      2  , XSR(t,x)     , x _, x        n+1x ,,

£--a2-----т— + £ 2p (t, x)------- + q (t, x) R (t, x) = O (e + ),  (t, x) eQ,(31)

(  dt          dx2   Jd

R(0,x) = O(£n+1),  xe[0,1], z(t,0) = z(t,1) = O(£n+1), t e[0,T], £^0.(32)

Применяя принцип максимума [7], получаем:

IR(t,x)| < max {q '(t,x)O(£n+1), O(£n+1)}. (t,x)eQ>

0<£<<1

Отсюда имеем:

R(t, x) = O(£n+1),  (t, x) e Q, £ ^ 0.

Теорема. Для решения задачи (1)-(3) при стремлении малого параметра к нулю в области Q справедливо асимптотическое разложение

z(t,x) = Un (t,x) + Vn (t, x) + П1,2n+1(t,n1) + П2,2n+1(t, n2) + W12n+1(^,П1) + W2,2 n+1(^,n2) + O(£   ), где функции Un (t,x) = Ek=0 £kuk (t,x), Vn (T, x) = Ek=0 £kvk (T,x), П j ,2 n+1( t ,nj ) = E 2=01 2knj, k (t, nj ), Wj ,2 n+1(т,П, ) = Ek =+1 ^wj, k (T^nj ) определены выше.

Заключение. Нами построено полное равномерное асимптотическое разложение по малому параметру решения первой краевой задачи для сингулярно возмущенного линейного неоднородного дифференциального уравнения в частных производных второго порядка с двумя независимыми переменными параболического типа. Доказано, что полученное разложение действительно является асимптотическим решением поставленной задачи на всем прямоугольнике. Данная работа для нас является началом исследования бисингулярно возмущенных задач параболического типа, в следующих работах, ссылаясь на эту статью, мы будем исследовать только бисингуляр-ные случаи.

  • 1.    Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. – М.: Наука, 1966. – 724 с.

  • 2.    Zauderer, E. Partial Differential Equations of Applied Mathematics / E. Zauderer. – New York etc.: John Wiley & Sons, Inc. – 891 p.

  • 3.    Алымкулов, К. Об одном методе построения асимптотических разложений бисингулярно возмущенных задач / К. Алымкулов, Д.А. Турсунов // Известия вузов. Математика. – 2016. – № 12. – С. 3–11.

  • 4.    Вишик, М.И. Регулярное вырождение и пограничный слой для линейных дифференциальных уравнений с малым параметром / М.И. Вишик, Л.А. Люстерник // УМН. – Т. 12, вып. 5(77). – С. 3–122.

  • 5.    Треногин, В.А. Развитие и приложение асимптотического метода Люстерника–Вишика / В.А. Треногин // УМН. – 1970, Т. 25, вып. 4(154). – С. 123–156.

  • 6.    Полянин А.Д. Справочник по линейным уравнениям математической физики. – М.: Физ.-мат. лит., 2001. – 575 с.

  • 7.    Protter, M.H. Maximum Principles in Differential Equations / M.H. Protter, H.F. Weinberger. – Prentice Hall, Englewood Cliffs, New Jersey, 1967, 261 p.

Кожобеков К.Г., Шооруков А.А.,                   Асимптотика решения первой краевой задачи

Турсунов Д.А.                 для сингулярно возмущенного дифференциального уравнения…

Поступила в редакцию 27 декабря 2021 г.

Турсунов Дилмурат Абдиллажанович – доктор физико-математических наук, директор ВШМОП, Ошский государственный университет, г. Ош, Кыргызская Республика, е-mail: , ORCID iD:

Bulletin of the South Ural State University Series “Mathematics. Mechanics. Physics” 2022, vol. 14, no. 1, pp. 27–34

Список литературы Асимптотика решения первой краевой задачи для сингулярно возмущенного дифференциального уравнения в частных производных второго порядка параболического типа

  • Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. - М.: Наука, 1966. - 724 с.
  • Zauderer, E. Partial Differential Equations of Applied Mathematics / E. Zauderer. - New York etc.: John Wiley & Sons, Inc. - 891 p.
  • Алымкулов, К. Об одном методе построения асимптотических разложений бисингулярно возмущенных задач / К. Алымкулов, Д.А. Турсунов // Известия вузов. Математика. - 2016. - № 12. - С. 3-11.
  • Вишик, М.И. Регулярное вырождение и пограничный слой для линейных дифференциальных уравнений с малым параметром / М.И. Вишик, Л.А. Люстерник // УМН. - Т. 12, вып. 5(77). - С. 3-122.
  • Треногин, В.А. Развитие и приложение асимптотического метода Люстерника-Вишика / В.А. Треногин // УМН. - 1970, Т. 25, вып. 4(154). - С. 123-156.
  • Полянин А.Д. Справочник по линейным уравнениям математической физики. - М.: Физ.-мат. лит., 2001. - 575 с.
  • Protter, M.H. Maximum Principles in Differential Equations / M.H. Protter, H.F. Weinberger. - Prentice Hall, Englewood Cliffs, New Jersey, 1967, 261 p.
Статья научная