Влияние магнитного поля и кривизны на турбулентную структуру течения

Вычислительные технологии Том 9, № 3, 2004

ВЛИЯНИЕ МАГНИТНОГО ПОЛЯ И КРИВИЗНЫ КАНАЛА НА ТУРБУЛЕНТНУЮ СТРУКТУРУ

ТЕЧЕНИЯ

У. С. АБДИБЕКОВ, Г. С. МАКАНАЛИНА Казахский национальный университет им. аль-Фараби, Алматы e-mail: uali@kazsu.kz, kvt@kazsu.kz

Turbulent flow of conducting fluid in a curved channel is considered under the action of a horizontal magnetic field. A mathematical model is constructed. Problem on the flow of a conducting fluid from a vertical slot into a channel of rectangular section, with longitudinal curvature has been solved numerically.

Введение

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

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

Пространственное течение проводящей жидкости в канале прямоугольного сечения в поперечном магнитном поле исследуется в течение многих лет [1, 2]. Пространственным течениям в условиях воздействия сильного магнитного поля посвящены работы [3, 4]. Большой интерес вызывает более сложное течение, когда на поток проводящей жидкости действуют кривизна канала и поперечное магнитное поле [5-10]. Данный вопрос мало исследован, так как экспериментальные исследования таких течений в основном посвящены определению сопротивления трубы.

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

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2004.

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

При изучении магнитогидродинамических машин принимается допущение о малости магнитного числа Рейнольдса (И.ет << 1), которое характеризует соотношение между процессами диффузии магнитного поля и его конвективным переносом при движении среды.

1. Постановка задачи и основные уравнения

В настоящей работе построена полуэмпирическая модель турбулентного сдвигового течения проводящей жидкости в канале со слабой кривизной H/R << 1, где H — высота канала, R — радиус кривизны канала.

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

dUi тт dUi тт dUi тт dUi 1 dP

+ Ui^- + U2—± + U3J-— = —— + dt dx1 ox2 dx3 p dx1

+ixS + dx;(-UiU3)+St0 (Uo - Ui) - ^R1'

dUo ÖUo ÖUo dU2 1 dP d , 2\ m

üt + Ui dxi + Uo äxx; + U3 dxx; =- + ^'

dUs , r, dU3 dU3 ÖU3 1 dP

-77Г + --г --+ U3—- =--тт--+

dt dx1 dx2 dx3 p dx3

d d U2

+— (-U3U1) + — (u.) + StoUs + -1, dx1 0x3 ' R

dUi + m + m = 0

dx1 dx2 dx3

Здесь x1 — продольная координата, x2 — поперечная координата, x3 — вертикальная координата; Ui, Ui — компоненты осредненной и пульсационной скоростей вдоль осей xi (i = 1, 2, 3); U0 — средняя расходная скорость по сечению канала; t — время; P — давление; p — плотность жидкости; St0 = (aB0L)/(pU0) — число Стюарта потока, характеризующее порядок отношения электромагнитной силы к силе инерции, где а — электрическая проводимость среды; B0 — величина магнитной индукции внешнего поля; L — характерный размер. Магнитное поле направлено параллельно оси x2, поэтому в системе уравнений (1) корреляции (u1u2) и (u2u3) отсутствуют [11].

Граничные условия задаются на твердой стенке для продольной составляющей в виде трения со скольжением

U = К (ycG),

U* к

где yc = xcU*/v; к = 0.4 — постоянная Кармана; G = 9 — коэффициент шероховатости стенки; U* — динамическая скорость на стенке; Ug — пристенная скорость; xc задается

так, чтобы выполнялось условие 30 < ус < 100. Для поперечной составляющей скорости задается условие прилипания и2 = и3 = 0.

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

и пи-

диг

+ пщ-

ди,

+ к

у/Ё (.

-у

2

п п — — бг] Е ) + - об,

3

Ч]

2

Е 3/2

3

гз

I

+

дхк дхк

+ и1 бзгПТп] — и бцп3п~ + [2п~п~Б3 Б3 — п~й;В] Б3 — u~п¡BiBs) = 0, К К р

где Е — энергия пульсационного движения; I = кх31 — х3/Н — масштаб турбулентности; а* = 1 — в, в учитывает условия замыкания пульсационных токов [1, 12], Бг, Б], Б3 — компоненты вектора магнитной индукции [г,], в = 1, 2, 3); б] — символ Кронекера.

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

2

о1

П2 3{ к/ о2/31 V дх3

ди 1

Пь П1 = ф2

П3

2! ди Л'

— П1П3 = I -- ^2,

дх

П

ф3 [1 + Бв)

2=

+ З^ЯЬф + 2^ — 6пБв + 4пБв2

п3 = - 1 - -

1

к) о2/3 \ дх.

ди 1

П3

ф

32

+ З^ЯЬф + 2^ — 8цБв + 2пБе

+ ^Ь) (ф2 + З^ЯЬф + 2^ — 6пБв + 4пБв2

п

1 + 2-

1

-х

ди 1

Па

Пл

ЗУ ' к) о2/3 \ дх3

ф2 (тф + ш^/цЗь) Ф — 2'рпф [1 + Б в) т(ф + у/цЯЪ) Ф

Ф = ф2 + З^ЯЬф + 2пЯЬ2 — 6пБе + 4^Бе2,

Е

о2/3

ди 1

дх3

+

ф=-

тЗЬ 66с2/3

ди 2

дх3

+ — в)1/3 — (^ +

ф2

1/3

2

2

2

2

3

2

2

2

2

1

2

6 = П 3 - ^ + Ф,

3/ 6 2'

Ф = 62 +

2\ з

и = + 2^3

2П + ЗД Я!2 - (6п + 3) °е + - 1 ) - 1,

Ф

ШГ)

Су3

Я! (Я!2 - 3Бе + 2Бе2) .

Безразмерные коэффициенты в (3) определены следующим образом:

/К _ а*аБ2/р

Бе

ои1 \

дхз )

Я!

ди 1

дхз

1

3/4

С \3/2 к)

1

3/4

2 Л с Р = 3 V - к

г Р у/ц = , Ш

2 - р, г = ш - р,

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

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

с

3

с

к

к

к

С

С

С

2. Численный алгоритм

Для численного решения системы (1) замкнутой моделью турбулентности (3) приведем следующий алгоритм, подробно описанный в работе [14]. Применяется метод расщепления по физическим процессам. На первом этапе предполагаем, что перенос осуществляется за счет конвекции и диффузии [15]:

/п+ з _ /п 1 _

/ / :(Ь1/п+ з + ^/п) + Ь2 /п + Ь3Г

т 2

2 1 /п+2 _ /п+ 3 1

т = 2

2

2

-2 - ь2/п),

/п+1 _ /п+ 3 1

1-1-= Ъ(Ь3/п+1 - Ь3/п),

т2

где Ь1, Ь2, Ь3 — операторы, включающие все производные по х1, х2, х3 соответственно; / = {и1*, и2*, и3*}. Из первого шага (4) находим /п+ з и подставляем во второй шаг системы (4). Подставляя полученное значение /п+ з в третий шаг, получаем /п+1. Каждое из уравнений (4) представляем в виде линейного алгебраического уравнения

А /+1 + Б/ + Сг/г-1 = Бг,

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

На втором этапе производится расчет давления, представленного в виде Р = Р1 + Р2

др

[16]. Продольная составляющая скорости определяется как и1 = а + 6——, где а, в —

дх1

рекуррентные коэффициенты, которые определяются из третьего шага (4). Давление Р1 находится из условия сохранения расхода в каждом поперечном сечении

S = JJ U1dx2dx3 п

в виде

N2-1

dPi дх1

S -Т 1 ^ 2

j=i

N3-1 , Ч N3-1 , ч

+ ai,j,k+1 ) д + у^ (ai,j+1,k + ^i,j + 1,k + 1) д

у^ (ai,j,k + ai,j,k+1) дХз + k=1 2

k= 1

2

ДХ2

N2-1

Е

j=1

2

N3-1

Е

k=1

(в ,j,k + ßi ,j,k+1 ) Л ™ , ^ (ßi ,j+1,k

2

-дхз +

+ в

i,j+1,k+1)

k= 1

2

Дхз

ДХ2

Две поперечные составляющие скорости представляются с поправкой давления в безразмерном виде:

U2 = U2 + т

дР2

U3 = U* + т

дР2

дх2 дх3

Подставляя соотношение (7) и найденное значение и1 в уравнение неразрывности, получаем двумерное уравнение Пуассона для каждого сечения по координате х\:

д2Р2 д2Р2 + 2

д^ дЩ дЩ

1 + ^ + 3

0x2, дх| V 0x1 ' дх2 ' дх3 которое решаем матричной прогонкой. Далее по (7) пересчитываем скорости и2, и3.

1

3. Результаты численного моделирования

На основе изложенной модели численно решается задача течения проводящей жидкости, вытекающей из вертикальной щели в канал прямоугольного сечения с продольной кривизной (рис. 1). Проводящая жидкость, вытекая из щели, попадает в искривленный заполненный канал, на который действует поперечное магнитное поле. Течение имеет следующие параметры: высота канала H = 90 мм, ширина щели d = 2 мм, число Рей-нольдса Re = 9300, параметр кривизны H/R = 0.025. Число Гартмана Ha = 82.5. Число Ha = B0L(a/pv) характеризует порядок отношения электромагнитной силы к силе вязкого трения, где v — коэффициент кинематической вязкости среды.

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

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

Рис. 1. Схема течения: I — вертикальная щель, II — расчетная область канала, III — выходное сечение канала.

Хз/Н3

х2/Н2

Рис. 2. Поле векторов скорости в поперечном сечении на расстоянии х\/И = 1.5 от начала канала.

По динамическим характеристикам можно судить, что поток во второй половине канала переходит в плавное течение, которое можно наблюдать по пространственным профилям продольной скорости, приведенным на рис. 3.

Проводящая жидкость, вытекая из вертикальной щели, попадает в заполненный криволинейный канал, где имеется поперечное магнитное поле. Плоское течение деформируется, как уже говорилось, центральная часть потока тормозится. Постепенно течение переходит в плавное движение, и верхняя часть потока отстает в силу торможения центробежными силами, вызванными кривизной канала. Все происходящие изменения средних динамических характеристик течения сильно отражаются на турбулентных характеристиках потока. Интерес представляет распределение кинетической энергии турбулентного движения (рис. 4). На рис. 5 приведены изолинии распределения кинетической энергии турбулентного потока в поперечном сечении. Представленные данные турбулентных характеристик свидетельствуют о взаимной корреляции пульсационных характеристик потока, динамике изменения пульсационных характеристик по всей длине канала, связанных с присутствием поперечного магнитного поля и кривизной канала.

Рис. 3. Пространственный профиль продольной скорости на расстоянии хх/И = 0.5 (а), хх/И = 1.5 (б) от начала канала.

Рис. 4. Распределение кинетической энергии турбулентного потока в поперечном сечении. х\/И = 0.5 (а), хх/И = 1.5 (б) от начала канала.

хз/нз

х2/Н2

Рис. 5. Изолинии распределения кинетической энергии турбулентного потока в поперечном сечении хх/И = 1.0 от начала канала.

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

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

[1] Бранновер Г.Г., ЦиноБЕр А.Б. Магнитная гидродинамика несжимаемых сред. М.: Наука, 1970. 380 с.

[2] Тананаев А.В. Течения в каналах МГД-устройств. М.: Атомиздат, 1979. 368 с.

[3] Hunt J.C.R., Leibovich S. Magnetohydrodynamic flow in channels of variable cross-section with strong transverse magnetic fields //J. Fluid Mech. 1967. Vol. 28. Pt 2. P. 241-260.

[4] Walker J.S., Ludford G.S.S., Hunt J.C.R. Three-dimensional MHD duct flows with strong transverse magnetic fields. Pt 2. Variable area rectangular ducts with conducting sides // J. Fluid Mech. 1971. Vol. 46. Pt 4. P. 657-670.

[5] Аитов Т.Н., КАлютик А.И., Тананаев А.В. Численное исследование динамики течения электропроводящей жидкости в области сложной геометрии при воздействии сильного магнитного поля // Магнитная гидродинамика. 1977. № 1. C. 73-78.

[6] БочЕнинский В.П., Тананаев А.В., Яковлев В.В. Экспериментальное изучение течения электропроводящей жидкости в изогнутых трубах кругового сечения в сильном магнитном поле // Магнитная гидродинамика. 1977. № 4. C. 62-65.

[7] Иевлев В.М. Численное моделирование турбулентных течений. М.: Наука, 1990. 273 c.

[8] Moser R.D., Moin P. The effects of curvature in wall bounded turbulent flows // J. Fluid Mech. 1987. Vol. 175. P. 479-510.

[9] Holloway A.G.L., TAVOULARIS S. The effect of curvature on sheared tubulence // J. Fluid Mech. 1992. Vol. 237. P. 569-603.

[10] ДжороговА Е.В., Новожилов В.В. Расчет установившегося турбулентного течения в криволинейном канале // Докл. АН СССР. 1983. T. 270, № 4. C. 819-824.

[11] Ковнер Д.С., Левин В.Б. О турбулентном течении электропроводной жидкости в трубе в продольном магнитном поле // Теплофизика высоких температур. 1964. Т. 2, № 5. C. 742-748.

[12] Воропаев Г.А., ПтухА Ю.А. Моделирование турбулентных сложных течений. Киев: Наукова думка, 1991. 166 c.

[13] АБДИБЕКОВ У.С. Влияние кривизны и поперечного магнитного поля на пульсаци-онную структуру турбулентного течения // Изв. НАН РК. Сер. физ.-мат. 2001. № 5. C. 60-63.

[14] АБДИБЕКОВ У.С., АмЕТОВ О.А., Смлгулов Ш.С. Численное моделирование турбулентного сдвигового течения в криволинейном канале // Вычисл. технологии. 2003. Т. 8 (спецвыпуск): Тр. Сов. рос.-каз. рабочей группы по вычисл. и информ. технологиям. C. 18-24.

[15] ЯНЕНКО Н.Н. Метод дробных шагов решения многомерных задач математической физики. Новосибирск: Наука, 1967. 196 с.

[16] ФЛЕТЧЕР К. Вычислительные методы в динамике жидкости. М.: Мир, 1991. Т. 2. 552 с.

Поступила в редакцию 22 июля 2003 г., в переработанном виде — 30 марта 2004 г.