CC BY
48
С.О. Ширяева, А.И. Григорьев, П.В. Мокшеев
ОБ УСТОЙЧИВОСТИ ВРАЩАЮЩЕЙСЯ ЗАРЯЖЕННОЙ КАПЛИ
Ярославский государственный университет им. П.Г. Демидова, ул. Советская,14, г. Ярославль, 150000, Россия
Введение. Расчет равновесных форм и устойчивости вращающихся заряженных капель представляет интерес в связи с изучением процессов в грозовых облаках, воронках смерчей и в других заряженных жидко-капельных системах естественного и искусственного происхождения [1-4]. Исследование устойчивости капель в различных силовых полях необходимо для понимания физических закономерностей возникновения и существования смерчей и грозовых облаков. Ввиду многочисленных технических и технологических приложений, в которых вращающаяся заряженная капля встречается как основной физический объект, исследование устойчивости капель непосредственно связано с изучением временной эволюции нелинейно осциллирующих капель, механизмов реализации их неустойчивости и закономерностей распада при силовых воздействиях [4-7]. Так, Е. Тринч и Т. Ванг [7], исследовавшие возбуждаемые акустическим полем осцилляции большой амплитуды капель, подвешенных в акустическом подвесе, выяснили, что осцилляции большой амплитуды весьма трудно возбудить из-за появления на поверхности капли неосесимметричной бегущей волны, которая в конце концов приводила к вращению капли как целого. Такой же эффект проявлялся и в экспериментах Н. Якоби и др. [8] со свободно висящими в условиях невесомости каплями, осцилляции большой амплитуды которых также генерировались акустическим полем. Корректное истолкование подобных экспериментов должно отталкиваться от известной равновесной формы заряженной вращающейся капли.
1. Постановка задачи. Пусть капля идеальной несжимаемой идеально проводящей жидкости с коэффициентом поверхностного натяжения с, массовой плотностью р, имеющая заряд Q, вращает-
ся с постоянной угловой скоростью Q в вакууме. Будем исследовать устойчивость капли в первом порядке малости по отношению амплитуды равновесной деформации исходной сферической формы капли, вызванной ее вращением, к радиусу R. Все рассмотрение проведем в неинерциальной системе отсчета, вращающейся вместе с каплей, в сферической системе координат, начало координат которой совпадает с центром масс капли, а полярный угол 3 отсчитывается от положительного направления оси вращения. Всю описанную систему, состоящую из вращающейся капли вместе с окружающим ее электрическим полем собственного заряда, будем считать замкнутой.
Как показано в [9], равновесная форма вращающейся заряженной капли в линейном приближении по амплитуде равновесной деформации является сплюснутым вдоль оси вращения сфероидом. Уравнение свободной поверхности капли, возмущенной осесимметричным капиллярным волновым движением теплового происхождения (т.е. порождаемым тепловыми движениями молекул жидкости [10]), будет иметь вид:
r = r (3) + 4(3, t); r (3) =R(l - e2 ■ P2 (p)/3); p = cos(3);
e
2
p Q2 ■ R3
—------------^ 1;
4a ■(l - W) ;
W
Q2 .
16rcaR3 ’
14 << ^
где P2 (p)- полином Лежандра; функция 4 (3, t) описывает деформацию равновесной сфероидаль-
ной формы вращающейся заряженной капли, связанную с осцилляциями капли под действием сил поверхностного натяжения. Амплитуда деформации |4 (3, t)| имеет величину ~VkT7 a , где k - ко-
эффициент Больцмана; T - абсолютная температура. Безразмерный параметр W характеризует устойчивость изолированной сферической капли электропроводной жидкости по отношению к собственному заряду [10]: капля претерпевает неустойчивость при W > 1.
© Ширяева С.О., Григорьев А.И., Мокшеев П.В., Электронная обработка материалов, 2007, № 4, С.42-45.
42
Движение жидкости в капле, связанное с осцилляциями ее поверхности, будем полагать потенциальным с потенциалом поля скоростей y (r, t), а электрическое поле собственного заряда
капли - электростатическим с потенциалом Ф(r, t) . Причем потенциальность электрического поля имеет место как в окрестности невозмущенной сфероидальной поверхности капли (при E, (9, t) * 0), так и при наличии зависящей от времени деформации (при E, (9, t) ^ 0), поскольку характерная скорость изменения рельефа поверхности капли (~ Vy (r, t) много меньше скорости света в вакууме
с: Vy (r, t) << с. Тогда математическая формулировка задачи об исследовании устойчивости заряженной вращающейся капли будет иметь вид
Ay ( r, t); ДФ( r, t ) = 0; r ^ 0: y (r, t) ^ 0; r : Ф(r, t) ^ 0; r = r(9) + E(9,t): Ф(r,t) = Ф5 (r,t);
j V / 9 л dy 1 dy dE
— = -V-VE (9, t)« — — ——; p + pa+ pQ = pc, dt sv ' dr r2 d9 89’ ° Q с
где p = — (dy/ p ■ dt) - гидродинамическое давление внутри капли, связанное с движением
жидкости; pQ — давление на свободную поверхность капли центробежных сил, связанных с вращением капли:
r (9)+E (9,t) 1
pQ = j pQ2xdx = 2pQ2[r(9)+E(9,t)]2;
pQ = -(-VO)2/8n — давление электрического поля собственного заряда на свободную поверхность капли; pa = с ■ divSn — давление сил поверхностного натяжения, n — единичный вектор нормали к свободной поверхности капли; ФS (t) — постоянное вдоль поверхности капли значение электрического потенциала.
Кроме того, форма поверхности капли должна удовлетворять условиям неизменности объема и неподвижности центра масс:
j r2-dr-dp-d y = ——R3; V = [0 < r < r (9) + E(9, t), 0 <9< п, 0 <ф< 2п];
V 3
j er ■ r3 - dr -dp- d ф = 0,
V
которые позволяют найти ограничения снизу на спектр мод, принимающих участие в формировании формы свободной поверхности капли, условием n > 2 [10, 11]. При осцилляциях капли также должен оставаться неизменным ее полный электрический заряд:
ф (П *VO)r2 - dp - dф = — 4пQ; S = [r = r (9) + E(9, t), 0 < 9 < п,0 < ф < 2n].
S
2. Решение сформулированной задачи. Дальнейший анализ проведем в соответствии со стандартной процедурой, подробно описанной в [1, 4, 10, 12], сводящейся к линеаризации задачи в окрестности сферической формы капли, отысканию неизвестных функций E, (9, t), y (r, t), Ф(r, t) в
виде бесконечных рядов по полиномам Лежандра и получению дисперсионного уравнения задачи из условия совместности гидродинамических граничных условий на свободной поверхности жидкости.
В итоге решение задачи будем искать в виде
43
^(3,t)=Z An • Pn (ц) ■ exp(iOt);
n=2
да
W, 3 t) = X Cn ■rn ■ Pn (ц) ■ exp(iOt);
n=2
Ф(г, 3, t) = X Dn ■ r (n+1) ■ Pn (ц) ■ exp(iot).
n=2
При записи этих выражений было учтено, что гидродинамический потенциал поля скоростей должен быть ограничен в центре капли, а потенциал электростатического поля, связанного с осцилляциями поверхности заряженной капли, должен стремиться к нулю на бесконечности. В выписанных проектах решений A Cn, Dn - неизвестные коэффициенты.
Удовлетворяя гидродинамическим граничным условиям на свободной поверхности жидкости, найдем связь между коэффициентами An, Cn, Dn, выражая Cn и Dn через An, и выпишем выражения
для давленийp, pn , pQ и ра для сферической поверхности капли:
r = R : p = -р ду(д3) = ° ■R ■ р ■ exp(at)[ X n{An +
3(n+1)(n2+n-4) (n +1)(n2 - n - 3)
+ A
+ —{A +2 3 2(3+2n)(5+2n)
(2n - 1)(2n + 3)
+ A
3(-1)4(2~n)n(n - 1У 2(2n - 3)(2n -1)
■}];
2 ^ „2 .. чг “ / л / n + n—1 ,
Pq=- 3 R ■ р ■&■ ■ exp(iOt)[ X(An(-3^—^^ +
3 n=2
+ An+2 (
3 (n +1)( n + 2)
2 (2n + 5)(2n + 3)
= Q2 + Q2
+ An-2 (
(2n-1)(2n+3)' 3(-1)^4n n(n -1)(n + 2) 2 (2n - 1)(2n - 3)
) Pn (ц)];
Pq 8nR4 ' 4nR5 n=2
{ X (n-1) An ■Pn (ц) -
2 ,(n+1)( n+2)(n-7K n(n +1)( n - 4) .
-e2 ■ exp(iot)[ X (-----------A„+7 + —^^An +
n o 2(3+2n)(5+2n) +2 3(3 + 2n)(2n -1) ”
(-1)4(2-n)n(n - 1)(n - 5)
+
2(n - 3)(2n -1)
An-2) Pn (ц)]};
2a 4a ^ а да
Pa = 17 - 71 e2P2(^ + [ -2 X (n-1)(n+2) AnPn (ц)]+ R 3R R n=2
+7 e2 exp(iot )£ (Д,„ 3<" + +10> +
t! 2(2n + 3)(2n + 5)
3R
. n(n + 1)(n2 +n- 4) . 3(-1)4(2 n) n(n - 1)(n2 - 3n + 6)4n,
+A„—-------------------- + An 7—— ------------------------)P„ (ц)].
(2n + 3)(2n -1) n-2 2(2n - 3)(2n -1)
Наконец, подставив эти выражения в условие баланса давлений на свободной поверхности капли (динамическое граничное условие) и воспользовавшись ортогональностью полиномов Лежандра, получим дисперсионное уравнение задачи, которое в безразмерных переменных, где R = a = р = 1 (то есть за основные единицы выбраны радиус капли R, коэффициент поверхностного натяжения с и плотность жидкости р), имеет вид
о2 = n( n - 1)[(n + 2) - W ] +
w n4 +17 n2 + 21n-21 (n+2)( n4 -5n3 + 4n2+25n-9)1
+e [W-----------------------—--------------------].
3(3+2n)(2n-1) 3(3+2n)(2n-1)
Условием реализации неустойчивости капли является прохождение квадрата частоты о2 через ноль в область отрицательных значений. Тогда для самой частоты получим два комплексно сопряженных решения (о12 = ±iy), одно из которых будет соответствовать затуханию движений с
44
декрементом —у, а второе - их нарастанию с инкрементом +у. Приравнивая в полученном дисперсионном уравнении квадрат частоты нулю, найдем критическую для реализации неустойчивости капли связь между параметрами W, п, в2. В приближении в2 << 1 обсуждаемую зависимость можно записать в простом виде:
Wcr = (п + 2) + в
5п4 + 23п3 + 47п2 + 30п - 24 3п(4п3 - 7п + 3)
Несложно отметить, что для произвольных номеров мод критическое значение параметра W с ростом в2 увеличивается.
Сказанное означает, что наличие вращения у капли приводит к увеличению ее устойчивости по отношению к осесимметричным осцилляциям, что совпадает с ранее полученными результатами численных [13] и аналитических исследований [12,14], в которых сплюснутая сфероидальная форма капли не обеспечивалась физическим воздействием, как в настоящей работе, а задавалась виртуально.
Заключение. Равновесная форма заряженной вращающейся капли в широком диапазоне угловых скоростей и зарядов капли может быть аппроксимирована сфероидом, сплюснутым вдоль оси вращения.
Работа выполнена при поддержке гранта РФФИ №06-01-00066-а.
ЛИТЕРАТУРА
1. Григорьев А.И., Григорьев О.А. О влиянии вращения и электрического заряда на устойчивость сферической капсулы // Электронная обработка материалов. 1991. № 3. С. 41-44.
2. Григорьев А.И., Синкевич О.А. О природе электрических явлений в воронке смерча // ЖТФ. 1986. Т.56. Вып.10. С.1985-1987.
3. Grigor’ev A.I., Shiryaeva S.O. The possible physical mechanism of initiation and growth of lightning // Physica Scripta. 1996. V.54. P. 660-666.
4. Brown R.A., Scriven L.E. The shape and stability of rotating liquid drop // Proc. R. Soc., London, 1980. V.A371. P.331-357.
5. Natarajan R., Brown R.A. Third-order resonance effects and the nonlinear stability of drops oscillations // J. Fluid Mech. 1987. V.183. P.95-121.
6. Natarajan R., Brown R.A. The role of three-dimensional shapes in the break-up charged drops // Proc. R. Soc., London. 1987. V.A410. P.209-227.
7. Trinch E., Wang T.G. Large amplitude free and driven drop-shape oscillations: experimental observations // J. Fluid Mech. 1982. V.122. P.315-338.
8. Jakobi N., Croonquist A.P., Elleman D.D. Wang T.G. Acoustically induced oscillations and rotation of a large drop in Space // Proc. 2-nd Int. Colloq. on Drop and Bubbles. Pasadena: 1982. JPL Publication 82-7. P.31.
9. Ширяева С.О., Григорьев А.И., Мокшеев П.В. Равновесная форма заряженной капли, вращающейся вокруг своей оси симметрии // Электронная обработка материалов. 2006. № 4. С. 46-52.
10. Schweizer J.W., Hanson D.N. Stability limit of charged drops // J. Coll. Int. Sci. 1971. V.35. № 3. P.417-423.
11. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. М.: Наука, 1986. 733 с.
12. Григорьев А.И., Фирстов А.А. Критические условия неустойчивости заряженной капли, имеющей форму сплюснутого сфероида // Электронная обработка материалов. 1992. № 6. С. 20-23.
13. Basaran O.A., Scriven L.E. Axisymmetric shapes and stability of isolated charged drops // Phys.Fluids A. 1989. V.1. № 5. P.795-798.
14.Щукин С.И., Григорьев А.И. Устойчивость заряженной капли, имеющей форму трехосного эллипсоида // ЖТФ. 1998. Т.68. Вып.11. С.48-51.
Summary
Поступила 06.03.07
It is derived the disperse equation for rotating about its own axis of symmetry charged drop. It is shown that rotation increased the stability of a charged drop to its axisymmetric oscillations.
45
