close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
ФИО соискателя: Ряжских Александр Викторович Шифр научной специальности: 05.13.18 - математическое моделирование, численные методы и комплексы программ Шифр диссертационного совета: Д 212.038.20 Название организации: Воронежский государственный унив
На правах рукописи
РЯЖСКИХ АЛЕКСАНДР ВИКТОРОВИЧ
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЛАМИНАРНОГО
ИЗОТЕРМИЧЕСКОГО ТЕЧЕНИЯ СТЕПЕННОЙ ЖИДКОСТИ НА
НАЧАЛЬНОМ УЧАСТКЕ ОСЕСИММЕТРИЧНЫХ
ГОРИЗОНТАЛЬНЫХ КАНАЛОВ
Специальность: 05.13.18 – Математическое моделирование,
численные методы и комплексы программ
Автореферат
диссертации на соискание ученой степени
кандидата физико-математических наук
Воронеж – 2012
Работа выполнена на кафедре технической механики в ФГБОУ ВПО «Воронежский
государственный университет инженерных технологий»
Научный руководитель:
Чертов Евгений Дмитриевич,
доктор технических наук, профессор
Официальные оппоненты: Артемов Михаил Анатольевич,
доктор физико-математических наук, профессор,
Воронежский государственный университет,
заведующий кафедрой программного обеспечения и
администрирования информационных систем
Семенов Михаил Евгеньевич,
доктор физико-математических наук, профессор,
Военный учебно-научный центр Военно-воздушных
сил “Военно-Воздушная Академия им. проф. Н.Е. Жуковского и Ю.А.Гагарина”, профессор кафедры №11
Ведущая организация:
Тамбовский государственный технический
университет
Защита состоится 28 ноября 2012 г. в 15 часов 10 минут на заседании диссертационного совета Д 212.038.20 в ФГБОУ ВПО «Воронежский государственный университет» по адресу: 394006, г. Воронеж, Университетская площадь, 1, ауд. 335.
Отзывы на автореферат ( в 2-х экземплярах), заверенные гербовой печатью учреждения, просим направлять по адресу: 394006, г. Воронеж, Университетская площадь, 1, ФГБОУ ВПО ВГУ, ученому секретарю диссертационного совета Д 212.038.20.
С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Воронежский государственный университет».
Автореферат разослан 24 октября 2012 г.
Ученый секретарь диссертационного совета,
кандидат физико-математических наук, доцент
2
С.А. Шабров
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность. Установление закономерностей явлений переноса при ламинарном изотермическом течении реологически сложных сред на начальном участке
горизонтальных прямолинейных каналов в различных компактных системах (тепломасообменные аппараты, экструдерные устройства, энергетические установки и
т.д.) химической и пищевой промышленности, а также в ракетно-космической и
криогенной технике имеет определяющее значение в достижимости расчетных характеристик при проектировании нового и модернизации существующего оборудования. В контексте этой проблемы идентификация длины гидродинамического начального участка в зависимости от реологических свойств жидкостей является одной из фундаментальных задач, результаты решения которой существенным образом влияют на точность определения тепломасообменных параметров.
Наибольший вклад в решение этой задачи внесли Boussinesq J., Langhaar H. L.,
Schiller L., Лейбензон А.С., Тарг С. М., Слезкин В.А., Жукаускас А., Matras Z., Gupta
R.C., Shlichting H. и др., которые использовали погранслойные представления и доказали прямопропорциональную зависимость между длиной гидродинамического
начального участка и числом Рейнольдса, причем основные результаты касались в
основном ньютоновских жидкостей.
Среди широкого спектра неньютоновских жидкостей особое значение имеет
модель степенной жидкости (жидкость Оствальда – де Виля), которой подчиняется
поведение большинства псевдопластиков и дилатантов, используемых в различных
прикладных областях.
Однако экспериментально-теоретический анализ задачи о начальном участке в
каналах при течении степенной жидкости затруднен из-за объективных трудностей
непосредственного измерения поля скоростей и давления, а также ввиду сложной
нелинейности уравнений гидродинамики, не позволяющей корректно применять современный аппарат вычислительной математики для их интегрирования. Поэтому
метод математического моделирования остается пока одним из эффективных
средств решения этой задачи.
Диссертационная работа выполнялась в соответствии с планом научноисследовательских работ Воронежского государственного университета инженерных технологий по теме «Создание инновационных технологических процессов и
оборудования для специализированных машиностроительных производств» (№ГР
01.2.011.14702)
Цель работы. Синтез и анализ математических моделей для идентификации
длины гидродинамического начального участка при изотермическом ламинарном
течении степенной жидкости в горизонтальных осесимметричных каналах.
Для достижения цели поставлены задачи:
1) на основе уравнений Навье-Стокса, записанных для степенного реологического
закона, разработать математические модели однонаправленного изотермического
ламинарного течения в плоском, круглого и кольцевого сечений прямолинейных горизонтальных каналах с проверкой их адекватности для случая ньютоновской жидкости;
2) с использованием представлений о ползущем течении разработать математические модели однонаправленного течения высоковязких ньютоновских жидкостей в
3
плоском, круглого, прямоугольного, а также кольцевого сечения каналах и получить
их аналитические решения;
3) разработать маршевый алгоритм, базирующийся на кончено-разностной дискретизации уравнений моделей, для их численного интегрирования и провести вычислительные эксперименты;
4) на примере плоского канала предложить класс приближенных аналитических решений, на основе которых идентифицировать длину гидродинамического начального участка для степенной жидкости.
При выполнении исследования в качестве основного инструментария был
применен метод математического моделирования совместно с методами теоретической гидродинамики, положениями теории дифференциальных уравнений математической физики, вычислительной математики и реологии. Достоверность и обоснованность полученных результатов основывается на использовании законов явлений
переноса, на проведении вычислительных экспериментов и сравнительном анализе с
классическими данными.
Научная новизна диссертации заключается в следующем:
1) предложена математическая модель, которая позволила свести задачу идентификации длины гидродинамического начального участка к решению начально-краевой
задачи с последующим масштабным переходом от времени к продольной координате на основании коэффициента трения установившегося течения по Фаннингу; в
случае ньютоновской жидкости получены приближенные аналитические решения
для плоского, круглого и кольцевого сечений каналов;
2) для численного интегрирования уравнений моделей, когда индекс течения степенной жидкости отличен от единицы, разработана маршевая по времени конечноразностная схема, отличающаяся модификацией начального условия к квазиоднородному виду;
3) обосновано применение идеализации ползущего течения к решению задачи об
идентификации гидродинамического начального участка высоковязких жидкостей,
формулировка которой сведена к интегрированию эллиптических уравнений Пуассоновского типа; получены аналитические решения для каналов плоского, круглого
и прямоугольного, а также кольцевого сечений;
4) на примере плоского канала показано, что определение длины начального участка
при течении степенной жидкости может быть сведено к установлению длины разгонного течения; получено приближенное аналитические решение с использованием
структуры решения для ньютоновской жидкости.
Практическая значимость состоит в разработке инструментария в виде
предметно-ориентированного программного комплекса, который позволяет оценивать длину гидродинамического начального участка при течении степенной жидкости в осесимметрчиных горизонтальных каналах в зависимости от индекса течения,
что создает базу для уточнения тепломассобменных характеристик и рационального
выбора режимных и конструкционных параметров при функционировании различных технических и технологических систем.
Апробация. Основные результаты диссертационного исследования доложены
и обсуждены на международных научных конференциях «Современные проблемы
прикладной математики и математического моделирования» (Воронеж, 2005);
4
«Авиакосмические технологии - VI» (Воронеж, 2005); на VII Всероссийском симпозиуме по прикладной и промышленной математике – весенняя сессия» (Кисловодск,
2006); «Математические методы в технике и технологиях – ММТТ-19» (Воронеж,
2006); «Современные проблемы прикладной математики, теории управления и математического моделирования – IV» (Воронеж, 2011); «Физико-математическое моделирование систем - VII» (Воронеж, 2011); на зимней математической школе «Современные методы теории функций и смежные проблемы» (Воронеж, 2011).
Публикации. Основное содержание диссертации изложено в 11 работах, из
них 4 в реферируемых журналах из списка ВАК РФ.
Структура и объем. Диссертация состоит из введения, четырех глав, основных выводов, списка литературы и приложения. Материал изложен на 150 страницах и содержит 50 рисунков и 5 таблиц.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертации, сформулированы
цель и задачи исследования, научная новизна и ее практическая значимость.
В первой главе проанализировано современное состояние проблемы математического моделирования течения реологически различных жидкостей на начальном
участке каналов.
Отмечено, что существует два подхода к понятию длины гидродинамического
начального участка: расстояние, отсчитываемое от входа до сечения, соответствующего слиянию пограничных слоев; сечение, при котором максимальное значение
скорости отличается с наперед заданной точностью от его установившегося значения.
Как правило, в основе построения математических моделей гидродинамики
при ламинарном изотермическом течении жидкости лежат уравнения, вытекающие
из законов сохранения вещества и количества движения:
Dv
D
p g ;
v ;
2 p ,
(1)
Dt
Dt
где – плотность жидкости; v – вектор скорости; p – давление; g – вектор ускорения силы тяжести; – тензор напряжений; D Dt t v – оператор субстанциональной производной по скалярному полю, t (время); – градиент скалярного поля.
Несмотря на различные подходы в определении начального участка, (1) трансформируется в приближенные уравнения пограничного слоя, в которых, как правило, инерционные члены учитываются через среднюю скорость по сечению. Решение
задачи в такой постановке дало возможность показать, что относительная длина начального участка прямо пропорциональна числу Рейнольдса с коэффициентом пропорциональности, зависящим только от геометрии проходного сечения канала. Тем
не менее, основные результаты полученные этим путем, относятся в основном к
ньютоновским жидкостям, причем известные данные различных авторов отличаются друг от друга существенным образом.
Многообразие спектра неньютоновских жидкостей привело к тому, что исследования разделились по типам жидкостей. В связи с этим в данной диссертации ог5
раничились анализом течения неньютоновских жидкостей, для которых напряжение
сдвига является степенной функцией от скорости сдвига
n ,
(2)
где – мера консистенции жидкости; n – индекс течения ( 0 n 1 – псевдопластики ; n 1 – ньютоновская жидкость; n 1 – дилатантная среда) . Но даже в этом
случае математическая формулировка задачи столь сложна в силу своей существенной нелинейности, что получение численного, а тем более аналитического ее решения в обозримом будущем пока проблематично.
Тем не менее в последнее время вновь с использованием погранслойного приближения, в научной литературе стал разрабатываться «трансформационный» метод
для получения оценочных решений, в основу которого положена аналогия с характеристиками ньютоновской жидкости. Одновременно с этим предпринимались попытки численного интегрирования задачи в общей постановке, но полученные различными авторами результаты практически мало коррелируются.
Сделан вывод о необходимости продолжения исследований в направлении
поиска новых решений, основанных на синтезе математических моделей, полученных из физически обоснованных допущений.
Вторая глава посвящена синтезу математических моделей гидродинамики
степенной и высоковязкой жидкостей на начальном участке каналов.
Вначале система (1), (2) в декартовой системе координат представлена для
плоского течения, когда остаются только поперечная x и продольная z координаты
и соответственно компоненты вектора скорости vx и vz при отсутствии действия
объемных сил.
Так как анализируется ламинарный режим течения несжимаемой жидкости, то
сделано допущение об однонаправленности движения на всем протяжении канала,
тогда vx 0 и из (1) следует, что vz z 0 , p x 0 , p z const . Т.о., система
(1) – (2) трансформируется в следующий вид
n 1 v z
p
vz 1 vz
(3)
t
z
x x 0 x
v 0, t с граничными условиями
vz x,0 v0 , vz h, t z
0
(4)
x
1h
и требованием выполнения балансового соотношения
(5)
vz x, t dx v0 ,
h0
где v0 const – скорость на входе в плоский канал шириной 2h ; 0 – скорость
сдвига в приведенном состоянии, обычно равная 1. По аналогии с (3) – (5) получены
уравнения в цилиндрической системе координат
n 1 vz
p vz 1 vz
,
r
(6)
t
z r r r 0 r
с граничными условиями для круглой трубы
6
vz r ,0 v0 ,
vz r0 , t vz 0, t 0,
r
(7)
где r0 – радиус трубы; и для кольцевого канала
vz r ,0 v0 ,
vz r1, t vz r2 , t 0 ,
(8)
где r1 , r2 – радиусы внутренней и внешней аксиальных труб; а также балансовые
соотношения
2 r2
2 r0
rv
r
,
t
dr
v
,
(9)
rvz r , t dr v0 .
z 0
r22 r12 r1
r02 0
Отметим, что (5) и (9) интерпретируются как критерии точности моделей (3), (4) и
(6) – (8). Если tl время установления стационарного профиля, то длина начального
участка zl может быть вычислена по соотношению
zl v0tl .
(10)
Далее рассмотрен случай стационарного однонаправленного движения высоковязкой жидкости ( достаточно большое, n близко к 1) в рамках представлений о
ползущем течении, т.е. когда силы вязкости существенно преобладают над силами
инерции, что позволяет пренебречь конвективной составляющей в операторе субстанциональной производной от скорости, поэтому (1) , (2), в предположении отсутствия градиента давления в поперечном направлении и его постоянства вдоль по потоку принимает вид
1 p
,
(11)
2v z z
где 2 – оператор Лапласа. Уравнение (11) дополняется условиями «прилипания»
на стенках канала, постоянства скорости на входе и отсутствием изменения скорости на .
В этой же главе приведена запись уравнений (1) , (2) в переменных Гельмгольца. Показано, что для однонаправленного течения преимущества формулировки
задачи в переменных Гельмгольца нивелируются по сравнению с задачей в действительный переменных.
В главе 3 рассмотрен анализ математических моделей (3), (4) и (6) – (8), который начат с получения аналитических выражений скоростей на основе одностороннего преобразования Лапласа по t для частного случая n 1 :
для плоского канала
m
1
4 1 2
V X ,T 1 X 2 3
4 2 2m 1 Br 3
2 Br
Br m 0 2m 1
2
2
cos 2m 1 X exp 2m 1 BrT ;
2
4
для круглой трубы
1 R 2 2 1 BrV0 n2 J 0 n R V R, T exp Br n2T ;
3
4 Br
Br n 1
J1 n n
7
(12)
(13)
для кольцевого канала
1 2 2
2 R 2 ln R
V R ,T R
R
R
1
2
1
4 Br R1
2 1 n2
Br n 1
J 0 R2n J 0
n3 R1J1 R1n
R J R N
R ln 2 R1 N R N R J R R N R exp BrT 0
1 n
1 n
0
0
n
2 n
0
n
2
n
N0 R2n R2 J 0 R1n N1 R2n 1 2 n 0 R1n R1J 0 R2 n N1 R1n ,
2
(14)
где X x h , R r r0 , R r r2 r1 , d hVr0 r2 r1 , T t t , V vz v , V0 v0 v ,
V , R 2 R 2 v 0 d , t v p z , Br 0 d p z , Br
0
2
1
R 22 R12 ln R 2 R1 8 , m , n – корни уравнений J 0 R1n N 0 R 2n J 0 R2n N 0 R1n 0 и J 0 n 0 ; J 0,1 , N 0,1 – функции Бесселя 1-го и 2-го родов
нулевого и первого порядков.
Характер
ZL
изменения осеV
V
Re
вой и средней
по сечению ско1
0,5 ростей течения 0,30
1,25
на
примере
плоского канала
0
1,0
(рис.1) показал, 0,20
0
1
3 T
2
что в области,
Рис. 1. К проверке адекватно- где происходит
сти математической модели: 1
0,10
установление
– осевая скорость; 2 – средняя профиля скоро0
0,25
R1 R2
скорость
сти, балансовое
Рис. 2. Относительная длина
соотношение (5) выполняется практически точно. Это
начального участка в кольцепозволило сделать вывод о принципиальной возможвом канале: — - расчет по моности использования синтезированных моделей для
дели; ○ – данные Heaton H.S.,
идентификации длины начального участка. Для конReynolds W.C., Kays W.M
кретизации (10) использовано определение коэффици(1964).
ента трения по Фаннингу, из которого найдено,
что в показателе экспоненты (12) вместо BrT необходимо подставить 2Z Re Br ,
где Z z h . Это позволило найти длину начального участка из условия отклонения
величины осевой скорости от ее стабилизированного значения с заданной наперед
точностью (обычно принимают 1%), причем показано, что это можно сделать с использованием лишь одного члена ряда (12). Найденные значения Z L Re для плоского канала и круглой трубы 0,196 и 0,156 коррелируют с классическими данными.
2
8
Для кольцевого канала результаты приведены на рис. 2. Различие в области
R R 0 (круглая труба) объяснено некорректностью применения авторами ли1
2
неаризации Лангхаара уравнений (1), (2) для кольцевого канала.
В общем случае, когда n 1 , для численного интегрирования уравнений моделей, записанных в безразмерном виде, предложены маршевые конечно-разностные
схемы. Непрерывная область, в которой искались решения, заменялась дискретной в
соответствии с центральным трехточечным шаблоном по геометрической координате и двухточечным по временной координате. Это позволило аппроксимировать
уравнения (3) и (6) дискретными аналогами:
1 n 1
2
1
1
1 n 2 2 Vi k 1 Vi k 1 Br ik
ik
fGi1 ik
n 1 ik
ik
ik T , (15)
где
ik Vi k1 Vik1
1
ik Vi k1 2Vik Vi k1
2
2 G ,
G X R R , X R 1 N , R R2 R1
Vi k V Gi , k T ,
G 2 ,
N , Gi iX iR R1 iR ; f 0 или f 1 соответственно для декартовой и цилиндрической сис-
тем координат; для плоского канала и круглой трубы 1 , для кольцевого канала
1 при R1 iR R1, R и 1 при R1 iR R , R 2 , где R удовлетворяет
интегральному уравнению
R R2
2
1
n
R2 2
R
R d
R
R
R
R1 R Дискретные аналоги граничных условий таковы:
для плоского канала
Vi0 1 2 X i V0 ,
R2
1
n
.
d R
(16)
3;
(17)
3;
(18)
VNk 0 ,
V0k 4V1k V2k
VNk 0 ,
V0k 4V1k V2
для круглой трубы
Vi0 2 3Ri V0 ,
для кольцевого канала
V0k
где
Ri , при R1 Ri R ,
0;
,
(19)
Ri , при R Ri R2 ,
0 0 V1 Ri R R1 V0 0 ;
V2 Ri R2 R V0 0 ;
R R
R R
V1
Vi V2
VNk
0
1
2
1
1
1
1 R22 R12 R3 R1R2 R13 R R1 2
6 3
2
1
1 1
R23 R 2 R2 R3 R2 R 1 0 ,
2
3 6
0 – точность аппроксимации начальных условий.
9
n
Безразмерный параметр Br равен: для плоского канала Br n V0 2n 1 ;
n
для круглой трубы Br 0,5 n V0 3n 1 ; для кольцевого канала
1
Br 2
R 2
R R R R22 R12 V0 R1 R1 V
а
5
1
2 n
R
R R R R .
dRdR
dRdR
R R R1 R1 V
4
3
2
5
4
3
2
R1
0
4
3
2
1
1
R2
в
1
1
1
0
V
б
5
1
n
1
2 n
R
0
R2
R1
R2
R1
Рис. 3. Профили скорости в кольцевом канале: а – n 0,75 ; б – n 1,0 ; в –
n 1,25 ; в моменты времени T : 1 0 ; 2 0,1 ; 3 0,4 ; 4 1 ; 5 2 .
Найдены условия устойчиZ L Re
4
вости и сходимости схемы
(15) – (19) ( T 0,01G и
3
0,3
0,2 n 3 ), которые оказа2
лись инварианты к геомет1
рии канала. Тестирование
0,2
схемы при T показало
ее работоспособность при
n 1 . Установлено, что коор0,1
дината максимальной скорости R с увеличением индекса
течения n смещается к
n
0
1
Рис. 4. Длина начального участка в кольцевом ка- R1 R2 2 , а Br 0 . Кроме
нале при R1 R2 : 1 0,5 (○ – Hartnett J.P. (1989), ■ – того, точность алгоритма поGupta R.C. (1990)); 2 0,25 ; 3 0,1; 4 0,05 (▲ – зволяет регистрировать разMehrotra A. K. (1990), ● – Pooley, Redley (2007), □ – личие в физической картине
течения при достаточно
Matras, Nowak (1982)) — - расчет по модели.
близких реологических па-
10
раметрах (рис. 3). Подтвержден факт, как и в случае других геометрий, что длина
начального гидродинамического участка с уменьшением n увеличивается (рис. 4),
причем максимальное ее увеличение наблюдается в области n 0, 2;0,4 .
Анализ математической модели (11) приведен в главе 4, в которой с помощью
интегральных синус-преобразования и преобразования Ханкеля получены аналитические выражения, определяющие поле скоростей на начальном гидродинамическом участке:
для плоского канала
1
1
1
1 U Z , X 2K (20)
cos m 1 2 exp m Z 2 sin m X ,
m m 1 m
m K где m m ; X x h , Z z h , V vz h , Re v0h ; K P Z Re ;
P p
v ; h - ширина канала; 2
0
– коэффициент кинематической вязкости жид-
кости; v0 const – скорость жидкости на входе канал;
для круглой трубы
J q R K
K U R, Z 2 0 m 2 1 2 exp qm Z ,
(21)
q
q
m 1 qm J1 qm m m
где qm – корни уравнения J 0 qm 0 , R r r0 , Z z r0 , V v z r0 , Re v0r0 ;
для кольцевого канала
2 2
2 pm J 0 pm R 2 pm U R, Z 2 m 1 J 2 p R J 2 p R
0
m 1
0
m 2
J 0 pm R N 0 pm R1 J 0 pm R1 N 0 pm R 1 K K 1 2 exp pm Z 3 ,
pm pm pm pm R2 J1 pm R2 R1J1 pm R1 N 0 pm R1 R2 N1 pm R 2 R1 N1 pm R1 J 0 pm R1 ;
где
(22)
pm – корни уравнения J 0 pm R2 N 0 pm R1 J 0 pm R1 N 0 pm R2 0 .
R r d ; R1 r1 d ; R2 r2 d ; d r2 r1 ; V v z d ; Re v0d ; r1 , r2 – радиусы
внутренней и внешней труб;
для канала прямоугольного сечения
1
U X ,Y , Z 4 cos n 1 cos m 1 n 1 m 1 n m
2 2 n2 m2 n2 m2 n
m
1 K 2 2 exp 2 2 Z K 2 2 B B A
B A
A
11
Y
X sin n sin m ,
(23)
A
B
где n n ; m m ; X x h , Y y h , Z z h , h h1h2 h1 h2 ; V vz h ;
Re v0h ; A h1 h2 h2 ; B h1 h2 h1 . Коэффициент K , найденный из условия, что при Z безразмерная средняя скорость по сечению есть 1, равен для
плоского канала 12 , для трубы 8 , для каналов кольцевого и прямоугольного сечений соответственно
1
J 02 pm R2 2 pm 2
K 2
,
2
2
2
2
R m 1 p J p R J p R R
1
2
m 0
m 1
0
m 2 (24)
1
1
n2 m2 2
2
K 4 2 2 cos n 1 cos m 1 2 2 .
(25)
A
B
n 1 m 1 n m
Получено, что результаты расчетов длины начальных участков по (20) – (25) коррелируют с результатами главы 3 при n 1 .
В этой же главе проведена приближенная аналитическая оценка длины начального участка в плоском горизонтальном канале при течении степенной жидкости в рамках модельных представлениях о слабодеформируемых течениях в приложении к ламинарному движению с использованием структуры уравнения (11), в результате чего получена краевая задача
n 1
n 1
U
U U
U (26)
K* ,
X X
X Z X
Z U X , U X ,0 1 ,
0,
(27)
Z
U 1, Z U 0, Z 0 ,
(28)
где X x h , Z z h , V vz 0 h ; Re v0 h ; K Re 1 dP dZ ; P p
v ;
2
0
K * Re1 n K ; U V Re .
Для n 1 показано, что длина начального участка может быть определена на
основе решения задачи о разгонном течении. Это позволило перейти к решению (26)
– (28) в формулировке разгонного течения, при этом путем выбора вспомогательной
функции W Z , X по правилу
U
X
сформулирована система
n 1
U W
,
X X
2W X , Z 2
U
X
n 1
2W X , Z 2
U W
Z Z
K* ,
X
Z
W X , W X ,0 W 0, Z W 1, Z 0 ,
Z
12
решение которой
W X , Z 2K * cos m 1 exp
m Z 1 sin m X ,
m3
где m m . Поэтому для разгонного течения поле скоростей
m 1
1
X
W
n
dX ,
при 0 X .5
0 X U Z , X ;
1
1
X
0,5 W n
W n
dX
dX , при .5 X 1
0 X 0,5 X (29)
1
1
0,5 X
0,5
W X , n
1 W X , n
K* dXdX
X
dX X
2
0
0
0
1
n
4Z L
hRe
n
1 X W dXdX .
(30)
X
0,5 0,5 Расхождение с известными данными в области «сильной» дилатантности объяснено выбором W , что накладывает ограничения на условия U X U Z и
W X W Z , которые с увеличением K * выполняются менее точно, о чем свидетельствует ве0,10
личина среднеинтегральной ошибки
11
W W
dXdZ .
0
1
n
Z
0 0 X
Рис. 5. Длина начального гид- Длина начального участка (рис. 5) рассчитывалась
родинамического участка в по модифицированному соотношению
плоском канале: — расчет по
(31)
U .5, Z 1 1 1 U .5, U .5, Z .
модели; □ – данные Gupta
В этой же главе показано, что длина гидродиR.C. (1989)
намического начального участка при разгонном течении эквивалентна длине инерционного течения.
Описание предметно-ориентированного программного комплекса (рис. 6), позволяющего вычислять длину начального участка, приведено в приложении.
Программный комплекс на алгоритмическом языке C++ построен по блочному типу для расчета длины начального участка при течении степенной жидкости в
плоском, круглом и кольцевом сечений каналов.
0,20
13
Begin
Reg
pipe
flat
circular
FL PI
n, T , X , N , K
FL CI
Br , Vi0 , V0k , VNk
k0
i0
Vik 1 Vik FLik
iN
i i 1
Vik Vik 1
k k 1
kK
INT
Scale
L
end
Рис. 6. Обобщенная блок-схема функционирования предметно-ориентированного
программного комплекса “ENTRY”.
14
ОСНОВНЫЕ ВЫВОДЫ
1. Показано, что задача идентификации длины гидродинамического участка
при ламинарном изотермическом течении степенной жидкости может быть сформулирована в рамках допущения об однонаправленности движения либо с помощью
понятия «ползущего» режима без использования погранслойных представлений.
2. Разработана маршевая конечно-разностная схема, позволяющая определять
динамику изменения профиля скорости в зависимости от индекса течения степенной
жидкости и геометрии канала, и предложен масштабный переход к длине гидродинамического участка, основанный на определении коэффициента трения по Фаннингу. Указан диапазон ее устойчивости и сходимости: шаг по безразмерному времени должен быть меньше на два порядка, чем шаг по координате, а индекс течения
изменяется в диапазоне от 0,2 до 3.
3. Полученные аналитические решения уравнений моделей для ньютоновской
жидкости, а также результаты вычислительных экспериментов подтвердили правомочность предложенного подхода для определения длины начального участка при
сравнительном анализе с известными данными.
4. Подтверждено, что длина начального участка с увеличением индекса течения уменьшается. Найдено, что его максимальная его величина находится по индексу течения в промежутке от 0,2 до 0,4 .
5. На примере плоского канала разработан алгоритм аналитического оценивания длины начального участка с использованием структуры решения для разгонного
течения степенной жидкости, который может быть применен для других геометрий.
6. Разработанный предметно-ориентированный программный комплекс позволяет автоматизировать процедуру оценивания длины гидродинамического начального участка в зависимости от индекса течения и геометрии канала.
Основные результаты диссертации опубликованы в следующих работах:
публикации в изданиях, рекомендованных ВАК РФ
1. Ряжских А.В. Гидродинамический начальный участок при течении высоковязкой
ньютоновской жидкости в круглой трубе // Вестник Санкт-Петербургского университета. Сер.10. – 2012. – вып. 3. – с. 98-102.
2. Ряжских А.В. Гидродинамический начальный участок при течении высоковязкой
ньютоновской жидкости в канале прямоугольного сечения // Вестник ВГТУ. – 2011.
– Т.7. – №11.3. – с. 39-42.
3. Ряжских А.В. Начальный участок в кольцевом канале при напорном течении ньютоновской среды // Обозрение прикладной и промышленной математики. – 2006. –
Т.13. – №2. – с. 350-351.
4. Ряжских А.В. Инерционное ламинарное течение в плоском канале // Вестник
ВГТУ. – 2006. – Т.2. – №6. – с.102-104.
статьи и материалы конференций
5. Чертов Е.Д., Ряжских А.В. Решение начально-краевой задачи для неоднородного
эллиптического уравнения в полуограниченной полосе / Современные методы теории функций и смежные проблемы: материалы Воронежской зимней математической школы. – Воронеж: Изд-во ВГУ. – 2011. – с. 356-357.
15
6. Ряжских А.В. Гидродинамический начальный участок при течении высоковязкой
жидкости в плоском канале / Физико-математическое моделирование систем: материалы VII Международного семинара. – Воронеж: Изд-во ВГТУ. – 2011. – с.199-202.
7. Ряжских А.В. Формулировка задачи в переменных Гельмгольца о гидродинамическом начальном участке при «медленном» течении ньютоновской жидкости в
плоском канале / Современные проблемы прикладной математики, теории управления и математического моделирования: материалы VI международной научной
конференции. – Воронеж: Изд-во ВГУ. – 2011. – с. 255-256.
8. Чертов Е.Д., Ряжских А.В. Идентификация ламинарного гидродинамического начального участка жидкости Оствальда-де Виля в цилиндрическом канале // Актуальные проблемы математики и информатики. – 2008. - №4. – с. 79-87.
9. Ряжских А.В. Гидродинамический начальный участок в кольцевом канале при
напорном течении ньютоновской среды (труды международной научной конференции «ММТТ-19»). Т.9. Секция 10. – Воронеж: Изд-во ВГТА. – 2006. – с. 47-48.
10. Ряжских А.В. К расчету ламинарного изотермического течения жидкости Оствальда-де Виля во входном участке плоского канала / Авиакосмические технологии:
труды XI международной научно-технической конференции. Ч.2. – Воронеж: Изд-во
ВГТУ. – 2005. – с. 226-230.
11. Ряжских А.В. Моделирование изотермической гидродинамики ньютоновской
жидкости на начальном участке цилиндрического канала / Современные проблемы
математики и математического моделирования: материалы конференции. – Воронеж: Изд-во ВГТА. – 2005. – с.200.
16
Документ
Категория
Физико-математические науки
Просмотров
142
Размер файла
336 Кб
Теги
кандидатская
1/--страниц
Пожаловаться на содержимое документа