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

Вычислительные технологии

Том 11, № 1, 2006

ТЕЧЕНИЕ ЭЛЕКТРОПРОВОДЯЩЕЙ ЖИДКОСТИ В КАНАЛЕ С ЧАСТИЧНО ПРОВОДЯЩИМИ

СТЕНКАМИ

А.М. Бубенчиков, В. С. Попонин, А. В. Колесникова Томский государственный университет, Россия e-mail: [email protected], [email protected]

The simplest algorithm that solves a self-similar problem of the flow of conductive fluid in a channel with partially conductive walls embedded in a transversal magnetic field. Calculations show that the velocity field has two maximums in this case.

Теоретическому и экспериментальному изучению установившихся ламинарных течений в круглых трубах посвящены многие работы. Аналитические решения задач о распределении скорости и индуцированного магнитного поля представляются в виде тригонометрических рядов [1], которые довольно громоздки и иногда плохо сходятся даже при больших значениях чисел Гартмана [2, 3]. Имеются также теоретические работы, базирующиеся на приближенном подходе, основанном на выделении в потоке характерных областей, для каждой из которых выписываются свои приближенные уравнения. Такой метод впервые был применен Дж. Шерклифом для трубы с непроводящими стенками [4] и затем обобщен на случай проводящих стенок [5]. Этот метод впоследствии был развит и усовершенствован А.Г. Куликовским [6]. Численное решение внутренних автомодельных задач магнитной гидродинамики выполнено в [7, 8]. Причем в [7] для реализации математической постановки задачи использован метод установления с расщеплением вычислительного процесса по координатным направлениям. В [8] интегрируются непосредственно стационарные уравнения магнитной гидродинамики и построенная численная процедура опирается на использование метода простой итерации.

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

Будем рассматривать течение электропроводящей жидкости на участке стабилизированного движения в канале круглого поперечного сечения, находящемся в однородном по длине канала электромагнитном поле. В работе [9] Дж. Хант показал, что в трубах

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

постоянного сечения при воздействии на них поперечным магнитным полем устанавливается течение с параллельными оси канала линиями тока. Это дает возможность исключить нелинейные конвективные члены в уравнении движения, а градиент давления считать функцией только продольной координаты. Поэтому в математической постановке задачи останется единственное из трех возможных скалярных уравнений импульсов, определяющее продольную компоненту скорости. Кроме того, в постановку задачи будет входить уравнение для продольной компоненты индуцированного магнитного поля. Если рассматривать течение с постоянным расходом, то определяющие уравнения лучше обезразмерить с помощью следующих масштабных параметров: г* = x* = R — радиус канала; U* = V0 — средняя скорость потока; B* = a^B0V0R; p* = цV0/R. Здесь B0 — величина магнитной индукции внешнего однородного магнитного поля; р, ц, a, ц0 — соответственно плотность, вязкость, проводимость и магнитная проницаемость жидкости. Тогда безразмерный вид уравнений будет следующим [10]:

ATT тт2( ndB sin Q8B\ '

△U=Ha\c°s ^ - —~£ш) + pz; (1)

BU sin Q8U

△ =c°s elT (2)

причем угловая координата в отсчитывается от направления вектора B0 оси x по стрелке

часов. Здесь г — радиальная координата; Ha = B0R\J a/ц — число Гартмана; p'z = dp/dz —

Л 1 д ( д\ 1 д2 перепад давления, приходящийся на единицу длины; △ = — —- г— +—— плоский

г дг \ дг J г2 дв2

оператор Лапласа, записанный в полярных координатах.

Уравнение импульсов (1) интегрируется с использованием условия прилипания жидкости на стенки канала:

U |г=1 = 0, (3)

а уравнение индукции (2) с условиями, характеризующими проводящие свойства стенок, имеет вид

дB

дг =0; b

дг Гьз

= 0. (4)

Г2.4

Г2

в= 0 г4

Рис. 1. Контур поперечного сечения с указанными на нем зонами изоляции и проводимости (жирные линии).

Как видно из рис. 1, стенки рассматриваемого канала выполнены в виде рельсовой камеры, в которой два рельса (Г1? Г3) являются участками идеальной проводимости, а два других (Г2, Г4) — участками изоляции канала.

Введем сетку, покрывающую область интегрирования с равномерными шагами к и П соответственно по г- и ^-координатам: к = 1/Ж, п = 2п/М. Здесь N и М — заданные числа, определяющие частоту разбиения исходных координатных отрезков г £ [0,1], в £ [0, 2п].

Аппроксимируем дифференциальные члены, входящие в (1), следующим образом:

1 ^ [ ди\

г дг \ дг )

г,]

1

кг.

иг,]+1 — иг, ] ]1-к--г]-2

иг,] иг,]—1 к

+ 0(к2);

1 д!и

г2 ~д¥

1 иг+1,] - 2иг,] + Цг-1,] + 0(^)

г,]

г]

п2

(5)

(6)

Тогда, если придерживаться обозначений [11], с учетом (5) и (6) разностный аналог (1) будет иметь вид

ариг] = ае11г-1,] + иг+1,] + asUг,j_1 + апиг,]+1 — Ь — р'г.

(7)

где

ае ат 2 2, as

г] п2

г-1 г] к2

ап =

г] к2

ар — ае + + + ап,

вт вдВ

(8)

Ь = На2^,], ^ = ( сов в— —

дг г дв

Задачи о течении жидкости в трубах и каналах с заданным расходом целесообразно решать с использованием подхода Л.М. Симуни [12]. Согласно этому подходу для сеточного значения скорости делаем представление

иг,] = ^ + ]. (9)

Тогда из уравнений (7) получим следующие две системы уравнений:

Ь; (10)

а] ^г,] = ае^г_1,] + ат Zг+l,j + asZг,j_l + ап^г,]+1 — 1- (11)

Проинтегрируем численно обе части уравнения (9), тогда для градиента давления можно получить соотношение

(1 — Л,)

Рг =

3,

(12)

Здесь и Зг — разностные представления следующих интегралов:

1 и - 1 /

s s

zds,

где Б — площадь поперечного сечения канала.

1

Таким образом, вместо системы (7), содержащей неизвестный градиент давления р/ , мы разрешаем две системы — (10) и (11), в которые указанный градиент не входит, после чего по соотношению (12) находим р'г и далее по (9) — сеточные значения продольной компоненты скорости.

Разностный аналог уравнения для В будет таким же, как и для Ш, т. е. иметь вид (10), за исключением коэффициентов Ь, которые в данном случае будут определяться соотношением

Ь = сое в-- _ ап в _ и*-1,з г 2к Т] 2п

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

Вц = («вВг-1] + аадВг+1] + asBij-l + апВг,]+1 _ Ь')/ар, (13)

где коэффициенты ав,а3,ап, аад определяются формулами (8). Формула (13) вместе с соотношением (14) для В0 используется для простого пересчета значений В во всех внутренних точках области, т. е. применяется многократно до сходимости.

Чтобы найти решение систем разностных уравнений (10), (11) и (13), опираясь на использование полярных координат (т, в), необходимо определить значения всех искомых функций в начале используемой системы координат, т. е. при т = 0. В работе [10] найдены соотношения, определяющие значения искомых функций в начале полярной системы

отсчета. Наиболее простая форма этих соотношений определяется как

1 М

Фо = мЕ ф1'к, ф = (14)

к=0

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

Соотношение (14) является приближенным и обладает погрешностью аппроксимации порядка к2. В случае, если Ь в (10) равно нулю (или Ь' = 0 в (11)), то оно переходит в точное разностное соотношение.

Как видно из (10), (11) и (13), дискретный аналог определяющих уравнений может быть представлен следующим образом:

арФр ^ ^ апЬФпЬ + Ь,

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

Ф

фр

(Е апьФПь + Ь)/ар, (15)

а в начале полярной системы отсчета — по соотношению

1

фо = ф1,к, (16)

где звездочкой отмечены соседние значения рассчитываемой величины, находящиеся в памяти вычислительной машины. Для соседних точек, к которым уже обращались в ходе текущей итерации, это новые рассчитанные значения, для остальных — значения с предыдущей итерации. Когда перебор всех внутренних точек по соотношениям (15), (16) закончен, а также перевычислены новые значения расчетных величин на границе (если это необходимо), одна глобальная итерация метода Гаусса — Зейделя считается законченной.

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

= 0, г = !^М. (17)

Здесь М — полное число точек по окружной координате. Аналогичные условия мы используем для Ш и Z. Так как участки проводимости расположены симметричным образом относительно в, значение М должно быть четным. Пусть М\ — число точек по в в секторе проводимости, тогда граничные условия (4) в разностном виде можно представить следующим образом:

_, , M M2

Вм,г = BN-м, i = 1 + M! U i = -у + 1, -гт + Mi; (18)

MM

BN,i = 0, i = Mi + 1, — [J i =y + Mi + 1,M. (19)

Реализация построенного алгоритма выполнена на языке программирования Visual C++.

Первоначально нами было рассмотрено течение проводящей жидкости в канале с непроводящими стенками и рассчитано сопротивление магнитно-гидродинамического потока для круглой трубы при различных числах Гартмана. На рис. 2 представлены результаты сопоставления расчетных величин с данными опытов Дж. Гартмана [13].

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

Рис. 2. Результаты расчета (сплошная линия) и экспериментальные данные Дж. Гартмана [13] (символы о). Все расчеты выполнены при а = 0. Здесь Л* = 64/Ие.

X 1

Рис. 3. Поверхность скоростей в канале с непроводящими стенками (На = 7, а = 0).

Рис. 4. Поверхность скоростей в канале с частично проводящими стенками (На = 7, а = п/2).

Так что в результате распределение скоростей остается симметричным лишь относительно двух центральных сечений (сечений, проходящих через ось канала), одно из которых параллельно, а другое перпендикулярно В0. При малых же числах Гартмана (На < 3) поверхность скоростей является практически осесимметричной. В то же время с ростом На описанные выше характерные особенности в распределении скоростей проявляются все более отчетливо.

Поверхность скоростей в канале с частично проводящими стенками в случае, когда угловой размер окна проводимости а равен п/2, приведена на рис. 4.

Как видим, в этом случае течение не является осесимметричным, но имеет две центральные плоскости симметрии, переходящие через ось трубы параллельно и перпендикулярно вектору напряженности внешнего однородного магнитного поля. Указанная симметрия поверхности скоростей есть отражение физического содержания задачи, поскольку при проведении вычислений мы не выделяли прямых или поверхностей, на которых выставлялись бы условия вида дГ/дп = 0, где п — нормаль к упомянутой прямой или поверхности. В принципе в данном случае можно было бы построить численные решения на одной четверти круга и далее дополнить его по соображениям симметрии. Однако при проведении настоящих расчетов мы использовали универсальную сквозную (по всей физической области) технологию проведения вычислений, опирающуюся на использование соотношений (14). Такая технология пригодна и для случая, когда нельзя выделить даже одной плоскости (линии) симметрии или антисимметрии.

Из рис. 4 следует, что в этом случае поверхность скоростей имеет два максимума, расположенных симметрично относительно сечения, которое проходит через ось трубки и вектор магнитной индукции. Причем эти максимумы выражены тем значительнее, чем выше число Гартмана. Безразмерный градиент давления в рассматриваемом случае совпадает с безразмерной величиной сопротивления единичного участка канала, т. е. с коэффициентом сопротивления, поэтому Cf = р/. Число Гартмана как безразмерный комплекс вида На = В0К (а/ц)0'5 может изменяться при вариации В0 либо геометрического размера задачи К, либо с изменением отношения физических параметров а/ц. Однако в настоящей работе для простоты рассуждений мы соотносим рост или уменьшение числа Гартмана с изменением напряженности внешнего однородного магнитного поля. Продолжая даль-

нейший анализ влияния числа Гартмана, следует отметить, что с его ростом существенно увеличивается Cf как для каналов с а = 0, так и для каналов с проводящими вставками.

При течении проводящей жидкости в канале с изолированными стенками возникает структура поперечных токов, показанная на рис. 5. Как видим, у стенок канала траектории зарядов загущены (эффект Дж. Гартмана для контуров тока). Расчеты показывают, что движением жидкости индуцируется такое электромагнитное поле, которое препятствует перемещению жидкости в направлении среднего потока. Поэтому, если среднее движение жидкости осуществляется от нас перпендикулярно к плоскости чертежа (рис. 5 или 6), то в ядре потока заряды должны перемещаться сверху вниз. Сила Лоренца будет направлена здесь против среднего движения. Так как траектории движения заряженных частиц все замкнутые, частицы, двигавшиеся изначально сверху вниз, рано или поздно приобретают обратное направление, т. е. снизу вверх. Обратное направление на этом участке движения имеет и сила Лоренца, которая теперь уже не сдерживает жидкие частицы, а, напротив, их ускоряет. В результате этого поверхность скоростей имеет заполненный вид (большие радиальные градиенты скорости у поверхности) и значительный плоский участок в приосевой зоне (см. рис. 3).

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

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

Таким образом, интегрально в случае канала с проводящими вставками заряженные

Рис. 5. Траектории движения заряженных частиц (а = 0).

Рис. 6. Траектории движения заряженных частиц (а = 2п/3).

я 2

S ~2

1.8-

1.6:

1.4-

1.2:

О 0.1 0.2 0.3 0.4 0.5 р

Рис. 7. Распределения относительного коэффициента трения q = Л/Ло в зависимости от углового размера окон проводимости p = а/2п. Здесь Ло отвечает расчету при а = 0; кривые 1-3 соответствуют Ha = 4, Т, 10.

частицы, так же как и в реперном варианте, движутся по замкнутым траекториям, только при этом в области течения отсутствуют участки возвратного движения зарядов, на которых сила Лоренца способствует перемещению жидких частиц в направлении основного потока. Поэтому с введением участков проводимости сопротивление движению жидкости всегда увеличивается. Если а = 0, то участки возвратного движения зарядов в какой-то мере компенсируют тормозящее действие поперечных токов в центральной части канала, в то время как при а = п/2 выстраивается целая полоса с направлением перемещения зарядов, перпендикулярным к вектору BQ. Эта полоса оказывает сдерживающее влияние на движение жидких частиц, осуществляемое в направлении среднего потока.

На рис. Т представлены распределения коэффициентов трения А, рассчитанные для канала с проводящими вставками, отнесенные к соответствующей величине А0, найденной при а = 0. Как видно из представленных обобщающих графиков, сопротивление только растет с увеличением как чисел Гартмана, так и углового размера окон проводимости.

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

[1] Ватажин А.Б., Любимов Г.А., Регирер С.А. Магнитодинамические течения в каналах. М.: Наука, 197G. 672 с.

[2] Ihara S., Tajima K., Matsushima A. The flow of conducting fluids in circular pipes with finite conductivity under uniform traverse magnetic fields jj Bull. Sci. and Engng. Res. Lab. Waseda Univ. 1963. Vol. 22. P. 1-9.

[3] Ihara S., Tajima K., Matsushima A. The flow of conducting fluids in circular pipes with finite conductivity under uniform traverse magnetic fields jj Trans. ASME. 1967. Vol. 89, N 1.

[4] Sheroliff J.A. Steady motion of conducting fluids in pipes under traverse magnetic fields // Proc. Cambrige. Philos. Soc. 1953. Vol. 49, N 1. P. 136-144.

P. 29-36.

Sheroliff J.A. Steady motion of conducting fluids in circular pipes under traverse magnetic fields // J. Fluid Mech. 1956. Vol. 1, pt 6. P. 644-666.

[6] Куликовский А.Г. О модельных стационарных течениях проводящей жидкости при больших числах Гартмана // Изв. АН СССР. Механика жидкости и газа. 1968. № 2. С. 3-10.

[7] Бубенчиков А.М., Ливаев Р.З. Некоторые автомодельные задачи магнитной гидродинамики // Вест. ТГУ. Бюл. оперативной науч. информации. 2001. № 4. С. 32-52.

[8] Бубенчиков А.М., Колесникова А.В. Стабилизированное течение электропроводящей жидкости в круглой трубе // Вычисл. технологии. 2005. Т 10, № 5. C. 5-12.

[9] Hunt J.C.R. Magnetohydrodynamic flow in rectangular ducts // Fluid Mech. 1965. Vol. 21, pt 4. P. 577-590.

[10] Gold R.R. Magnetohydrodynamic pipe flow. Pt. 1 // Fluid Mech. 1962. Vol. 13, pt 4. P. 505-512.

[11] Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. 152 с.

[12] Симуни Л.М. Численное решение задач теплообмена при неизотермическом движении вязкой жидкости в плоской трубе // Инж.-физ. журн. 1966. Т. 10, № 1. С. 89-91.

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

Поступила в редакцию 24 марта 2005 г., в переработанном виде — 14 июля 2005 г.