ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2018, том 28, № 4, c. 127-134
МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ
УДК 539.19 + 532.5.032 © Б. П. Шарфарец
ПРИМЕНЕНИЕ СИСТЕМЫ УРАВНЕНИЙ ЭЛЕКТРОГИДРОДИНАМИКИ ДЛЯ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ НОВОГО СПОСОБА ЭЛЕКТРОАКУСТИЧЕСКОГО ПРЕОБРАЗОВАНИЯ
Проведен анализ системы уравнений электрогидродинамики (ЭГД) в контексте описания физических процессов, протекающих при возбуждении акустической энергии в электроакустическом преобразователе нового типа. Выявлено, что при реализации устройства преобразователя необходимо придерживаться ряда ограничений и рекомендаций. Во избежание возбуждения кратных частот необходимо прикладывать однородное электрическое поле, жидкая среда также должна быть однородной и без примесей. Рабочая жидкость в преобразователе должна обладать малой удельной проводимостью, в противном случае происходит ее перегрев, возникают другие паразитные явления. Для повышения уровня прикладываемого электрического поля необходимо повышать электрическую прочность рабочей жидкости. Использование математического моделирования при помощи системы ЭГД-уравнений позволит оптимизировать устройство электроакустического преобразователя нового типа.
Кл. сл.: электрогидродинамика, пондеромоторные силы, электродинамика, уравнение Навье—Стокса, уравнение сохранения энергии, уравнения непрерывности, электрическая прочность рабочей жидкости
ВВЕДЕНИЕ
В патенте [1] был предложен новый способ электроакустического преобразования, а в работах [2, 3] приводятся первичные результаты по его теоретическому обоснованию. Характерной особенностью этого преобразования является необходимость приложения к некой пористой среде мощного постоянного электрического поля наряду с умеренным переменным электрическим полем, несущим акустическую информацию. Изучение гидродинамических процессов в условиях доминирования либо магнитного, либо электрического полей занимаются соответственно пограничные дисциплины: магнитная гидродинамика (МГД) и электрогидродинамика (ЭГД) [4, с. 650].
В случае МГД при движении электропроводящей среды (газа, жидкости), находящейся в магнитном поле, в ней индуцируются электрические поля и токи, на которые действует магнитное поле и которые сами могут повлиять на магнитное поле. Таким образом, возникает сложная картина взаимодействия магнитного и гидродинамического явлений, и которая должна рассматриваться на основе совместных уравнений гидродинамики и электромагнитного поля.
В средах с очень малой электропроводностью и без приложенного извне большого магнитного
поля при скоростях жидкости v, много меньших скорости света |у| ■ с, определяющим во взаимодействии электромагнитного поля со средой является не магнитное, а электрическое поле. Эту область со своим кругом интересных задач и изучает ЭГД [4, с. 652, 653]. Электрическое поле описывается в ЭГД законами электростатики, а его воздействие на среду часто учитывают, ограничиваясь электрической частью силы Лоренца реЕ, где Е — вектор электрической напряженности стороннего электрического поля; ре — объемная плотность заряда в жидкости. Именно в таком виде была получена Максвеллом сила1 (ее еще называют пондеромоторной силой), действующая на несжимаемый диэлектрик в произвольном неоднородном электрическом поле при неизменной диэлектрической проницаемости [5, с. 138] (см., например, также [6, с. 309]). Такой же вид имела объемная сила в уравнении движения Навье—
1 Максвеллом было получено следующее выражение
для пондеромоторной силы (см. [5, с. 138], а также выЕ 2
ражение (1)): f = рЕ--Уе (обозначения см. в тек-
8л
сте). Недостающий здесь третий член из (1) был позднее добавлен Гельмгольцем. Приведенное выражение записано в системе СГСЭ.
Стокса в работах [2, 3], где рассматривалось движение жидкости под воздействием постоянного или переменного электрического поля при описании нового принципа электроакустического преобразования.
ПОСТАНОВКА ПРОБЛЕМЫ
Целью настоящей работы является описание условий, при которых справедливо представление объемной силы f в уравнении Навье—Стокса движения жидкости под воздействием электрического поля в укороченном виде f = рЕ, а также запись замкнутой системы уравнений ЭГД применительно к задачам, поставленным в работах [2, 3].
ПОНДЕРОМОТОРНЫЕ СИЛЫ
Известно, что электрическая часть силы Лоренца р Е является лишь одной из составляющих пондеромоторной силы, действующей на среду в произвольном неоднородном электрическом поле. Пондеромоторные силы, действующие в идеальном диэлектрике, описаны, например, в работах [5, § 32 и 34], [7, § 15], [8, глава I, § 1], [9, с. 86], [10, с. 128-130], [11, § 32] и др. Запишем выражение для этой силы в системе СИ [10, с. 130]:
f = р E - е E 2Ve + e V ( E 2р—Л
2
2
dp
(1)
f = p E + £°^ ^VE2. 2
(1а)
Выражение (1) упрощается в случае, если
известна зависимость е = е(р). Например, в [8, с. 24] и в [10, с. 130] используется формула Клаузиаса—Моссотти для однородных ди-
электриков
е-1
(е + 2 ) р
= C
(2)
где CK-M = const для данной среды (ее значение см. [4, с. 373]). Из (2) следует равенство [1°, с. 13°] и в [8, с. 24]:
Р
de (е- 1)(е+ 2 )
dp ~ 3
(3)
Здесь f — пондеромоторная сила; Е = Е ; е0 — электрическая постоянная; е — относительная диэлектрическая проницаемость среды; р — ее плотность. Здесь и далее для перехода в выражениях для удельной силы f в гауссовую систему необходимо второй и третий (нелинейные) слагаемые для f в СИ домножать
на коэффициент —1— (см., например, выражение 4яе0
для f в гауссовой системе в [5, с. 138] или [8, с. 24]).
Здесь отметим один предельный случай формулы (1), который касается слабо поляризуемых диэлектриков, для которых справедлива оценка е «1 (например, воздух) [11, с. 153]:
С учетом (2) выражение (1) принимает вид [9, с. 130]
f = рЕ - еЕ2Уе + еУ(Е2 (е - 1)(е + 2)) . (1б)
Формула Клаузиаса—Моссотти справедлива для неполярных диэлектриков,2 подчиняющихся определенному условию [4, с. 373, выражение (2)]. Для жидкостей, состоящих из неполярных молекул, формула (2) выполняется с высокой точностью. Для газообразных полярных3 диэлектриков и сильно разбавленных полярных жидкостей в неполярном растворителе формула (2) также верна [4, с. 374].
Обобщением формулы Клаузиаса—Моссотти на случай полярных молекул является формула Ланжевена—Дебая [4, с. 575], в которой принимается вместо постоянной Ск-М своя постоянная С
В области оптических частот видимого и УФ-диапазонов формула Клаузиаса—Моссотти (2) переходит в формулу Лоренца—Лоренца [4, с. 611]. В этом случае в формуле вида (2) фигурирует не относительная диэлектрическая проницаемость е , а квадрат показателя преломления оптически прозрачной среды п2, т. к. для оптического диапазона относительная магнитная проницаемость ц = 1 и е = п2. Формула (2) в этом случае преобразуется к виду (там же)
Из (1а) видно, что при е ^ 1 нелинейный член становится пренебрежимо малым даже при умеренном УЕ2 (е0 » 8.85 -1012).
Неполярные диэлектрики — диэлектрики с неполярными молекулами. В отсутствие внешнего электрического поля дипольный момент в молекулах отсутствует. Примеры: М2, Н2, С02, СН4. 3 Полярные диэлектрики — диэлектрики с полярными молекулами. В отсутствие внешнего электрического поля в молекулах присутствует дипольный момент. Примеры: Н2О, МН SO2, СО.
n2 -1 (n2 + 2) р
= C
(4)
Для газов формула Лоренца—Лоренца (4) преобразуется к виду [12, с. 7]
n -1 Р
= C,
где, например для воздуха, C = 0.22635 см3/г.
Полезно привести следующее соображение о частотной привязке приведенных формул. Формула Ланжевена—Дебая выполняется и для диэлектриков, находящихся в переменном электрическом поле, т. к. дипольные моменты молекул еще успевают ориентироваться вслед за полем вплоть до СВЧ-диапазона. В более высокочастотных микроволновом и ИК-диапазонах формула Ланжевена—Дебая переходит в формулу Клаузиу-са—Моссотти, а в видимом и УФ-диапазоне, когда уже и ионы не успевают смещаться вслед за полем, формула Ланжевена—Дебая переходит в формулу Лоренца—Лоренца (см. URL: http ://knowledge. su/l/lanzhevena-debaya-formula).
Анализ формул (1), (1б).
Как видно из выражений (1), (1б), объемная пондеромоторная сила f является нелинейной по амплитуде вектора электрической напряженности. Рассмотрим последствия этого на примере, актуальном для работ [2, 3]. Пусть вектор E представляет собой сумму векторов E = E0 + E1, где E0 = const, а E1 — коллинеарный ему вектор, совершающий гармонические колебания частоты со : E
E1 =—0 E1 cos cot (временная зависимость принята
E0
не в экспоненциальной форме ввиду нелинейности задачи). Здесь E0 и E1 — модули (амплитуды) соответствующих векторов. Тогда имеем для модуля вектора E :
E = E0 + E cos cot.
(5)
Отсюда
E2 = (E0 + E1 cos cot)2 =
E 2
= E02 + 2E0E1 cos cot + (1 + cos2ct) =
E 2 + E + 2
2
E2
+ 2E0E1 cos cot + -^cos2ct.
(6)
Из (6) видно, что если не равны нулю вторые и третьи слагаемые в (1), (1б), содержащие Е2, то кроме постоянной составляющей реЕ0 и гармони-
ческой составляющей peE1 cos cot частоты со возникают дополнительные постоянные составляющие и составляющие кратной частоты 2с . И если вторичные постоянные составляющие не искажают исходный информативный сигнал частоты с, то кратная гармоника является паразитной. Впрочем, такова специфика нелинейных процессов. В случае, если второе и третье слагаемые в (1), (1б) равны нулю или пренебрежимо малы по сравнению с линейным членом peE, то паразитные составляющие отсутствуют, либо пренебрежимо малы.
Второе и третье слагаемые в (1), (1б) не равны нулю, если стоящие под знаком градиента величины неоднородны по пространству. Второй член неоднороден по пространству, если непостоянна по пространству диэлектрическая проницаемость s . Третий член отличен от нуля либо когда неоднородно поле E , либо когда непостоянна по пространству диэлектрическая проницаемость s , либо когда они неоднородны одновременно. Что касается выражения (1а) для пондеромоторной силы в слабополяризуемых диэлектриках, то при малой неоднородности поля E она также вырождается к случаю f = peE .
На неоднородность величины s влияет пространственное возмущение плотности среды p. Возмутителями могут быть, в частности, вариации температуры и акустического давления. В первом случае причиной может быть, например, неоднородное температурное поле T . Пусть температурное поле является суммой T = T0 + T1. Ему соответствует плотность среды p = p0 + p1, где индексу 0 соответствует невозмущенное состояние поля, а индексу 1 — его возмущенное состояние. Тогда, согласно [13, с. 306], имеем для возмущения плотности p1:
Фо
Р =1^0 I T =-Ро^.
1 Г Фо
Здесь 3 =--1 —0 I — температурный коэффи-
Ро I ^ ) р циент расширения жидкости (газа).
Далее будет рассматриваться ситуация, когда электрическая проводимость среды отлична от нуля. Следствием этого будет отличие от нуля и тока проводимости в среде. Это в свою очередь вызывает в общем случае неоднородный нагрев среды, что и приводит к неоднородности ее плотности (см., например, в работе [14] процесс на примере кругового капилляра), а следовательно, и диэлектрической проницаемости е .
Другим возмутителем плотности является аку-
стическое давление p1 (общее давление p вновь обозначаем в виде суммы p = p° + p1, где p° — равновесное давление в среде). Хорошо известна связь между возмущением давления p1 и возмущением плотности р1 (акустический процесс обычно рассматривается как адиабатический, т. е. 5 = const, где 5 — удельная энтропия) [13, с. 351]
Р1 =
dp dP
Р.
° /
VxE = -dB,
dt
V-D = p
dD
Vx H = j +—, V-B = °, dt
дополняемые материальными уравнениями D = ee°E, B = Mi°H, j = ctE. В гауссовой системе:
(7а)
Система уравнений электрогидродинамики.
Электрогидродинамика (ЭГД) представляет собой уже устоявшийся междисциплинарный раздел науки, в котором тесно взаимодействуют такие направления, как гидродинамика, электродинамика в основном в своей электростатической части, электрохимия и т. д. [15]. Объектом исследований в ЭГД являются жидкости (газы), обладающие слабой проводимостью 7 (< < 10-7 ^ 10-12 Ом-1см-1). Библиография ЭГД весьма обширна, здесь же, кроме упомянутых выше монографии [8] и обзора [15], в котором подробно представлена современная библиография по теме, сошлемся только на ключевые работы [16-24]. Приведем далее два коротких, но емких определения ЭГД-системы, следуя [16, с. 15, 16]:
- ЭГД-система — система, в которой действующее электрическое поле много больше индуцированного электрического поля и действующее магнитное поле порядка индуцированного магнитного поля;
- ЭГД-системой называется такая система, в которой определяющими являются электростатические и гидродинамические силы.
Основным инструментом ЭГД является замкнутая система уравнений, приведенная с различными вариациями в упомянутых выше работах об ЭГД-системах. В систему, как правило, входит уравнение сохранения импульса (Навье—Стокса), гидродинамическое и электродинамическое уравнения непрерывности, уравнение сохранения энергии и ряд уравнений электродинамики. Коротко остановимся на них и приведем накладываемые ЭГД ограничения на физические процессы, подпадающие под его изучение.
Электродинамические уравнения.
В качестве электродинамической основы берется система уравнений Максвелла для подвижной среды. Вначале, однако, выпишем исходную систему уравнений Максвелла для неподвижной среды в обеих системах: в СИ и гауссовой системе ([25] и др.).
В системе СИ:
VxE = -1dB, с dt
^ 4я. 1 dD
Vx H = — j +--,
с с dt
V- D = 4жре, V-B = °,
(8)
дополняемые материальными уравнениями
D = еE, B = ¡Ж, \ = стE. (8а)
Здесь, помимо уже принятых выше, введены следующие обозначения: D — вектор электрической индукции; H — вектор магнитной напряженности; B — вектор магнитной индукции; ] — вектор плотности тока проводимости4; 7 — удельная проводимость (например, в системе СИ [См/м], или [Ом-1м-1]); с — скорость света в ва-
кууме; ¡л°
магнитная постоянная, равная
¡и0 = 4^-10 7 Гн/м; и — относительная магнитная проницаемость.
Следуя, например, [20, § 1], выпишем адаптированную к ЭГД систему уравнений Максвелла и ограничения, накладываемые на исходные уравнения Максвелла при получении уравнений ЭГД.
Если поле меняется во времени, то, согласно (7), (8), в среде будет индуцироваться магнитное поле H за счет токов проводимости ] = стE и то-
дБ
ков смещения ]см = — .
д1
Ограничение на величину магнитного поля.
В электрогидродинамике обычно полагают, что магнитное поле равно нулю H = B = 0, что, согласно уравнениям в (7), (8), где фигурирует Ух E , означает, что Ух E = 0 , что говорит о потенциальности вектора E, т. е. потенциальности электрического поля. В [20, § 1] установлены границы такого приближения.
Мерой пренебрежения магнитными эффектами может служить отношение плотностей энергий магнитного и электрического полей
Плотность тока — векторная физическая величина, имеющая смысл силы электрического тока, протекающего через элемент поверхности единичной площади.
5
uu H 2
ss0 E2
«1.
(9)
Ограничение на ток проводимости.
. dD
Согласно уравнению в (7) УхН = ]н--, на-
дt
пряженность магнитного поля, индуцированного током проводимости, по порядку величины составляет
H ~ а El.
(10)
где I — характерный пространственный масштаб, где существенно изменяется Н. Подстановка (10) в (9) дает ограничение на величину удельной проводимости а, при которой можно пренебречь магнитным полем
1
а «l
ss,
\
1/2
UU
(11)
-о у
При принятых в [20] значениях (l = 10 2 м, s = № = 1) знак равенства в (11) достигается при а « 0.3 Ом-1 м-1. Для примера, удельная проводимость а составляет (информация по адресу в Википедии "Удельная электропроводность" https://slovar.wikireading.ru/2618966): для морской воды 3 Ом-1м-1, для воды дистиллированной 10-4 Ом-1м-1, для воздуха у поверхности земли а = 2.2 -10"18 Ом-1м-1. Для такой жидкости, как трансформаторное масло, а = 10"13 -^10"12 Ом_1м_1 (ссылка по адресу
https://studopedia.ru/10_13 0433_kompensatsionnaya-obmotka.html).
Как видно из этих данных, использование, например, морской воды в качестве рабочей жидкости чревато возникновением больших токов проводимости и, как следствие, перегревом электроакустического преобразователя.
n dD
Ограничения по токам смещения ]см = —.
. dt
Следуют из второго уравнения Максвелла в (7). Как видно из этого уравнения, ток смещения индуцирует магнитное поле напряженностью H по-ss El
рядка [20, с. 13] H--0--ss0Elf, где t0 = 1/ f
t0
— период колебаний частоты f. Подстановка полученной оценки в (9) дает ограничение по частоте колебаний [20, с. 13]
показатель преломления электромагнитной волны в среде. При принятых в [20] значениях (п «3, I = 10"2) получена оценка f ■ 10 МГц. Например, для дистиллированной воды (п «1.33) ограничение (12) почти вдвое смягчится — f ■ 22.5 МГц, а для воздуха (п«1) смягчится втрое — f « 30 МГц.
Отметим, что частота ультразвуковых колебаний, применяемых в промышленности и биологии, лежит в диапазоне от нескольких десятков кГц до единиц МГц (ссылка по адресу: https://ru.wikipedia.org/wiki/%D0%A3%D0%BB%D 1 %8С%01 %82%01 %80%®0%В0%Л0%В7%®0% В2%®1%83%®0%ВА). Это означает отсутствие рассматриваемого ограничения для УЗ-излучателей, построенных на принципе, описанном в [13].
Ограничение на скорость движения жидкости.
Еще одно ограничение, допускаемое в ЭГД, касается соотношения между скоростью движения жидкой среды |у| и скоростью света с . Уравнения ЭГД выводятся, исходя из разумного допущения
v « c
(13)
Очевидно, что это неравенство, безусловно, выполняется в рассматриваемом здесь случае.
Уравнения электрогидродинамики.
Получаются синтетическим объединением системы уравнений гидродинамики и системы уравнений Максвелла для движущихся сред с учетом принятых выше ограничений и замечаний относительно системы уравнений Максвелла. Отметим, что материальные уравнения (7а), (8а) в движущейся среде видоизменяются и определяются общими уравнениями Минковского для движущихся сред [7, с. 364] (см. также [16, с. 15]). Например, (8а) имеет вид [7, с. 364]
D = sE + ^^ v х H, c
B = jH + ^^ E х v.
(14)
Системы ЭГД-уравнений выписаны в упомянутых выше работах [8, 15-24], однако здесь выпишем их версию для несжимаемой жидкости, представленную в работах [20, 23, 24], которая имеет вид [20, с. 14]5:
г С f « — , nl
где c — скорость света в вакууме; n
(12) ; n = yJsjU —
Аналогичная система ЭГД-уравнений применительно к сжимаемой жидкости выписана в работе [16, § 1.1, 1.2].
c
+ ( V V V 1 = + + Рg + Г
дТ + V 'УГ = х^Т + ^^ , дт рсР
V' V = о, Е =
д(еЕ)
] = стЕ + р V + £0 +Ух( Р х V):
дт
V' ] = 0,
f = рЕ - £ Е2Уе+ £ VI Е2 р —
2
2
др
(15)
(16)
(17)
(18)
(19)
(20)
(21)
УТ
Здесь г) — динамическая вязкость жидкости; g — вектор ускорения силы тяжести; f — объемная внешняя сила; Т — абсолютная температура в жидкости; х — коэффициент температуропроводности жидкости; р — скалярный потенциал поля Е; Р — вектор поляризации. Как видно из (21), объемная сила f принимается равной пондеромоторной силе, действующей в идеальном
£ де ^
диэлектрике (1).6 Слагаемое — V Е2р— в (21)
2 I дР)т
называется электрострикционной силой. Вследствие потенциальности этой силы она может быть включена в градиент давления. Это означает, что электрострикция приводит к понижению давления
£о Г7 Гг72 д£ ^ на величину — V Е р —
2
др
Т
Замечания по приведенным уравнениям.
В уравнении теплопроводности (16) справа стоит тепловой источник, связанный с мощностью выделения джоулева тепла в единице объема ^ = ]' Е при плотности тока ] = стЕ, т. е. по умолчанию принимается, что остальные слагаемые плотности тока (19) пренебрежимо малы. Впрочем, в [20] по этому поводу даны комментарии. Однако этот вопрос в [20] не затронут для случая электроосмотических процессов, когда слагаемое р V явно отлично от нуля в двойном электрическом слое. Поэтому этот случай требует отдельного специального рассмотрения (техника подобных
6 Такое же выражение присутствует для силы Г в большинстве цитированных работ по ЭГД при наличии ряда отличий. Например, в обзоре [15] вместо силы Г присутствует тензор с такой же структурой слагаемых (но только тензорных), что и в (21).
выкладок приведена, например, в [14]).
В обзоре [15] в ЭГД-системе присутствует электродинамическое уравнение непрерывности др
—- + V' ] = 0, которого нет в системе (21), что не
дт
принципиально при обеспечении замкнутости ЭГД-системы.
Пробивная напряженность (электрическая прочность) диэлектриков (по материалам [24]).
При напряженности электрического поля, превосходящей предел электрической прочности диэлектрика, наступает пробой. Пробой представляет собой процесс разрушения диэлектрика, в результате чего диэлектрик теряет электроизоляционные свойства в месте пробоя. Величину напряжения, при котором происходит пробой диэлектрика, называют пробивным напряжением ипр, а
соответствующее значение напряженности электрического поля называется электрической прочностью диэлектрика Епр .
Для равномерного электрического поля электрическая прочность (пробивная напряженность) диэлектрика Епр определяется по формуле
и
Е
пр d '
где d — толщина диэлектрика в месте пробоя.
Пробой воздуха, как и прочих газов, является следствием ударной ионизации. Содержащиеся в газе положительные и отрицательные ионы и электроны, находящиеся в хаотическом тепловом движении, при наложении поля получают некоторую добавочную скорость и начинают перемещаться в соответствующем направлении. Если эта энергия достаточно велика, будет происходить ионизация молекул. В результате ионизации создаются положительный ион и электрон. Этот электрон может, в свою очередь, ионизировать, создавать новые электроны и положительные ионы. Возникшие положительные ионы наряду с электронами также участвуют в ионизации. В результате лавинного процесса образования свободных зарядов резко возрастает сила тока и наступает пробой газа. Пробивная напряженность воздуха в нормальных условиях (давление 760 мм рт. ст. и температура 20 °С) составляет 3.2 кВ/мм. Пробивная напряженность трансформаторного масла 50-180 кВ/мм.
Заметим также, что электрическая прочность диэлектрика зависит от окружающего давления. Существуют также различные присадки для повышения электрической прочности диэлектриков.
При проектировании электроакустического излучателя необходимо иметь в виду величину электрической прочности диэлектриков.
ВЫВОДЫ
В результате проведенного анализа системы уравнений ЭГД в контексте описания физических процессов, протекающих при возбуждении акустической энергии в электроакустическом преобразователе нового типа [1-3], выявлено, что при реализации устройства преобразователя необходимо придерживаться ряда ограничений и рекомендаций:
- для избежания возбуждения кратных частот необходимо прикладывать однородное электрическое поле; жидкая среда также должна быть однородной и без примесей;
- рабочая жидкость в преобразователе должна обладать малой удельной проводимостью, в противном случае происходит ее перегрев, возникают другие паразитные явления;
- для повышения уровня прикладываемого электрического поля необходимо учитывать рекомендации по повышению электрической прочности рабочей жидкости;
- использование математического моделирования при помощи системы ЭГД-уравнений позволяет оптимизировать устройство электроакустического преобразователя нового типа.
СПИСОК ЛИТЕРАТУРЫ
1. Shishov S. V., Andrianov S.A., Dmitriev S.P., Ruchkin D. V. Method of converting electric signals into acoustics oscillations and an electric gas-kinetic transducer. Patent US N 8,085,957,B2. Dec. 27, 2011.
2. Сергеев В.А., Шарфарец Б.П. Об одном новом методе электро-акустического преобразования. Теория, основанная на электрокинетических явлениях. Ч. I. Гидродинамический аспект // Научное приборостроение. 2018. Т. 28, № 2. С. 25-35.
URL: http://213. 170.69.26/mag/2018/abst2.php#abst4.
3. Сергеев В.А., Шарфарец Б.П. Об одном новом методе электроакустического преобразования. Теория, основанная на электрокинетических явлениях. Ч. II. Акустический аспект // Научное приборостроение. 2018. Т. 28, № 2. С. 36-44.
URL: http://213.170.69.26/mag/2018/abst2.php#abst5.
4. Физическая энциклопедия. Т. 2 / Под ред. А.М. Прохорова. М.: Сов.энциклопедия, 1990. 703 с.
5. Сивухин Д.В. Общий курс физики. Т. 3. Электричество. М.: Физматлит, 2004. 656 с.
6. Ньюмен Дж. Электрохимические системы. М.: Мир, 1977. 464 с.
7. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика.
Т. VIII. Электродинамика сплошных сред. М.: Наука, 1982. 621 с.
8. Остроумов Г.А. Взаимодействие электрических и гидродинамических полей: физические основы электрогидродинамики. М.: Наука, 1979. 320 с.
9. Физическая энциклопедия. Т. 4 / Под ред. А.М. Прохорова. М.: Сов. Энциклопедия, 1994. 703 с.
10. Стрэттон Дж.А. Теория электромагнетизма. М.: ГИТТЛ, 1948. 540 с.
11. Тамм И.Е. Основы теории электричества. М.: Наука, 1976. 616 с.
12. Герасимов С.И., Ерофеев В.И., Сперанский А.В., То-тышев К.В. Алгоритмы определения параметров схем теневого фотографирования. Саровский физико-технический институт НИЯУ МИФИ, 2013. 95 с.
13. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика. М.: Наука, 1986. 736 с.
14. Шарфарец Б.П., Курочкин В.Е. К вопросу об определении стационарного температурного поля в капилляре при прохождении в нем электрического тока, а также об изменении полей концентрации примесей в этом температурном поле // Научное приборостроение. 2015. Т. 25, № 2. С. 53-60.
URL: http://213. 170.69.26/mag/2015/abst2.php#abst6.
15. Жакин А.И. Электрогидродинамика. Обзор // УФН. 2012. Т. 182, № 5. С. 495-520.
16. Рубашов И.Б., Бортников Ю.С. Эектрогазодинамика. М.: Атомиздат, 1971. 168 с.
17. Мелчер Дж., Тейлор Дж. Электрогидродинамика: обзор роли межфазных касательных напряжений // Механика: Сб. переводов. 1971. № 5. С. 66-99.
18. Мелчер Дж.Р. Электрогидродинамика // Магнитная гидродинамика. 1974. Т. 2. С. 3-30.
19. Гогосов В.В., Полянский В.А. Электрогидродинамика // Итоги науки и техники: Механика жидкости и газа. М.: ВИНИТИ, 1976. Т. 10. С. 5-85.
20. Болога М.Н., Гросу Ф.П., Кожухарь И.А. Электроконвекция и теплообмен. Кишинев: Штиинца, 1977. 320 с.
21. Стишков Ю.К., Остапенко А.А. Электрогидродинамические течения в жидких диэлектриках. Л.: ЛГУ, 1989. 176 с.
22. Electrohydrodynamics / Ed. by Antonio Castellanos. Vien: Springer-Verlag, 1998. 362 p.
23. Остроумов Г.А. К вопросу о гидродинамике электрических разрядов // ЖТФ. 1954. Т. 24, № 10. С. 1915.
24. Остроумов Г.А. Электрическая конвекция (обзор) // ИФЖ. 1966. Т. 10, № 5. C. 683-695.
25. Яворский Б.М., Детлаф А.А. Справочник по физике для инженеров и студентов. М.: Наука, 1971. 939 с.
26. Богородицкий Н.П., Пасынков В.В., Тареев Б.М. Электротехнические материалы. Учебник для вузов. Л.: Энергоатомиздат, 1985. 304 с.
ISSN0868-5886 NAUCHNOE PRIBOROSTROENIE, 2018, Vol. 28, No. 4, pp. 127-134
Институт аналитического приборостроения РАН, Контакты: Шарфарец Борис Пинкусович, Санкт-Петербург [email protected]
Материал поступил в редакцию 13.08.2018
APPLICATION OF THE SYSTEM OF ELECTROHYDRODYNAMICS EQUATIONS FOR MATHEMATICAL MODELING OF A NEW METHOD OF ELECTRO-ACOUSTIC TRANSFORMATION
B. P. Sharfarets
Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia
An analysis of the system of EHD equations in the context of describing the physical processes occurring during the excitation of acoustic energy in an electroacoustic transducer of a new type is carried out. It is revealed that when implementing a converter device, it is necessary to adhere to a number of limitations and recommendations. To avoid excitation of multiple frequencies, it is necessary to apply a uniform electric field, the liquid medium must also be homogeneous and without impurities. The working fluid in the converter must have a low specific conductivity, otherwise it overheats, other parasitic phenomena occur. To increase the level of the applied electric field, it is necessary to increase the electrical strength of the working fluid. The use of mathematical modeling with the help of the system of EHD equations will allow us to optimize the device of an electroacoustic transducer of a new type.
Keywords: electrohydrodynamics, ponderomotive forces, electrodynamics, the Navier-Stokes equation, the equation of conservation of energy, equations of continuity, electric strength of working fluid
REFERENCES
1. Shishov S.V., Andrianov S.A., Dmitriev S.P., Ruch-kin D.V. Method of converting electric signals into acoustics oscillations and an electric gas-kinetic transducer. United States Patent N US 8,085,957,B2. Dec. 27, 2011.
2. Sergeev V.A., Sharfarets B.P. [About one new method of electroacoustic transformation. A theory based on electro-kinetic phenomena. Part I. The hydrodynamic aspect]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2018, vol. 28, no. 2, pp. 25-35. Doi: 10.18358/np-28-2-i2535. (In Russ.).
3. Sergeev V.A., Sharfarets B.P. [About one new method of electroacoustic transformation. A theory based on electro-kinetic phenomena. Part II. The acoustic aspect]. Nauch-noe Priborostroenie [Scientific Instrumentation], 2018, vol. 28, no. 2, pp. 36-44. Doi: 10.18358/np-28-2-i3644. (In Russ.).
4. Prohorov A.M., ed. Fizicheskaya enziklopediya [Physical encyclopedia]. Vol. 2. Moscow, Soviet encyclopedia Publ., 1990. 703 p. (In Russ.).
5. Sivuhin D.V. Obshchij kurs fiziki. T. 3. Elektrichestvo [General course of physics. Vol. 3. Electricity]. Moscow, Fizmatlit Publ., 2004. 656 p. (In Russ.).
6. Newman J. Elektrochimicheskie sistemy [Electrochemical
Systems]. Moscow, Mir Publ., 1977. 464 p. (In Russ.).
7. Landau L.D., Lifshiz E.M. Teoreticheskaya fizika. T. VIII. Elektrodinamika sploshnyh sred [Theoretical physics. Vol. VIII. Electrodynamics of continuous environments]. Moscow, Nauka Publ., 1982. 621 p. (In Russ.).
8. Ostroumov G.A. Vzaimodejstvie ehlektricheskih i gidro-dinamicheskih polej: fizicheskie osnovy ehlektrogidrodi-namiki [Interaction of electric and hydrodynamic fields: physical fundamentals of electrohydrodynamics]. Moscow, Nauka Publ., 1979. 320 p. (In Russ.).
9. Prohorov A.M. ed. Fizicheskaya enziklopediya [Physical encyclopedia]. Vol. 4. Moscow, Soviet encyclopedia Publ., 1994. 703 p. (In Russ.).
10. Stretton J.A. Teoriya elektromagnetizma [Theory of elec-tromagnetism]. Moscow, GITTL Publ., 1948. 540 p. (In Russ.).
11. Tamm I.E. Osnovy teorii elektrichestva [Bases of the theory of electricity]. Moscow, Nauka Publ., 1976. 616 p. (In Russ.).
12. Gerasimov S.I., Erofeev V.I., Speranskij A.V., Toty-shev K.V. Algoritmy opredeleniya parametrov skhem te-nevogo fotografirovaniya [Algorithms of determination of parameters of schemes of shadow photography]. Sarov physics and technology institute of NIYaU MEPhI, 2013. 95 p. (In Russ.).
13. Landau L.D., Lifshiz E.M. Teoreticheskaya fizika. T. 6. Gidrodinamika [Theoretical physics. Vol. 6. Hydrodynamics]. Moscow, Nauka Publ., 1988. 736 p. (In Russ.).
14. Sharfarets B.P., Kurochkin V.E. [To the question about the definition of a stationary temperature field in the capillary during the passage of electric current and the change in water concentration of impurities in the temperature field]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2015, vol. 25, no. 2, pp. 53-60. (In Russ.).
URL: http://213. 170.69.26/mag/2015/abst2.php#abst6.
15. Zhakin A.I. [Electrohydrodynamics. Review]. Uspekhi Fi-zicheskikh Nauk [Achievements of physical sciences], 2012, vol. 182, no 5, pp. 495-520.
Doi: 10.3367/UFNr.0182.201205b.0495. (In Russ.).
16. Rubashov I.B., Bortnikov Yu.S. Elektrogazodinamika [Electrogas dynamics]. Moscow, Atomizdat Publ., 1971. 168 p. (In Russ.).
17. Melcher G., Tejlor G. [Electrohydrodynamics: review of a role of interphase tangent tension]. Mekhanika: Sbornik perevodov [Mechanics: Collection of translations], 1971, no. 5, pp. 66-99. (In Russ.).
18. Melcher G.R. [Electrohydrodynamics]. Magnitnaya gidrodinamika [Magnetohydrodynamics], 1974, vol. 2, pp. 3-30. (In Russ.).
19. Gogosov V.V., Polyanskij V.A. [Electrohydrodynamics]. Itogi nauki i tekhniki: Mekhanika Zhidkosti i Gaza [Results of science and technology: Fluid Dynamics], 1976,
Contacts: Sharfarets Boris Pinkusovich, [email protected]
vol. 10, pp. 5-85. (In Russ.).
20. Bologa M.N., Grosu F.P., Kozhuhar' I.A. Elektrokonvek-ciya i teploobmen [Electroconvection and heat exchange]. Chisinau, Shtiintsa, 1977. 320 p. (In Russ.).
21. Stishkov Yu.K., Ostapenko A.A. Elektro-gidrodinamicheskie techeniya v zhidkih diehlektrikah [Electro-hydrodynamic currents in liquid dielectrics]. Leningrad, LGU, 1989. 176 p. (In Russ.).
22. Castellanos Antonio, ed. Electrohydrodynamics. Vien: Springer-Verlag, 1998. 362 p.
23. Ostroumov G.A. [To a question of hydrodynamics of electric discharges]. Zhurnal tekhnicheskoj fiziki [Journal of technical physics], 1954, vol. 24, no. 10, pp. 1915. (In Russ.).
24. Ostroumov G.A. [Electric convection (review)]. Inzhener-no-Fizicheskii Zhurnal [Journal of Engineering Physics and Thermophysics], 1966, vol. 10, no. 5, pp. 683-695. (In Russ.).
25. Yavorskij B.M., Detlaf A.A. Spravochnik po fizike dlya inzhenerov i studentov [The reference book on physics for engineers and students]. Moscow, Nauka Publ., 1971. 939 p. (In Russ.).
26. Bogorodickij N.P., Pasynkov V.V., Tareev B.M. Elektrotekhnicheskie materialy. Uchebnik dlya vuzov [Electrotechnical materials. The textbook for higher education institutions]. Leningrad, Energoatomizdat Publ., 1985. 304 p. (In Russ.).
Article received in edition 13.08.2018
HAYHHOE nPHEOPOCTPOEHHE, 2018, tom 28, № 4