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

УДК 532.529.5

И. В. Моренко, В. Л. Федяев

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

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

теплоотдача.

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

Keywords: viscous incompressible fluid flow with particles, flow past the circular cylinder, drag, heat transfer.

A viscous incompressible fluid flow with particles past the circular cylinder in a channel and heat transfer was considered. A numerical simulation of unsteady laminar flow of the mixture was carried out using Euler continuum approach. The corresponding boundary value problem is solved numerically using the finite volume method. The influence of density, specific heat, thermal conductivity of the particle material and particle diameter on the drag and heat transfer of the cylinder was estimated.

Введение

Разработка новых технологий и оборудования разделения, очистки гетерогенных сред [1], перемешивания их [2], охлаждения либо нагрева тел, контактирующих с этими средами, и многих других процессов требует проведения исследований динамики многофазных сред в областях сложной формы.

Результаты численного моделирования обтекания кругового цилиндра потоком газа с частицами с использованием уравнений Навье-Стокса, когда число Рейнольдса Re=100, число Стокса Stk=0.01, 0.1, 1, 10, 100, приводятся в статье [3]. При описании движения частиц учитывается уточненная сила сопротивления Стокса, силы гравитации, подъемная сила Саффмана.

Влияние твердой дисперсной примеси на структуру двухфазного следа за поперечно обтекаемым круговым цилиндром (Re=100) рассматривается в работе [4]. При этом двумерное нестационарное течение сжимаемого газа описывается уравнениями Навье-Стокса с дополнительными источниковыми членами, моделирующими воздействие примеси. На входе в расчетную область частицы располагаются неравномерно, со сгущением к середине потока. Объемная концентрация их в натекающем потоке варьируется в диапазоне 0 + 3 -10"6. При этом массовая концентрация изменяется в пределах 0 + 0.58 . Отмечается, что, несмотря на очень малую инерционность частиц, малое значение числа Стокса (Stk и 0.05), характер движения частиц примеси существенно отличается от течения несущего газа. В частности, частицы выпадают на поверхность цилиндра. Это происходит вследствие того, что частицы примеси движутся в потоке несущего газа в свободномолекулярном режиме, поэтому сила аэродинамического сопротивления, действующая на них со стороны газа, сильно отличается от силы

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

В работе [5] акцентируется внимание на особенностях движения частиц в вихревом следе за обтекаемым круговым цилиндром при Re=30^150. Траектории частиц рассчитываются с использованием Лагранжева подхода. Подробно анализируется дисперсия частиц, характеризующая отклонение частиц в вертикальном направлении от среднего положения. Показывается, что дисперсия зависит от времени, числа Рейнольдса потока, числа Стокса. При малых числах Стокса частицы движутся вдоль линий тока несущей среды. С увеличением числа Стокса частицы концентрируются в сравнительно узких областях вблизи границ вихрей. Поведение частиц в следе за цилиндром зависит, главным образом, от процесса образования областей больших градиентов скоростей между смежными вихрями с противоположной циркуляцией.

Результаты экспериментального и численного исследования турбулентного обтекания водной суспензией кругового цилиндра, находящегося в канале над слоем осадка, приводятся в работе [6]. Концентрация частиц в воде достигает 50%, их диаметр составляет 3.6-10-4 м. Оценивается толщина слоя осадка, изменение формы верхней поверхности этого слоя под воздействием потока суспензии.

Исследование дисперсии частиц в следе за круговым цилиндром при высоких числах Рейнольдса потока с использованием метода дискретных вихрей и подхода Лагранжа для расчета траекторий частиц выполнено в [7]. Числа Стокса варьируются от 0.01 до 1000.

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

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

Имеющиеся в литературе выводы о влиянии многофазности среды на конвективный теплообмен, зачастую, противоречивы.

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

Однако наблюдается и обратное [10]. Более предметное исследование влияния размера частиц примеси на вынужденную конвекцию при обтекании цилиндра вязкой жидкостью показало, что без учета локального теплового равновесия между несущей и дисперсной фазами при изменении отношения диаметра цилиндра к диаметру частиц от 10 до 100, числа Рейнольдса от 1 до 250, отношения коэффициентов теплопроводности материала частиц и несущей среды в пределах 0.01-1000 с увеличением радиуса частиц число Нуссельта может как увеличиваться, так и уменьшаться.

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

Комплексное влияние наночастиц (Н100 нм) на гидродинамику и теплообмен поперечно обтекаемого кругового цилиндра исследуется в [12]. При моделировании ламинарного течения жидкости используются уравнения Навье-Стокса с дополнительными членами, учитывающими влияние наночастиц, а также уравнение энергии. Расчеты проводятся при числах Рейнольдса 80^180. Объемная концентрация частиц варьируется в интервале от 0 до 25%. Установлено, что с ростом концентрации частиц в потоке увеличивается коэффициент сопротивления цилиндра, число Струхаля, число Нуссельта, уменьшается амплитуда подъемной силы цилиндра.

Следует заметить, что большинство работ по численному моделированию обтекания тел гетерогенной смесью выполнено с использованием Лагранжева подхода [4, 12-14]. При этом объемная концентрация частиц примеси не может превышать 10%. Вместе с тем, не практике часто встречаются

более густые смеси. Поэтому для расчета движения их, теплообмена с обтекаемыми телами следует применять Эйлеров континуальный подход [6, 15 и др.].

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

Определяющие соотношения

Математическое описание

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

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

При проведении расчетов основной несущей фазой является вода; характерная скорость движения смеси порядка 10"5 ^ 2.5 -10~3 м/с. Следовательно, числа Маха малы, сжимаемостью среды можно пренебречь.

Расчетная область представляет собой прямоугольник длиной 6.0 м, шириной 0.4 м, в котором на удалении 1 м от входного сечения на равном расстоянии от боковых границ располагается обтекаемый цилиндр с радиусом К =0.05 м. Данная область покрывается сеткой треугольных элементов, сгущающихся к поверхности обтекаемого тела. На поверхности тела находится 100 узлов, общее количество узлов сетки 8 тыс., треугольных элементов 16 тыс. Начало системы координат х10х2 располагается в центре цилиндра, ось 0х; направляется вдоль набегающего потока.

Основные уравнения. Уравнения сохранения массы записываются в виде [15]:

д -

—(а,.р,.)+У,- ■У(а1р1) = 0, I = 1,2; (1)

„ - д - д - т -

где т - время, У = I--ь ]-, VI = ул I + уп ] -

дх1 дх2

вектор скорости I -ой фазы; ах, а2 - объемная

концентрация, соответственно, дисперсной фаз (а1 + а2 = 1);

д

несущей

— {а2)+ -У(а2 ) = 0. (2)

Уравнение импульсов для , -ой фазы в отсутствие гравитационных сил имеет вид

рV)+ V,- ^(р^) = -а,.Чр, + У-П, + Я,. (3)

Здесь р1 - среднее давление , -ой фазы;

П =а1^1(уу, +(у7у,) ^) + а1 ^ст -—Уу,1 - тензор

напряжений в , -ой фазе; ¿и,, - сдвиговая, объемная вязкость; I - единичный тензор; (-1)'

Ri = -

V

-Jy((2 -Vl)

- vi (dV - осредненная сила

взаимодействия между фазами, V - контрольный

объем, у = а2р2fj8, 8 = 2р2а2 ¡9ц, f =

CD Re2

24

коэффициент сопротивления частицы

124 (l + 0.15Re2f87 )/Re21, Re21 < 1000

0.44,

Re21 > 1000

Re2i =-

2р1 |v2 - vi а

число

Рейнольдса

относительного движения фаз.

В отсутствие соударения частиц коэффициент динамической вязкости ¿2 определяется по формуле [6]:

¿2 =¿2/ + ¿2к , где и2 / - коэффициент трения,

H2k =

аа2р2

(1 + 0.4(1 + Cd )(3CD - 1)а2g0 ]

3 (3 - сд )

слагаемое, обусловленное кинетической энергией частиц, £0 - радиальная функция распределения. Объемная вязкость ст2 полагается равной нулю.

Уравнение переноса тепла в несущей и дисперсной фазах имеет вид:

д_

~дт

(pfij,)+ v, ■V(pfJi) = V.(VTi) + Qt: (-1)

(4)

где Q = i^-JЬпага2р21 ( -T)dV ,

V

Р2Х = Л/Nu2//2a - коэффициент теплообмена между несущей и дисперсной фазами, число Нуссельта рассчитывается по формуле:

№21 = 2 + 0^е1212Рг2/13, (5)

число Прандтля Рг21 = р/cplv/|\ .

Уравнения (1)-(4) дополняются граничными и начальными условиями. Во входном сечении, перпендикулярном оси 0 х1, примем: уп = у21 = ™, у12 = у22 = 0; операционное давление р0 =0.1 МПа; температура несущей и дисперсной

фаз одинакова: T1 = T2 = T0 = const; концентрация а2 = а0 = const.

На выходе из расчетной области воспользуемся следующими граничными

dv., dv. 2 др. дТ.

условиями: -1 = = = = 0 (i = 1,2),

дх1 дх1 дх1 дх1

да2

дх1

= 0 ; что означает выравнивание на выходе как

гидродинамических, так и тепловых, концентрационных характеристик сред.

На боковых границах - стенках канала реализуются условия прилипания несущей и

дисперсной фаз v1 = v2 = 0 ; температура:

T = T = T да2 = 0 11 2 10 > ~ u •

дх2

Считается, что на поверхности обтекаемого цилиндра имеет место прилипание несущей и дисперсной фаз v1 = v2 = 0 ; идеальный тепловой контакт: T1 = T2 = Tw. При этом удельный тепловой поток qt = дНГ^дп = а* (Tw - T*) (i = 1,2 ), где п -

т*

внешняя нормаль к поверхности цилиндра, Ij -температура i -ой фазы на удалении от цилиндра (среднемассовая температура), а* - искомый локальный коэффициент конвективного

теплообмена.

В начальный момент времени т = 0 среда мгновенно начинает двигаться, v(1 = w, vj2 = 0 ,

Р = Po, T = T0, а2 = а0 .

Интегрирование соответствующих

уравнений (1)-(4) производится с помощью программного комплекса. Выбирается неявный решатель второго порядка точности, используется процедура коррекции давления SIMPLE. Критерий сходимости по всем параметрам принимается равным 10-8, максимальное количество итераций на каждом временном шаге 40. Шаг по времени выбирается таким образом, что при изменении времени в пределах этого шага малый объем жидкости смещается на одну ячейку сетки.

Входные параметры

При проведении вычислительных экспериментов число Рейнольдса потока

2p1wR

Re = —1- принимает значения 1, 10, 40, 100, 250.

В результате движение среды является ламинарным, корректна плоская постановка задачи. Из этих данных находится скорость w несущей среды и частиц на входе в расчетную область. Во входном сечении фиксируется объемная концентрация частиц а0 = 0.5, температура жидкости и частиц T0=300K На боковых стенках канала поддерживается постоянная температура T0 , постоянна и температура поверхности обтекаемого цилиндра Tw =350К.

и

D

Для данного диапазона изменения температур полагается, что плотность рх, коэффициент динамической вязкости , удельная теплоёмкость ср1 несущей среды - воды зависят от

температуры линейно (см. табл. 1), а коэффициент теплопроводности Я1 =0.6 Вт/(м-К). Средние значения теплофизических характеристик материала, диаметр вбрасываемых частиц в ходе эксперимента варьируются.

Таблица 1 - Теплофизические характеристики воды

т, к рх, кг/м3 /и1, кг/(м-с) ср1, Дж/(кг-К)

300 996.59 0.000852 4164

325 987.09 0.000529 4127

350 973.46 0.000369 4082

Результаты расчетов

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

На первом этапе вычислительных экспериментов варьируется только плотность материала частиц р2 при фиксированных ср2 =4164 Дж/(кг-К), 12=0.6 Вт/(м-К), диаметре

частиц ёр = 103 м и скорости набегающего потока

^ = 10"4 м/с. Изменение продольных составляющих действующих на цилиндр сил трения , давления

¥р, силы сопротивления ¥в в зависимости от р2

показывается на рис. 1. Из приведенных данных следует, что на графиках кривых = (р2),

¥р = ¥р (р2) и ¥в = ¥в (р2) можно выделить три характерных участка. На первом и третьем участках, когда 1000 кг/м3 < р2 < 1500 кг/м3 ,

2000 кг/м3 < р2 < 3000 кг/м3 соответственно, данные функции близки к линейным, причем на третьем участке прямые более пологие. Второй участок (1500 кг/м3 < р2 < 2000 кг/м3 ) является переходным.

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

Например, если р2 =1000 кг/м3 , то коэффициент

а12 =23 Вт,

г/(м2 • К);

р =3000 кг/

м

а12=30 Вт^(м2 • К). Интенсификация теплообмена,

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

Р, Н

---------

, , , , , . , , , , , , -

1000

1500

2000

2500

рп. кг/м-

Рис. 1 - Зависимость продольной составляющей силы сопротивления трения (пунктирная

линия), давления ¥р (штрих-пунктирная линия)

и ¥в = + ¥р (сплошная линия)

Обнаруженную закономерность

подтверждают эксперименты, в которых скорость среды на входе в расчетную область меняется в названном интервале. Причем с ростом скорости среды влияние плотности материала частиц на теплоотдачу цилиндра увеличивается нелинейно (рис. 2).

Рис. 2 - Коэффициент теплоотдачи цилиндра: плотность материала частиц р2 =160 кг/м3 (нижняя кривая), р2 =8030 кг/м3 (верхняя кривая)

Изменение коэффициента теплоотдачи цилиндра в случае, когда изменяется скорость среды

на входе в расчетную область и характеристики материала частиц в соответствии с табл. 2 [8], иллюстрирует рис. 3. При этом диаметр частиц остается постоянным, равным 10-3 м.

Таблица 2 - Теплофизические характеристики материала частиц

Материал Муассанит Сталь Дерево

р2 , кг/м3 3.21 8030 160

Ср2 , Дж/(кг-К) 298 502.48 1200

Л,, Вт/(м-К) 420 16.27 0.055

о;:,Вт/(м2-К) 2500

2000

1500

1000

500

-

1 ч

/

; / /2

Г /4 о *

0.001 0.002 ^„/с

Рис. 3 - Коэффициент теплоотдачи цилиндра: 1 -муассанит; 2 - сталь; 3 - однородная среда (вода); 4 - дерево

Установлено, что для названных материалов доминирующее влияние на теплообмен оказывает величина коэффициента

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

Следует заметить, что изменение удельной теплоёмкости материала частиц ср2 в пределах

300^1200 Дж/(кг-К) не оказывает заметного влияния на теплоотдачу цилиндра и его сопротивление; рассчитанные значения коэффициента теплоотдачи цилиндра с учетом теплового взаимодействия фаз по формуле (5) несколько выше, чем без его учета.

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

(=420.0 Вт/(м-К)), теплообмен ухудшается (рис. 4). Этот эффект проявляется тем заметнее, чем больше скорость среды на входе в канал.

2500

2000

1500

1000

500

-

- X

/ /.'

/ /

//

; /

0.001 0.002 ^м/с

Рис. 4 - Коэффициент теплоотдачи цилиндра при обтекании потоком жидкости с частицами:

-3 .. __„_________„________ 1 щ-2

ёр = 10^ м -

сплошная линия;

ёр = 10" м -

пунктирная линия

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

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

При плотности материала частиц меньше или равной плотности несущей среды частицы, обладая сравнительно малой инерционностью, будут двигаться, в основном, непосредственно с потоком несущей среды. Поэтому объемная концентрация их вблизи поверхности обтекаемого тела будет близка к объемной концентрации частиц в набегающем потоке, особенно частиц малых размеров. Действительно, из результатов расчетов следует, что при плотности материала частиц р2=3.21 кг/м3, существенно меньшей плотности

воды, при диаметре частиц

10 -3 м

и скорости

набегающего потока = 10-5 ^2.5•Ю-3 м/с средняя объемная концентрация частиц на поверхности

" (1р = 10-2 м

цилиндра «2 = 0.492 + 0.493; при

«2 = 0.317 ^ 0.298.

Следовательно, определяя эффективный коэффициент теплопроводности смеси по формуле [19]

д = д 2Л1 + Л - 2а2 (Л -Л2 )

2Д + Л + а2 (Л - Л ) убеждаемся, что вблизи поверхности обтекаемого цилиндра при увеличении размера частиц примеси коэффициент Д2 уменьшается, а значит будет уменьшаться коэффициент теплоотдачи цилиндра, что и наблюдается.

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

Заключение

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

При уменьшении плотности,

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

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

Если частицы «тяжелые», плотность их материала больше плотности несущей среды, то увеличение размера частиц при хорошей теплопроводности материала теплообмен улучшает, а при плохой - ухудшает.

Таким образом, при фиксированной концентрации частиц примеси на входе в расчетную область для интенсификации теплообмена обтекаемого дисперсной смесью цилиндра за счет увеличения размера частиц примеси необходимо, чтобы материал примеси был бы «легким», с плохой

теплопроводностью, что характерно для теплоизоляционных материалов; либо «тяжелым» и хорошо проводящим тепло, таким, например, как материалы металлов.

Необходимо также иметь в виду, что при одинаковых размерах частиц в первом случае гидродинамическое сопротивление обтекаемого цилиндра будет меньше, чем во втором.

Литература

1. М.И. Фарахов, А.Г. Лаптев, Вестник Казанского технологического университета, 9, 152-158 (2011).

2. Л.Р. Минибаева, А.Г. Мухаметзянова, А.В. Клинов, Вестник Казанского технологического университета, 11, 201-210 (2010).

3. J. Yao, Y. Zhao, G. Hu, J. Fan, K. Cen, Aerosol Science and Technology, 43, 174-187 (2009).

4. А.Н. Волков, Ю.М. Циркунов, Математическое моделирование, 15, 7, 98-110 (2003).

5. L. Liu, F. Ji, J. Fan, K. Cen, Journal of thermal science, 13, 4, 344-349 (2004).

6. Zh. Zhao, H.J.S. Fernando, Environ Fluid Mech, 7, 121142 (2007).

7. Y. Huang, W. Wu, H. Zhang, Powder Technology, 162, 73-81 (2006).

8. G.F. Al-Sumaily, A. Nakayama, J. Sheridan, M.C. Thompson, International Journal of Heat and Mass Transfer, 55, 3366-3378 (2012).

9. S.W. Hsiao, P. Cheng, C.K. Chen, Int. J. Heat Mass Transfer, 35, 12, 3407-3418 (1992).

10. K. Nasr, S. Ramadhyani, R. Viskanta, J. Heat Transfer, 116, 1, 73-80 (1994).

11. P.X. Jiang, Z. Wang, Z.P. Ren, B.X. Wang, Ехр. Therm. Fluid Sci, 20, 1, 45-54 (1999).

12. S. Sarkar, S. Ganguly, G. Biswas, International Journal of Heat and Mass Transfer, 55, 4783-4799 (2012).

13. И.В. Моренко, В.Л. Федяев, Экологический вестник научных центров Черноморского экономического сотрудничества, 4, 52-58 (2010).

14. И.В. Моренко, В.Л. Федяев, Тепловые процессы в технике, 3, 6, 242-252 (2011).

15. Р.И. Нигматулин, Динамика многофазных сред. Ч. I. Наука, Москва, 1987. 464 с.

16. T. Farrant, M. Tan, W.G. Price, Computers & Fluids, 30, 211-236 (2001).

17. S. Mittal, A. Raghuvanshi, Int. J. Numer. Meth. Fluids, 35, 421-447 (2001).

18. D.L. Young, J.L. Huang, T.I. Eldho, Comput. Methods appl. Mech. Eng, 190, 5975-5998 (2001).

19. J.W. Smith, J.C. Paradi, J. Pipelines, 3, 43-52, (1982).

© И. В. Моренко - канд. техн. наук, с.н.с., ИММ КазНЦ РАН; В. Л. Федяев - д-р техн. наук, зав. лаб. МТП, ИММ КазНЦ РАН, morenko@mail.knc.ru.