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

УДК 532.5.01

Точные решения модельного уравнения динамики

вязкой среды в канале

Зайцев А.В., Кудашов В.Н. zai_@inbox.ru

Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики.

Институт холода и биотехнологий

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

Ключевые слова: уравнение Бюргерса, течение вязкой жидкости, точное решение.

Precise Solutions of a Model Dynamical Equation

of the Viscous Medium in the Channel

Zaitsev A.V., Kudashov V.N. zai_@inbox.ru

St. Petersburg National Research University of Information Technologies, Mechanics and Optics.

Institute of Refrigeration and Biotechnologies

Here is offered the use of Burgers equation as the model equation when developing a technique of numerical research of convergence and stability parameters of difference schemes at calculation of the viscous medium dynamics by the grids method. Here are given precise solutions of the presented equation.

Keywords: Burgers equation, current of viscous fluid, precise decision.

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

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

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

Запишем уравнение Бюргерса

ди ди д2и ^

д£ дх дх2

Здесь аналогами физических величин и функций процесса динамики вязких сред являются: и - скорость потока; £ - время; х - координата вдоль потока; ^ - вязкость.

Мы хотим построить решения уравнения Бюргерса при х е [0, Ь], £ е [0, то), с граничными условиями

и(0,£) = и(Ь,£) = 0, £ е [0, то). (2)

Для построения таких решений воспользуемся преобразованием Коула-Хопфа [2], сводящее нелинейное уравнение (1) к уравнению теплопроводности, являющееся линейным.

Пусть функция у(х,£) является решением уравнения

ду д2у

Тогда функция

^ 1 ду ,, ч

и=- 2|а-—. 4

у дх

удовлетворяет уравнению Бюргерса. Преобразование (4) называется преобразованием Коула-Хопфа.

Чтобы найти решения уравнения с граничными условиями (2), как видно из (4), достаточно найти функцию -и(х,г), при х е [0, Ь], г е [0, то), удовлетворяющую уравнению (3) и граничным условиям

= = (е[0,оо). (5)

При этом необходимо, чтобы выполнялось неравенство

у(х, г) Ф 0, х е [0,Ь], г е [0, то). (6)

Используя метод Фурье (см. [3]) нетрудно видеть, что существует счётное множество функций, удовлетворяющих уравнению (3) и граничным условиям (5):

уп{х,г) = е~^соз(Хпх), Хп = = 1,2,.... (7)

Ь

Так как |г>п(х, г)| < 1 при г > 0, то функции

г;п(х, г) = ^п(х,г) + С, (8)

являются решениями уравнения (3) с граничными условиями (5). Неравенство (6) очевидно выполнено, если константа С > 1.

Так как дуп/дх = -Хпе"^^вт(А,пх), то из (4) получаем требуемые функции

ип[Х' е~^соБ(Хпх) +С' Перепишем их в виде

2иАп в1п(Хпх) „ли Заметим, что и-ая функция имеет ровно и - 1 корней внутри интервала

(0,Ь).

На рис. 1, 2 приведены графики первых двух функций распределения величин щ и и2 вдоль координаты х в различные моменты времени при = 0,1, С = 1,5.

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

Рис. 1. Распределение П\ вдоль x: 1- t = 0,3; 2- t = 0,5; 3 - t = 0,7

1- t = 0,3; 2- t = 0,5; 3 - t = 0,7

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

[1] Зайцев А.В., Пеленко Ф.В. Моделирование течения вязкой жидкости в трубе // Научный журнал СПбГУНиПТ. Серия: Процессы и аппараты пищевых производств (электронный журнал) / СПбГУНиПТ. - № 1. - Март. 2012. Режим доступа к журн.: http://processes.open-mechanics.com/journals свободный.

[2] Уизем Дж. Линейные и нелинейные волны: Пер. с англ. / Под. ред А.Б. Шабата. - М.: Мир, 1977. - 622 с.

[3] Владимиров В.С. Уравнения математической физики: — М.: Наука, 1988. - 512 с.