Серия «Математика»
Том 1 (2007), № 1, С. 113-117
Онлайн-доступ к журналу: http://isu.ru/izvestia
ИЗВЕСТИЯ
Иркутского государственного университета
УДК 532.5.013, 517.948
Ом и
существовании непрерывных решении в одной модельной задаче теории пограничного слоя *
А. И. Дрегля ([email protected])
Иркутский государственный университет
Аннотация. Доказана теорема существования непрерывных решений краевой задачи для нелинейного дифференциального уравнения третьего порядка, возникающего в теории пограничного слоя.
Ключевые слова: краевая задача, нелинейные дифференциальные уравнения, принцип Шаудера, компактность.
Рассмотрим краевую задачу
х"(г) + м(х(г),г)х (г) = о, а <г < в, (1)
x(a) = a, x (а) = b, x(ß) = c. (2)
К такой задаче сводятся некоторые задачи теории пограничного слоя в
гидромеханике [5], [9] и теории полимеров [1].
Пусть
А. Функция M(x,t) определена и непрерывна в области
D = {x,t\ \x\ < \a\ + \b\\ß\ + \c — a — bß\, а < t < ß} ,
m = min M (x, t) M = max < max M (x, t), 0 > . x,teD [x,teD J
В частном случае (M(x,t) = x, а = 0, ß = 1) в книге В. А. Треногина ([12] стр. 412 — 413) было доказано существование решения этой задачи. В этой заметке доказано существование решения краевой задачи (1), (2) на конечном отрезке [а, ß] в общем случае, когда функция M(x, t) любая непрерывная функция, удовлетворяющая условию А.
* Работа частично поддержана Kluber Lubrication KG и Дублинским технологическим университетом.
114 А. И. ДРЕГЛЯ
Теорема 1. Пусть выполнено условие А. Тогда краевая .задача (1), (2) на конечном отрезке [а, в] имеет в классе С^,^ решение.
Доказательство:
Сведем задачу (1), (2) к нелинейному интегральному уравнению (6). Для этого введем обозначение и = х'' и перепишем уравнение (1) в виде
и' + М (х,Ь)и = 0
с условием
и(а) = х''(а), где х''(а) определим позже. Тогда
ь
- У м (х(з),з)аз
х''(Ь) = х''(а)е а . (3)
Путем интегрирования (3) от а до Ь с учетом условия х'(а) = Ь получим
с
г -/ м(х(з),з^з
х'(Ь) = Ь + х''(а) е а йа. (4)
Интегрируя (4) с учетом условия х(а) = а получим
с
} П -/ М(х(з),з^з
х(Ь) = а + ЬЬ + х''(а) е а йайп- (5)
а а
Полагая в (5) Ь = в и используя условие х(а) = а найдем
с — а — Ьв
х (а) = - _
в П -/ М(х(з),з)4з
/ / е а йайг/
а а
Итак, интегральное уравнение эквивалентное задаче (1), (2) имеет вид
с
гп -/м(х(з),з)аз / / е а йайг/ х(Ь) = а + ЬЬ + (с — а — Ьв) —-с-. (6)
в п -/ М(х(з),з)4з
/ / е а йайг/
а а
Введем обозначение
с
п -/ М(х(з),з)4з
Рх(ц) = е а йа, а <п < в- (7)
а
О СУЩЕСТВОВАНИИ НЕПРЕРЫВНЫХ РЕШЕНИЙ 115
Очевидно, функция Рх(п) положительная при ц > а. С учетом обозначения (7) интегральное уравнение (6) принимает вид
г
/ Рх(П)йП
х(Ь) = а + ЬЬ + (с — а — Ьв)а-= Рг(х), а<Ь<в- (8)
¡Рх(П)йП
а
При а < Ь < в и Ух € Сав] имеем неравенство
о < а Рх шч < г РхШп
Следовательно, для Ух € С[а,в]
\\Рг(х)\\< |а| + \Ь\\в\ + \с — а — Ьв\ = г. (9)
Тем более, при \\х\\ < г
\\Рг(х)\\< г,
т.е. интегральный оператор Рг : С[а в] ^ С[а в] переводит шар \\х\\ < г самого в себя.
Докажем, что оператор Рг вполне непрерывен. Для этого воспользуемся теоремой Арцелла ([12] стр. 207 - 209) о предкомпактности множеств в С[а в]. Ввиду оценки (9) образ множества \\х\\ < г при отображении Рг(х) равномерно ограничен постоянной г. Далее
¿2
\Ри (х) — а (х)\<\Ь\\Ь1 — Ь2\ + \С— а — Ьв\ \/ РхШп\- (10)
\ I Рх(п)йп\ ¿1
а
Отметим, что так как функция Рх(Ь) положительна, то
¿2
= У Рх(п)йп, и < Ь2.
¿2 ¿2 / Рх(п)йп
¿1 ¿1
Рассмотрим два случая: Случай I: (М > 0)
в П вема
У РхШп > УI е-м(а-а)йайп = I е— (е-Мп — е-Ма)йп
еМа
~М~
__1_ е-Мв + _!_ е-Ма
ме + ме
+ М(в — а)
в
116 А. И. ДРЕГЛЯ
{M(в — а) + e-M(в-а) — 1} > 0, УМ> 0.
M 2
Случай II: (M = 0)
в в п в
J Fx(n)dn — J J dadn = J(n — a)dn =
a a a a
= 2(в2 — a2) — а(в — a) = (в — а)(в + *»> = ^
Следовательно, существует const Ci > 0 такая, что при Ух € D, справедливо неравенство
в
J Fx(n)dn — Ci > 0. (11)
a
С другой стороны,
п в
Fx(n)dn < J e-m(a-a)da < J e-m(a-a)da =
aa
в — a, m = 0 -, m> 0
m (е\т\в — e]m]a), m < 0.
Следовательно, существует const C2 > 0 такая, что при Ух, t € D, справедливо неравенство
0 < Fx(n) < C2. (12)
В итоге неравенства (11), (12) дают такое уточнение оценки (10):
\Ptl (х) — Pt(x)\ < |6||ti — t2| + |С — a— Ьв 1 C2|ti — t2| = l\ti — t2|,
Ci
где l = \ь\+\с—а-Ьв]{ C2. Итак, существует константа l такая, что при У(х, ti, t2) € D1
Pti (x) — Pt2 (x^ < l^i — t2|.
При этом константа l не зависит от х, ti, t2 из D. Мы доказали, что образ Pt(х) шара S(0, r) равноограничен и равностепенно непрерывен в C^^] и следовательно, на основании теоремы Арцелла предкомпактен
в C[a,0]-
Следовательно, оператор P вполне непрерывен. Итак, мы доказали, что оператор P удовлетворяет условиям теоремы Шаудера [12]. Поэтому интегральное уравнение (6) имеет в Ca,в] решение.
О СУЩЕСТВОВАНИИ НЕПРЕРЫВНЫХ РЕШЕНИЙ 117
Отметим, что ввиду структуры оператора P производные dip
— ,% = 0 I, 2, 3 Ух € см
при Ух(Ь) € С^ав] будут непрерывны. Поэтому решение уравнения (6) будет три раза непрерывно дифференцируемо. Примеры численного решения подобных краевых задач теории погранслоя были рассмотрены нами в работах [8], [3], [1], [2].
Автор благодарит проф. Д. Гильберт и проф. Н. А. Сидорова за полезные советы и внимание.
Список литературы
1. Aliona Dreglea, Some analytical solutions in the theory of modeling polymers. Reports in Mathematics, School of Math. Sc., Dublin Inst. of Tech., No. 2, 2007, available online: http://www.maths.dit.ie/preprints/2007_4_1.html.
2. Aliona Dreglea, Model of melt spinning process. Ph.D. thesis. Dublin Inst. of Tech., Dublin.- 2007, 150 p.
3. Дрегля А.И., Пимшина Л.П., О применении сеток Шишкина в одной сингулярной задаче с малым параметром, Труда 3ей межвузовской занальной конференции, посвященной памяти проф. Б.А.Бельтюкова - Иркутск: Изд-во Иркут. гос. пед. ун-та, 2007.- C.43-46.
4. Glauert M.B., Lighthill M.J., The axisymmetric boundary layer on a long thin cylinder, Proc. R. Soc. London, 1955, 320.- P. 188-203.
5. Schlichting H. Boundary Layer Theory, 7 th edition, McGraw Hill.- 1951.
6. Farrel P. A., Hegarty A. F, Miller J. J. H, O'Riordan E, Shishkin G. I., Robust Computational Techniques for Boundary Layers, Chapman and hall CRC, Florida, USA.- 2000.
7. Петровский Г.И. Лекции об уравнениях с частными производными М: Наука.-1966.
8. Dreglea A.I., Shishkin G.I. Robust numerical method for a singularly perturbed equation with unboundedly growing convective term at infinity, Proceedings of International Conference on Computational Mathematics, 2004, June, 21-25, Novosibirsk.- 2004.- p. 83-87.
9. Ландау Л.Д., Лифшиц Е.М. Гидродинамика, Глава 4, М: Наука.- 1988.
10. Джозеф Д. Устойчивость движений жидкости, М. Мир.- 1981.
11. Дрегля А.И. Об одной нелинейной задаче, возникающей в плоской модели полимеров, Труды 13ой Байкальской международной школы-семинара "Методы оптимизации и приложения", ИСЭМ СО РАН, Иркутск. т.3.- 2005.
12. Треногин В. А. Функциональный анализ, М: Наука.- 1980, стр. 408-418.
A. I. Dreglea
Continuous solutions in boundary layer problem
Abstract. Existence theorem for nonlinear ordinary differential equation with boundary conditions is proved. Such equations appear in boundary layer theory applied to melt spinning process.