Сингулярно возмущенное уравнение в частных производных первого порядка в случае смены устойчивости Текст научной статьи по специальности «Математика»

УДК 517.956.226

СИНГУЛЯРНО ВОЗМУЩЕННОЕ УРАВНЕНИЕ В ЧАСТНЫХ ПРОИЗВОДНЫХ ПЕРВОГО ПОРЯДКА В СЛУЧАЕ СМЕНЫ УСТОЙЧИВОСТИ

E. A. Деркунова

SINGULAR PERTURBED PARTIAL DIFFERENTIAL EQUATION OF FIRST ORDER IN A CASE OF STABILITY CHANGE

E. A. Derkunova

Формулируется и доказывается теорема о дифференциальных неравенствах для уравнений в частных производных первого порядка. Иллюстрируется ее применение к начальной задаче с внутренним переходным слоем.

Ключевые слова: сингулярные возмущения, асимптотический метод дифференциальных неравенств, смена устойчивости

The theorem on differential unequalities for partial differantial equations of the first order is formulated and proofed. Its applicaition for initial problem with inner transition layer is illustrated.

Keywords: singular perturbations, asimptotic method of differential unequalities, stability change

Введение

После основополагающих работ академика А. Н. Тихонова [1] и последовавших за ними работ А. Б. Васильевой [2] ведется активное исследование сингулярно возмущенных задач асимптотическими методами. В частности, рассматриваются задачи, где нарушается условие изолированности корня вырожденного уравнения. Поясним это на примере начальной задачи в скалярном случае:

= Р(и0 <* - Т} и(°) = и°'

Пусть корни вырожденного уравнения

^(и, 0) = 0

неизолированы, например, корней два (и = и и = <^(£)), и графики их пересекаются во

внутренней точке рассматриваемого отрезка [0, Т]. Кроме того, пусть по прохождении точки пересечения корней они меняются ролями в отношении устойчивости (происходит смена знака производной Еи, взятой на каждом из корней, или, как говорят, происходит смена устойчивости). Возникает вопрос: как будет вести себя решение нашей задачи при е ^ 0? В работе [3], а затем в работе [4] для тихоновской системы доказана при определенных условиях теорема о предельном переходе при е ^ 0 от решения исходной задачи к решению

вырожденной задачи, которое строится с использованием устойчивого составного (вообще говоря, негладкого) корня вырожденного уравнения. Для доказательства существования решения и предельного перехода в работе [4] был применен метод дифференциальных неравенств. В последующие годы этот метод использовался для целого ряда других сингулярно возмущенных задач, в том числе задач со сменой устойчивости.

1. Метод дифференциальных неравенств для уравнений в частных производных первого порядка

Пусть рассматривается уравнение

£г(ди + Л(х,*)= / (и,х,*,£), (х,*) € Б (1.1)

с начальным условием

и(х, 0,£) = и°(х), 0 < х < 1, (1.2) £ > 0, г > 0. Решение ищется в области

Б = {(х, *) : х°(*) < х < Ж!(¿), 0 < * < Т}, (1.3)

где х = х°(*) и х = Ж!(¿) - характеристики, выходящие соответственно из точек (0, 0) и

¿Ж

(1, 0) и определяемые уравнением — = Л(х,*) (считаем, что Л(х,*) - гладкая функция и

все характеристики, выходящие из точек начального отрезка {* = 0, 0 < х < 1} существуют при 0 < * < Т).

1.1. Лемма о дифференциальных неравенствах

Лемма. Пусть функции и и v непрерывны и имеют кусочно непрерывные частные производные, а функция /(и, х, £) непрерывно дифференцируема, и, кроме того,

(Яи ди \

— + Л(х, *) — ) - /(и,х,*,£) < 0, (х,*) € Б;

(2) > 0, (х,*) € Б;

(3) и(х, 0) < v(x, 0), 0 < х < 1.

Тогда имеет место неравенство

и(х,*,£) < v(x,í,£) (х,*) € Б. (1.4)

Доказательство. В силу условий на функцию /(и, х, £) (непрерывной дифференцируе-мости) найдется постоянная N > 0, такая что

|/(и!, х, £) — /(и2, х, £)| < N|и! — и2|

при (х, ¿, £) € Б х (0; £°], |и| < М. Предположим, что

и(х,*,£) >v(x,t, £) (1.5)

в некоторой точке (х°,*°) € Б. Рассмотрим характеристику х = х(*), определяемую урав-дх

нением — = Л(х, *) и проходящую через точку (х°,*°). На этой характеристике найдется точка (х!,¿1), *! < такая что имеет место равенство и(х!,*!,£) = v(x!,*!,£), а в точках

(х,£) характеристики х = х(£) между точками (ж1, ^ 1) и (хо,£о) выполняется неравенство «(ж, е) > -и(х, ¿, е).

Выберем на рассматриваемой характеристике точку (х2,^2), ¿2 > ¿1, для которой

¿2

У ЖШ < е'.

¿1

Обозначим (хз, ¿3) точку, доставляющую максимум разности «(ж, е) — ^(х, ¿, е) = ¿(ж, е) > > 0 при ¿1 < £ < ¿2 и х = х(£). Отметим, что .(х^^е) = 0. Получаем следующую цепочку неравенств (интеграл вычисляется вдоль характеристики х = х(£)):

¿3

I

е'*<хз,<з,е) = /е'((| + Л(х,()-(| +А(х,()£))* <

¿1

¿3 ¿3

</(^(«,х(^),^,е) — / (V, х(£), ¿, е)^ ^ Жг(х(£), е)^ <

¿1 ¿1

¿2

¿1

Полученное противоречие доказывает ошибочность предположения (1.5), что приводит к утверждению (1.4).

1.2. Теорема о нижнем и верхнем решениях

Определение. Функции и(х, е) и и(х, ¿, е) называются нижним и верхним решениями задачи (1.1), (1.2), если выполнены неравенства

10. = е' (д| + Л(х, ¿)Ц) — /(и, х, е) < 0 < ^е(й),

_ (х, ¿) € Д

20. и(х, 0, е) < и0(х) < и(х, 0, е), 0 < х < 1.

Нижнее и верхнее решения называются упорядоченными, если

и(х, ¿, е) < и(х, ¿, е), (х,£) € Д.

Из Леммы следует, что нижнее и верхнее решения и и и начальной задачи (1.1), (1.2) являются упорядоченными.

Теорема 1. Если существуют нижнее и верхнее решения и и и задачи (1.1), (1.2), то эта задача имеет решение и(х,£,е), удовлетворяющее неравенствам

и(х, ¿, е) < и(х, е) < и(х,£,е), (х,£) € Д. (1.6)

Доказательство. Из Определения следует, что для функции и выполнено условие (1), а для функции и - условие (2) Леммы. Рассмотрим вспомогательное уравнение

и < и,

и < и < и, (1.7)

и > и

< ди+Л(х,«) а=

/ (и, х, ¿, е), / (и, х, ¿, е), / (и, х, ¿, е),

с начальным условием (1.2), удовлетворяющем неравенствам 20 Определения.

Уравнения (1.7) и (1.1) совпадают в области, где и < и < и. Правая часть вспомогательного уравнения ограничена, поэтому у задачи (1.7), (1.2) существует решение, и если будет доказано, что это решение не выйдет из полосы, задаваемой функциями и и и, то отсюда последует существование решения задачи (1.1), (1.2), удовлетворяющего неравенствам (1.6). Доказательство проведем от противного. Предположим, что

и(х,.,е) >Щх,.,е) (1.8)

в некоторой точке (хо,.о) € Б. Рассмотрим характеристику х = х(.), определяемую урав-Лх

нением — = Л(х, ¿) и проходящую через точку (хо,.о). На этой характеристике найдется

точка (х!,^), ¿1 < ¿о, такая, что имеет место равенство и(х1,.1,е) = и(х^.^е), а в точках (х,.) характеристики х = х(.) между точками (х^^) и (хо,.о) выполняется неравенство и(х,.,е) > и(х,., е). Поскольку из (1.7) имеем равенство

КI +Л(х(» I) =/ <и-х.(--).

а из определения верхнего решения - неравенство

еЛ вт+Л(х-() ах)>/(и 'х-(-г)

на характеристике х = х(.) между точками (х1, ¿1) и (хо, ¿о), то получаем (интеграл вычисляется вдоль характеристики):

ег(и - и) |(,0.10.«, = ег / (+ Л(х,.)) Л. < 0 ,

¿1

но последнее противоречит (1.8) в области Б. Следовательно, и < и. Аналогично можно доказать, что и > и в области Б. Таким образом, решение (1.7), (1.2) удовлетворяет неравенствам (1.6), что и доказывает теорему.

2. Начальная задача для уравнения с малым параметром при производных

2.1. Постановка задачи и условия

Рассмотрим следующую задачу (е > 0 - малый параметр):

+ Л(х,.)^ = /(и,х,е), (х,.) € Б, (2.1)

и(х, 0,е) = ио(х), 0 < х < 1. (2.2)

Решение задачи (2.1), (2.2) ищется в области вида (1.3).

Хорошо известно [5, с.135 - 139], что если вырожденное уравнение

/(и, х, ¿, 0)=0, (2.3)

имеет корень и = <^(х,£), который устойчив в области Б, т.е. выполнено неравенство

/„(<р(х,.), х, 0) < 0, (х, ¿) € Б, (2.4)

и если начальная функция u0(x) принадлежит области притяжения этого корня, то решение u(x,t, е) задачи (2.1), (2.2) существует и удовлетворяет предельному равенству

lim u(x,t,e) = <^(x,t), xo(t) < x < xi(t), 0 <t < T.

£—>0

Более сложная ситуация возникает тогда, когда уравнение (2.3) имеет корни u = <^i(x, t) и u = (x, t), пересекающиеся по некоторой кривой, проекция которой на плоскость (x, t) лежит в области D. В этом параграфе будет рассмотрен случай, когда проекция пересечения корней расположена выше начального отрезка {t = 0, 0 < x < 1}. Пусть выполнены следующие условия.

Условие 1. Функция f (u, x,t,e) достаточно гладкая (дважды непрерывно дифференцируемая) в области G = Iu х D х [0,ео], где Iu - некоторый интервал ео > 0 - некоторое число.

Условие 2. Уравнение (2.3) в области 1и х Д имеет относительно и два корня и = ^(х^), и = ^2(х, ¿), имеющие ту же гладкость, что и /, и удовлетворяющие соотношениям

(х^) = (х^) при £ = ^(х), (2.5)

где ^(х) - гладкая функция, 0 < ^(х) < Т. Для краткости будем обозначать кривую t = ^(х) буквой Г. Образуются две подобласти области Д: , ограниченная характеристиками х = х^) и х = х2^), начальным отрезком {(х,^ : t = 0, 0 < х < 1} и кривой Г, а также подобласть ^2 =

(х, ¿) > ^2(х, ¿) при (х, ¿) € А\Г; (2 6)

(х, ¿) < ^2(х, ¿) при (х, ¿) € ^2\Г. ( . )

И пусть

fu(pi(x,t),x,t, 0) < 0, fu(^2(x,t),x,t, 0) > 0 fu(^i(x,t),x,t, 0) > 0, fu(^2(x,t),x,t, 0) < 0

при (x, t) € Di\r, при (x, t) € D2\Г.

Условие 2 означает, что корни и = ^(х, ¿) и и = ^(х, ¿) уравнения (2.3) пересекаются по кривой, проекция которой Г на плоскость (х, ¿) лежит в области Д и описывается уравнением t = ^(х).

Не ограничивая общности, можно считать, что функция, для которой выполняются Условия 1,2, представима в виде

/(и, х, ¿, е) = —/0(и, х, ¿) (и — ^1(х, ¿)) (и — ^>2(х, ¿)) + е/1(и, х, ¿, е), (2.7)

где /0 > 0, , ^2, /1 - достаточно гладкие, и ^1, ^>2 удовлетворяют (2.5), (2.6). С использованием корней и ^2 определим составные решения вырожденного уравнения (2.3):

U(x, t) = U(x, t) =

^i(x, t), (x, t) € Di, ^2(x,t), (x,t) € D2;

^2(х, ¿), (х, ¿) € (х, ¿), (х, ¿) € ^2.

Отметим, что составные корни и и и являются непрерывными, но, вообще говоря, негладкими на кривой Г. Из Условия 2 следует, что

U(x,t) >U(x,t), (x, t) € D\r;

и(х,.) = и(х,.), (х,.) € Г;

/и(и(х, .), х, 0) < 0, /и(и(х,.), х, 0) > 0, (х, ¿) € Б\Г; (2.8)

/и(и(х, .), х, ¿, 0) = /и(и(х,.), х, ¿, 0) = 0, (х, ¿) € Г. (2.9)

На основании неравенств (2.8) можно назвать корень и(х, ¿) устойчивым, а корень и(х, ¿) неустойчивым (см. (2.4)). Следует, однако, заметить, что в точках кривой Г неравенство /и(и(х, .), х, ¿, 0) < 0 не выполнено (см. (2.9)), и это обстоятельство не позволяет однозначно ответить на вопрос о том, будет ли решение и(х,.,е) задачи (2.1), (2.2) стремиться при е ^ 0 к составному корню и(х,.) в области Б (за исключением начального отрезка). Как оказывается, при Условиях 1, 2 и некоторых дополнительных условиях этот предельный переход будет иметь место.

Потребуем выполнения еще одного условия, которое обеспечивает принадлежность начальной функции ио(х) области притяжения составного устойчивого корня и(х,.).

Условие 3. ио(х) > и(х, 0) при 0 < х < 1, ио(х) € /и, где 1и - интервал из Условий 1, 2.

Для доказательства предельного перехода от решения задачи (2.1), (2.2) к устойчивому составному решению вырожденного уравнения (при е ^ 0) потребуются еще два условия:

Условие 4. /ии(и(х, ¿), х, ¿, 0) |(ж,г)ег < 0.

Заметим, что для функции (2.7) это условие выполнено, так как /ии(и, х, ^(х), 0) = —2/о(и(х,'0(х)),х,'0(х)) < 0.

Условие 5.

'дй . , . дй \ т + Л(М) дХ

- /£(ü(x,i),x,i, 0)|i=^(x) < 0.

Ниже будет показано, что справедливо следующее утверждение Теорема 2. При выполнении Условий 1-5 решение задачи (2.1), (2.2) существует, и имеет место предельное равенство

lim u(x,t,e) = ü(x,t) (2-10)

£—>0

для всех (x, t) € D кроме начального отрезка {t = 0, 0 < x < 1}.

При этом разность u(x, t, е) — ü(x,t) имеет порядок О(у^е) в малой окрестности кривой Г и порядок O(e) в остальной части области D (за исключением малой окрестности начального отрезка).

2.2. Асимптотическое поведение решения

Рассмотрим задачу в области, ограниченной характеристиками x = xo(t) и x = xi(t), начальным отрезком и кривой t = ^(x) — 5, где 5 > 0 - сколь угодно малое, но фиксированное при е ^ 0 число. Здесь применима стандартная теория [5]: уравнение (2.1) с начальным условием (2.2) имеет единственное решение, и для него справедливо асимптотическое представление:

u(x, t, е) = u(x, t) + üo(x, т) + O(e), (2.11)

где ü(x,t) = ^>i(x,t) - функция регулярной части асимптотики, no(x,T) - пограничная

функция, т =--погранслойная переменная. Для функции no(x,T) имеем задачу:

е

дПо ^ . ^ ______ „ ^ о,

дт

= /(u(x, 0) + no(x, т), x, 0, 0), no(x, 0) = u0(x) — u(x, 0).

В силу Условий 2, 3 функция П0 экспоненциально убывает при т ^ то.

Замечание 1. Заметим, что пограничная функция определена в прямоугольнике (0 < х < 1) х (0 < t < Т), который не совпадает, вообще говоря, с областью Д. Чтобы определить ее во всей области Д, поступим так. Продолжим гладким образом функции /(и, х, 0, 0), и0(х) и и(х, 0), входящие в задачу для пограничной функции, за границы отрезка 0 < х < 1, и рассмотрим эту задачу при Х0 < х < Х1, т > 0, где X = пип^ х0^),

X = шах^ х^). Тогда пограничная функция будет определена во всей области Д.

На кривой t = ^(х) — 5 из представления (2.11) имеем:

и(х, ^(х) — 5, е) = и(х, ^(х) — 5) + О(е), (2.12)

так как пограничная функция экспоненциально мала.

Рассмотрим теперь задачу в области Д;, ограниченную характеристиками х = х^), х = х2(^ и кривыми t = ^(х) — 5, t = ^(х) + 5.

Замечание 2. Кривые t = ^(х) — 5 и t = ^(х) + 5 определены при тех же х, что и кривая t = ^(х). Гладко продолжим их за область определения, чтобы их концы лежали на характеристиках х = х0(^ и х = х^). Построим нижнее и верхнее решения в виде:

и = и(х, ¿) — Ае и = и(х, ¿) + А^е,

где число А > 0 будет выбрано ниже. Проверим выполнение условий 10, 20 Определения применительно к области Д; для указанных функций. Имеем:

^=<1+Л(х,^) И)—/ (и,х^е)=

= е ( ди + Л(х, ¿) дх) — / (х, ¿) + /и(х, ¿) Ае — / (х, ¿)е + о(е) = = е^й + Л(х, ¿) дх — /е(х,^ — / (х,*) + /и(х, ¿)Ае + о(е).

Здесь / (х^) = /(«(х, ^х^, 0), и такой же смысл имеют обозначения /(х^), /(х^). Так как /(х^) = 0 и /^х^) < 0, и так как неравенство из Условия 5 выполняется на кривой t = ^(х), а значит ив 5 - окрестности нее, если только 5 достаточно мало, то для достаточно малых е

< 0, (х^) € Д;.

«? = ^+л(М) Ц) -/ (и,х,(,е) =

/д« „ , „ д«\ у, „ х , „ , ^ х , „ 1 у . . ,2

-/ии (х, ¿) А2--

2 ии

(ди д«\ 1

— + Л(х, ¿) — ) — / (х, ¿) — /и(х, ¿)А\/е — /£(х, ¿)е — 2 /ии (х, ¿)А2е + о(е),

Так как / = 0, /и < 0 и выполняется Условие 4, то при достаточно большом А и достаточно

малых е слагаемое — ^/^(х^А^ будет доминирующим и обеспечит выполнение неравенства

> 0, (х^) € Дй.

Таким образом, условие 10 Определения выполнено. Проверим выполнение условия 20 Определения, имея в виду, что в Д; начальным множеством является кривая t = ^(х) — 5. При достаточно большом А и достаточно малых е из (2.12) имеем:

(и — «) = — Ае + О(е) < 0

(и — и) |4=^(Л)-Й = А^е + О(е) > 0.

Тем самым условия Определения применительно к области Б; выполнены. Уравнение (2.1) с начальным условием (2.12) имеет решение и(х,£,е) в области Б;, и справедливо представление

и(х, е) = и(х, £) + о(^ё), (х, £) € . (2.13)

5

Рассмотрим область между характеристиками х = хо(£), х = х1(£), кривой £ = ^(х) +— (см. Замечание 2) и линией £ = Т .В этой области снова работает стандартная теория. Корень и(х,£) = ^>2(х,£) в силу (2.6) и (2.8) является изолированным и устойчивым. Как следует из (2.13), начальное условие для решения и(х,£,е) в этой области представимо так:

и ^х,^(х) + 2, е^ = и ^х,^(х) + 0 + О(^ё). (2.14)

Пограничная функция, как это следует из (2.14), имеет порядок О(у'е), однако, уже в области выше кривой £ = ^(х) + 5 в силу экспоненциального затухания По для решения и(х, е) справедлива асимптотика:

и(х, е) = и(х, £) + О(е).

Таким образом, доказана Теорема 3. Если выполнены Условия 1-5, то при достаточно малых е .задача (2.1), (2.2) имеет единственное решение и(х,£,е), и для него справедливо асимптотическое представление

и(х, е) = и(х, £) + По(х, £/е) + -ш(х, е),

где -ш(х,£, е) = О(уе) в 5 - окрестности кривой £ = ^(х), и -ш(х,£, е) = О(е) в остальной части области Б.

Из этого асимптотического представления решения задачи (2.1), (2.2) непосредственно следует предельное равенство (2.10).

Литература

1. Тихонов, А. Н. О зависимости решений дифференциальных уравнений от малого параметра / A. Н. Тихонов // Математический сборник. - 1948.- Т. 22(64), № 2.— С. 193 — 204.

2. Васильева, А. Б. Асимптотика решений некоторых задач для обыкновенных нелинейных дифференциальных уравнений с малым параметром при старших производных / А. Б. Васильева // УМН. — 1963. — Т.18, № 3.— С. 15 — 86.

3. Lebovitz, N. R. Exchange of stabilities in autonomous system — II. Vertical bifurcaition / N. R. Lebovitz, R. J. Schaar // Stud. Appl. Math.— 1977: V. 56. — P. 1 — 50.

4. Nefedov, N. N. Singularly perturbed systems: Case of exchange of stability / N. N. Nefedov, K. R. Schneider// Weierstrafi - Institut fur Angewandte Analysis und Stochastik, Berlin.— Preprint № 158.

5. Васильева, А. Б. Асимптотические методы в теории сингулярных возмущений / А. Б. Васильева, В. Ф. Бутузов. — М.: Высш. шк., 1990. — 208 с.

Кафедра функционального анализа, Южно-Уральский государственный университет derk@math.susu.ac.ru

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