Метод последовательных приближений в параболической начально-краевой задаче Текст научной статьи по специальности «Математика»



Серия «Математика» 2013. Т. 6, № 2. С. 77—83

Онлайн-доступ к журналу: http://isu.ru/izvestia

УДК 517.956

Метод последовательных приближений в параболической начально-краевой задаче *

В. А. Терлецкий

Иркутский государственный университет

Е. В. Тучнолобова

Иркутский государственный университет

Н. Ю. Ульянова

Иркутский государственный университет

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

Ключевые слова: начально-краевая задача; параболическое уравнение; фундаментальное решение; метод последовательных приближений.

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

* Работа выполнена при финансовой поддержке РФФИ (проекты 11-01-00713, 1201-98009) и федеральной целевой программы "Научные и научно-педагогические кадры инновационной России" на 2009-2013 гг.

1. Введение

При исследовании задач оптимального управления системами и уравнениями с распределенными параметрами наиболее часто прибегают к обобщенным решениям по Соболеву. По-видимому, это объясняется, прежде всего, простотой и универсальностью такого подхода. Однако указанный способ имеет и существенные недостатки, отмеченные, в частности, в [2, с. 575]. Например, это необходимость отдельного, не связанного с обоснованием существования обобщенного решения, нахождения оценок его роста относительно входных данных; доказательства формулы интегрирования по частям; вывод сопряженной задачи; выяснение и обоснование аналитических свойств решений исходной и сопряженной задач.

Указанных недостатков удается избежать при конструировании обобщенного решения дифференциальной задачи как решения эквивалентной интегральной задачи. Именно такой подход является общепринятым в системах обыкновенных дифференциальных уравнений при построении решения по Каратеодори [3; 4]. Так же успешно он применяется и при обосновании решения в "широком смысле" (непрерывного, но с разрывными частными производными) в гиперболических системах одномерных полулинейных уравнений [1]. Эта идея определения обобщенного решения использовалась в гиперболических системах одномерных [5] и многомерных [6; 7] полулинейных уравнений, а также в различных начально-краевых задачах для волнового уравнения [8].

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

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

В прямоугольнике П = 5 х Т плоскости переменных (в,Ь), Б = (80,81), Т = (¿0,^1) рассмотрим уравнение теплопроводности

хг - a2(s)xss = f (x,xs,s,t) (1.1)

с начальными

х(8,Ь)= х0(8), 8 е Б (1.2)

и граничными

х^) = д0^), х(8,) = д1^), t е Т (1.3)

условиями.

Здесь t - время, 8 - пространственная переменная, х = х(,з,{) -искомое решение; а = а(,з) - гладкая на Б функция, причем 0 < а0 < а(,з) < ао < ,з е Б; f = (х, у, 8, ¿) непрерывна по Липшицу по

переменным (х,у) в Я2 при фиксированных е П и суммируема в

П при фиксированных значениях (х,у) е Я2; функции х0 = х0^), дг = д0^), г = 0,1 суммируемы в Б и Т соответственно.

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

3. Интегральный эквивалент

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

^ = а(в)

с начальным условием ,в(п0) = 80.

В силу положительности функции а = а(,в), е Б, функция = ,в(п) является монотонно возрастающей на интервале (^0, П1), где П1 выбирается из условия ,3(^1) = 81. Поэтому существует и обратное отображение П = П^), переводящее отрезок [80, в отрезок [щ,щ]. Кроме того, на интервале Б справедливо равенство

¿п 1

а(8)

и Ф0) = П0, = П1.

Нетрудно видеть, что

хп = х3а(8), хпп = xssa2(s) + х^^а^).

Следовательно, в терминах новой пространственной переменной п уравнение (1.1) принимает вид

хь ~ хт = / (х, ф^))'5^)'*) ~

Таким образом, при естественных предположениях о гладкости, неотрицательности и ограниченности функции а = а(.в) уравнение вида (1.1) всегда может быть сведено к уравнению того же вида, но с единичным коэффициентом при второй частной производной по пространственной переменной от решения в дифференциальном операторе. Это означает, что не ограничивая общности задачи (1.1)-(1.3), далее будем считать, что в уравнении (1.1) а= 1, .в е Б. Для функций х, имеющих непрерывные первую производную по времени и вторую производную по пространственной переменной, будем использовать обозначение Dx — хг xss.

Перейдем к построению интегрального эквивалента задачи (1.1)— (1.3), используя известное (см., например, [9, 10]) фундаментальное решение для уравнения теплопроводности:

р \4t~T)

Х(з-£,1-т) =

-т)"

Действуя по схеме [10, с. 370-373], можно установить, что на гладких решениях задачи (1.1)-(1.3) функции

/те

х(в - £,г - го)х0(£)с£,

■те г те

х(/; в,г) = / X(в - £,г - г)/[£,тЩсСт

Jto ■) — те

обладают свойствами

Ох(х°; в,г)=0, (в,г) € П, х(х0; в,г0) = х°(в), в € Б,

х(х°; в0,г) = х(х°; в1,1) = 0, г € Т,

Бх(/; в,г) = /(х(в,г),хэ(в,г),в,г), (в,г) € п, х(/; в,го) = 0, в € Б,

х(/; во,г) = х(/; в1,г) = 0, г € Т.

Здесь функции х0 = х°(£), / = /[£,т] = /(х(£,т),хс(£,т),£,т) строятся по функциям х0 и / как их антисимметрические продолжения по пространственной переменной на всю вещественную ось:

х°(£ ± Ы) = (-1)кх0(£),

/[£ ± к1, т] = (-1)к/(х(£, т)х(£, т),£, т),

£ € Б, т € Т, I = в1 - во, к = 0,1,2,...

Для завершения построения интегрального эквивалента задачи (1.1)-(1.3) остается указать ее решения, "отвечающие"за граничные условия (1.3). Пусть

уг(£ ± 2к1) = -(-1У2к, £ € Б, к = 0,1,..., г = 0,1. В плоскости переменных (в, г) определим функции

уг(в г) = {—те х(в - £,г - т)уг(£г > 0 [0, г < 0, г = 0,1.

Аналогично предыдущим устанавливаем справедливость равенств

'1, г = з, г> 0,

Буг = 0, (в,г) € п, уг(в,го) = 0, угв,г) =

0, г = 3, г,з = 0,1.

е X 2 \

По схеме [9, с. 47] можно доказать, что функции

[ д ■

х(дг]8,г) = 1^—уг(8,1-т)дг(т)<1т обладают свойствами

Ох(д1; в,г) = о, (в, г) е п,

х(дг; в,го) = 0, в е Б, х(дг; в,,г) = = _ . .

[0, г = з, г е Т, г,3 = 0,1.

Обозначим для краткости х(д; в, г) = х(д°; в, г) + х(д1; в, г). Тогда искомый вид интегрального эквивалента следующий:

х(в, г) = х(х°; в, г) + х(д; в, г) + х(/(х, х3, в, г), в, г). (2.1)

Интегральное уравнение (2.1) уже не содержит старших производных хг и х33, которые входят в дифференциальный оператор D уравнения (1.1), что позволяет сформулировать следующее

Определение 1. Обобщённым решением задачи (1.1)-(1.3) назовем функцию х = х(в,г), удовлетворяющую уравнению (2.1).

4. Метод последовательных приближений

Равенство (2.1) естественным образом может быть положено в основу метода последовательных приближений. Однако, в общем случае, когда правая часть / уравнения (1.1) зависит от частной производной х3 решения х, его следует дополнить уравнением

х3(в, г) = х3(х°; в, г) + х3(д; в, г) + х3(/(х, х3, в, г); в, г), (3.1)

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

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

х(к+1)(в, г) = х(х0; в, г) + х(д; в, г) + х(/(х(к), х8(к),в, г); в, г),

х.3{к+1)(в,г) = х3(х°; в,г)+ х3(д; в,г) + х(/(х(к),х3{к),в,г); в,г), х(0)(в,г) = х3{0)(в,г) = 0, (в,г) е п,к = о, 1,...

Для доказательства сходимости рядов

те те

^(х(к+1)(в,г) - х(к)(в,г)),^2(ха(к+1)(в,г) - х.3(к)(в,г)), к=0 к=0

а, следовательно, и существования обобщенного решения достаточно установить сходимость рядов

те те

к=0 к=0 где приращения 5(к) и 53(к определяются равенствами б(к)(в,г) = \х(к+1\в,г) - х(к\в,г)\, (к)(в,г) = \х/к+1)(в,г) - х/к)(в,г)\.

Тогда из (2.1), (3.1) и условия Липшица для функций f по фазовым переменным следуют неравенства

гЬ гте

5(к+1\в,г) < С / Х(в - - т)[6к(£,т) + 63(к)&т)№т,

¿Ьо ■)—те

{3.1)

гЬ рте

53{к+1)(в,1) < С / X (в - - т)[6(к) (£,т) + 63 (к\£,т)№т,

¿Ьо ■) — те

в которых константа С > 0 определяется только константой Липшица.

Система неравенств (3.2) обосновывает существование и единственность обобщенного решения задачи (1.1)—(1.3), которое вместе со своей частной производной х3 по пространственной переменной является точкой равновесия отображения (2.1) и (3.1). Справедливость интегральных уравнений (2.1), (3.1) позволяет проанализировать аналитические свойства обобщенного решения и получить оценки его роста относительно входных данных , д0, д1 и f задачи (1.1)-(1.3).

Список литературы

1. Рождественский Б. Л. Системы квазилинейных уравнений и их приложения к газовой динамике / Б. Л. Рождественский, Н. Н. Яненко. - М. : Наука, 1978. -688 с.

2. Васильев Ф. П. Методы оптимизации / Ф. П. Васильев. - М. : Факториал Пресс, 2002. - 824 с.

3. Сансоне Дж. Обыкновенные дифференциальные уравнения. Т. 2 / Дж. Сан-соне. - М. : Иностр. лит., 1954. - 415 с.

4. Коддингтон Э. А. Теория обыкновенных дифференциальных уравнений / Э. А. Коддингтон, Н. Левинсон. - М. : Иностр. лит., 1958. - 474 с.

5. Терлецкий В. А. Обобщенное решение одномерных полулинейных гиперболических систем со смешанными условиями / В. А. Терлецкий // Изв. вузов. Математика. - 2004. - № 12. - С. 75-83.

6. Терлецкий В. А. Обобщенное решение многомерных полулинейных гиперболических систем / В. А. Терлецкий // Изв. вузов. Математика. - 2001. - № 12. -С. 68-76.

7. Терлецкий В. А. Обобщенное решение в задачах оптимального управления гиперболическими системами / В. А. Терлецкий // Изв. вузов. Математика. -2004. - № 4. - С. 68-78.

8. Терлецкий В. А. Обобщенное решение нелинейного волнового уравнения с нелинейными граничными условиями первого, второго и третьего рода / В. А. Терлецкий, Е. А. Лутковская // Дифференц. уравнения. - 2009. - Т. 45, № 3. -С. 403-415.

9. Годунов С.К. Уравнения математической физики / С. К. Годунов. - М. : Наука, 1979. - 392 с.

10. Михайлов В. П. Диффернциальные уравнения в частных производных /В. П. Михайлов. - М. : Наука, 1983. - 424 с.

V. A. Terletsky, E. V. Tuchnolobova, N. Yu. Ulyanova Method of successive approximations of parabolic initial-boundary problem

Abstract. The existence and uniqueness of generalized solution of an initial-boundary problem of parabolic equation with discontinuous enterance data with respect to independent variables and nonlinear right part of equation with respect to spacial variables is received by method of successive approximations.

Keywords: initial-boundary problem, parabolic equation, fundamental solution, method of successive approximations

Терлецкий Виктор Анатольевич, кандидат физико-математических наук, доцент, Институт математики, экономики и информатики, Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1, тел.: (3952)52-12-82 ([email protected])

Тучнолобова Екатерина Владимировна, аспирант, Институт математики, экономики и информатики, Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1, тел.: (3952)52-12-82 ([email protected])

Ульянова Надежда Юрьевна, старший преподаватель кафедры, Институт математики, экономики и информатики, Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1, тел.: (3952)52-12-82 ([email protected])

Terletsky Viktor, Ph. D., Irkutsk State University, 1, K. Marks St., Irkutsk, 664000, Phone: (3952)52-12-82 ([email protected])

Tuchnolobova Ekaterina, Irkutsk State University, 1, K. Marks St., Irkutsk, 664000, Phone: (3952)52-12-82 ([email protected])

Ulyanova Nadegda, Irkutsk State University, 1, K. Marks St., Irkutsk, 664000, Phone: (3952)52-12-82 ([email protected])