Асимптотическое исследование нелинейных нелокальных моделей типа реакция - диффузия - адвекция с пограничными и внутренними слоями тема автореферата и диссертации по математике, 01.01.03 ВАК РФ

003466415

На правах рукописи

НИКИТИН АНДРЕИ ГЕННАДЬЕВИЧ

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

01.01.03 - математическая физика

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

О 9 АПР 2009

Москва 2008

003466415

Работа выполнена на кафедре математики физического факультета Московского государственного университета имени М.В.Ломоносова

Научный консультант доктор физико-математических наук профессор

Нефёдов Николай Николаевич

Официальные оппоненты: доктор физико-математических наук профессор

Галкин Валерий Алексеевич

доктор физико-математических наук профессор Дмитриев Михаил Геннадьевич

доктор физико-математических наук профессор Курина Галина Алексеевна

Ведущая организация: Ярославский государственный университет имени П.Г.Демидова

Защита состоится «23» ¿г^и-А^/Л 2009г. в « /-Г"» часов на заседании диссертационного совета Д 501.002.10 при Московском государственном университете имени М.В.Ломоносова по адресу: 119991, ГСП-1, Москва, Ленинские горы, МГУ имени М.В.Ломоносова, дом 1, строение 2, физический факультет, аудитория .

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ.

Автореферат разослан « » ^¿--¿-¿у^-чГ2009г.

Ученый секретарь диссертационного совета,

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

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

В теории сингулярных возмущений контрастные структуры ранее исследовались в нелинейных эллиптических краевых задачах с малыми параметрами при старших производных, рассматриваемых в ограниченных областях. Впервые существование контрастных структур в сингулярно возмущенных задачах для обыкновенных дифференциальных уравнений было доказано в работах А.Б. Васильевой и В.Ф. Бутузова [2-5], Результат по существованию двумерных контрастных структур типа ступеньки принадлежит П. Файфу (P. Fife) и У. Гринли (W. Greenlee) [6]. Асимптотические разложения решений типа контрастной структуры по малому параметру можно построить на основе метода пограничных функций [2, 32-35]. Для одномерных задач это сделано в [2-5,7, 3234] и ряде других работ, для некоторых двумерных задач - в [8,9, 35]. Обширная библиография по этой проблематике содержится в [7].

Важным вопросом, как с теоретической, так и с прикладной точки зрения, является вопрос об устойчивости контрастных структур как стационарных решений соответствующих параболических задач (в смысле Ляпунова). Для одномерных задач этот вопрос был решен в работе А.Б. Васильевой [10], В.Ф. Бутузова [10], работах [32-34, 36, 37], S. Angenent., J. Mallet-Paret и L. Peletier [12], J. Hale и К. Sakamoto [13] и др. Устойчивость периодической контрастной структуры типа ступеньки в пространственно двумерном случае была впервые получена в работе [35] путем исследования спектра сингулярно возмущенной двумерной задачи на собственные значения. Вопросы устойчивости и локальной единственности решений сингулярно возмущенных нелинейных эллиптических задач, а так же важная проблема формирования контрастных структур в сингулярно возмущенных параболических задачах, были решены В.Ф. Бутузо-

вым и И.В. Неделько с помощью предложенного ими метода параметрических барьеров [14, 15].

Наиболее эффективным методом доказательства существования контрастных структур и оценки остаточных членов асимптотических разложений является асимптотический метод дифференциальных неравенств H.H. Нефедова [7, 8]. Суть его состоит в том, что верхнее и нижнее решения конструируются путем модификации формальной асимптотики. Рассматривая эллиптическую задачу как стационарную задачу для соответствующего параболического уравнения, этим методом можно также доказать устойчивость по Ляпунову и локальную единственность решения исходной задачи.

В настоящее время большой интерес вызывают более сложные модели, которые включают эффекты обратной связи или нелокального взаимодействия. Различные направления теории нелокальных нелинейных моделей интенсивно разрабатываются как у нас в стране, так и за рубежом. Как правило, эти модели представлены сингулярно возмущенными интегродифференциальными уравнениями, описывающими важные для приложений процессы, в которых необходимо принять во внимание последствия или задержку, во многих областях естествознания, в частности, в задачах динамики реакторов, моделях генетики популяций, химической кинетике [16], теории фазовых переходов [17-19], социологии [31] и других областях [1, 20, 21]. Так модели, обладающие наследственными свойствами, описываются практически только интегродифференциальными уравнениями [22, 23]. В частности, например, в теории фазовых переходов при рассмотрении теоретической модели процесса разделения фаз в двойной полимерной смеси возникает следующая задача для уравнения Кана-Хилиарда (Cahn-Hilliard) [17-19]

где и(х,у,1) - концентрация одной из компонент' смеси, О - ограниченная область с гладкой границей, п - единичный вектор внешней нормали к границе области Ш, f{u) = W'{u), где IV(и) - двойная потенциальная яма, е - диапазон межмолекулярных сил. Простыми вычислениями [18] задача (1) сводится к задаче для нелокального уравнения реакция-диффузия

ди сп

и, = Д[/(и) - ¿?2Ли], (х,у) е Q, t > 0, ^г[/(м)-£2Ам] = 0, (х,у)едП, u(x,y,0) = g(x,y),

дп

(1)

гг Дм - f{u) + J/(M)c/Q, (х,у) е Q, О 0,

п

дп

>и

- = 0, (x,>')e3Q, u(x,y,Q)~ g{x,y).

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

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

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

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

Научная новизна. Все основные результаты работы являются новыми.

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

Апробация работы.

Основные результаты диссертации докладывались и обсуждались на семинаре МГУ по малому параметру (руководители: профессора А.Б. Васильева, В.Ф. Бутузов, H.H. Нефедов), на научных семинарах факультета ВМ и К (руководитель: профессор И.А. Шишмарев), кафедры математики физического факультета МГУ, НИВЦ МГУ (руководители: профессора А.Г. Ягола, А.Б. Бакушинский и A.B. Тихонравов), на международной конференции "Теория и приложения методов малого параметра", посвященной 90-летию со дня рождения академика А.Н. Тихонова (Обнинск, 1996), на международной конференции, посвященной 70-летию академика А.М. Ильина (Уфа, 2002), на международных конференциях "Nonlinear partial differential equations" (Алушта, 2003, 2005), на международной конференции, посвященной 100-летию A.A. Андронова (Нижний Новгород, 2001), на международных конференциях "Математические идеи П.Л.

Чебышева и их приложение к современным проблемам естествознания" (Обнинск, 2002, 2004, 2008), на седьмой Крымской международной математическая школе "Метод функций Ляпунова и его приложения" (Алушта 2004), на VI международном конгрессе по математическому моделированию (Нижний Новгород, 2004), на международной конференции "Tikhonov and contemporary mathematics" (Москва, 2006), на международной конференции "Дифференциальные уравнения и смежные вопросы", посвященной памяти И.Г. Петровского (Москва, 2007), на Тихоновских чтениях (Москва, 2002, 2006, 2008), на Ломоносовских чтениях (Москва, 2006, 2008), на вторых - шестых (1994-1998), восьмых (2000), десятых (2002), одиннадцатых (2003), пятнадцатых - семнадцатых (2006-2008) математических чтениях РГСУ (МГСУ), на международной конференции Science Links Moscow-Berlin-Paris Workshop 2008, посвященной 50-тию сотрудничества МГУ им. М.В. Ломоносова с Гумбольдтским университетом (Берлин, Германия).

Публикации Основные результаты, полученные автором и изложенные в диссертации, опубликованы в работах [32-46] (список литературы приведен в конце автореферата). По материалам диссертации опубликованы 15 научных работ и сделано 29 докладов на научных конференциях. Результаты, содержащиеся в работах, выполненных в соавторстве, и включенные в диссертацию, получены автором лично и включены в диссертацию с согласия и одобрения соавторов этих работ.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, разбитых на параграфы (некоторые параграфы, в свою очередь, разбиты на пункты), заключения и списка литературы, содержащего 70 наименований. Нумерация формул своя в каждом параграфе (пункте). В работе для формул принята двойная нумерация: первое число - номер параграфа (пункта), второе - порядковый номер формулы в параграфе (пункте). Объем диссертации составляет 198 страниц, включая 8 страниц цитированной литературы.

Содержание работы.

Во введении обосновывается актуальность темы диссертации, научная новизна полученных результатов, а также кратко изложено содержание и основные результаты работы.

В первой главе изучаются задачи Коши для сингулярно возмущенного ин-тегро-дифференциального уравнения типа Вольтерра и Фредгольма, которые описывают нелокальные модели, например, в биологии при описании распространения инфекционных болезней [22, 23]. Рассмотрен случай, когда вырожденное уравнение имеет изолированное решение,и случай, когда решение вырожденного уравнения имеет угловую точку (т.н. случай со сменой устойчивости). Результаты, представленные в данной главе, опубликованы в работах [43, 45]. При построении асимптотики решения используется метод пограничных

функций [2], в котором, в отличие от случая обыкновенных дифференциальных уравнений, вводятся дополнительные члены в разложение нелинейности уравнения, порождаемые интегральным оператором. Для обоснования результата используется развиваемый для нового класса задач асимптотический метод дифференциальных неравенств, предложенный в свое время Н.Н. Нефедовым для обоснования асимптотических разложений решений сингулярно возмущенных задач для обыкновенных дифференциальных уравнений и уравнений в частных производных [8,9].

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

£и' = 1(м,и,<,£), 0</<1, (3)

г

где ^g(u(í,£),v(s,s),t,s,£)ds + е>0-малый параметр,

о

и(0,£) = и°. (4)

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

Ци,и,(,0) = 0

имеет изолированное гладкое решение, и случай, когда есть пересекающиеся гладкие решения. В последнем случае рассматривается ситуация, когда происходит так называемая смена устойчивости, и решение вырожденного уравнения, приближающее решение исходной задачи, имеет, вообще говоря, угловую точку при некотором значении аргумента / — (0 £(0,1). Для случая со сменой устойчивости развиваются идеи изучения таких задач в различных классах дифференциальных уравнений (см. [24] и ссылки, приведенные в этом обзоре) на новый класс интегродифференциальных задач. Отметим, что в случае смены устойчивости решение имеет «слабый» внутренний слой (скачок производной), объясняющий скачок скорости химической реакции в интегро-дифференциальной модели реакции-диффузии аналогично тому, как это было в дифференциальных моделях. Первый параграф состоит из следующих разделов. В разд. 2 строится формальная асимптотика решения задачи (3), (4) в случае гладкого устойчивого решения, в разд. 3 проводится обоснование асимптотики для данного случая с помощью развиваемого для нового класса задач асимптотического метода дифференциальных неравенств [8, 9]. В разд. 4 § 1 рассмотрен случай со сменой устойчивости.

Для построения асимптотики решения применяется модифицированная схема метода пограничных функций [2], предложенная впервые в [38]. Асимптотика решения задачи строится в виде формального разложения

ОО

u(t,e) = Yl £"(«„(') + Пли(т)) = u(t,e) + Пи(т,г), (5)

(1=0

где йй(0 члены регулярной части асимптотики, П„и(т) — погранслойные члены асимптотики, г = t/e - погранслойная переменная. Модификация состоит в

I

том, что интегральный оператор I(u,t,e) & j g(u(t),u(s),t,s,e)ds нужно предо

ставить в виде, аналогичном (5)

1 = 7 + П/. (6)

Для этого подынтегральную функцию g нужно представить в виде более сложном, чем разложение (5). Это представление подробно описано в п. 2 § 1. Задача (3), (4) рассматривается при выполнении следующих условий:

Условие А1. Функции fug являются достаточно гладкими (При построении асимптотики с остаточным членом порядка 0(гл+1) достаточно потребовать, чтобы они были п + 2 раза непрерывно дифференцируемы).

Условие А2. Пусть вырожденное уравнение

L(u,u,t, 0) = 0 (7)

имеет изолированный корень u(t) = cp(t) на отрезке [0,1], и Lu((p{t\(p{s),t,Q)< 0 при t € [0,1] ^

I

где Lu(<p(t),<p(s),t,0) = fu(<p(t),t,0)+ jgu(<pit),<p(s),t,s,0)ds (условие устойчиво-

о

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

Условие A3. Начальное условие и0 принадлежит области влияния устойчивого корня (pit), т.е. задача

^ = /(й(г),0,0), г > 0, ат

имеет решение й = й(т), причем й(т) гр(0) при т -> so.

Условие А4. Пусть функция ?,.?,£) будет монотонно неубывающей

по переменной V в некоторой области изменения этой переменной при любом фиксированном значении переменной и из той же области изменения, что и для второй переменной, при (/,5) е [0,1] х [0,1] и достаточно малых в (так называемое условие квазимонотонности, область, о которой идет речь в данном условии, будет уточнена в Теореме 1.2).

Регулярный член асимптотики нулевого порядка является решением нелинейного интегрального уравнения

1(й0,г70,/,0) = 0 .

В силу условия А2 можно выбрать в качестве его решения г70 (/) = <р(0- По-гранслойный член нулевого порядка определяется из нелинейного дифференциального уравнения

Пои/ = Л9(0) + Пом(т),0,0), т>0,

с дополнительными условиями

П0и(0) = м° - <р(0), Псм(+°°) = 0.

Условие АЗ гарантирует разрешимость данной задачи Коши, а условие А2 -экспоненциальную оценку пограничной функции

|П0м(г)|<Сехр(-кт), г>0, (8)

где Сиг- некоторые положительные постоянные.

Регулярные члены первого и следующих порядков находятся из линейных интегральных уравнений

1

ьи((р,<р,1,0)щ (0 + 0)", 0)^ =Р, (0-

0

где /?.(?)определяются через уже известные члены асимптотики (¿ = 1,2,...). В силу свойств интегрального оператора Вольтерра эти уравнения однозначно разрешимы. Погранслойные члены первого и следующих порядков вблизи начальной точки ( = 0 определяются из линейных дифференциальных задач

П,-и' = /и(<р{0) + П0м(г),0,0)П,и +я"Дг), г>0, П(и(0) = -«((0),

где функции яг, (г) выражаются через уже определенные члены асимптотики (¿ = 1,2,...). Для погранфункций П,м(г)также справедлива экспоненциальная оценка, аналогичная оценке (8). Таким образом, построена формальная асим-

птотика решения задачи (3), (4) с остаточным членом порядка 0(£"n+1), где п -произвольное натуральное число.

Обоснование построенного асимптотического разложения решения задачи (3), (4) проведено развиваемым в диссертации асимптотическим методом дифференциальных неравенств. Приведем классическое определение верхнего и нижнего решения для задачи (3) ,(4) (см. [1]).

Определение 1.1. Функция fl(t,c) е С1 ((0,1]) П С([0,1]) называется верхним решением задачи (3), (4), если

£2p'-L{l3,p,t,£)> 0, Г €(0,1), J3(0,£)>u°.

Аналогично .функция a(t,e) называется нижним решением, если выполнены неравенства с обратным знаком.

При обосновании асимптотики используется следующая известная теорема о дифференциальных неравенствах [1].

Теорема 1.1. Пусть существуют функции a(t,e) и /?(/,£), являющиеся соответственно нижним и верхним решением задачи (3), (4), причем a{t,e)< P(t,e) при i с [0,1]. Пусть, кроме того, функция g непрерывна вместе со своими частными производными по первому и второму аргументу и удовлетворяет условию A4 на промежутке [а, 0]. Тогда супцествует единственное решение задачи (3), (4) u(t,s) такое, что

a(t,e)<u(t,s)<0(t,s).

На основе конструктивного метода построения нижнего и верхнего решений доказана следующая теорема:

Теорема 1.2. Пусть выполнены условия Л1 - A4. Тогда для достаточно малых £ существует единственное решение задачи (3), (4), удовлетворяющее неравенству

™х| u(t,e)-Un(t,e) |< Cf,+1, где U„(t,£) - частичная сумма порядка п асимптотического ряда (5). Приведем вид верхнего решения для случая п = 0

fi(t, £) = Ç(t) + Щ(1) + w(t)) + now(r) + m „ (т), где w(t) > 0 - некоторое положительное решение неравенства

Ьи{<р,<р,и 0М0+ <0, 0<г<1.

с

Показано, что положительные решения денного неравенства существуют при выполнении условия А4. Функция - модифицированная пограничная

функция.

В случае произвольного п соответствующие нижнее и верхнее решения о„(, и /?л+1 получаются из такой же модификации £/„м(г,е) - частичной суммы порядка п +1 асимптотического ряда (5).

Далее в первом параграфе рассмотрена задача (3), (4) в случае, когда вырожденное решение имеет на отрезке [0,1] два различных корня е/)(г) и ^(г), пересекающихся в некоторой внутренней точке отрезка (т.е. <р(гй) - ф{10), /0 е (0,1)).

Условие А5. Пусть

¥?(0><К') при 0</</0, </?(?)<0(0 при г0</<1, причем,выполнены следующие условия устойчивости корней

Определение 1.2. Решение й(1) называется устойчивым составным решением.

Таким образом, мы имеем вырожденное решение, имеющее, вообще говоря, угловую точку при / = . Наличие точки пересечения корней не позволяет построить асимптотику выше нулевого порядка для всех значений Г е[0,1]. Кроме того, нам понадобятся более сильные условия на интегральный оператор Ь в точке С0:

Условие А6. Пусть выполнено неравенство

4(<р,(р,1,0) + 4(<р,<р,1,0)<0при Ьи{(р,<р,(,0) + Ц((р,<р,(,0)>0 при (а<1<1, Ьи{ф,ф>и0) + Ц(ф,ф,1,0) > О при 0 < / < , 4 & Л 0)+4 ф, 1,0) < о при 10 < г < 1,

о»

о

Составим следующее решение вырожденного уравнения

[LJÛ(t,),û(s)Afi) + 2LJÛ(QMslt6,Q) + LJÛ(t0),û(s)J0M < 0 •

Условие A7. Пусть выполнено неравенство

«'(/„) -Le(û(Q,û(s),t0,0) - < 0,

о

= g(ç(t),m + nou(Ç),t,0,0)-g{ç{t),<P(0),f ДО).

Последнее неравенство следует понимать в смысле предела слева и справа при / -> г0. В качестве ujt) берется устойчивый вблизи начальной точки корень <p(t).

Теорема 13. Пусть выполнены условия Al —А7. Тогда для достаточно малых е существует единственное решение задачи (3), (4), удовлетворяющее неравенству

maxj u(t,£)-U0(t,£) I < С-Je,

m

где U0(t,£) = û(t) + П0и(г).

Приведем вид верхнего и нижнего решения для данного случая. Верхнее решение P{t, е) имеет следующий вид

P(t, £) = û(t) + П0«(г) + 4еу + sfell ( (г),

/

где у>0~ некоторая постоянная, а П, (т) - модифицированная пограничная

ïв

функция. Нижнее решение a{t,встроится в виде, отличном от верхнего решения

a(t, е) = й(0 + П 0«(г) - еу + £П1а. (г),

где П1а(т) - модифицированная пограничная функция . В §2 главы 1 рассмотрена задача Коши для сингулярно возмущенного обыкновенного нелинейного интегродифференциального уравнения Фредгольма

eu'= L(u,v,t,e), 0</< 1, (9)

и(0,е) = и\ (10)

]

где L(u,v,t,e) = ]#(«(*,£), v(i,г),¿.s,г)ds.

о

Асимптотика решения задачи строится методом, аналогичным примененному в §1 данной главы, в виде разложения (5). Будем считать, что выполнены условия В1 - В4, аналогичные условиям А1 - А4 §1. Кроме того, в случае интегрального оператора Фредгольмовского типа потребуется следующее дополнительное условие:

Условие В5. Пусть следующее неравенство имеет положительное решение

Л О-

I I

y{t)\gu(fP(f)Ms),t,sfi)ds + \gv(<p(t),<P(.slt,s$)y(s)ds < 0,0 < t < 1.

о о

Замечание. Достаточным условием для выполнения условия В5 является, например, требование

I I

+ J£Л<Р(!),<РШ^№ < О,

о о

в этом случае требуемым решением будет положительная постоянная, или требование на собственные значения следующего интегрального уравнения Фредгольма с иеотрш{ательиым (в силу условий В2 и В4) ядром K(t,s)

I

о

где

K(t j)-t

\gu{<P(t),(p{.s),t,s,0)ds

о

а именно, требуется, чтобы главное позитивное собственное значение А0 было меньше единицы (см. [25], стр.83, теорема 2.1).

Заметим также, что в §1 условие В5 для интегрального оператора Вольтерра было следствием условий А1 - А4.

В п. 2 §2 строится формальная асимптотика решения задачи (9), (10). Дня построения асимптотики решения применяется модификация метода пограничных функций, аналогичная предложенной в п. 2 § 1 данной главы. В п. 3 проводится обоснование асимптотики с помощью развиваемого для данного нового класса задач асимптотического метода дифференциальных неравенств [8, 9]. В п. 4 рассматривается случай, когда вырожденное уравнение (7) имеет непрерывное негладкое (имеющее скачок производной в некоторой точке t0 е (0,1)) решение. В отличие от случая оператора Вольтерра,вводятся дополнительные ограничения на нелинейность, и дается другое определение негладкого решения вырожденной задачи.

Как и в § 1,для обоснования построенного асимптотического решения ведущую роль играет развиваемый для нового класса задач асимптотический метод дифференциальных неравенств. Определение верхнего и нижнего решения для задачи (9), (10) аналогично определению 1 § 1. Также в рассматриваемом случае справедлива теорема о дифференциальных неравенствах, аналогичная теореме 1.1 § 1. На основе конструктивного метода построения нижнего и верхнего решений доказана следующая теорема:

Теорема 1.5. Пусть выполнены условия В1 - В5. Тогда для достаточно малых е существует единственное решение задачи (9), (10), удовлетворяющее неравенству

тах| ии,£)-ип(1,£) |<С£"+1,

где Цп(1,£) - частичная сумма порядка п асимптотического ряда (5). Приведем вид верхнего решения для случая п — 0

Р(1,£) = ср(1) + £(щ(1) + у({)) + П0м(т) + £П ^(г),

где >>(/) > 0 - некоторое положительное решение неравенства 1 1

о о

которое существует с силу требования В5.

В п. 4. § 2 рассмотрена задача (9), (10) в случае, когда вырожденное уравнение имеет на отрезке [0,1] два различных корня $о(?) и фи), пересекающихся в некоторой внутренней точке отрезка (т.е. ) = /0 е (0,1)).

Условие В2а . Пусть вырожденное уравнение 1

_[ £(м0 (0> «о С* )> ^ я,0)ск = 0

0

имеет составное решение следующего вида

Ш, о^, [<«/),/0</<1,

1

удовлетворяющее условию в,0)<1х<(!1 при 0<t<t0, г0 < ? < 1 (та-

0

кое решение мы будем называть устойчивым составным решением). Заметим, что составить решение й(1) из (р(1) и ф{1) возможно не для всякого вырожденного уравнения. В § 2 приведен пример уравнения, для которого это

возможно сделать. Таким образом, мы имеем вырожденное решение, имеющее, вообще говоря, угловую точку при t = /0. Наличие точки пересечения корней не дает возможности построить асимптотику выше нулевого порядка для всех значений / е [0,1], так мы не можем построить регулярные члены асимптотики выше нулевого порядка из-за неразрешимости интегрального уравнения для регулярных членов выше нулевого порядка при t>ta. Заменим требование В5 на следующее:

Условие П5а. Пусть при 0<t<t0u t0<t< 1 выполнено следующее неравенство

1 I

Jgu(fiQ)Msy,s,0)ds + jg„(u(t),ü(s),t,s,0)ds < 0.

о о

Кроме того, нам понадобятся более сильные условия на нелинейный интегральный оператор в точке ta:

Условие В6. Пусть выполнено неравенство

1 I

[ j 8ии ("<Л ),ü(s),ta,s,0)ds + 2 J gm (ü(t0), ü(s),t0,s,0)ds +

о о

1

о

Условие В7. Пусть выполнено неравенство

I »

и ('ob ¡Q0g(t0,£)dZ<0.

о о

Последнее неравенство следует понимать в смысле предела слева и справа при t~*t0-

Теорема 1.6. Пусть выполнены условия В1, В2а, ВЗ, В4, В5а, В6 иВ7. Тогда для достаточно малых е существует единственное решение задачи (9), (10), удовлетворяющее неравенству

maxi u(t,e) - Ua(t,e) | < С-Je, [0,1] 1 u 1

где Uü(t,s) = ü(t) + ГТ0м(т).

Верхнее решение и нижнее решения имеют вид, аналогичный случаю со сменой устойчивости в § 1.

Во второй главе диссертации изучаются краевые задачи для сингулярно возмущенных интефодифференциальных уравнений. В § 1 рассмотрены задачи Дирихле и Неймана для обыкновенного интегродифференциального урав-

нения с пограничными слоями [38]. Сначала рассмотрена вторая краевая задача для нелинейного обыкновенного сингулярно возмущенного интегродифферен-циального уравнения:

е2и" = Ци,и,х,е), 0<х<\, (11)

и'(0,е) = а, и\1,е) = Ь. (12)

Асимптотика задачи строится с помощью модификации метода пограничных функций [2], предложенной в § 1 гл.1 и обусловленной наличием интегрального члена в уравнении (11) в виде разложения

и(*,*) = ъ/ (ад + П>( г) + П 1и(т.)) ^

з и(х, е) + П'м(т, £•) + П2и(г,, £■),

где ип(х) - члены регулярной части асимптотики, П);2м - погранслойные члены асимптотики соответственно на левом и правом концах отрезка [0,1], где Т = х! £,т, - (1 - х)!е - погранслойные переменные. Сформулируем требования на нелинейности, входящие в интегродифференциальный оператор, которые понадобятся для построения асимптотического разложения решения и его обоснования.

Условие С1. Функция g является достаточно гладкой (При построении асимптотики с остаточным членом порядка 0(£'"+1) достаточно потребовать, чтобы она была п + 2 раза непрерывно дифференцируема). Условие С2. Пусть вырожденное уравнение

Ь(и,и,х,0) = 0, 0<х<1, имеет изолированный корень и = <р(х) на отрезке [0,1], причем

Ьи{ф,(р,х,0)>0 при хе[0,1],

(условие разрешимости и устойчивости).

Условие СЗ будет сформулировано ниже, так как оно не требуется для построения и обоснования асимптотики задачи Неймана.

Условие С4. Пусть функция g{u,v,x,s,e) будет монотонно невозрас-тающей по переменной v в некоторой области изменения этой переменной при любом фиксированном значении переменной и из той же области изменения, что и для второй переменной, при (х,х)е[0,1]х[0,1] и достаточно малых е (область, о которой идет речь в данном условии, будет уточнена в п. 3 данного параграфа при обосновании асимптотики).

Условие С5. Пусть следующее интегральное неравенство имеет положительное решение у(х):

где K(x,s) = -gv((p(x),<p(s),x,s,0)/Lu(<p,<p,x,0).

Заметим, что ядро данного интегрального неравенства неотрицательно в силу условий С2 и С4. Достаточным условием для выполнения условия С5 является, например, требование выполнения неравенства

в этом случае требуемым решением будет положительная постоянная, или требование на собственные значения следующего интегрального уравнения Фред-гольма с неотрицательным ядром K(t,s)

а именно, требуется, чтобы главное позитивное собственное значение Л 0 было меньше единицы (см. [25] стр.83, теорема 2.1).

В качестве нулевого приближения регулярной части выбирается щ = ср(х). Особенностью асимптотики задачи Неймана является то, что она не содержит пограничных функций нулевого порядка, то есть П'0,2м = 0.

Регулярный член первого порядка находится из линейного интегрального уравнения

Регулярные члены разложения более высокого порядка определяются из линейных интегральных уравнений с таким же интегральным оператором, как в левой части данного уравнения, а неоднородность в правой части выражается через уже известные члены асимптотики. В силу требования С5 и неотрицательности ядра К(х,з) это уравнение имеет единственное решение.

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

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

min[Lu(<p,(p,x,0) + Lv(<p,<p,x,0)] > О,

о

V

/

венств. Приведем определение верхнего и нижнего решения для задачи (11), (12) (см. [1]).

Определение 2.1. Функция Р(х,е) 6 С2 ((0,1)) п С" ([0,1]) называется верхним решением задачи (11), (12), если

£2р"<Ь{р.,р,х,£), *6(0,1) Р'(0,£)<а, Р'{\,£)>Ъ. Аналогично,функция Ос(х,£) называется нижним решением, если выполнены неравенства с обратным знаком.

Использована теорема о дифференциальных неравенствах [1]:

Теорема 2.1. Пусть существуют функции а(х,£) и р(х,е), являющиеся соответственно нижним и верхним решением задачи (11), (12), причем а(х,£)< Р(х,£) при хе[0,1]. Пусть, кроме того, функцш % непрерывна вместе со своими частными производными по первому и второму аргументу и удовлетворяет условию С4 на промежутке [а, Р\. Тогда существует решение задачи (1.1), (1.2) и(х,Е) такое, что а(х,£) < и(х,£) < Р(х,£).

Применяя конструктивный метод построения нижнего и верхнего решений, доказана следующая теорема:

Теорема 2.2. Пусть выполнены условия С1, С2, С4 и С5. Тогда для достаточно малых £ существует единственное решение задачи (11), (12), удовлетворяющее неравенству

тах| и(х,£) - и„(х,£) | < Сея+\

где 1]п(х,£) - частичная сумма порядка п асимптотического ряда (13). Приведем вид верхнего и нижнего решения для случая п = 0

а(х,е) - <р(х) - е(А(х) + ехр (~кх/£) + схр(-к(\ - х)/£)), Р(х,£) = <р(х) + £{А(х) + ехр (~кх/£) + ехр(-£(1 - х)!е)),

где к- некоторая положительная постоянная, А(х) - положительное на отрезке [0,1] решение интегрального уравнения

1

1и(<р,<р,х,0)у(х) + ^у{(р(х),<р(х),х,з,0)у(х)с1з = В> 0, 0 < * < 1.

о

Здесь В ~ некоторая положительная постоянная. Существование положительного решения этого уравнения гарантируется выполнением условия С5 и неотрицательностью ядра К(х,в), которая есть следствие требований С2 и С4.

Далее в § 1 главы 2 рассмотрено уравнение (11) с краевыми условиями Дирихле

и(0,е) = и°, и(1,£) = м1. (14)

Для построения и обоснования асимптотики требуется условие принадлежности граничных значений области влияния устойчивого корня:

Условие СЗ. Пусть р

| Ци,(р,i,ö)du>0 при всех р е (<p(i), и'], ¿ = 0,1.

т

Теорема 2.3. Пусть выполнены условия CI - С5. Тогда для достаточно малых е существует решение и(х,е) задачи (11), (14), для которого Öп{х,е) является равномерным на отрезке [0,1] асимптотическим приближением, причем

maxj и(х,£) - U„(x,£) [ = 0(е"+1).

Доказательство теоремы 2 основано на теореме 2.1 о дифференциальных неравенствах. Краевые условия для верхнего и нижнего решения в случае задачи Дирихле имеют следующий вид

а(0,£)<и° < ß(0,£), а(\,£)<и{ < ß{\,e).

Верхнее и нижнее решения задачи (11), (14) представляются в виде

ß(x,£) = U„(x,£) + £"(А(х) + П >(г) + T\2nßu{r,)\ а(х,£) = U„(x,£) + £"Ы(х) + П>(г) + П2ми(г,)),

где А(х) > 0 - определяется так же, как в предыдущем случае, а П\ß(T) - является модифицированной пограничной функцией.

В § 2 главы 2 изучаются решения с внутренним переходным слоем (контрастной структурой) [39] для задачи (11), (14). Контрастной структурой типа ступеньки будем называть решение, удовлетворяющее следующему предельному соотношению

г , л г л У-\х,у),0<х<у<\,

|y+W). 0<;у<х<1.

Асимптотическое разложение решения задачи (11), (14) ищется в следующем виде

и(х,с) = |>" (*С (х,у) + Рп1)и(1])+и\и(т) + П 2„и(т,)) ^

п=о (15)

В и1' (х,у,Б) +• Рыи(т],е) + п'гг(г,£) + П2и(т„£),

где ип(х,у) - члены регулярной части асимптотики, Рп и(т]) - члены асимптотики, описывающие внутренний переходный слой, т) = (х-у)/£ - переменная переходного слоя, у - точка перехода. Знак " - " означает, что х<у (г] < 0), а знак " + ", что х>у (> 0). Точка перехода х = у определяется условием

и(у,£) = срт{у), (16)

где (рт(у) - некоторая гладкая функция, такая, что

(р'~\у,у)«р(\у)<^\у,у\ УФЛ

Если решение с одним переходным слоем типа ступеньки существует, то уравнение (16) имеет единственное решение у = у.. Асимптотику точки перехода также будем искать в виде ряда по степеням малого параметра е:

у, = х0 + £Х, + £2х2 +... . (17)

Функции Пявляются погранслойными членами асимптотики соответственно на левом и правом концах отрезка [0, 1], а т = х/е, т,-(\-х)1 £ - погранслойными переменными. Следуя схеме, предложенной в §1 представим интегральный оператор I в виде, аналогичном (15)

1 = Т+Р1 + П>1 + П21.

Для определения коэффициентов разложения (17) потребуем выполнения условия гладкого сопряжения членов асимптотики в точке х = у,

еи\у,-0,£) = £и'(у, + 0,е) (18)

(множитель е введен здесь для удобства при дальнейшем разложении производных решения в асимптотические ряды).

Потребуем выполнения условия гладкости Ш, аналогичного условию С1. Для члена нулевого порядка регулярной части асимптотики получим систему двух связанных нелинейных интегральных уравнений. Следующее условие гарантирует ее разрешимость и нужные для построения и обоснования асимптотики свойства решения.

Условие 1)2. Пусть существуют две достаточно гладких функции ср^(х,у) и (р^\х,у), определенные соответственно в областях = {(х,у):0< х< у <1} и ~ {(х, у):0 < у < х <Ц, которые для любого за-

данного у е (0,1) удовлетворяют системе двух связанных интегральных уравнений

У I

}(Х>У)>Р1-'С*)+ = 0, 0<х<у,

о ^

У 1

}ё(?(,>(х,у),<рн(х,у),хлО)Ж+ = 0, у<х<1,

О у

и неравенству <р<_>(.у,.у)<<РЬ){у,у)- Пусть также для всех х,уе[0,1] выполняются неравенства Ьи((р['\х,у),<р(з,у),х,0)>0 (/ = -,+), где

¡Р^у),^)^,

В качестве разделяющей функции <р(0\у) удобно выбрать функцию, являющуюся решением уравнения

I

В силу условия 1)2 функция /*"(«, у) = ¿(и, <р(я, у), у, 0) является гладкой функцией и для любого заданного у е [0,1] достигает нулевого значения хотя бы в одной точке и е (<р( )(у,у), У>У)) • Поэтому существует хотя бы одна гладкая функция /р(ю(у), являющаяся решением данного уравнения.

Условие Б2 существенно отличает алгоритм построения асимптотики от случая, когда решения типа ступеньки ищутся в сингулярно возмущенных задачах без интегрального члена (см., например, [35]), где его аналогом является условие наличия двух изолированных корней, и составное решение слева от точки разрыва равно одном)' корню, справа - другому. При этом значение решения в любой точке отрезка слева и справа от точки разрыва х = у не зависит от у. Ясно, что условие 02 является более общим и включает в себя условие существования составного решения для обыкновенных дифференциальных уравнений без интегрального члена. Для построения пограничных функций требуется условие БЗ, аналогичное условию СЗ § 1. Алгоритм построения пограничных функций мало отличается от алгоритма § 1.

Регулярный член асимптотики нулевого порядка выбирается в следующем виде

Щ(х,у.) = <р(х,у,) .

Вводя в рассмотрение следующую функцию (интеграл энергии)

<рм(у,у)

■Ду) = / Ци,р(.ч,у),у,0)(1и,

определим точку перехода в нулевом приближении из следующего условия:

Условие 04. Пусть существует такая точка х0 е (О, 1), что

•/(*„) = О и /(х0)<0.

Условие 04 обеспечивает построение гладкой функции переходного слоя в нулевом порядке асимптотики. Для регулярного члена первого порядка щ(±)(х,у) разложения (15) получим систему двух связанных линейных интегральных уравнений, разрешимость которой обеспечивается условием квазимонотонности Г)5, аналогичным условию С5 § 1, и следующим условием

Условие 1)6. Пусть система связанных линейных интегральных неравенств

у \

м>1'\х,у) + {к'-^хлу)™'-^,}>)Ж + > О, 0 <*<>',

О у

У 1

м/+)(*,зО + 'Хх,$,у)-м{~\ь,у)с1& + р^^х^.уУЧя,^ > о, у <* < 1,

О у

где К"\х,з,у) н 0) = любого

у е (0,1) положительное решение.

Коэффициент х1 разложения (17) определяется из условия (18) С1 -сшивания асимптотики решения в точке перехода, которое в первом порядке приводит к линейному уравнению относительно х,. Коэффициент при х1 в этом линейном уравнении с точностью до постоянного множителя равен -У'Сх,,). Следовательно, в силу требования 04 неизвестное х, определяется однозначно. Следующие коэффициенты разложения (17) определяются из линейных уравнений, аналогичных уравнению для х,. Алгоритм построения формальной асимптотики определен для произвольного порядка малости по степеням малого параметра £.

При обосновании построенного асимптотического решения ведущую роль играет асимптотический метод дифференциальных неравенств, модифицированный для рассматриваемого во воротом параграфе главы 2 нового класса задач. Приведем известное классическое определение верхнего и нижнего решения для задачи (11), (14) (см. [1])- В отличие от случая, рассмотренного в §1 главы 2, нам понадобится определение верхнего и нижнего решения, которое может иметь скачок первой производной некоторой точке х отрезка [0,1].

Определение 2.2. Функции а(х,£) и ß(х, £) называются соответственно нижним и верхним решением задачи (11), (14), если

а{х,£),ß(x,е) е С([0,1]) п С1 ([0,1] \ х) о С2((0,х]) о С2([х, 1)), где х е (0,1), а(х,£)< ß{x,£), при X G [0,1], £2ß"< L(ß, ß, x, £), ега" > Ца, а, х, е), х е (0,1), а'(х + 0,£)-а'(х-0,е)>0, ß'(x + 0)-ß'(x-0)<0, .îe(0,l) ß(0,£) <u°< a( 0,£), ß(l,£) >u'> a(l,£).

Для определенных таким образом верхнего и нижнего решений справедлива теорема 2.4 о дифференциальных неравенствах [1], аналогичная теореме 2.1, опираясь на которую, на основе конструктивного метода построения нижнего и верхнего решения доказана следующая теорема:

Теорема 2.5. Пусть выполнены условия Dl -D6. Тогда для достаточно малых е существует решение задачи (11), (14), удовлетворяющее неравенству

maxi u(x,£)-U„(x,£) I <С£п. [0,1] 1 1

Здесь t/„(x,у„+,,£) - частичная сумма « -ною порядка асимптотического представления (15)

= W„+l) + P^ui-n) + П\и(т) + П]и(т,)),

/=о

л+1

где y = y„+i='^£'xr Для случая п = 0 верхнее и нижнее решения ß(x,£) и

i=0

а(х,£) определены выражениями

ß(x,£) = Т41\х,у0) + еЦ(±)(х,У/1) + sA{±\x,yp) + P^u{Vß) + +eF§XVp) + ПЬи(г) + е(П\ßu(r) + П>(т.)),

a(x,s) = ü0(±Wa) + £И^(х,уа)- £^-\х,ха) + PfMlJ +

+£Р^и(Ла) + ni«(r) + П20и(т.) + £(П\аи(т) + Il 1и(т,)). где v„ = _v, + t'S, yß - }) - eS (S> 0 - некоторая положительная постоянная), 7]ß = (x - уß)/£, tja=(x-ya)/£, A(x,y) - некоторое положительное на множестве [0,1] \ {х = у} решение системы связанных интегральных неравенств из требования D6. Функции tffuiflß), Р^и(т]а) - модифицированные функции переходного слоя, а функции П^и, П\,2и - модифицированные пограничные

функции. В п.4 § 2 приведен конкретный пример интегродифференциальной задачи, для которой выполнены все условия теоремы 2.5.

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

Условие Б4а. Пусть для любого у е (0,1) выполняется равенство J{y) = 0.

Нелинейность, удовлетворяющая условию Б4а, называется сбалансированной нелинейностью. В этом случае в нулевом приближении можно построить гладкую функцию переходного слоя для любой точки отрезка [0,1]. Тем самым, коэффициент ха разложения (17) не определяется в нулевом порядке асимптотики. Такой случай в теории сингулярных возмущений получил название критического случая. В силу требования Б4а У(у) = 0 для любого у е (0,1), то есть коэффициент при ас, в линейном уравнении для определения лг, будет нулевым, а ха определяется из условия разрешимости этого уравнения. Условие 1)1 требует, чтобы у нелинейного уравнения для определения х0 существовало решение, принадлежащее интервалу (0,1). Неизвестное же х1 определяется в следующем порядке асимптотики из условия разрешимости линейного уравнения для определения х2, множитель при котором такой же, как при х, в первом порядке асимптотики, и, следовательно, тоже равен нулю. При этом из условия разрешимости для определения х, получается уже линейное уравнение. Условие Б8 требует, чтобы коэффициент при неизвестном х1 в этом уравнении был отрицательным. Продолжая этот алгоритм далее, можно построить асимптотику точки перехода любого порядка.

В отличие от "грубого" случая верхнее и нижнее решения для п- 0 в случае сбалансированной нелинейности имеют следующий вид

Р{х,е) = П^{х,У/!) + ей/^х,^) + е2й'(*,>>) + £~Ащ(х,ур) +

+вР^Члр)++п'ы(г)+п^(0+

+ я(П^и(г) + П^н(г.)) + е2(П' ¿и(т) + П^и(г.)), а(х,е) = й^(х,уа) + ещ(1)(х,уа) + в2й^(х,уа)- £2А(±\х,ха) + ) + еР^и(т]а) + е2Р£и{Ла) + П\и(т) + П2«(г.)

+ £{П\аи{т)+П?аи(г.)) + е2(П'в«(г) + П2а«(г,)).

Проверив выполнение условий теоремы о дифференциальных неравенствах, получим для случая сбалансированной нелинейности теорему:

Теорема 2.6. Пусть выполнены условия -03,04а, В5~Э8. Тогда для достаточно малых в существует решение задачи (11), (12), удовлетворяющее неравенству

max u(x,e)-U„(x,e) \ <C£n.

[0,1] 1

В третьей главе рассмотрены решения с внутренними переходными слоями для сингулярно возмущённых эллиптических интегродифференциальных уравнений [40, 41, 44]. А именно, изучается следующая интегродифференциальная задача

s2 Аи = L(u,u,x,s),x = efi, (19)

ди

on

= 0, (20)

где L(u,v,x,s) = |j,g(M(x),v(s),x,i)<ii , Q - замкнутая односвязная область с

Q

гладкой границей 6Q, п - внутренняя нормаль к границе 8Q. Под решением типа ступеньки мы будем понимать решение, имеющее следующий предельный вид

lmu(x,£) = <

где Q, с Q, граница области Q, - простая гладкая замкнутая кривая С, целиком лежащая внутри области Q. Введем в окрестности кривой локальную систему координат (г,1), где | г | - расстояние от точки М(х) до кривой С по той нормали к кривой, которая проходит через точку М, а / - координата той точки на С, из которой эта нормаль выходит. При этом будем выбирать знак г положительным, если точка М находится в области Q,, ограниченной кривой

С, и отрицательным, если в области П\0,.. Если кривая С достаточно гладкая, а окрестность кривой достаточно мала, то нормали, выпущенные из различных точек кривой, не пересекаются в этой окрестности. Следовательно, существует взаимно однозначное соответствие между координатами (х,,х2) и (г,1) точки М. Кривая С в локальной системе координат задается уравнением r-0, а уравнение неизвестной пока кривой перехода С, в этой системе координат будем искать в виде

Г = /(/,£) = £Г1(/) + £2Г2(/) + ..., (21)

где г (Г) - достаточно гладкие периодические функции. Переменная переходного слоя 7] представляет собой растянутое с коэффициентом б расстояние от кривой С.: r/ = (r-f(l,s))ls. Аналогичным образом в окрестности границы cQ вводится локальная система координат (р,т), где р - расстояние от точки М е Q до границы SQ вдоль той внутренней нормали к границе, которая про-

ходит через точку М, а т - координата точки границы, из которой выпущена нормаль. Погранслойная переменная т = р!е является растянутым с коэффициентом 1 /е расстоянием от границы области SQ.

Асимптотическое разложение решения задачи (19), (20) будем искать с помощью модификации алгоритма метода пограничных функций, использованной в § 2, в следующем виде

и(х,е)= I г"(й„(х,С) + РяЧ^О + Пли(г,/и)) =

п"° (22) = t¡(x,C,E) + Pu(rj,l,s) + Tlu(T,m,£),

где йДх) - регулярные члены асимптотики, P¡u(ri,l) - функции внутреннего переходного слоя, П¡и(т,т) - функции пограничного слоя. Потребуем выполнение следующих условий:

Условие El. Функция g является достаточно гладкой.

Условие Е2. Пусть существуют две достаточно гладких функции (р{){х,С) и (р>Л\х,С), определенные соответственно в областях х е Q,( ) = Q\Q¡ и х е Q<+) s Q.r которые для любой простой замкнутой кривой С = 8Q¡ удовлетворяют системе двух связанных интегральных уравнений

J g((p{-\x,C)M~\s,C),x,s,Q¡)ds + J gi^-'ix^l^is.Clx^Wds = 0,

xeQf'\

J g(^^x,C),(p(-\S,C),x,s,0)ds + J g(<pM(x,C),<p(L)(s,C),x,s,0)ds = 0,

*6QW,

и неравенству (p{'\x,C)| c < <£>(+)(x,C)| . Пусть также для всех х е íi и любой кривой С s 5Q. выполняются неравенства

L.A(puXx,C),y(s,C\x,0)>0 (/ = -,+), где

coi П = К>(5'С)'

n

Представляя интегральный оператор в виде, аналогичном представлению (22),

L = L + PL + UL,

и, переходя в операторе Лапласа там, где он действует на Р-функции и

П-функции, соответственно к переменным (77,/) и (г, т), построим асимптотическое разложение решения в областях и Г2,+). Для члена нулевого порядка регулярной части асимптотики получим систему двух связанных нелинейных интегральных уравнений. В качестве ее решения выберем

"о (х,С,) = <р(х,С,).

В качестве дополнительного условия, позволяющего определить коэффициенты г(1) разложения (1.7) для кривой С,, будем использовать условие непрерывности производной по г решения и на кривой С,. В нулевом порядке асимптотики это условие гарантируется требованием ЕЗ. Введем в рассмотрение следующую функцию (интеграл энергии)

р'^.С'о)

Дх0) = | Ь(и,(р(х0,С0),х0,0)с1и,

где х0 е С0. Рассмотрим интеграл энергии для кривых С, лежащих в окрестности С0,

(И*\х,С)

J(x) = | 1(и,ф(х, С),х,0)Ли,

с)

где х е С, и перейдем в функции J(x) к локальным координатам (г,/), как это было описано в § 1. В этих координатах функцию /(*) можно представить в следующем виде

где г = г(1) - уравнение кривой С в локальных координатах. Приближение нулевого порядка кривой перехода С, определим из следующего условия.

Условие ЕЗ- Пусть на некоторой простой замкнутой кривой С0 е П

<0,

причем для всех и е ^<_)(л:0,С0),<р(+)(л:0,С0)^ выполняется неравенство

У(х0) = 0 и^(г{1))

•=о

| Цк,р(5,С0),х0,0)<йг > 0 при и е (^"'(^СД^

<р1~Чх0,Со)

Как и в одномерном случае^ для регулярного члена первого порядка и^~}(х, С) разложения (22) получим систему двух связанных линейных интегральных уравнений, разрешимость которой обеспечивается условием квазимонотонно-

сти Е4, аналогичным условию С4 § 1 гл. 2, и следующим условием:

Условие Е5. Пусть система связанных линейных интегральных неравенств

iv

С1<~> 0м

гэе К^,у) ^ = -.+). имеет для любой

кривой С ^ дС11 положительное решение.

Из условия С1 -сшивания асимптотики на кривой С, получается для определения г,(/) линейное уравнение, в котором коэффициент при неизвестном гх(1)

равен с точностью до положительного множителя — (х0,С0). В силу требова-

<1г

ния ЕЗ неизвестное г,(/) определяется однозначно. Члены переходного слоя более высоких порядков определяются из систем обыкновенных линейных дифференциальных уравнений аналогичных уравнению для />(±>и(77,/), а члены разложения (21) из уравнений, аналогичных уравнению для г,(/). Таким образом, асимптотика может быть построена для произвольного порядка точности по степеням малого параметра е.

Введем определение верхнего и нижнего решения для задачи (19), (20), имеющего скачок нормальной производной на некоторой замкнутой кривой.

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

а(х,б)<р(х,б), при хе£1,

а(х,е),Р(х,е) е С0 (О) п С'(П \ С) о С2(Пе/дП) п С2(Й,).

б2 А р<1(р,р,х,е), £ Лаг > Цех

,сс,х,б), х е О ^

8г+

-а

хеС дг

-а

хе С

>0,

дг+р

хеС дг

Р

хеС

<0,

дп

Р

л

<0<—а хе дО дп

хедО.

Используя теорему 3.1 о дифференциальных неравенствах [1], аналогичную теореме 2.1, на основе конструктивного метода построения нижнего и верхнего решений доказана следующая теорема.

Теорема 3.2. Пусть выполнены условия Е1 - Е5. Тогда для достаточно малых е существует решение задачи (19), (20), удовлетворяющее неравенству

т$х\и(х,£)-ип(х,е)\<С£п*\ (23)

Здесь ип(х,£) - частичная сумма порядка п асимптотического ряда (22), в котором кривая С, заменена кривой, уравнение которой является уравнением (21)

п+1

с правой частью, равной частичной сумме е'г:;(/). Вид верхнего и нижнего

¿=1

решений аналогичен случаю одномерной ступеньки, так как переменные / и т входят в асимптотическое разложение решения как параметры. Заметим, что теорема 3.2 и аналогичные ей теоремы их § 2 гл. 2 не обеспечивают локальной единственности решения задачи вблизи построенной асимптотики. Локальную единственность решения задачи (19), (20) можно получить, рассмотрев полученное решение как стационарное решение соответствующей интегропараболической краевой задачи

Зы

— = £2Ди-Ци,и,х,£), хеР, />0, (24)

3/

ди дп

= 0, (25)

яг

с начальными условиями, лежащими в окрестности построенного асимптотического решения. Доказательство асимптотической устойчивости по Ляпунову построенного решения задачи (19), (20), лежащего в малой окрестности и„(х,е), как стационарного решения интегропараболической задачи (24), (25) проводится также методом асимптотических дифференциальных неравенств. Для решения интегропараболической задачи (24), (25) с начальными условиями, лежащими между верхним и нижним решениями исходной задачи (19), (20), предлагаются следующие барьерные функции

Р(х^,е) = и, + (Д, - «,)ехр(-А/),

а(х, г, £) = щ + (а„ - и$) ехр(~Л(),

где и( - стационарное решение задачи (24), (25), существующее вследствие Теоремы 3.2 и удовлетворяющее неравенству (23), ап и Д, - соответственно нижнее и верхнее решения п -ного порядка точности задачи (19), (20), Л > 0 -достаточно малый, зависящий от £ параметр. Непосредственной проверкой условий для барьерных функций получим следующую теорему:

Теорема 3.3. При выполнении условий Е1-Е5 и достаточно малых £ существует локально единственное и асимптотически устойчивое стагщонар-

ное решение и „задачи (24), (25), удовлетворяющее неравенству

шах|м, -ип{х,е)\<Сс".

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

Ъи д2и Эй

М[и]з£— = е2——-еА{х,£)—--1(и,и,х,£), 0<х<1,/>0, (26) т Эх Эх

-^(0,Г,£) = 0, -^(и,£) = 0, и(х,0,£) = и°(х,£), (27)

Эх Эх

1

где £(и,у,х,с)= , £•>0- малый параметр, и°(х,£) - не-

0

которая начальная функция, имеющая вид контрастной структуры типа ступеньки. В § 2 главы 2 рассматривались стационарные решения задачи (26), (27) в случае А = 0. Целью данной главы является развитие методов, предложенных во второй главе, на новый более сложный класс задач, получение асимптотического представления для движущейся контрастной структуры типа ступеньки (фронта) и исследование вопроса существования такого решения. Результаты данной главы могут быть также рассмотрены как расширение результатов работы [26] на нелокальные краевые задачи. Некоторые относящиеся к данной задаче вопросы рассматривались в работе [21], где изучались решения типа движущихся волн для нелокальной задачи с бистабильной нелинейностью и линейным интегральным членом. Задача (26), (27) описывает многие важные практические приложения. В частности;является одномерным аналогом задачи вида (1), (2), возникающей в теории фазовых переходов [18].

Для построения и обоснования асимптотики требуется выполнение ряда условий. Так условия И и ¥2 аналогичны условиям 01 и 02. Задача состоит в том, чтобы показать, что при некоторых дополнительных предположениях, решение задачи (26), (27) имеет вид движущегося внутреннего слоя (фронта), имеющего следующий предельный вид:

Нш и(х,1,£) = <р(х,х0{0), хе[0,1]\{х0(0).

где х0(?) - нулевое приближение точки перехода >;(?,£), которая ищется в виде асимптотического разложения

у,Ц,£) = х0(,<) + Ex.it) + £2хг(1) +...,

и определяется условием пересечения решения и с некоторой гладкой функцией (р{0)(х), такой что (р{'](х,х) < <рт(х) < (р'+'(х,х):

иЫ(,£),ие) = <рт[у.(!,е)). (28)

В качестве такой разделяющей функции удобно выбрать функцию^ являющуюся решением уравнения

В § 2 гл. 4 строится асимптотические разложение решения с внутренним переходным слоем задачи (26), (27) в предположении, что начальная функция и°(х,е) имеет вид фронта, то есть такой, что для малых е она близка к некоторому решению ср(х,х00) из семейства разрывных решений вырожденной задачи для 0 < х < х00 - А, х00 + А < х < 1, где хю - некоторая точка из интервала (0,1), в которой в начальный момент времени локализован фронт, а А>0 - некоторое малое не зависящее от £ число. Тем самым предполагается, что

Для построения формального асимптотического разложения начально-краевая задача (26), (27) с условием (28) представляется в виде краевой задачи для системы двух связанных интегродифференциальных уравнений

Э и 2д2и ,, чЭ и т/ . . . . .

е-—- = е —~т-еА(х,е)— -¿(и,и,х,е), 0 < х <у.0,е), г>0, (29) Э/ дх ах

|^(0,/,£0 = 0, и(уЛ,е),ие) = 9т(у.(?,е)) =0, ах (30)

и(х, 0, е) = м° (х, е), 0 < х < (/, е),

ди 2д2и . Эн . , . „ .„,.

£- = £2^-£А{х,£)-—Ци,и,х,е), у.ц,е)<х< 1, ¿>0, (Л) ш ах ах

и(у.((,ф,£) = <р(0)(уМ,£))^ 0, ~(1,(,£) = 0,

дх (32)

и(х,0,£) = и\х,£), >',(?,£) <х< 1.

Асимптотика решений левой и правой задачи строится при помощи алгоритма, являющегося модификацией алгоритма, предложенного в § 2 гл.2 диссертации. Будем использовать сведущее асимптотическое представление решения для задачи (29) - (32):

- 1У("Г + Р?и(т,,1) + П\и(т,г) + п>(г„0) - (33)

л=0

= и(х, , е) + Рыи(т},1, е) + П'и(г,/, е) + П2м(г,, г, г),

где ип (х, >'.(/, с),/) - члены регулярной части асимптотики, Рп и(г],{) - члены асимптотики, описывающие внутренний переходный слой и существенные в малой окрестности точки перехода, Т] = (х-у, &,£)) 1е - переменная переходного слоя. Знак " - " означает, что х (/,£■) (77 <0), а знак " + ", что х>у, (/,£■) (?7 > 0). Функции П]'2и являются погранслойными членами асимптотики^ существенными соответственно на левом и правом концах отрезка [0, 1], а г = х!е, г, = (1 - х)!е - погранслойными переменными.

Члены регулярной части асимптотики определяются способом, аналогичным использованному в § 2 гл.2. Для построения членов первого и следующих порядков регулярной части асимптотики предполагается выполнение условий Р4 и Е5, аналогичных условиям Б5 и Б6 соответственно.

Так как мы не рассматриваем процесс выхода фронта на границу отрезка [0,1], то пограничные функции не будут оказывать никакого влияния на функции переходного слоя. Поэтому их построение ничем не будет отличаться от стационарного случая, рассмотренного в §§ 1, 2 второй главы. А аргумент I будет входить в них только через коэффициенты разложения у,((,е). Краевые условия Неймана приведут к тому, что в нулевом приближении пограничный слой будет отсутствовать.

Переходя в дифференциальном операторе

_ ди 2 д2 ,, > Э

И = -£-+ в -г- - еА{х, £)-

9/ дх ох

к переменным (т/,1), получим

Используя данный оператор, и, вводя в рассмотрение следующую непрерывную функцию, описывающую переходный слой в нулевом приближении,

<Р(~)(У,У) + Ро"(т},У), П<0, <Рщ{у), 77 = 0, <Р^ХУ>У) + РАЛ>У)> до-

получим для и задачу

д2й э-г 1

. си + vl

оТ~= ¡ё(й(п,У),Ф>У)>У^,^, 1] е (-со,+со),

<У "дч - . V , . (34)

где у0 . Уравнение в задаче (34) является обыкновенным

V ей ^

дифференциальным уравнением второго порядка. Из условия Р2 следует, что (р{±\у,у) - седловые точки. Такая задача хорошо изучена (см., например, [26, 27]). Основные результаты исследования задачи (34) можно сформулировать в виде следующей леммы:

Лемма 4.1. Для любого у е (0,1) существует единственное значение у(| такое, что задача (34) имеет единственное гладкое монотонное решение й, удовлетворяющее оценке

| й07>У)~<Р™(У,У)\£ Сехр(-ф|),

где С и к - некоторые положительные постоянные. А функция У0(у) удовлетворяет равенству

9{*\у,у)

| Ци,<р(з,у),у,0)с1и

Из леммы 4.1 следуют экспоненциальная оценка для Р0и(^,у):

р0и(г?,у)^Сехр(-ф\),

и то, что нулевое приближение хй({) точки перехода у = у,(*>£) может быть определено как решение начальной задачи

ск т

где х00 определяет начальное положение точки перехода. Как показано в § 2 гл. 2, корни функции (у0(х0)-М(х0,0)) определяют положение стационарных точек перехода внутреннего слоя. В данной главе изучается только движение фронта, поэтому предполагается выполнение следующего требования:

Условие ГЗ. Пусть (у0(х0) + Л(х0,0))>0 для любого у е [0,1].

Для определения я, (7) используем первый порядок условия С -сшивания асимптотики:

ЭР{~"'и г)Рми Эй><+)

^(0,0 +^-(х0(О,*ь(О) =—^(0-0 (35)

ОТ/ дх 67] ох

дР(±)и

Из явного интегрального представления для —^—(0,?) и условия (35) полуду

чим для определения х, (/) линейное обыкновенное дифференциальное уравнение

Иг

где функции В(х0) и /^(х0) выражаются через уже построенные члены асимптотики. Если входные данные задачи достаточно гладкие, то асимптотика решения может быть построена до любого порядка п, и формальная асимптотика

удовлетворяет задаче с невязкой порядка е"+1.

Так же, как в предыдущих главах, для обоснования полученного выше формального асимптотического разложения применим модификацию метода асимптотических дифференциальных неравенств [8, 9] на новый класс задач -интегропараболические уравнения с движущимся фронтом.

Определение 4.1. Функция /?(х,?,с) называется верхним решением задачи (26), (27), если

Ах, е) е С([0,1] к [О,Г]) о С2']((0,х(0] х [0,Г]) о С2'([х(0,1) * [0, где х(?) е (0,1) при I е [0,7"] и является гладкой функцией,

М\Р\<0 для всех (х,Ое{(0,х(0]х[0,Г]}п{[х(0,0)х[0,Г]}, (36)

^(х + 0,(,е)-^-(х-0,(,е)<0, (37)

ах йх

^(0,г,£)<0 и ^(М,£)>0, (38)

дх дх

/?(х,0 ,£)>«°(х,£), (39)

где М - интегродифференциалъный оператор уравнения (26). Аналогично, нижним решением называется функция ог(х,/,£), принадлежащая тому же классу гладкости и удовлетворяющая всем неравенствам с противоположным знаком.

Доказательство существования и единственности решения задачи (26), (27) основывается на следующей теореме о дифференциальных неравенствах (см. [1], гл. 2. разд. 2.7, см. также [28-30]).

Теорема 4.1. Если существуют функции /?(х,гг) и а(х,1,е), такие, что:

(a) Р(х,1,с) и а(х,1,£)естъ верхнее и нижнее решение задачи (26), (27) соответственно;

(b) Р(х,[,с) > а(х,I,в) для любых е (1.0,1] х [0,Г|};

(c) еС'СОД], gu(u,v,x,s) и ёХи^,х,$)еС{[а(х,1,-),Р(х,1,-)} х[а(5Л-),А*,г,-)]х[0,1]2);

(ф £„(...)^0 длявсех g(u,v,x,s)e[a(x,t,■),P(x,t,■)]x[a(s,t,•), /?(5,г,-)]х[0,1]2; то задача (26), (27) имеет единственное классическое решение и{х,1,с) такое, что Р(х,^в)>и(х,(,в)>а(х,1,в) для (лг,*)е{[0,1]х[0,Г]}.

Верхнее решение задачи (26), (27) строится в следующем виде

Д,(х,/,гг) = 110(х,ур(1,е),Е) + £й^(х,у^, £)) +

+£А{±\х,у^,£)) + вР^и(Т]р,х^),х^)) + с(е " + е'"'),

Здесь ио(х,у,е) - частичная сумма нулевого порядка разложения (33) с

уДЦ,е) = х0Ц) + Ех,/)( 0,

Растянутая переменная переходного слоя в верхнем решении определяется как Т]р = (х-ур^^))/£. Коэффициент х1/;(?) будет определен ниже. Функции А(±)(х,у) являются некоторым положительным решением системы связанных интегральных уравнений, существующим в силу условия Г5. Функции Р/^и -

модифицированные функции переходного слоя. Вместо условия С' -сшивания асимптотики в точке т]р = 0 для верхнего решения требуется выполнение условия скачка производной (37) из определения 4.1

Э Р$и д<р(+)

—г (0,дг0,х,„) Н — (Х^Хд)" д}]„ н дх

(дР^и Э^-»

<-£Т,

У

где а - некоторая положительная постоянная. Из последнего условия получим уравнение для определения х1р{1)

^ - В(х0(0)х1/г = Р\р (х0(О) - -^-|^(0,х0),

й?/ СО (X„) 07]

где функция /^(х^)) получается некоторым преобразованием из функции

^ (х0(0) ■ Выбирая начальное условие для х|/((?) в виде

где 5 - некоторая положительная постоянная, получим, что при достаточно больших значениях 5

х[р(1)<0 для всех /е[0,Г].

Итак, верхнее решение /?0(х,г,£) определено. Условия (36), (38) определения 4.1 проверяются аналогично тому, как это было сделано в п. 3 § 2 гл. 2. Условию (39) всегда можно удовлетворить с помощью выбора подходящей начальной функции и0(х,е).

Нижнее решение аа(х,1,£) определяется аналогично. Условие скачка производной для а0(х,г,£) будет следующим

дР^и д<р{+)

„ (и,,» ) + (х„,х0)

дт}а р дх

f'AP^U, -лЛ-) Л

>(У,

Эр д<р{

dij ,X<),X[ß' М)

а х1ц(г) определяется из задачи

Выбирая начальное условие для xIa(t) в виде

х|а(0) = £,

получим, что при достаточно больших значениях ö

xla(t) > 0 для всех t е [0,Г].

Выполнение условий теоремы 4.1 для ß0(x,t,c) и a0(x,t,e) доказывает, что задача (26), (27) имеет единственное классическое решение, имеющее предельный вид

lim u(x,t,s) = <p(x,x0(t)), xe[O,l]\{x0(f)},

для любой начальной функции такой, что

а0 (х, t,a)<u(x, е)< ß0 (х, t,s). Проверка случая и> 0 для ßn(x,t,£) и an(x,t,s), где

ß„(x,t,e) = U„(x,yß{t,e),e) + £л+,й„+,(х,уДг,£)) + £"'Аи>(х,у^, с)) + + £"+lP^ßu(nß,x0(t), ...,xn{t),x,n^ß{t)) + ir»V" + ),

ап(х,1,в) = и„(х,уа(1,Е),£) + е"\+1(х,уа(м)) ~ £^А(±\х,уа(1,с)) +

+ 0, *<„.)«('>)- «"'(е-" + ^),

проводится аналогично.

Теорема 4.2. Пусть выполнены условия Р1 - Р5. Тогда для достаточно малых е задача (26), (27) с начальной функцией и°(х,в), удовлетворяющей неравенству

а0 0,0.Е) ¿и(х, £)< Д, (х, 0, с), имеет при / е [0,7"] единственное классическое решение и(х,1,£) такое, что Нш и(х,?,г) = г/;(л-,х0(/)), хе[0,1]\{х0(/)}.

Если же начальная функция и°(х,£) удовлетворяет неравенству

аЙ{х,Ъ,Б)йи\х,е)йРп{хМ,

то задача (26), (27) имеет при ? е [0,Г] единственное классическое решение и(х,1,е) такое, что справедлива 01(енка

шах| и(х,Г,е) -II„(х,>„,,(/,£),в) | <Св",

п+1

где у^(1,£) = ^Г£'х,(Г).

¡=о

Заключение

Основные результаты работы, полученные лично автором:

1. Построены асимптотические приближения решений для следующих новых классов нелинейных сингулярно возмущенных задач:

- Начальные задачи с нелинейными интегральными операторами типа Воль-терра и Фредгольма, в том числе в случае смены устойчивости корня вырожденного уравнения.

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

- Краевые задачи для эллиптических интегродифференциальных уравнений с пограничными и внутренними слоями (двумерными контрастными структурами типа ступеньки).

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

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

В заключение автор выражает глубокую признательность профессору Николаю Николаевичу Нефёдову за постоянное внимание к работе и полезные обсуждения.

СПИСОК ЛИТЕРАТУРЫ

1. Рао C.V. Nonlinear Parabolic and Elliptic Equations. New York: Plenum, 1992.

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

3. Васильева А.Б. К вопросу о близких к разрывным решениях в системе с малым параметром при производных условно устойчивого типа // Дифференциальные уравнения. 1972. Т. 8. N 9. С. 1560-1568.

4. Бутузов В.Ф., Васильева А.Б. Об асимптотике решения типа контрастной структуры // Математические заметки. 1987. Т. 42. N 6. С. 831-841.

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

6. Файф П., Гринли В. Внутренние переходные слои для эллиптических краевых задач с малым параметром // Успехи мат. наук. 1974. Т. 29. N 4. С. 103— 131.

7. Бутузов В.Ф., Васильева А.Б., Нефедов H.H. Контрастные структуры в сингулярно возмущенных задачах // Фундаментальная и прикладная математи-ка.1998. T.4.N3.C. 799-851.

8. Нефедов H.H. Метод дифференциальных неравенств для некоторых классов нелинейных сингулярно возмущенных задач с внутренними слоями // Дифференциальные уравнения. 1995. Т. 31. N 7. С. 1132-1139.

9. Нефедов H.H. Метод дифференциальных неравенств для некоторых сингулярно возмущенных задач в частных производных // Дифференциальные уравнения. 1995. Т. 31. N 4. С. 719-722.

Ю.Васильева А.Б. Об устойчивости контрастных структур // Математическое моделирование. 1991. Т. 3. N 4. С. 114-123.

11.Бутузов В.Ф. О неустойчивости контрастных структур типа всплеска // Математические модели и методы в социальных науках. (Труды вторых математических чтений МГСУ 26 января - 2 февраля 1994). М.: МГСУ. 1994. С. 14-18.

12.Angenent S., Mallet-Paret J., Peletier L. Stable transition layers in a semilinear boundary value problems //J. Diff. Equations. 1987. V. 67. N 2. P. 212-242.

13.Hale J. К., Sakamoto К. Existence and stability of transition layers// Japan J. of Appl. Math. 1988. V. 5. N 3. P. 367-405.

14. Бутузов В.Ф., Неделько И.В. О глобальной области влияния устойчивых решений с внутренними слоями в двумерном случае // Известия РАН (серия математическая). 2002. Т. 66. N 1. С. 3-42.

15.Бутузов В.Ф., Неделько И.В. О формировании контрастной структуры типа ступеньки в параболической системе с разными степенями малого параметра//Доклады РАН. 2003. Т. 390. N 1. С. 15-18.

16. Raquepas J., Dockery J. Dynamics of a reaction-diffusion equation with nonlocal inhibition // Physica D. 1999. V. 134. P. 94-110.

17.Novick-Cohen A. The Cahn-Hilliard equation: Mathematical and Modelling Perspectives // Advances in Math. Sci. and Appl. 1998. V. 8, 965-985.

18. Rubinstein J., Sternberg P. Nonlocal reaction-diffusion equations and nucleation // IMA J. Appl. Math. 1992. V. 48. P. 249-264.

19. Okada K. Intermediate dynamics of internal layers for a nonlocal reaction-diffusion equation // Hiroshima Math. J. 2005. V. 35. P. 263-308.

20.Bates, P., Zhao, G. Existence, uniqueness and stability of the stationary solution to a nonlocal evolution equation arising in population dispersal // J. Math. Anal. Appl. 2007. V. 332. N 1, P. 428-440.

21. Bates, P., Chen, F. Spectral analysis of traveling waves for nonlocal evolution equations // SIAM J. Math. Anal. 2006. V. 38. N. 1. P. 116-126.

22. Kot M., Lewis M., Driessche P. Dispersal data and the spread of invading organisms // Ecology. 1996. V. 77. N 7. P. 2027-2042.

23. Medlock J., Kot M. Spreading disease: Integro-differential equations old and new // Mathematical Biosciences. 2003. V. 184. N 2. P. 201-222.

24. Butuzov V.F., Nefedov N.N., Schneider K.R. Singularly perturbed boundary value problems in case of exchange of stabilities // J. Math. Analys. and Appl. 1999. V. 229. P. 543-562.

25. Забрейко П.П., Кошелев А.И., Красносельский M.A. и др. Интегральные уравнения. М.: Наука, 1968.

26. Fife P., Hsiao L. Generation and Propagation of Internal Layers //' Nonlinear Anal. 1988. V. 12. N l.P. 19-41.

27. Perko L. Differential Equations and Dynamical Systems. New York: Springer, 2001.

28.Amann H. Periodic Solutions of Semilinear Parabolic Equations, Nonlinear Analysis: a Collection of Papers in Honor of Erich Rothe. New York: Academic, 1978, pp. 1-29.

29. Sattinger D. Monotone Methods in Elliptic and Parabolic Boundary Value Problems // Indiana Univ. Math. J. 1972. V. 21. N 11. P. 979-1001.

30. Fife P., Tang M. Comparision Principles for .Reaction-Diffusion systems: Irregular Comparision Functions and Applications to Question of Stability and Speed Propagation of Disturbances // J. Diff. Equations. 1981. V. 40, P. 168-185.

31. Михайлов А.П. Моделирование системы "власть-общество". М.: ФИЗМАТЛИТ, 2006.

СПИСОК ОСНОВНЫХ ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

32.Никитин А.Г. Неустойчивость контрастных пространственных структур типа "всплеска" в системе реакции-диффузии // Журнал вычислительной математики и математической физики. 1992. Т. 31, N 3. С. 443-452.

33.Васильева А.Б., Никитин А.Г., Петров А.П. Асимптотический метод исследования контрастных структур и его приложения к теории гидромагнитного динамо // Математическое моделирование, т. 7, 1995, № 2, с. 61 - 71.

34.Vasil'eva A., Nikitin A., Petrov A. Stability of contrasting solutions of nonlinear hydromagnetic dynamo equations and magnetic fields reversals in galaxies // Geophys. Astrophys. Fluid Dynamics, v. 78,1995, pp. 261 - 279.

35.Васильева А.Б., Никитин А.Г. К вопросу об устойчивости периодических контрастных структур в пространственно двумерном случае // Дифференциальные уравнения, 1996. Т. 32. № 10. С. 1355-1361.

36.Никитин А.Г. О главной собственной функции одной сингулярно возмущенной задачи Штурма-Лиувилля // Журнал вычислительной математики и математической физики, 1999, т. 39, № 4, с. 558-591.

37.Никитин А.Г., Петров А.П. О предельном переходе по малому параметру для собственных значений сингулярно возмущенной задачи Штурма-Лиувилля // Дифференциальные уравнения, 1999, т. 35, №6, с. 843-845

38.Нефедов Н.Н., Никитин А.Г. Асимптотический метод дифференциальных неравенств для сингулярно возмущенных интегро-дифференциальных уравнений //Дифференциальныеуравнения, 2000,т.36, № 10, с. 1398-1404.

39.Нефедов Н.Н., Никитин А.Г. Развитие асимптотического метода дифференциальных неравенств для решений типа ступеньки в сингулярно возмущенных интегро-дифференциальных уравнениях // Журнал вычислительной математики и математической физики, 2001, т. 41, № 7, с. 1057 -1066.

40.Нефедов Н.Н., Никитин А.Г., Асимптотическая устойчивость контрастных структур типа ступеньки в сингулярно возмущённых интегро-дифференциальных уравнениях в двумерном случае //Математическое, моделирование, 2001,т. 13, № 12, с. 65-74.

41.Nikitin A.G. Contrast structures in the integro-differential equations // Progress of nonlinear science. Proceeding of international conference dedicated to the 100th Anniversary of A.A Andronov Vol. I. Mathematical Problems of Nonlinear Dynamics,University of Nizhny Novgorod, 2002,323-326

42.Нефедов H.H., Никитин А.Г. Сингулярно возмущенные интегро-дифференциальные уравнения в случае сбалансированной нелинейности // Труды второй международной конференции "Нелинейные дифференциаль-

ные уравнения в частных производных» (Алушта, 2005), В сб.: Нелинейные граничные задачи, Институт прикладной математики HAH Украины,2006, с. 186-192.

43.Нефедов H.H., Никитин А.Г., Уразгильдипа Т.А. Задача Коши для интегро-дифференциального уравнения Вольтерра. Журнал вычислительной математики и математической физики, 2006, т. 46, №5, с. 805-812.

44.Нефедов H.H., Никитин А.Г. Метод дифференциальных неравенств для контрастных структур типа ступеньки в сингулярно возмущенных интегро-дифференциальных уравнениях в пространственно двумерном случае. Дифференциальные уравнения, 2006, т. 42, №5, с. 690-700.

45.Нефедов H.H., Никитин А.Г. Задача Коши для интегро-дифференциального уравнения Фредгольма. Журнал вычислительной математики и математической физики, 2007, т. 47, №4, с. 655-664.

46.Nefedov N.N., Nikitin A.G., Recke L-. Moving Internal Layers in the Singular Perturbed Integro-Parabolic Reaction-Diffusion-Advection Equations. Preprint Nr. 2007-22. Humboldt University of Berlin, Institute of Mathematic, pp. 1-17.

Подписано к печати i9.0i.D9 Тираж j2.fi_Заказ 4

Отпечатано в отделе оперативной печати физического факультета МГУ

Введение

Глава 1 Начальные задачи для сингулярно возмущенных интегродифференциальных уравнений.

§ 1 Задача Коши для сингулярно возмущенного интегродифференциального уравнения типа Вольтерра

§ 2 Задача Коши для сингулярно возмущенного интегродифференциального уравнения Фредгольма.

Глава 2 Краевые задачи для обыкновенных сингулярно возмущенных интегро диф ф еренциальных уравнений типа реакция-диффузия.

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

§ 2 Внутренние слои в нелинейной краевой задаче для обыкновенного сингулярно возмущенного интегро-дифференциального уравнения типа реакциядиффузия

Глава 3 Краевые задачи для сингулярно возмущенных интег-родифференциальных уравнений эллиптического типа

§ 1 Постановка задачи.

§ 2 Формальная асимптотика решения.

§ 3 Существование и асимптотическая устойчивость решения типа ступеньки.

Глава 4 Движущиеся фронты в интегропараболическом уравнении реакция-диффузия-адвекция.

§ 1 Постановка задачи.

§ 2 Формальное асимптотическое разложение решения с внутренним переходным слоем.

§ 3 Обоснование формального асимптотического разложения решения с движущимся внутренним переходным слоем (фронтом)

§ 4 Пример.

Актуальность темы.

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

В теории сингулярных возмущений контрастные структуры ранее исследовались в нелинейных эллиптических краевых задачах с малыми параметрами при старших производных, рассматриваемых в ограниченных областях. Впервые существование контрастных структур в сингулярно возмущенных задачах для обыкновенных дифференциальных уравнений было доказано в работах А.Б. Васильевой и В.Ф. Бутузова [2-5]. Результат по существованию двумерных контрастных структур типа ступеньки принадлежит П. Файфу (Р. Fife) и У. Гринли (W. Greenlee) [6]. Асимптотические разложения решений типа контрастной структуры по малому параметру можно построить на основе метода пограничных функций [2, 34-36]. Для одномерных задач это сделано в [2-5, 7, 34-35] и ряде других работ, для некоторых двумерных задач -в [8, 9, 36]. Обширная библиография по этой проблематике содержится в [7].

Важным вопросом, как с теоретической, так и с прикладной точки зрения, является вопрос об устойчивости контрастных структур как стационарных решений соответствующих параболических задач (в смысле Ляпунова). Для одномерных задач этот вопрос был решен в работе. А.Б. Васильевой [10], В.Ф. Бутузова [10], работах [33-35, 37, 38], S. Angenent., J. Mallet-Paret и L. Peletier [12], J. Hale и К. Sakamoto [13] и др. Устойчивость периодической контрастной структуры типа ступеньки в пространственно двумерном случае была впервые получена в работе [36] путем исследования спектра сингулярно возмущенной двумерной задачи на? собственные значения. Вопросы устойчивости и локальной единственности решений сингулярно возмущенных нелинейных эллиптических задач, а так же важная проблема формирования контрастных структур в сингулярно возмущенных параболических задачах, были решены В.Ф. Кутузовым и И.В. Неделько с помощью предложенного ими метода параметрических барьеров [14, 15]:

Наиболее эффективным методом доказательства существования контрастных структур и оценки.остаточных членов асимптотических разложений является асимптотический метод дифференциальных неравенств Н.Н. Нефедова [7, 8]. Суть его состоит в том, что верхнее и нижнее решения конструируются путем модификации формальной асимптотики. Рассматривая эллиптическую задачу как стационарную задачу для соответствующего параболического уравнения, этим методом можно также доказать устойчивость по Ляпунову и локальную единственность решения исходной задачи.

В настоящее время большой интерес вызывают более сложные модели, которые включают эффекты обратной связи или нелокального взаимодействия. Различные направления теории нелокальных нелинейных моделей интенсивно разрабатываются как у нас в стране, так и за рубежом. Как правило, эти модели представлены сингулярно возмущенными интегродифференци-альными уравнениями, описывающими важные для. приложений процессы, в которых необходимо принять во внимание последствия или задержку, во многих областях естествознания, в частности, в задачах динамики реакторов, моделях генетики популяций, химической кинетике [16], теории фазовых переходов [17-19], социологии [20] и многих других областях [1, 21-24]. Так модели, обладающие наследственными свойствами, описываются практически только интегродифференциальными уравнениями [23, 24]. В частности, например, в теории фазовых переходов при рассмотрении теоретической модели процесса разделения фаз в двойной полимерной смеси возникает следующая задача для уравнения Кана-Хилиарда (Cahn-Hilliard) [17-19] ut = A [/(w)-£2 Aw], />0,

1) f^ = 0> = u(x,y,0) = g(x,y), где u(x,y,t) - концентрация одной из компонент смеси, Q — ограниченная область с гладкой границей, п - единичный вектор внешней нормали к границе области ЭП, /(и) = W'(u), где W(u) - двойная потенциальная яма, е -диапазон межмолекулярных сил. Простыми вычислениями [18] эта задача сводится к задаче для нелокального уравнения реакция-диффузия ut = £2Ди - f{u) + \f(u)<m, {х,у)е Q, />0, Q

2) ди = 0, О, у) е 8Q, и(х, у, 0) = g(x, у). дп

Изучение новых моделей, подобных приведенной выше, требует развития соответствующих методов математической физики, адекватных сложности задач.

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

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

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

Научная новизна; Все основные результаты работы являются новыми.

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

Апробация работы. Основные результаты диссертации докладывались и обсуждались на семинаре МГУ по малому параметру (руководители: профессора А.Б. Васильева, В.Ф. Бутузов, Н.Н. Нефедов), на научных семинарах факультета ВМ и К (руководитель: профессор И.А. Шишмарев), кафедры математики физического факультета МГУ, НИВЦ МГУ (руководители: профессора А.Б. Бакушинский, А.В. Тихонравов и А.Г. Ягола), на международной конференции "Теория и приложения методов малого параметра", посвященной 90-летию со дня рождения академика А.Н. Тихонова (Обнинск, 1996) [48], на международной конференции "Асимптотики решений дифференциальных уравнений", посвященной 70-летию академика A.M. Ильина (Уфа, 2002) [49], на международных конференциях "Nonlinear partial differential equations" (Алушта, 2003, 2005) [50, 51], на международной конференции, посвященной 100-летию А.А. Андронова (Нижний Новгород, 2001) [42], на международных конференциях "Математические идеи П.Л. Чебышева и их приложение к современным проблемам естествознания" (Обнинск, 2002, 2004, 2008) [52-54], на седьмой Крымской международной математической школе "Метод функций Ляпунова и его приложения" (Алушта 2004) [55], на VI международном конгрессе по математическому моделированию (Нижний Новгород, 2004) [57], на международной конференции "Tikhonov and contemporary mathematics" (Москва, 2006) [57], на международной конференции "Дифференциальные уравнения и смежные вопросы", посвященной памяти И.Г. Петровского (Москва, 2007) [58], на-Тихоновских чтениях (Москва; 2002, 2006, 2008), на Ломоносовских чтениях (Москва, 2006, 2008) [59, 60], на вторых - шестых (1994-1998) [61-65], восьмых (2000) [66], десятых (2002) [67], одиннадцатых (2003) [68], пятнадцатых - семнадцатых (2005-2008) [69, 70] математических чтениях РГСУ (МГСУ), на международной конференции Science Links Moscow-Berlin-Paris Workshop 2008, посвященной 50-тию сотрудничества МГУ им. М.В. Ломоносова с Гумбольдтским университетом (Берлин, Германия).

Публикации Основные результаты, полученные автором и изложенные в диссертации, опубликованы в работах [33-47]. По материалам диссертации опубликованы 15 научных работ и сделано 29 докладов на научных конференциях. Результаты, содержащиеся в работах, выполненных в соавторстве, и включенные в диссертацию, получены автором лично и включены в диссертацию с согласия и одобрения соавторов этих работ.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, разбитых на параграфы (некоторые параграфы, в свою очередь, разбиты на пункты), заключения и списка литературы, содержащего 70 наименований. Нумерация формул своя в каждом параграфе (пункте). В работе для формул принята двойная нумерация: первое число - номер параграфа (пункта параграфа), второе - порядковый номер формулы в параграфе (пункте параграфа). Объем диссертации составляет 198 страниц, включая 8 страниц цитированной литературы.

Основные результаты работы, полученные лично автором:

1. Построены асимптотические приближения решений для следующих новых классов нелинейных сингулярно возмущенных задач:

- Начальные задачи с нелинейными интегральными операторами типа Вольтерра и Фредгольма, в том числе в случае смены устойчивости корня вырожденного уравнения.

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

- Краевые задачи для эллиптических интегродифференциальных уравнений с пограничными и внутренними слоями (двумерными контрастными структурами типа ступеньки).

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

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

В заключение автор выражает глубокую признательность профессору Николаю Николаевичу Нефёдову за постоянное внимание к работе и полезные обсуждения.

191

ЗАКЛЮЧЕНИЕ

1. Рао С. V. Nonlinear Parabolic and Elliptic Equations. New York: Plenum, 1992.

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

3. Васильева А.Б. К вопросу о близких к разрывным решениях в системе с малым параметром при производных условно устойчивого типа // Дифференциальные уравнения. 1972. Т. 8. N 9. С. 1560-1568.

4. Бутузов В.Ф., Васильева А.Б. Об асимптотике решения типа контрастной структуры // Математические заметки. 1987. Т. 42. N 6. С. 831-841.

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

6. Файф П., Гринли В. Внутренние переходные слои для эллиптических краевых задач с малым параметром // Успехи мат. наук. 1974. Т. 29. N 4. С. 103-131.

7. Бутузов В.Ф., Васильева А.Б., Нефедов Н.Н. Контрастные структуры в сингулярно возмущенных задачах // Фундаментальная и прикладная математика. 1998. Т.4. N 3. С. 799-851.

8. Нефедов Н.Н. Метод дифференциальных неравенств для некоторых классов нелинейных сингулярно возмущенных задач с внутренними слоями // Дифференциальные уравнения. 1995. Т. 31. N 7. С. 1132-1139.

9. Нефедов Н.Н. Метод дифференциальных неравенств для некоторых сингулярно возмущенных задач в частных производных // Дифференциальные уравнения. 1995. Т. 31. N 4. С. 719-722.

10. Васильева А.Б. Об устойчивости контрастных структур // Математическое моделирование. 1991. Т. 3. N4. С. 114—123.

11. Бутузов В.Ф. О неустойчивости контрастных структур типа всплеска // Математические модели и методы в социальных науках. Труды вторых математических чтений МГСУ. М.: МГСУ, 1994. С. 14-18.

12. Angenent S., Mallet-Paret J., Peletier L. Stable transition layers in a semilinear boundary value problems // J. Diff. Equations. 1987. V. 67. N 2. P. 212-242.

13. Hale J. K., Sakamoto K. Existence and stability of transition layers// Japan J. of Appl. Math. 1988. V. 5. N 3. P. 367-405.

14. Бутузов В.Ф., Неделько И.В. О глобальной области влияния устойчивых решений с внутренними слоями в двумерном случае // Известия РАН (серия математическая). 2002. Т. 66. N 1. С. 3-42.

15. Бутузов В.Ф., Неделько И.В. О формировании контрастной структуры типа ступеньки в параболической системе с разными степенями малого параметра // Доклады РАН. 2003. Т. 390. N 1. С. 15-18.

16. Raquepas J., Dockery J. Dynamics of a reaction-diffusion equation with nonlocal inhibition // Physica D. 1999. V. 134. P. 94-110.

17. Novick-Cohen A. The Cahn-Hilliard equation: Mathematical and Modelling Perspectives // Advances in Math. Sci. and Appl. 1998. V. 8, 965-985.

18. Rubinstein J., Sternberg P. Nonlocal reaction-diffusion equations and nuclea-tion // IMA J. Appl. Math. 1992. V. 48. P. 249-264.

19. Okada K. Intermediate dynamics of internal layers for a nonlocal reaction-diffusion equation // Hiroshima Math. J. 2005. V. 35. P. 263-308.

20. Михайлов А.П. Моделирование системы "власть-общество". М.: ФИЗМАТЛИТ, 2006.

21. Bates, P., Zhao, G. Existence, uniqueness and stability of the stationary solution to a nonlocal evolution equation arising in population dispersal // J. Math. Anal. Appl. 2007. Y. 332. N 1, P. 428-440.

22. Bates, P., Chen, F. Spectral analysis of traveling waves for nonlocal evolution equations // SIAM J. Math. Anal. 2006. V. 38. N. 1. P. 116-126.

23. Kot M., Lewis M., Driessche P. Dispersal data and the spread of invading organisms // Ecology. 1996. V. 77. N 7. P. 2027-2042.

24. Medlock J., Kot M. Spreading disease: Integro-differential equations old and new // Mathematical Biosciences. 2003. V. 184. N 2. P. 201-222.

25. Butuzov V.F., Nefedov N.N., Schneider K.R. Singularly perturbed boundary value problems in case of exchange of stabilities // J. Math. Analys. and Appl. 1999. V. 229. P. 543-562.

26. Забрейко П.П., Кошелев А.И., Красносельский M.A. и др. Интегральные уравнения. М.: Наука, 1968.

27. Fife P., Hsiao L. Generation and Propagation of Internal Layers // Nonlinear Anal. 1988. Y. 12. N 1. P. 19-41.

28. Perko L. Differential Equations and Dynamical Systems. New York: Springer, 2001.

29. Amann H. Periodic Solutions of Semilinear Parabolic Equations, Nonlinear Analysis: a Collection of Papers in Honor of Erich Rothe. New York: Academic, 1978, pp. 1-29.

30. Sattinger D. Monotone Methods in Elliptic and Parabolic Boundary Value Problems //Indiana Univ. Math. J. 1972. V. 21. N 11. P. 979-1001.

31. Fife P., Tang M. Comparision Principles for Reaction-Diffusion systems: Irregular Comparision Functions and Applications to Question of Stability and Speed Propagation of Disturbances // J. Diff. Equations. 1981. V. 40, P. 168185.

32. Нефедов H.H., Омельченко O.E., Рекке Л. Стационарные внутренние слои в интегро-дифференциальной системе реакция-адвекция-диффузия // Журнал вычислительной математики и математической физики, 2006, т. 46. N 4. с. 623-645.

33. СПИСОК ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

34. Никитин А.Г. Неустойчивость контрастных пространственных структур типа "всплеска" в системе реакции-диффузии // Журнал вычислительной математики и математической физики. 1992. Т. 31. N 3. С. 443-452.

35. Васильева А.Б., Никитин А.Г., Петров А.П. Асимптотический метод исследования контрастных структур и его приложения к теории гидромагнитного динамо // Математическое моделирование. 1995. Т. 7. N 2, С. 61 -71.

36. Vasil'eva, A. Nikitin and A. Petrov Stability of contrasting solutions of nonlinear hydromagnetic dynamo equations and magnetic fields reversals in galaxies // Geophys. Astrophys. Fluid Dynamics. 1995. V. 78. P. 261 279.

37. Васильева А.Б., Никитин А.Г. К вопросу об устойчивости периодических контрастных структур в пространственно двумерном случае // Дифференциальные уравнения. 1996. Т. 32. N 10. С. 1355-1361.

38. Никитин А.Г. О главной собственной функции одной сингулярно возмущенной задачи Штурма-Лиувилля // Журнал вычислительной математики и математической физики. 1999. Т. 39. N 4. С. 558-591.

39. Никитин А.Г., Петров А.П. О предельном переходе по малому параметру для собственных значений сингулярно возмущенной задачи Штурма-Лиувилля // Дифференциальные уравнения. 1999. Т. 35, N 6. с. 843-845

40. Нефедов Н.Н., Никитин А.Г. Асимптотический метод дифференциальных неравенств для сингулярно возмущенных интегро-дифференциальных уравнений // Дифференциальные уравнения. 2000. Т. 36. N 10. С. 13981404.

41. Нефедов Н.Н., Никитин А.Г. Асимптотическая устойчивость контрастных структур типа ступеньки в сингулярно возмущённых интегродифференциальных уравнениях в двумерном случае // Математическое моделирование. 2001. Т. 13. N 12. С. 65-74.

42. Нефедов Н.Н., Никитин А.Г., Ур аз гиль дина Т. А. Задача Коши для интег-ро-дифференциального уравнения Вольтерра // Журнал вычислительной математики и математической физики. 2006. Т. 46. N 5. С. 805-812.

43. Нефедов Н.Н., Никитин А.Г. Задача Коши для интегро-дифференциального уравнения Фредгольма // Журнал вычислительной математики и математической физики. 2007. Т. 47. N 4. С. 655-664.

44. Nefedov N.N., Nikitin A.G., Recke L. Moving Internal Layers in the Singular Perturbed Integro-Parabolic Reaction-Diffusion-Advection Equations // Preprint Nr. 2007-22. Humboldt University of Berlin, Institute of Mathematic. P. 1-17.

45. Nikitin A. Boundary and internal layers in the integro-differential equations // Nonlinear Partial Differential Equations. International Conference. (Book of Abstracts). P. 148. Alushta, 2003.

46. Nefedov N.N., Nikitin A. Singular perturbed integro-differential equations with balanced nonlinearly // Nonlinear Partial Differential Equations. International Conference (Book of Abstract). P. 67. Alushta, 2005.

47. Никитин А.Г. Контрастные структуры в сингулярно возмущенных интегро-дифференциальных уравнениях // Седьмая Крымская международная математическая школа. Метод функций Ляпунова и его приложения (Тезисы докладов). С. 114. Алушта, 2004.

48. Nikitin A.G. Contrast structures in the integro-differential equations with balanced nonlinearity // VI International Congress on Mathematical Modeling. Book of Abstract. P. 523. Nizhny Novgorod, 2004.

49. Никитин А.Г. Нелокальные сингулярно возмущенные уравнения "реакция-диффузия-адвекция" // Ломоносовские чтения 2008. Подсекция "Теоретическая и математическая физика" (Сборник расширенных тезисов докладов). С. 156-161. Москва, 2008.

50. Васильева А.Б., Никитин А. Г. Об устойчивости контрастных структур в пространственно двумерном случае // Математические модели и методы в социальных науках (Труды вторых математических чтений МГСУ). С. 2425, М.: МГСУ, 1994.

51. Никитин А. Г., Петров А.П. Структура спектра сингулярно возмущенной задачи Штурма-Лиувилля // Математические методы и приложения (Труды третьих математических чтений МГСУ). С. 60-62. М.: МГСУ, 1995.

52. Никитин А. Г. О главной собственной функции сингулярно возмущенной задачи Штурма-Лиувилля // Математические методы и приложения (Труды четвертых математических чтений МГСУ). С. 25-28. М.: МГСУ, 1996.

53. Никитин А. Г., Давыдова М.А. Устойчивость контрастной структуры типа "ступеньки" в случае слабой зависимости правой части от первой производной // Математические методы и приложения (Труды пятых математических чтений МГСУ). С. 50-54. М.: МГСУ, 1997.

54. Нефедов Н. Н., Никитин А. Г. Асимптотика решения сингулярно возмущенной интегро-дифференциальной краевой задачи // Математические методы и приложения (Труды шестых математических чтений МГСУ). С. 22-23. М.: МГСУ, 1999.

55. Никитин А.Г. Сингулярно возмущенные интегро-дифференциальные уравнения в критическом случае // Математические методы и приложения (Труды десятых математических чтений МГСУ). С. 20-23. М.: МГСУ, 2003.

56. Никитин А.Г. О корректности вырожденных задач для сингулярно возмущенных интегро-дифференциальных уравнений // Математические методы и приложения (Труды одиннадцатых математических чтений МГСУ). С. 37. М.: МГСУ, 2004.

57. Нефедов Н. Н., Никитин А. Г. Начальная задача для интегродифференциального уравнения Фредгольма с малым параметром при производной // Математические методы и приложения (Труды пятнадцатых математических чтений РГСУ). С. 92-96. М.: РГСУ, 2007.

58. Нефедов Н. Н., Никитин А. Г., Рекке Л. Движущиеся фронты в сингулярно возмущенной интегро-параболической задаче реакция-диффузия-адвекция // Математические методы и приложения (Труды шестнадцатых математических чтений РГСУ). С. 123-132. М.: РГСУ, 2008.