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

УДК 539

РАСПРЕДЕЛЕНИЕ ТЕМПЕРАТУРЫ ВНУТРИ КОЛЬЦА НЕКЛАССИЧЕСКОЙ ЖИДКОСТИ С ДВУМЯ СВОБОДНЫМИ ГРАНИЦАМИ

В. О. Бытев

Е. А. Гербер

Тюменский государственный университет (Институт математики и компьютерных наук), кафедра математического моделирования E-mail: [email protected], [email protected]

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

Ключевые слова: гидродинамика, уравнения Навье - Стокса, распределение температуры, жидкое кольцо, свободные границы.

ВВЕДЕНИЕ

Temperature Distribution Inside the Ring of Liquid with Two Free Borders in the Nonclassical Model of Hydrodynamics

V. O. Bytev

E. A. Gerber

Tyumen State University (Institute of Mathematics and Computer Sciences),

Chair of Mathematical Modeling

E-mail: [email protected], [email protected]

The article deals with the task to determine the velocity field and temperature field inside the insulated ring of incompressible viscous fluid with two free boundaries in the non-classical model of hydrodynamics. The solution to the Navier - Stokes equations and heat equation obtained by numerical methods. Analysis of the results of numerical experiments reveal the effect of nondissipative viscosity on the temperature distribution inside the ring.

Keywords: hydrodynamics, Navier - Stokes equations, temperature distribution, liquid ring, free borders.

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

В работе Л. В. Овсянникова [2] впервые определяются основы подхода к изучению задач со свободными границами. Однозначная разрешимость задачи по описанию динамики кольца в рамках классической системы уравнений Навье - Стокса была доказана в [3], там же была установлена и асимптотика его поведения при £ ^ го. Динамика свободного кольца при стремлении коэффициента кинематической вязкости к нулю (V ^ 0 ) рассматривалась в работе В. В. Пухначева [4]. В работах О. М. Лавреньевой [5, 6] изучено поведение жидкого кольца при наличии сил поверхностного натяжения.

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

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

1. ПОСТАНОВКА ЗАДАЧИ

В классической модели стоксовой жидкости тензор напряжений имеет вид Т = —р1+ 2рБ, где I — шаровой единичный тензор, р — обычная динамическая вязкость, Б — тензор скорости деформации.

© Бытев В. О. , Гербер Е. А, 2011

81

Впервые в работе [9] авторами был рассмотрен вариант обобщения модели стоксовой жидкости, в которой вместо обычной динамической вязкости используется тензор динамической вязкости в следующем виде:

V Vo\

~Vo V /

M =

где р,0 — недиссипативная составляющая общей вязкости, которая формально может иметь любой знак. Таким образом, рассматривается своеобразный вариант анизотропной вязкости. Система уравнений для описания движения кольца и потока тепла в нем с учетом вида тензора вязкости представлена в виде

du ,, . . „ р— — МАи + Vp = О, dt

div u = 0, dT

РСР— = div (fcVT) + 2(MD) : D,

(1)

где и(иг(Ь,т); иф(Ь,т)) — вектор скорости, а для уравнения теплопроводности использовано обозначение скалярного произведения тензоров. В предположении осесимметричности плоского движения свободного жидкого кольца выпишем в явном виде уравнения, связанные с переносом импульса и теплоты в полярной системе координат. После проведения несложных преобразований вид слагаемого, связанного с диссипацией энергии в уравнении теплопроводности, совпал с таковым в классической трактовке, при этом все слагаемые связанные с недиссипативной вязкостью сократились. Таким образом, систему уравнений (1) можно записать в следующей форме:

диг + _и( _

dt

% r

\rdr

dt

+

, диф ,

+ Щ-^г + Vo dr

„ (dT dT\

= к

1д_

r dr

+ d2ur\ dr2 | + Vo

dur + d2urs

rdr dr2

1 д (rur)

r dr

dT Y

V

Г dr )

+ V

= 0,

2 ( dur I dr

диф д2иф \ 1 др J p dr

rdr dr2

^ф диф _

r2 rdr dr2

+ 2 ( —

r

+

dr

(2)

(3)

(4)

(5)

Если все теплофизические величины жидкости являются константами, то система уравнений распадается на две независимые группы. Первые три уравнения не содержат слагаемых, связанных с температурой, а решение последнего уравнения можно найти только в том случае, если известно решение первых трех уравнений для поля скоростей. Конвективная составляющая температурного поля зависит от радиальной составляющей скорости, а диссипативная — от обеих составляющих скорости. Граничные условия для уравнений (2)-(4) можно определить из условия равенства нулю тензора напряжений на свободных границах при г = Я12а именно из следующих равенств:

ГТ dur f диф

Trr = -р + 2 ¡л—- + fio —— dr \ dr

(диф иф\ диг

Тгф = ц, —---- 2ц,о—— = 0.

dr r dr

Пф

r

= 0,

(6)

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

2

r

u

u

2

2

2

d

эквивалентной форме записи системы уравнений для скоростей в виде системы уравнений в переменных скорость - давление или эквивалентной ей в переменных функция тока и вихрь. Приведем для справки систему уравнений в переменных скорость - вихрь:

Г Ж

~dt dp,

+ u

d О

+

dr д 2pr

rdr dr2

= P

1Ж dHl

r dr dr2

Vo

16>0 d2tt\ J

r dr dr2 J | r

1 dur d2ur

r dr

+

dr2

+ ( ^U \ dr

связано со слагаемым v0

где О = — + —— — вспомогательная переменная для описания вихря. Основное отличие полученной г дг

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

1 дП д20\

-——Ь -г-о" в уравнении для определения давления (в классическом г дг дг2 )

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

Главной проблемой поиска численного решения системы уравнений (2)-(5) является то, что область интегрирования меняется в зависимости от времени. Для решения этой проблемы авторами [7,10] использовался прием, связанный с введением новых переменных, позволяющих осуществить отображение изменяющейся области интегрирования на фиксированную область. В качестве переменной, отслеживающей изменение положения внутренней границы кольца, взято отношение

R2(t)

r

£(£) = —■ Можно показать, что новая переменная г\ = — С + определяющая положение любой

Д2

точки жидкого кольца, будет изменяться в фиксированных пределах, а именно 0 < г] < а = — 1.

Що

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

иф =

Vu(r, t)r

R2

r20

T =

vt

r20

где V — обычная кинематическая вязкость, а Ф(£), £) — неизвестные функции, определяющие закономерности изменения радиальной и угловой составляющих скорости, то система уравнений (2)-(5) вместе с начальными и краевыми условиями перепишется в следующем виде:

¿Ф

dT

аФ(Ф - 4)

С(а + С) ln(1+ С/а) 2Фи

+

1

ln(1+ С/а)

a du 2 4е—— +uj2 dr/

dn, Ф(0) = Фо,

dT

= 2Ф, С(0) = 1,

du

57 +

_ gduj ч д2ио

(С + v) дг) дг]2'

/ дсо е

\дг] (С + г])2

Ф

= 0, u(n, 0)= uo(n), (7)

П=0,а

Ф2

. d © „ d , d ©,

+ (С + n)2

duj_\ dn

2

дв dn

= 0, ©(n, 0) = ©o(n),

П=0,а

где использована безразмерная величина е = v0/v, равная отношению соответствующих вязкостей, A, B — комбинация теплофизических характеристик жидкости, © = T(t,n)/Tbaz — относительная температура.

Для проведения экспериментов по численному моделированию использовалась программа «Ring v 1.1» [10], написанная на языке Delphi на основе конечно разностной аппроксимации системы уравнений (7). Тестирование программы осуществлялось путем сравнения численного решения с известным аналитическим решением при начальном условии Ф0 = 4 для первого уравнения этой системы. Анализ результатов численного моделирования показал достаточно хорошее согласие между этими решениями. Использование этой программы позволило исследовать характер влияния недис-сипативной вязкости на динамику жидкого кольца в зависимости от ее величины и знака отношения

= V

r

a

вязкостей ^ путем варьирования ее значений. Для осуществления сквозного счета программа «Ring v1.1», была доработана, и к блоку численного решения уравнений движения и неразрывности был добавлен блок для численного расчета температурного поля.

2. ИТОГИ РАСЧЕТА

На первом этапе численного моделирования распределения температуры во вращающемся жидком кольце проводились для воды, т. е. использовались ее теплофизические характеристики. Значения параметра г выбиралось из интервала [0.001 ^ 40]. В нулевой момент времени относительная температура кольца принималась одинаковой, что соответствует значению для относительной температуры 0(0, п) = 1. Этот вариант начального распределения был выбран для того, чтобы процесс диссипации энергии движения жидкого кольца был доминирующим при условии, что границы кольца теплоизолированы. Несмотря на то что поле скоростей менялось достаточно существенно, никаких эффектов влияния недиссипативной составляющей вязкости обнаружено не было. Результаты численного эксперимента с большой точностью выдавали 0(т, п) = 1, для любого набора значений параметра г, т. е. первоначальное распределение температуры в жидком вращающемся кольце практически оставалось неизменным.

Для объяснения подобного поведения температурного поля была проведена оценка величины дис-сипативной части энергии, которая перейдет в тепло для реального вещества, например в нашем случае для воды. Характерные значения рассчитанных составляющих скорости равны соответственно: иг = Ъ • 10-6^ и иф = 0.5 • 10-6^. Абсолютное значение модуля скорости \и\ рз 5 • 10-6^. Если вся кинетическая энергия перейдет в тепло, то справедливо равенство

и2

шу = тСрАТ. (8)

Подставляя соответствующие значения в (8), получим оценку для разности температур АТ =

25 • 10

-12

2Ср 2 • 2.6 • 103

5 • 10 15 град. Величина абсолютно неуловимая на фоне ошибок, получаю-

щихся в ходе процесса численного решения системы (7). Таким образом, показано, что для ньютоновской жидкости, которой является обычная вода, в зоне ламинарного течения (число Рейнольдса 5 • 10-6 • 0.5

Не = —-—— ^ 3)), влияние недиссипативнои вязкости ничтожно мало.

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

Результаты численного моделирования в этом случае показывают заметное влияние величины недиссипативной вязкости на поведение температурного поля. Иллюстрация поведения относительной температуры представлена на рис. 1, где показан характер ее изменения в радиальном направлении в разные моменты относительного времени. Поле скоростей было построено для начального радиального распределения скорости при Ф0 = 4, что соответствует тому, что радиальная составляющая скорости направлена от центра кольца в начальный момент времени.

© т 1200 1000 800 600

2

0 5 10 15

♦ Зависимость температуры от параметра s;

20 25

Полиномиальный ряд

Рис. 1. Изменение распределения температуры вдоль радиуса в разные моменты

времени

2

u

0

Относительная температура кольца меняется достаточно динамично. Вначале идет разогрев наружной части кольца, а затем и внутренней. Подобное поведение можно объяснить разностью абсолютных величин скоростей на границах кольца. По истечении некоторого времени температура кольца достигает своего максимального значения и становится практически постоянной по всей площади кольца, так как кольцо практически останавливается. Величина максимальной относительной температуры ©max = const напрямую определяется лишь абсолютным значением е, т. е. величиной недиссипативной вязкости, так как по условию в нашем случае, v = 1.

Используем данные по составляющим скорости для этого случая, например при е = 10 получились следующие средние по радиусу составляющие скорости: радиальная — ur_sr = 5^, угловая — иф_зг =80^. Оценка числа Рейнольдса дает величину Re ~ 80, что позволяет сделать вывод о том, что движение жидкости для этого случая все еще осуществляется в ламинарном режиме.

Так как температура кольца стремится к насыщению, то становится интересным характер зависимости ©шах(е) — максимального значения относительной температуры от абсолютной величины отношения вязкостей. На рис. 2 изображена зависимость ©шах(е) при прочих фиксированных характеристиках. Набор точек, полученных в результате численного моделирования для величины нагрева кольца А©тах(е) = ©шах(е) — 1 достаточно хорошо описывается зависимостью

А©тах(е) = 2е2. (9)

Рис. 2. Зависимость ©max(s)

Объяснение подобного поведения, скорее всего, возможно на основании соотношения (8). Действительно, при всех имеющихся допущениях прирост относительной температуры за счет диссипации будет пропорционален квадрату средней скорости вращения. Поскольку начальное распределение угловой составляющей скорости должно быть согласовано с граничным условием для третьего уравнения системы (7), то в качестве начального распределения угловой скорости, соответствующего этому условию, в данной программе использовалось выражение <(0, п) = еФо/(1 + п). Очевидно, что в данном случае начальное распределение угловой части скорости, являющейся основной составляющей полной скорости элемента жидкости (рис. 3), прямо пропорционально отношению вязкостей, что и объясняет квадратичную зависимость прироста относительной температуры. Таким образом, можно считать, что нет прямого влияния недиссипативной составляющей вязкости на распределение относительной температуры. Что касается характера влияния недиссипативной вязкости на поведение составляющих полной скорости, то следует отметить, что с ростом ее величины заметным образом изменяется характер поведения радиальной составляющей скорости движения жидкого кольца.

,•; • 1

s — ±0. Ь

1

- У ..—,— £ = ±0, 2

1 /

If

щ

1

а б

Рис. 3. Изменение вектора скорости скорости при различных значениях е и Ф0 = -4: а — изменение радиальной составляющей скорости, б — изменение угловой составляющей скорости

ЗАКЛЮЧЕНИЕ

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

Библиографический список

1. Андреев В. К., Бублик В. В., Бытев В. О. Симметрии неклассических моделей гидродинамики. Новосибирск, 2003. 352 с.

2. Овсянников Л. В. Общие уравнения и примеры // Задача о неустановившемся движении жидкости со свободной границей: сб. работ теор. отдела ИГ СО АН СССР. Новосибирск, 1967. С. 5-75.

3. Бытев В. О. Неустановившиеся движения кольца вязкой несжимаемой жидкости со свободными границами // ПМТФ. 1970. № 3. С. 88-98.

4. Пухначев В. В. Неклассические задачи теории пограничного слоя. Новосибирск, 1980.

5. .Лаврентьева О.М. Неустановившееся движение вращающегося кольца вязкой капиллярной жидкости // Динамика сплошной среды: сб. науч. тр. ИГ СО АН СССР. Новосибирск, 1978. Вып. 31. С. 52-60.

6. Лаврентьева О. М. Предельные режимы движения вращающегося вязкого кольца // Динамика сплошной

УДК 539.374

Д. В. Гоцев

Военный авиационный инженерный университет, Воронеж, кафедра математики E-mail: [email protected]

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

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

среды: сб. науч. тр. ИГ СО АН СССР. Новосибирск, 1980. Вып. 44. С. 15-34.

7. Бытев В. О., Гербер Е.А. Об одной задаче с двумя свободными границами // Современные проблемы математики и её прикладные аспекты: материалы Всерос. науч.-практ. конф. Пермь, 2010. С. 100.

8. Бытев В. О., Гербер Е.А. Численное моделирование динамики жидкого кольца // Наукоемкие информационные технологии: тр. XIV молодежной науч.-практ. конф. Переславль-Залесский, 2010. С. 109-114.

9. Bytev V. O. Building of Mathematical Models of continuum media on the basis of invariante principle // Acta Appl. Math. 1989. Vol. 16. P. 117-142

10. Бытев В. О., Гербер Е. А. О восстановлении точного решения и о распространении температуры в кольце жидкости // Математическое и информационное моделирование: сб. науч. тр. Тюмень, 2009. Вып. 11.

Stability of the Cylindrical Cover with the Elastic-Viscous-Plastic Filler at Axial Compression

D. V. Gotsev

Military Aviation Engineering University, Voronezh, Chair of Mathematics E-mail: [email protected]

Within the limits of the exact three-dimensional equations of stability of an equilibrium state of a cylindrical cover with a filler is investigated at axial compression. Calculations were spent for a case when the cover material was modelled by an elastic body, and a filler material - environment with difficult rheological properties - elastic-viscous-plastic. The estimation of influence on size of critical pressure of parameters of a cover and a filler is given.

Key words: exact three-dimensional equations, stability of balance, cylindrical cover, filler, difficult rheology, axial compression.

УСТОЙЧИВОСТЬ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ С УПРУГОВЯЗКОПЛАСТИЧЕСКИМ ЗАПОЛНИТЕЛЕМ ПРИ ОСЕВОМ НАГРУЖЕНИИ

( Гоцев Д. В., 2G11