CC BY
32
КРАТКИЕ СООБЩЕНИЯ
УДК 517.9
МОДЕЛЬ ТЕРМОКОНВЕКЦИИ НЕСЖИМАЕМОЙ ВЯЗКОУПРУГОЙ ЖИДКОСТИ НЕНУЛЕВОГО ПОРЯДКА. ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ
О.П. Матвеева
Целью статьи является численное исследование решения начально-краевой задачи для модели термоконвекции ненулевого порядка. Рассматривается система, которая моделирует эволюцию скорости, градиента давления и температуры несжимаемой вязкоупругой жидкости Кельвина - Фойгта ненулевого порядка. Используя метод Галер-кина, разработан алгоритм численного решения начально-краевой задачи для системы, моделирующей плоскопараллельную термоконвекцию несжимаемой жидкости ненулевого порядка, и реализована программа для персональных компьютеров нахождения численного решения указанной задачи. Получена графическая иллюстрация численного решения системы при заданных параметрах. Проведенное исследование основано на результатах теории полулинейных уравнений соболевского типа, поскольку начальнокраевая задача для соответствующей системы дифференциальных уравнений в частных производных сводится к абстрактной задаче Коши для уравнения соболевского типа.
Ключевые слова: уравнение соболевского типа, термоконвекция, несжимаемая вязкоупругая жидкость.
Система уравнений
/ к (1 - ХЧ2)уг = VЧ2у — (V ■ в1 У2—1 - 9$ — Ур + ¡,
1=1
0 = У (У ■ V),
дъи
= V + а—г, а е к_, I = 1, к,
дЬ
„ 01 = жУ20 — V ■ У0 + V ■ q
моделирует эволюцию скорости V = ^1,...^п), VI = Vi(x,t), градиента давления Ур = (р\,...,рп), pi = р^х, Ь) и температуры 0 = 0(х,Ь) несжимаемой вязкоупругой жидкости Кельвипа-Фойгта порядка к > 0 х е П С Кп, п = 2, 3, 4 — ограниченная область с границей д П класс а С[1]. Парамет ры Л е К, V е К+ и ж е К+ характеризуют упругость, вязкость и теплопроводность жидкости соответственно; 9 е К+ — ускорение свободного падения; вектор q = (0,..., 0,1) — орт в Кп. Параметры вг е К+ , I = 1, к определяют время ретардации (запаздывания) давления. Свободный член f = (/\,...,/п), ^ = ^(х) отвечает внешнему воздействию па жидкость. Начально-краевые задачи для моделей тер-мокопвекции были изучены в [2-4].
В области П = [0, п] х [0, п] рассмотрим систему (1) в гаде (к = 1)
(1 — ЛУ2^1 = VУ2v — (V ■ У^—вУ2 — — 9q0 — Ур,
0 = У ■ V,
д— (2) —— = V + а—, а е К_,
(1)
дЬ
01 = жУ20 — V ■ У0 + V ■ q.
Введем функцию тока, определенную уравнениями VI (2)
д(ф, У2ф)
= Ц ,v2 = - дХ, ГДвф = ф(X,y,t).
d (x,y)
dy
dx
dw1 дф
~wr = тг + awi , dt dy
dw2 дф _
= - ^—+ aw2 , a € M-, dt dx
et = V + §M + ?X.
d(x,y) dx
Для системы (3) поставим задачу Коши - Бенара
ф^, 0, t) = ф^, п, t) = V2ф(x, 0, t) = V2ф(x, п, t) = 0, ф(о, y, t) = ф(п, y, t); V2ф(0, y, t) = V2ф(п, y, t),
9(x, 0, t) = 9(x, п, t) = 0; 9(0, y, t) = 9(п, y, t),
9(x,y, 0) = 9o(x, y); ф^,у, 0) = фо^,у), wi(x, 0, t) = wi(x, п, t) = 0; wi(0, y, t) = wi(n, y, t), к wi(x,y, 0) = wio(x,y); i = 1, 2.
(3)
(4)
Целыо данной статьи является проведение вычислительного эксперимента по исследованию решения задачи (3), (4).
Вычислительный эксперимент
На основе теоретических результатов [5] для системы (3), моделирующей эволюцию скорости, градиента давления pi температуры несжимаемой вязкоупругой жидкости Кельвипа-Фойгта, в системе компьютерной математики Maple разработана программа [6], которая позволяет:
1. По заданным коэффициентам а, в, Х, ж, v па основе метода Галеркина численно находить решение системы.
2. Получить графическое изображение решения системы.
Для реализации вычислительных алгоритмов программы использовались встроенные функции и стандартные операторы языка программирования Maple. Для получения графического изображения подключен пакет plots.
(4) (3)
целыо выберем в качестве базисных функций метода Галеркина собственные функции следующей задачи
-V2^ = Х^,
<p(x, 0) = <^(x, п) = 0, <p(0,y) = р(п,у).
Нетрудно получить akl = П sin(2kx)sin(ly), вы = П cos(2kx)sin(ly), il = П sin(ly) -
ортонормированное в смысле L2 множество собственных функций. Галеркинское приближение к решению задачи (4) для системы (3) возьмем в виде ф = a(t)aii, 9 = b(t)eii + c(t)j2-, wi = d(t)eii + f (t)Y2, w2 = g(t)eii + h(t)Y2.
На следующем этапе, умножив скалярно уравнения системы (3) на функции aii, ви, 12-получим систему обыкновенных дифференциальных уравнений. Зададим начальные условия из окрестности точки пуль. Затем численно решим задачу Копти для системы обыкновенных дифференциальных уравнений с заданными начальными условиями.
Пример 1. Требуется найти численное решение задачи (3), (4) при заданных коэффициентах V = 2 а = —1 в = 2, Л = 1, ж = 1, а также получить графическое изображение этого решения.
Умножим скалярно уравнения системы (3) на собственные функции ац,вц,^2- Получим систему дифференциальных уравнений
' —5а(Ь)(1 + 5Л) — 2^а(Ь) + 10вд(Ь) + 19.6Ь(Ь) = 0 ¿(Ь) — ай(Ь) = 0, д(Ь) + 2а(Ь) — ад(Ь) = 0
< Ь(Ь) + 5жЬ(Ь) + а(Ь)с(Ь) — 2а(Ь) = 0
/ (Ь) — а/(Ь) = 0, Н(Ь) — аН(Ь) = 0
ч с(Ь) + 4жс(Ь) — а(Ь)Ь(Ь) = 0.
Зададим начальные условия из окрестности точки нуль. Пусть а(0) = 0,2, Ь(0) = 0,1,с(0) = 0,2,й(0) = 0,1,/(0) = 0, 3,д(0) = 0,2,Н(0) = 0, 2. Решим задачу Коши для данной системы уравнений. Графическая иллюстрация решения системы представлена па рисунке. Результаты численного решения частично приведены в таблице.
Решение системы при а = — 1, в = 2,А = 1,и = 2, ж = 1
Автор выражает признательность профессорам Т.Г. Сукачевой и Г.А. Свиридюку за внимание к данным исследованиям и обсуждение результатов.
Литература
1. Осколков, А.П. Начально-краевые задачи для уравнений движения жидкостей Кельвипа-Фойгта и жидкостей Олдройта / А.П. Осколков // Тр. мат. итт-та АН СССР. - 1988. - № 179. - С. 126-164.
2. Свиридток, Г.А. Разрешимость задачи термокоттвекции вязкоупругой несжимаемой жидкости / Г.А. Свиридток // Изв. вузов. Математика. - 1990. - № 12. - С. 65-70.
Таблица
Численное решение системы с начальными условиями а(о) = 0, 2, 6(0) = 0,1, с(0) = 0, 2, ¿(0) = 0,1,/(0) = 0, 3, д(0) = 0, 2, к(0) = 0, 2
ь а(Ь) ь(ь) С(ь) ¿(ь)
од 0,1854240039 0,0886898129 0,1354017426 0,0904837418
0,3 0,1525995343 0,0730387274 0,0624894731 0,0740818220
0,4 0,1357883884 0,0663610716 0,0426287976 0,0670320046
0,6 0,1035998576 0,0535727678 0,0200176562 0,0548811636
0,8 0,0749935012 0,0414719893 0,0095013049 0,0449328964
1,0 0,0509756644 0,0305398288 0,0045377105 0,0367879441
1,2 0,0317854606 0,0211874342 0,0021642916 0,0301194211
1,4 0,0171784838 0,0135878727 0,0010216194 0,0246596925
1,6 0,0066299269 0,0077117905 0,0004733233 0,0201189652
ь / (ь) д(ь) Чь)
0,1 0,2714512254 0,1442806994 0,1809674836
0,3 0,2222454661 0,0569715237 0,1481636441
0,4 0,2010960138 0,0241357331 0,1340640092
0,6 0,1646434907 -0,023355431 0,1097623271
0,8 0,1347986891 -0,051193908 0,0898657927
1,0 0,1103638322 -0,064457824 0,0735758881
1,2 0,0903582633 -0,067515694 0,0602388422
1,4 0,0739790888 -0,063933144 0,0493193925
1,6 0,0605689548 -0,056483131 0,0403793032
3. Свиридток, Г.А. Некоторые математические задачи динамики вязкоупругих несжимаемых сред / Г.А Стшридток, Т.Г. Сукачева /'/ Весттт. МаГУ . - 2005. - № 8. - С. 5-33.
4. Сукачева, Т.Г. Исследование математических моделей несжимаемых вязкоупругих жидкостей: дис. ... д-ра физ.-мат. паук /Т.Г. Сукачева; Новгород, гос. утт-т. - Великий Новгород, 2004. - 249 с.
5. Сукачева, Т.Г. Задача термокоттвекции несжимаемой вязкоупругой жидкости Кельвитта-Фойгта ненулевого порядка /Т.Г. Сукачева, О.П. Матвеева /'/ Изв. вузов. Математика. - 2001. - .№ 11 (474). - С. 46-53.
6. Численное решение ттачальтто-краевой задачи для модели термокоттвекции несжимаемой вязкоупругой жидкости ненулевого порядка / Матвеева О.П. (Ш.-); правообладатель: Матвеева О.П. (Б,и). - №2012612862, зарегистр. 22.03.2012, Реестр программ для ЭВМ.
Ольга Павловна Матвеева, кафедра алгебры и геометрии, Новгородский государственный университет им. Ярослава Мудрого, (г. Великий Новгород, Российская Федерация), oltan.72@mail.ru
Bulletin of the South Ural State University. Series «Mathematical Modelling, Programming & Computer Software:»,
2013, vol. 6, no. 1, pp. 134-138.
MSC 35R20, 35G25, 35Q35, 35Q72, 76A05
Model of Thermoconvection of Incompressible Viscoelastic Fluid of Nonzero Order. Computational Experiment
O.P. Matveeva, Novgorod State University, Velikiy Novgorod, Russian Federation, olt. an .72® mail .r u
The purpose of this paper is the numerical investigation of the solution of the initial-boundary value problem for the model of thermal convection of the nonzero order. We consider the system that models the evolution of the velocity gradient of the pressure and temperature of the incompressible viscoelastic Kelvin-Voigt fluid of nonzero order. Using the Galerkin method, the algorithm of the numerical solution of the initial-boundary value problem for the system modeling plane-parallel thermal convection of the incompressible fluid of the nonzero order is created, and the program for personal computers to find numerical solutions of this problem is implemented. A graphical illustration of the numerical solution with the given parameters is obtained. The study was based on the results of the theory of semi-linear Sobolev type equations, because the initial boundary value problem for the corresponding system of differential equations in partial derivatives is reduced to the abstract Cauchy problem for the Sobolev type equation.
Keywords: sobolev type equation, thermal convection, incompressible viscoelastic fluid.
References
1. Oskolkov A.P. Initial-Boundary Value Problems for Equations of Motion Kelvin-Voight and Oldroyd Fluids [Nacharno-kraevye zadachi dlya uravneniy dvizheniya zhidkostey Kel’vina-Foygta i zhidkostey OldroytaJ. Trudy Mat. In-ta AN SSSR, 1988, no. 179, pp. 126-164.
2. Sviridyuk G.A. Solubility of the Thermal Convection of Viscoelastic Incompressible Fluid [Razreshimost’ zadachi termokonvektsii vyazkouprugoy neszhiinaeinoy zhidkostij. Russian Mathematics (Izvestiya VUZ. Matematika), 1990, no. 12, pp. 65-70.
3. Sviridyuk G.A., Sukacheva T.G. Some Mathematical Problems of the Dynamics of Viscoelastic Incompressible Media [Nekotorye inatematicheskie zadachi dinamiki vjazkouprugih neszhiinaeinyh sredj. Vestnik MaGU, 2005, no. 8, pp. 5-33.
4. Sukacheva T.G. Issledovanie rnaternaticheskikh modeley neszhimaemykh vyazkouprugikh zhidkostey: dis. ... d-ra fiz.-mat. na.uk [The Study of Mathematical Models of Incompressible Viscoelastic Fluids: dis. Dr. Science]. Velikiy Novgorod, 2004. 249 p.
5. Sukacheva T.G., Matveeva O.P. The Problem of Thermal Convection of an Incompressible Viscoelastic Kelvin-Voigt Fluid of Nonzero Order [Zadacha termokonvekcii neszhiinaeinoj vjazkouprugoj zhidkosti Kel’vina-Fojgta nenulevogo porjadkaj. Russian Mathematics (Izvestiya VUZ. Matematika), 2001, no. 11, pp. 46-53.
Поступила в редакцию 81 августа 2012 г.
