Effect of Charge Relaxation Effect on Electromagnetic Radiation Intensity of Oscillating Viscous Liquid Drop
- Autores: Grigoryev A.I.1, Kolbneva N.Y.2, Shiryaeva S.O.2
-
Afiliações:
- Ishlinsky Institute for Problems in Mechanics of the RAS
- Demidov Yaroslavl State University
- Edição: Volume 88, Nº 5 (2024)
- Páginas: 692-721
- Seção: Articles
- URL: https://journals.rcsi.science/0032-8235/article/view/280963
- DOI: https://doi.org/10.31857/S0032823524050046
- EDN: https://elibrary.ru/JPXANL
- ID: 280963
Citar
Texto integral
Resumo
Theoretical asymptotic methods have shown that the electric charge relaxation effect affects the physical characteristics of the electromagnetic radiation of the oscillating charged droplet. Analytical expressions are obtained for frequencies, decrements of attenuation of capillary oscillations of droplets due to viscous attenuation and energy losses on radiation. It has been shown that the frequencies of electromagnetic radiation of cloud droplets, realized in the ranges of hundreds of kilohertz and megahertz units, decrease with an increase in the radius and charge parameter of the emitting droplet, as well as attenuation decrements associated with radiation. Intensity of emission of electromagnetic waves decreases with decrease of electrical conductivity and mobility of charges in liquid.
Palavras-chave
Texto integral
- Введение. Впервые проблема электромагнитного излучения от осциллирующей заряженной капли в облаке была рассмотрена в [1] и затем развита в [2–5]. В [5], в частности было указано на ошибку в расчетах [1] (неправильно была определена асимптотика функции Ханкеля), и эта ошибка перекочевала в [2–4]. В [6] было указано на необходимость учета кривизны заряженной поверхности жидкости при записи уравнения баланса поверхностной плотности заряда в жидкости с конечной электропроводностью. Все это сделало актуальным повторение расчетов [3–4].
Сама задача о расчете интенсивности радиоизлучения (конвективного) грозового облака актуальна в связи с широким спектром разнообразных приложений в геофизике, гидрометеорологии, технике и химической технологии. Большое количество работ посвящено теоретическому и экспериментальному изучению капиллярных осцилляций и электростатической неустойчивости жидкой капли, заряженной собственным или индуцированным внешним электрическим полем зарядами (см. [1–5] и цитируемую там литературу).
В частности, исследование излучения осциллирующей каплей электромагнитных волн привлекает внимание из-за проблем радиолокационного зондирования облаков и туманов [7–10], проблем радиопомех от огней Св. Эльма, появляющихся на обшивке самолетов, пролетающих через облака [11].
Что касается эффекта релаксации заряда, то он возможен только при отличной от нуля вязкости жидкости. Это связано тем, что на свободной поверхности электрические напряжения, возникающие из-за релаксационных движений, могут быть скомпенсированы только вязкими напряжениями.
- Постановка задачи. Рассмотрим уединенную сферическую каплю радиуса R вязкой несжимаемой жидкости, характеризующейся массовой плотностью ρ, коэффи-
циентами поверхностного натяжения σ и кинематической вязкости ν, удельной проводимостью γ1 ≡ γ, диэлектрической проницаемостью ε1. Примем, что капля несет заряд и находится в вакууме с удельной проводимостью γ2 = 0 и диэлектрической проницаемостью ε2 = 1. Будем считать проводимость внутренней среды достаточной, чтобы обеспечить отсутствие объемного заряда в капле. Индекс j, принимает значения 1 и 2.
В сферической системе координат (r, θ, φ) с началом в центре масс капли уравнение ее поверхности, совершающей осесимметричные капиллярные осцилляции, в произвольный момент времени представим в виде:
, (2.1)
где ξ (θ, t) – малое возмущение равновесной сферической формы, генерируемое капиллярным волновым движением тепловой природы. При наличии внешних неконтролируемых силовых воздействий (коагуляция, дробление, столкновение, трение о воздух и т.д.) тепловая амплитуда |ξ| может возрастать до десятков процентов от радиуса капли [12, 13].
Отношение |ξ| / R, имеющее смысл безразмерной амплитуды капиллярного волнового возмущения, примем в качестве малого параметра задачи: ε ≡ |ξ| / R << 1. Так как движение жидкости вызывается осцилляциями капли, будем считать, что поле скоростей течения внутренней среды имеет тот же порядок малости, что и ξ (θ, t): .
В нижеследующих рассуждениях учтем, что собственный заряд капли равномерно распределен на равновесной поверхности с постоянной поверхностной плотностью κ(0). В виду конечности электропроводности жидкости перенос электрического заряда при волновом искажении формы капли приводит к эффекту релаксации заряда.
Выпишем систему электрогидродинамических уравнений, описывающих капиллярные осцилляции заряженной вязкой капли конечной проводимости и возникающего при этом электромагнитного излучения.
Система состоит из уравнения Навье–Стокса:
, (2.2)
уравнения непрерывности:
,
уравнений Максвелла:
условий ограниченности:
, (2.3)
и граничных условий на свободной поверхности капли F (r, θ, t) ≡ r − R − ξ (θ, t) = 0:
кинематического:
,
динамического для касательной компоненты тензора напряжений:
и динамического для нормальной компоненты тензора напряжений:
Чтобы учесть эффект релаксации электрического заряда, введем дополнительные граничные условия на поверхности капли F (r, θ, t) ≡ r − R − ξ (θ, t) = 0:
скачка нормальной компоненты электрической индукции:
, (2.4)
непрерывности касательных компонент напряженности электрического поля:
,
и баланса поверхностной плотности заряда [6]:
(2.5)
Кроме того, потребуем выполнения естественных интегральных условий сохранения объема капли при ее осцилляциях, неизменности положения центра масс:
а также постоянства суммарного заряда:
(2.6)
В приведенной математической формулировке задачи нижний индекс 1 относится к капле, индекс 2 – к внешней среде, μj – магнитная проницаемость j-й среды (ниже принимается μ1 = μ2 = 1), κ (r, θ, t) – поверхностная плотность электрического заряда на капле, c – скорость распространения электромагнитных волн в вакууме, b – поверхностная подвижность носителей заряда, , – векторы напряженности и индукции электрического поля, создаваемого собственным зарядом, Ejn (r, θ, t), Ejτ (r, θ, t) – компоненты вектора вдоль ортов нормали и касательной к поверхности капли (2.1), P – давление внутри капли, – давление сил поверхностного натяжения, – давление электрических сил, Πinτ, Πexτ – электрические составляющие касательной компоненты тензора напряжений, относящиеся к капле (с нижним индексом «in») и внешней среде (с нижним индексом «ex»), – оператор поверхностной дивергенции (символом Σ обозначена поверхность капли (2.1)), , – компоненты вектора скорости движения жидкости в капле и напряженности электрического поля внутри капли, лежащие в плоскости касательной к поверхности капли (2.1), и вычисляемые по формулам:
Заметим, что в граничном условии (2.5) скорость изменения поверхностной плотности заряда (первое слагаемое в (2.5)) на возмущенной поверхности капли обусловливается несколькими физическими факторами. В этом случае второе слагаемое представляет собой нормальный к свободной поверхности капли ток проводимости за характерное время максвелловской релаксации τr ~ ε1 / γ. Появление третьего слагаемого обусловлено изменением площади свободной поверхности капли в процессе ее капиллярных осцилляций. Четвертый член имеет смысл конвективного потока заряда касательными составляющими течения жидкости за время порядка периода осцилляций капли τ ~ (R3ρ / σ)1/2. Пятый член связан с направленным переносом носителей заряда вдоль свободной поверхности капли под действием касательной составляющей электрического поля.
В уравнении (2.5) отсутствует вклад диффузионной компоненты, так как влияние механизма диффузии на конечность скорости переноса носителей заряда существенно лишь для плохо проводящих жидкостей типа жидкого гелия.
Отыскание решения приведенной выше системы уравнений проведем методом прямого разложения в рамках теории возмущений [14]. Ограничиваясь линейным приближением по малому параметру ε ≡ |ξ| / R << 1, все искомые функции представим в виде асимптотических разложений по :
(2.7)
Номер верхнего индекса указывает на порядок малости по ε соответствующих компонент. Индексом «0» обозначены величины, характеризующие равновесное состояние системы.
Обратим внимание, что в граничных условиях задействованы орты и . В виду этого возникает необходимость в определении явного вида этих функций.
Вычислим единичный вектор внешней нормали к свободной поверхности капли вида F (r, θ, t) = 0 по общей формуле: . В линейном по ε приближении получим:
(2.8)
В плоскости, касательной к возмущенной поверхности капли в любой ее точке, выберем ортогональные орты , : орт касательной к параллелям и орт касательной к меридианам .
Представим касательную к параллелям как , где в качестве постоянного вектора выбран орт оси симметрии . Несложно убедиться, что вектор , имеющий единичную длину, совпадает с азимутальным ортом сферической системы координат:
(2.9)
Определим орт второй касательной при помощи векторного произведения: . Используя (2.8), (2.9) и учитывая, что орт должен иметь единичную длину, найдем единичный вектор касательной в меридиальном направлении в виде:
, (2.10)
, , – орты сферической системы координат.
- Решение задачи нулевого порядка малости по ε. Подставляя разложения (2.7) в систему уравнений (2.2)–(2.6) и собирая вместе слагаемые ~ ε0, получим задачу нулевого порядка малости по безразмерной амплитуде осцилляций ε для невозмущенного состояния системы:
где , , – орты нормали и касательных к равновесной сфере, P0 – равновесное давление внутри капли.
Решение вышеприведенной системы уравнений нетрудно найти:
(3.1)
В итоге приходим к тому, что равновесная форма поверхности капли совпадает с исходной сферической.
- Решение задачи первого порядка малости по ε. Вывод дисперсионного уравнения. Подставляя асимптотические выражения (2.7) в полную математическую формулировку задачи (2.2)–(2.6) и выделяя слагаемые ~ ε1, выпишем гидродинамическую часть задачи в первом порядке малости:
(4.1)
(4.2)
(4.3)
(4.4)
(4.5)
(4.6)
(4.7)
(4.8)
(4.9)
(4.10)
для (4.11)
для (4.12)
(4.13)
(4.14)
, (4.15)
где , , – проекции вектора на орты сферической системы координат.
Скаляризуем уравнения (4.1), (4.2) по методике, подробно изложенной в [15]. Для этого разложим искомое поле скоростей на сумму трех ортогональных векторных полей:
, (4.16)
где ψm – скалярные функции, а дифференциальные векторные операторы в сферической системе координат определяются соотношениями:
, (4.17)
и удовлетворяют условиям ортогональности:
, (4.18)
и условиям коммутативности с оператором Лапласа Δ:
(4.19)
В выражениях (4.17), (4.18) – радиус-вектор, – операторы, эрмитово-сопряженные к , имеющие вид:
(4.20)
Уравнение (4.16) означает, что поле скоростей течения вязкой жидкости представлено суперпозицией потенциальной и вихревой (обусловленных вязкостью) составляющих: оператор выделяет потенциальную компоненту, а операторы – вихревую тороидальную и – вихревую полоидальную части движения реальной жидкости относительно оси симметрии z капли.
Учитывая разложение (4.16) и свойства операторов (4.18), (4.19), приведем систему векторных выражений (4.1), (4.2) к системе трех скалярных уравнений для функций ψm:
, (4.21)
и добавки к давлению внутренней среды капли, вызванной искажением ξ (θ, t) сферической поверхности:
, (4.22)
δ1m – дельта-символ Кронекера.
Прежде чем записать граничные условия (4.6)–(4.9) в терминах скалярных функций ψm, перепишем систему (4.6)–(4.9) через проекции вектора скорости на орты сферической системы координат:
(4.23)
(4.24)
для (4.25)
для (4.26)
Входящие в (4.23)–(4.26) компоненты Vr, Vθ, Vφ поля скоростей выражаются в виде:
где – угловая часть оператора Лапласа в сферических координатах.
После подстановки этих соотношений в (4.23)–(4.26) система скаляризованных граничных условий примет удобный для дальнейшего решения вид:
(4.27)
(4.28)
(4.29)
(4.30)
В (4.28) выражения для линейных по ε поправок к давлениям электрического поля и капиллярных сил определяются в виде:
(4.31)
(4.32)
В случае стоячих капиллярных волн положим зависимость малых величин ξ (θ, t), ψm (r, θ, t) от времени периодической:
,
где i – мнимая единица, ωn – частота n-ой моды капиллярных осцилляций капли, в общем случае комплексная: ωn = Reωn ± iImωn, где вещественная часть Reωn определяет собственную частоту осцилляций, а мнимая часть – декремент затухания в случае отрицательной мнимой компоненты: Imωn < 0 или инкремент капиллярной неустойчивости капли по отношению к ее собственному заряду в случае положительной мнимой компоненты Imωn > 0.
Решения уравнения (4.21) для скалярных функций ψm (r, θ, t) с учетом условий ограниченности (4.4) представим в виде разложений [16]:
(4.33)
, (4.34)
где μ ≡ cosθ, Pn (μ) – полином Лежандра n-ого порядка, n – целое положительное число, In (x) – модифицированная функция Бесселя первого рода [17]. В (4.34) нижний индекс m в амплитудном коэффициенте Bmn указывает на вихревую тороидальную (m = 2) или полоидальную (m = 3) компоненту поля скоростей движения жидкости в капле.
Отклонение поверхности осциллирующей капли от равновесной сферической запишем в форме ряда по осесимметричным полиномам Лежандра:
(4.35)
Здесь суммирование начинается с n = 2, так как в первом порядке малости невозможно возбуждение нулевой (n = 0) и первой (n = 1) мод, отвечающих за радиальные пульсации капли и ее движение как целого соответственно, в силу дополнительных условий (4.14).
В (4.33)–(4.35) неизвестные амплитудные коэффициенты An, Bnm, Mn имеют первый порядок малости по и вычисляются из системы граничных условий (4.27), (4.29), (4.30).
Сосредоточим свое внимание на вычислении в динамическом граничном условии для нормальной компоненты тензора напряжений (4.28) линейных по ε компонент давлений P(1), , , вызванных деформацией ξ (θ, t) равновесной формы поверхности капли.
Используя решение (4.33) для функции ψ1 (r, θ, t), из формулы (4.22) при r = R выпишем добавку к гидродинамическому давлению внутри капли P(1):
(4.36)
Для поправки к давлению сил поверхностного натяжения из (4.32) с учетом разложения (4.35) для возмущения ξ (θ, t) получается соотношение:
(4.37)
Чтобы найти выражения для добавки к давлению электрических сил и электрических составляющих касательной компоненты тензора напряжений , , необходимо решить краевую задачу (4.3)–(4.5), (4.10)–(4.13), (4.15) для напряженности поля .
Для решения векторных уравнений (4.3) в рамках метода операторной скаляризации представим напряженность электрического поля в виде суперпозиции:
, (4.38)
где – неизвестные скалярные функции, – векторные дифференциальные операторы, описываемые выражениями (4.17) и подчиняющиеся свойствам (4.18), (4.19).
Подставляя (4.38) в уравнение непрерывности электрического поля (4.3) и применяя условие ортогональности, придем к уравнению Лапласа для скалярной функции :
(4.39)
В силу свойства коммутативности операторов-проекторов с оператором Лапласа (4.19) подстановка разложения (4.38) в волновое уравнение (4.3) позволяет выписать систему трех скалярных уравнений:
Положим временную зависимость функций экспоненциальной: . Тогда подействовав слева на выше записанное уравнение поочередно эрмитовосопряженными операторами , при использовании свойства (4.18) перейдем к уравнениям Гельмгольца:
, (4.40)
где k – волновое число.
Из уравнения (4.40) при значении индекса m = 1 с учетом (4.39) легко прийти к решению: , что соответствует отсутствию потенциальной составляющей электрического поля. В результате поправка к напряженности поля , связанная с возмущением ξ (θ, t) равновесной формы капли, представляется только вихревой частью:
Учитывая явный вид векторных операторов , , запишем составляющие вектора в сферических координатах (r, θ, φ), выраженные через неизвестные функции , :
(4.41)
(4.42)
Решение уравнения Гельмгольца (4.40) при j = 1 (внутри капли), ограниченное в начале координат, представим в виде суперпозиции стоячих волн:
(4.43)
Решение уравнения (4.40) при j = 2 (вне капли), удовлетворяющее условию ограниченности на бесконечном удалении от капли, естественно искать в виде суперпозиции бегущих волн:
(4.44)
Здесь jn (z) – сферическая функция Бесселя первого рода, – сферическая функция Бесселя третьего рода или функция Ханкеля [18], , – неизвестные коэффициенты первого порядка малости по ε, вычисляемые из граничных условий (4.10)–(4.13).
Несложно показать, что для выполнения условия равенства касательных компонент при необходимо потребовать, чтобы все константы , в решениях (4.43), (4.44) для функций (j = 1, 2) обращались в ноль. Из сказанного следует отсутствие тороидальной части поля: . Таким образом, напряженность электрического поля внутри и вне капли описывается лишь полоидальной составляющей:
(4.45)
(4.46)
Исходя из условия непрерывности тангенциальных компонент электрического поля для орта касательной , найдем соотношение, связывающее амплитуды Mn, , разложений (4.35), (4.45), (4.46):
(4.47)
Условие скачка нормальной компоненты вектора электрической индукции на заряженной границе жидкости (4.10) после подстановки в него решений (4.45), (4.46) при использовании амплитудного коэффициента (4.47) позволяет определить поправку первого порядка малости к равновесной плотности заряда:
(4.48)
Применяя (4.23), перепишем векторное условие баланса поверхностной плотности заряда в скалярном виде:
из которого с учетом уравнений (4.33), (4.34) при m = 3, (4.45), (4.48) выразим коэффициент в (4.46) через амплитуды Mn, An, B3n в виде:
(4.49)
Заметим, что интегральное условие сохранения полного заряда капли (4.15) выполняется тождественно.
Наконец, подставляя в (4.31) полученные решения (3.1), (4.46) для напряженности поля , с учетом амплитуды (4.49), запишем при r = R выражение для добавки к давлению электрических сил :
(4.50)
Далее перейдем к вычислению электрических частей касательных компонент тензора напряжений в (4.29), (4.30). Воспользовавшись соотношениями (2.8)–(3.1), (4.45), (4.46), убедимся, что в первом порядке малости по ε компоненты , а рассчитывается по формуле , из которой при r = R несложно найти:
Рассмотрим теперь граничные условия (4.27), (4.29), (4.30). Подставим найденные электрические составляющие в (4.29), (4.30) и подействуем на эти уравнения слева оператором . При использовании свойства ортогональности операторов (4.18) и равенства преобразованное уравнение (4.29) и условие (4.30) с учетом разложений (4.33), (4.34) при m = 3 позволяют определить связь неизвестных констант An, B3n с амплитудами капиллярных осцилляций Mn:
а из граничного условия (4.30) для неизвестной скалярной функции ψ2 (r, θ, t) после подстановки в него представления (4.34) при m = 2 несложно получить уравнение вида:
Легко заметить, что последнее равенство может выполняться, когда все коэффициенты B2n положить равными нулю: B2n = 0, либо приравнять нулю выражение в квадратных скобках. Из первого случая следует обнуление скалярной функции ψ2 (r, θ, t), что равносильно отсутствию вихревых тороидальных движений внутренней среды капли при ее капиллярных осцилляциях. Во втором случае получаем дисперсионное уравнение, определяющее апериодически затухающие вихревые тороидальные движения жидкости. При этом краевая задача для отыскания ψ2 (r, θ, t) полностью автономна и никак не связана с остальной частью анализируемой системы. В виду этого в нижеследующем изложении функция ψ2 (r, θ, t) рассматриваться не будет.
Для дальнейшего анализа решения удобней рассмотреть динамическое граничное условие для нормальной составляющей тензора напряжений (4.28) относительно амплитуд возмущения Mn, которые легко оценить в экспериментальных и естественных условиях [11]. Для этого подставим в (4.28) выражения (4.36), (4.37), (4.50) для компонент давлений P(1), , и разложения (4.33), (4.34) для скалярных функций ψ1 (r, θ, t), ψ3 (r, θ, t) с учетом полученных амплитудных коэффициентов An, B3n. Применяя известные рекуррентные соотношения для функции Бесселя In (x) (см. [18], стр. 263):
,
и пользуясь ортогональностью полиномов Лежандра, получим общий вид дисперсионного уравнения для полоидальных движений жидкости в излучающей капле:
(4.51)
В (4.51) параметр Релея W ≡ Q2 / 4πσR3 характеризует устойчивость заряженной поверхности проводящей капли по отношению к собственному заряду (т.е. в случае, когда силы поверхностного натяжения преобладают над электрическими силами). При этом условие устойчивости капиллярных осцилляций n-ой колебательной моды представляется в виде W < (n + 2) [19]. Для основной моды n = 2 критическая для начала реализации электростатической неустойчивости проводящей капли величина ее заряда выписывается в виде: .
Несложно заметить, что учет конечности скорости перераспределения заряда капли приводит к увеличению порядка алгебраического уравнения (4.51) по сравнению с дисперсионным уравнением при отсутствии эффекта релаксации заряда. Это объясняется усложнением структуры спектра реализующихся движений внутренней среды капли за счет появления релаксационных апериодических движений жидкости, связанных с перераспределяющимся по возмущенной поверхности электрическим зарядом. Возникающие приповерхностные электрические токи, обусловленные сдвигом фаз между капиллярным волновым движением жидкости и течением, связанным с перемещением носителей заряда, приводят к диссипации кинетической энергии в виде тепла.
Обратим внимание, что влияние вязкости на капиллярные осцилляции капли (на основной (второй) колебательной моде n = 2) оценивается величиной безразмерного коэффициента кинематической вязкости [20]. Как видно, безразмерная комбинация μ зависит не только от самого коэффициента кинематической вязкости ν, но и от плотности жидкости ρ, коэффициента поверхностного натяжения σ и радиуса капли R. В связи с этим при неизменных величинах ν, ρ, σ жидкость может проявлять себя как маловязкая (μ << 1) и сильновязкая (μ >> 1) в зависимости от размеров капель. Таким образом, выполнение критерия сильной вязкости (μ >> 1) возможно для очень мелких капель, когда их осцилляции подавляются вязкостью. Однако осцилляции капель имеют место для малых значений вязкости (μ << 1), соответствующих большим радиусам.
Оценим величину безразмерного параметра μ, принимая во внимание, что линейные размеры R туманных, облачных и дождевых капель варьируются в пределах от R = 1 мкм до R = 3.5 мм [21]. Используя средние характеристики водяных капель σ = 73 дин/см, ρ = 1 г/см3, ν = 0.01 см2/с и принятые радиусы R, легко показать, что значение параметра μ изменяется в диапазоне от μ = 2 · 10−3 при R = 3.5 мм до μ ≈ 0.117 при R = 1 мкм. Это означает, что во всех случаях величина безразмерной вязкости удовлетворяет условию μ << 1, что применительно к осцилляциям рассматриваемых капель.
Исходя из выше сказанного, для перехода к пределу малой вязкости (μ << 1) воспользуемся асимптотическим выражением функции Бесселя In (x) при больших значениях аргумента [18, стр. 262]:
,
с учетом которого имеем:
, (4.52)
Ο – символ порядка [14].
Учитывая асимптотическое поведение функции Бесселя jn (z) при малых значениях аргумента [18, стр. 257]:
,
для которой справедливо рекуррентное соотношение (см. [18, стр. 258]):
,
найдем асимптотическое выражение:
(4.53)
Используя разложение сферической функции в степенной ряд [18]:
,
для случая z = z1 << 1 получим решение:
(4.54)
После подстановки асимптотических разложений (4.52)–(4.54) в (4.51), выпишем кубическое дисперсионное уравнение в относительно простом виде:
(4.55)
,
где ωn0 – собственная частота капиллярных осцилляций заряженной капли проводящей жидкости.
Как видно из (4.55), учет конечности скорости переноса заряда описывается соотношением , которое при ε1 → ∞ обращается в ноль. В этом случае отсутствие эффекта релаксации заряда обуславливается мгновенным перераспределением носителей заряда. В результате (4.55) превращается в дисперсионное уравнение для затухающих капиллярных осцилляций заряженной капли идеально проводящей маловязкой жидкости с учетом энергопотерь на электромагнитное излучение:
(4.56)
где – волновое число.
- Решение дисперсионного уравнения. Обратимся теперь к решению дисперсионного уравнения (4.55). Данные натурных измерений средних зарядов на каплях естественного происхождения с типичными размерами от R = 1 мкм до R = 3.5 мм показывают, что в реальных условиях большая часть капель несет малые заряды, не превышающие одной десятой Рэлеевского предела [21]. Тогда в случае устойчивой поверхности капли по отношению к собственному заряду общий вид комплексной частоты задается в виде: ωn = Reωn − iImωn. Вещественная компонента Reωn соответствует собственной частоте осцилляций, а мнимая отрицательная часть определяет полный декремент затухания Imωn = η1 + η2, связанный с потерями первоначально запасенной энергии капиллярных осцилляций капли на диссипацию в виде тепла за счет вязкости (η1) и на генерацию электромагнитных волн (η2).
При помощи пакета аналитических вычислений, учитывая условие η2 << Reωn (это следует из того, что формула для η2, определяемая слагаемыми в (4.55), зависящими от параметра Релея W, будет содержать малый параметр z1 = kR << 1), получим комплексные решения кубического уравнения (4.55). Их реальные и мнимые (j = 1−3) составляющие для разных диапазонов радиусов капель представляются в виде:
для R < 42 мкм:
(5.1)
(5.2)
(5.3)
(5.4)
(5.5)
(5.6)
для R ≥ 42 мкм:
(5.7)
(5.8)
(5.9)
(5.10)
(5.11)
(5.12)
В выражениях (5.1)–(5.12) первый корень соответствует зарядово-релаксационному апериодическому движению жидкости, второй и третий корни определяют затухающие капиллярные осцилляции капли (, ).
В (5.1)–(5.12) верхний индекс в круглых скобках (j = 1−3) указывает на порядковый номер корня дисперсионного уравнения (4.55).
Вспомним, что волновое число и, принимая во внимание отсутствие реальной части (5.1), (5.7) мнимой частоты , из (5.4), (5.10) найдем декремент затухания зарядово-релаксационного апериодического движения жидкости:
для R < 42 мкм:
(5.13)
для R ≥ 42 мкм:
Из полного декремента затухания (j = 2, 3), определяемого выражениями (5.5), (5.6) при R < 42 мкм и (5.11), (5.12) при R ≥ 42 мкм, выделим отдельные решения для декремента затухания (j = 2, 3) полоидального движения жидкости в капле (при вязкости ν ≠ 0):
для R < 42 мкм:
(5.14)
для R ≥ 42 мкм:
Заметим, что полученные соотношения , определяются двумя первыми слагаемыми в (5.5), (5.6), (5.11), (5.12), не зависящими от волнового числа (в полном декременте затухания принимается c → ∞).
Подставляя в третье слагаемое в (5.5), (5.6), (5.11), (5.12) волновое число (j = 2, 3) и выражения (5.2), (5.3) при R < 42 мкм, (5.8), (5.9) при R ≥ 42 мкм для реальных компонент (j = 2, 3), из полного декремента затухания (в отсутствие вязкости ν = 0) выразим поправки (j = 2, 3), связанные с энергопотерями осциллирующей заряженной капли на излучение электромагнитных волн:
для R < 42 мкм:
(5.15)
для R ≥ 42 мкм:
где коэффициенты y0, y зависят от физических параметров задачи.
Отметим, что поправки (j = 2, 3) к декременту затухания на электромагнитное излучение зависят от волнового числа k и определяются конечным значением скорости распространения электромагнитных волн в вакууме c.
Для отыскания численных решений дисперсионного уравнения (4.55), положим, что осцилляции капель связаны с возбуждением основной (второй) моды n = 2, соответствующей сфероидальным деформациям [20]. Примем для расчетов согласно справочным данным [21]: σ = 73 дин/см, ρ = 1 г/см3, ν = 0.01 см2/с, ε1 = 80, γ = 1 · 106 СГСЭ, b = 300 СГСЭ. Тогда для внутриоблачных капель, имеющих средние заряды Q = 2 · 10−5 СГСЭ (~ 0.06Qкр при R = 3 мкм и ~ 0.002Qкр при R = 30 мкм), при R = 3 мкм справедливо решение:
Для капли размером R = 30 мкм найдем:
Используя средние заряды дождевых капель Q = 7 · 10−4 СГСЭ (~ 2 · 10−3Qкр при R = 0.25 мм и ~ 3 · 10−5Qкр при R = 3.5 мм), для мелкой облачной капли с радиусом R = 0.25 мм получим следующие результаты:
Для крупной дождевой капли с характерным линейным размером R = 3.5 мм имеем:
Анализ численных расчетов показывает, что первый корень является чисто мнимым и соответствует релаксационному апериодическому затухающему движению жидкости, порожденному перетеканием заряда по возмущенной поверхности капли. Второй и третий корни содержат одинаковые мнимые части и разные по знаку реальные компоненты. При этом , ответственны за быстро затухающие капиллярные осцилляции облачных и дождевых капель. В этом случае величина декремента затухания (j = 2, 3) полоидального движения жидкости в капле на один-два порядка слабее частоты собственных осцилляций (j = 2, 3).
Кроме того, заметим, что среди решений , физически значимым является второй корень , имеющий положительную вещественную часть > 0. В связи с этим для дальнейшего анализа будем использовать решение .
В рассматриваемой ситуации декремент затухания проявляется в подавлении осцилляций капли и, следовательно, на прекращении радиоизлучения на соответствующей частоте. Так, осцилляции (электромагнитное излучение) на основной моде n = 2 облачной капли размером R = 3 мкм исчезают при значении кинематической вязкости ν = 0.08 см2/с, а для капли R = 30 мкм – при ν = 0.26 см2/с.
Численные расчеты показывают, что по сравнению с каплей идеально проводящей жидкости учет конечности скорости перераспределения заряда приводит к небольшому увеличению собственной частоты осцилляций при снижении радиуса капли.
Из рис. 1, на котором представлены расчетные зависимости (R, W) по формуле (5.2), прослеживается весьма быстрое уменьшение частоты осцилляций при возрастании размера R и заряда Q капли (величины параметра Релея W). Заметим, что кривые построены для основной моды n = 2 при значениях параметра Релея, достаточно близких к критическому Wкр = 4. Однако в реальных условиях грозового облака заряды на каплях такие, что величина W << 1, поэтому приведенные зависимости имеют лишь качественный характер.
Рис. 1. Зависимость собственной частоты капиллярных осцилляций вязкой заряженной облачной капли:
а: от ее радиуса R; кривые 1–4 соответствуют значениям W = 0.1, 2, 3, 3.9
б: от величины параметра Релея W; кривые 1–4 соответствует значениям R = 5, 10, 15, 20 мкм
Расчеты проведены по формуле (5.2) при n = 2, σ = 73 дин/см, ρ = 1 г/см3, ν = 0.01 см2/с, ε1 = 80, γ = 1 · 106 СГСЭ, b = 300 СГСЭ.
Графики зависимостей (R, ν), изображенные в соответствие с (5.14) на рис. 2, показывают, что декремент затухания , при возрастании радиуса R – снижается по обратной квадратичной зависимости, а с ростом кинематической вязкости ν увеличивается по линейному закону.
Рис. 2. Зависимость декремента затухания капиллярных осцилляций вязкой облачной капли c зарядом Q = 2 · 10−5 СГСЭ (~ 0.06Qкр при R = 3 мкм и ~ 0.002Qкр при R = 30 мкм) от:
а: радиуса R (а) кривые 1–4 соответствуют значениям ν = 0.01, 0.03, 0.05, 0.1 см2/с
б: коэффициента кинематической вязкости ν, кривые 1–4 соответствует значениям R = 5, 10, 15, 20 мкм
Расчеты проведены по формуле (5.14) при тех же значениях параметров, что на рис. 1.
Согласно численным расчетам по (5.13) электропроводность жидкости γ и поверхностная подвижность носителей заряда b не оказывают влияние на декремент затухания капиллярных осцилляций. При этом зависимость собственной частоты от величин γ, b весьма мала.
Из общефизических соображений очевидно, что периодически меняющаяся амплитуда Mn (t) возмущения ξ (θ, t) равновесной формы капли за счет вязкости затухает со временем по экспоненциальному закону:
(5.16)
,
где ε0 = 0.1 – безразмерная амплитуда возмущения.
На рис. 3, а иллюстрируются численно рассчитанные по соотношению (5.16) зависимости от времени t амплитуды M2 (t) второй моды n = 2 капиллярных осцилляций мелкой облачной капли при различных радиусах. На приведенных графиках можно наблюдать, что внутриоблачная капля наименьшего размера R = 3 мкм, осциллирующая на частоте = 4.61 · 106 рад/с с декрементом затухания = 0.56 · 106 рад/с и периодом затухающих осцилляций T = 1.37 мкс, возвращается к равновесной форме, совершив ≈ 8 колебаний. Для капли большего радиуса R = 4 мкм с характеристиками осцилляций = 3 · 106 рад/с, = 0.31 · 106 рад/с, T = 2.1 мкс число колебаний до полного их исчезновения увеличивается до ≈ 10. Капля размером R = 5 мкм при = 2.15 · 106, = 0.2 · 106 рад/с, T = 2.93 мкс совершает уже ≈ 11 колебаний за время вязкого затухания, при котором искажение равновесной формы снижается до нуля.
Рис. 3. Зависимость амплитуд Mn (t) различных мод возмущения равновесной формы вязкой облачной капли c зарядом Q = 2 · 10−5 СГСЭ (~ 0.06Qкр при R = 3 мкм и ~ 0.03Qкр при R = 5 мкм) от времени t.
а – n = 2, б – n = 5. Кривые 1–3 соответствуют значениям R = 3, 4, 5 мкм.
Расчеты проведены по формуле (5.16) при тех же значениях параметров, что и на рис. 1.
Из результатов численных расчетов и рис. 1, а можно сделать заключение, что снижение радиуса капли приводит к уменьшению периода T затухающих осцилляций так, что чем меньше размер капли, тем ее осцилляции быстрее гасятся вязкостью.
Для сравнения на рис. 3, б приведены временные зависимости амплитуды M5 (t) пятой моды n = 5, построенные по (5.16) при тех же значениях физических параметров, что на рис. 3, а. Несложно видеть, что скорость снижения амплитуды отклонения поверхности капли от равновесной сферы увеличивается с ростом номера моды. Указанное обстоятельство объясняется повышением величины декремента затухания на высоких колебательных модах (см. рис. 4): увеличение номера моды в три раза приводит к возрастанию на порядок.
Рис. 4. Зависимость декремента затухания капиллярных осцилляций капли радиуса R = 5 мкм c зарядом Q = 2 · 10−5 СГСЭ (~ 0.03Qкр) от номера моды n.
Расчеты проведены по формуле (5.14) при тех же значениях параметров, что на рис. 1.
Для иллюстрации изобразим временную зависимость экспоненциально убывающей со временем амплитуды Mn (t) возмущения ξ (θ, t) в случае релаксационного движения внутренней среды капли, порождаемого переносом заряда по искривленной сферической поверхности, в соответствие с законом:
, (5.17)
где безразмерная амплитуда возмущения принимается ε0 = 0.1.
Из рис. 5, полученному по выражению (5.17) при n = 2, видно, что релаксационное движение жидкости имеет апериодический быстро затухающий характер. В соответствие с численными расчетами для декремента затухания зарядово-релаксационного движения жидкости прослеживается линейная зависимость от проводимости γ и подвижности зарядов b. Причем, увеличение γ сказывается на весьма слабом возрастании декремента затухания , а увеличение подвижности в два раза приводит к удвоению величины b.
Рис. 5. Временная эволюция амплитуды M2 (t) основной моды n = 2 возмущения равновесной формы капли радиуса R = 5 мкм с зарядом Q = 2 · 10−5 СГСЭ (~ 0.03Qкр), соответствующая зарядово-релаксационному апериодическому движению жидкости.
Расчеты проведены по формуле (5.17) при тех же значениях параметров, что на рис. 1.
На рис. 6 и 7 изображены рассчитанные по (5.15) зависимости декремента затухания , связанного с излучением электромагнитных волн, от удельной электропроводности жидкости γ и подвижности носителей заряда b. Приведенные кривые качественно ожидаемы из общефизических тенденций: поправка к декременту затухания почти линейно возрастает при изменении величин γ, b.
Рис. 6. Зависимость величины поправки декремента затухания капиллярных осцилляций, связанного с излучением электромагнитных волн капли радиуса R = 10 мкм и заряда Q = 2 · 10−5 СГСЭ (~ 0.01Qкр), от удельной проводимости γ.
Расчеты проведены по формуле (5.15) при тех же значениях параметров, что на рис. 1.
Рис. 7. Зависимость величины поправки декремента затухания капиллярных осцилляций, связанного с излучением электромагнитных волн капли радиуса R = 10 мкм и заряда Q = 2 · 10−5 СГСЭ (~ 0.01Qкр), от подвижности носителей заряда b.
Расчеты проведены по формуле (5.15) при тех же значениях параметров, что на рис. 1.
Что касается малой добавки к декременту затухания капиллярных осцилляций, определяемой энергопотерями капли на электромагнитное излучение, рассмотрим ее поведение при конечных значениях электропроводности жидкости (рис. 8, а), а также в двух предельных случаях: для капли идеально проводящей жидкости (ε1 → ∞) (рис. 8, б) и капли идеально диэлектрической жидкости (ε1 → 0) (рис. 8, в). Из графиков рис. 8, построенных в соответствие с (5.15), можно заметить, что величина поправки к декременту затухания, связанного с излучением электромагнитных волн, сильно заряженных капель в случае конечной удельной проводимости жидкости на порядок ниже по сравнению с идеальным проводником и на два порядка выше по сравнению с идеальным диэлектриком. Из этого факта следует, что чисто диэлектрическая капля, для которой заряд равномерно распределен по поверхности и вморожен в нее, при осцилляциях излучает электромагнитные волны с наименьшей интенсивностью.
Рис. 8. Зависимость поправки декремента затухания капиллярных осцилляций, связанного с излучением электромагнитных волн заряженной облачной капли, от величины параметра Релея W.
а – ε1 = 80, б – ε1 = ∞, в – ε1 = 0; кривые 1–3 соответствуют значениям R = 26, 28, 30 мкм.
Расчеты проведены по формуле (5.15) при тех же значениях параметров, что на рис. 1.
Как видно из рис. 8, увеличение параметра Релея W приводит к весьма быстрому нарастанию добавки к декременту затухания капиллярных осцилляций капли, связанного с радиоизлучением. Достигнув наибольшего значения, величина снижается до нуля в широком диапазоне изменения параметра Релея W (вплоть до критических значений). Кроме того, рис. 8 показывает, что возрастание размера капли сказывается на быстром уменьшении малой поправки .
- Электромагнитное излучение. Известно, что вихревая полоидальная часть течения связана с вязкостью реальной жидкости. Тогда механическая энергия капиллярных осцилляций капли, определяемая кинетической энергией полоидальной составляющей поля скоростей движения реальной жидкости, теряется на диссипацию в виде тепла. В случае маловязкой жидкости положим диссипацию энергии капиллярных осцилляций пренебрежимо малой. Учитывая это обстоятельство, в приближении идеальной жидкости (ν = 0) вязкое затухание отсутствует. Тогда согласно закону сохранения энергии потери механической энергии капиллярных осцилляций, определяемые кинетической энергией потенциальной составляющей поля скоростей течения идеальной жидкости, связываются с излучением электромагнитных волн, вызванным наличием ускоренно движущегося заряда. Принимая во внимание, что энергия поверхностных колебаний n-й моды ϑn снижается по экспоненциальному закону: ϑn ~ exp (−2η2), запишем общее выражение для мощности излучения [1]:
, (6.1)
где величина η2 находится из дисперсионного уравнения (4.56) в виде [22]:
(6.2)
При использовании теоремы вириала представим η2 в виде удвоенной средней за период кинетической энергии потенциального движения жидкости в капле:
, (6.3)
где Mn = εR – размерная амплитуда осцилляций n-ой колебательной моды.
В итоге, подставляя (6.2), (6.3) в (6.1), получим окончательное аналитическое выражение для мощности излучения единичной заряженной капли идеальной проводящей жидкости, осциллирующей на фиксированной частоте ωn0:
(6.4)
Согласно данным натурных измерений большая часть облачной воды сосредоточена в каплях с типичными радиусами от 3 до 30 мкм и концентрацией ~ 103 см–3. Максимум функции распределения таких капель по размерам приходится на диапазон 3–7 мкм [21].
Подставляя в (6.4) приведенные ранее средние характеристики внутриоблачных капель, осциллирующих на основной (второй) моде n = 2 с амплитудой M2 = 0.1R, несложно оценить интенсивность электромагнитного излучения, равную I ~ 1 · 10−39 эрг/с, для капли размером R = 3 мкм и зарядом Q = 5 · 10−5 СГСЭ (~ 0.16Qкр), осциллирующей на частоте ω20 = 4.6 · 106 рад/с с поправкой к декременту невязкого затухания η2 = 2 · 10−33 рад/с, идущего на генерацию электромагнитного излучения. При возрастании радиуса капли на порядок R = 30 мкм, имеющей заряд Q = 1.1 · 10−3 СГСЭ (~ 0.11Qкр), мощность излучения снижается на два порядка величины: I ~ 7 · 10−42 эрг/с при ω20 = 146 · 103 рад/с и η2 = 1 · 10−37 рад/с.
На рис. 9 приведены результаты расчетов по формуле (6.4) в виде зависимостей I (R, W). Видно, что интенсивность излучения I единичной облачной капли идеальной идеально проводящей жидкости, осциллирующей на основной (второй) моде n = 2, при возрастании радиуса R убывает по обратному квадратичному закону. Несложно заметить, что график функции I (W) (рис. 9, а) качественно повторяет зависимость поправки (W) к декременту затухания капиллярных осцилляций капли маловязкой идеально проводящей жидкости, определяемого потерями на электромагнитное излучение (рис. 8, б). Интересно отметить, что учет малой вязкости по сравнению с идеальной жидкостью приводит лишь к весьма малым количественным изменениям.
Рис. 9. Зависимость интенсивности электромагнитного излучения I заряженной капли идеальной идеально проводящей жидкости от параметра Релея W (а) и радиуса R (б).
а: Кривые 1–3 соответствуют значениям R = 26, 28, 30 мкм.
б: Кривые 1–3 соответствуют значениям W = 2, 2.5, 3.
Расчеты проведены по формуле (5.4) при n = 2, σ = 73 дин/см, ρ = 1 г/см3.
Сравним теперь между собой интенсивности излучения заряженной невязкой идеально проводящей капли в двух случаях: учитывая в степенном асимптотическом разложении сферической функции Бесселя (z) при малых z (см. [18, стр. 257, ф-ла 10.1.17]) все члены ряда по методике Н.А. Богатова [7] и сохраняя только старшие слагаемые как это делалось в [1].
Зависимость I (R, W) (см. рис. 10), построенная для второй моды осцилляций n = 2 по выражению [1]:
, (6.5)
указывает на качественные и количественные расхождения по сравнению с численными расчетами по формуле (6.4) при n = 2 в методике Н.А. Богатова. Из выражения (6.5) при амплитуде возмущения M2 = 0.1R и рис. 10 видно, что мощность излучения не зависит от размера капли и на 31–33 порядка величины выше (при изменении радиуса облачной капли от R = 3 мкм до R = 30 мкм).
Рис. 10. Зависимость интенсивности электромагнитного излучения I заряженной капли идеальной идеально проводящей жидкости от параметра Релея W, построенная по [1].
Следует отметить, что использование корректного асимптотического представления сферической функции Бесселя (z) (а именно: учет всех членов степенного ряда) приводит к качественному и количественному согласованию полученных результатов интенсивности излучения капли идеальной идеально проводящей жидкости, осциллирующей на основной (второй) моде n = 2, выполненных на основе физико-химической гидродинамики двумя методами: на основе закона сохранения энергии и в рамках теории излучения [11].
Ясно, что при снижении электропроводности жидкости уже сказывается конечность скорости переноса электрического заряда при деформациях равновесной формы капли. Данное обстоятельство приводит к снижению характеристик излучения, полученных без учета эффекта релаксации заряда. Примем для оценки γ = 1 · 106 СГСЭ, b = 300 СГСЭ. Так для наименьшей облачной капли R = 3 мкм с зарядом Q = 5 · 10−5 СГСЭ (~ 0.16Qкр), осциллирующей на частоте = 4.56 · 106 рад/с, получена интенсивность излучения I ~ 8 · 10−40 эрг/с при добавке к декременту затухания капиллярных осцилляций, равной = 1 · 10−33 рад/с. Для наиболее крупной внутриоблачной капли R = 30 мкм, несущей заряд Q = 1.1 · 10−3 СГСЭ (~ 0.11Qкр), интенсивность излучения составила I ~ 5 · 10−42 эрг/с при = 146 · 103 рад/с, = 7 · 10−38 рад/с.
Из проведенных численных расчетов, проиллюстрированных на рис. 11, для заряженной капли вязкой жидкости с конечной электропроводностью, осциллирующей на второй моде n = 2, выясняется, что интенсивность излучения I при изменении размера капли R снижается по закону ~ R−2. При этом зависимость мощности излучения от величины параметра Релея W (рис. 11, а) качественно схожа с графиком зависимости малой поправки (W) к декременту затухания, связанного с излучением электромагнитных волн (рис. 8, а).
Рис. 11. Зависимость интенсивности электромагнитного излучения I заряженной вязкой капли конечной электропроводности от параметра Релея W (а) и радиуса R (б).
Расчеты проведены при тех же значениях параметров, что на рис. 1.
Из сравнения рис. 9 и 11 можно сделать заключение, что учет конечности скорости переноса электрического заряда при капиллярных осцилляциях капли маловязкой жидкости приводит снижению (примерно на порядок величины или в пределах порядка величины в зависимости от параметра Релея W) интенсивности излучения.
Заключение. Согласно вышеизложенному, учет конечности скорости переноса заряда приводит к уменьшению интенсивности радиоизлучения от осциллирующей заряженной капли, и в первую очередь этому способствует вязкость жидкости, снижающая частоты осцилляций.
Работа выполнена в ИПМех РАН в рамках Государственного задания АААА-А20-120011690131-7.
Sobre autores
A. Grigoryev
Ishlinsky Institute for Problems in Mechanics of the RAS
Autor responsável pela correspondência
Email: grigorai@mail.ru
Rússia, Moscow
N. Kolbneva
Demidov Yaroslavl State University
Email: kolbneva-nata@yandex.ru
Rússia, Yaroslavl
S. Shiryaeva
Demidov Yaroslavl State University
Email: shir@uniyar.ac.ru
Rússia, Yaroslavl
Bibliografia
- Kalechits V.I., Nakhutin IE, Poluektov P.P. On the possible mechanism of radio emission of convective clouds // Dokl. AN USSR, 1982, vol. 262, no.6, pp. 1344–1347.
- Shiryaeva S.O., Grigoriev A.I., Golovanov A.S., Rybakova M.V. Electromagnetic radiation generated by linear vibrations of a charged drop // ZhTF, 2002, vol.72, no. 1S, pp. 8–14.
- Grigor’ev A.I., Shiryaeva S.O. Electromagnetic radiation of an oscillating charged drop of finite conductivity // Izv. RAS. Fluid&Gas Mech., 2002, no. 5, pp.74–80.
- Shiryaeva S.O., Grigoriev A.I., Kryuchkov O.S. On oscillations of a charged drop of viscous liquid with final conductivity // ZHTF, 2007, vol. 77, no.6, pp.13–21.
- Bogatov N.A. Electromagnetic field generated by capillary oscillations of drops // VI Int. Conf. “Solar-Earth Connections and Physics of Earthquake Precursors.” Petropavlovsk-Kamchatsky, FEB RAS. Abstracts. 2013. pp.22–26.
- Belonozhko D.F., Grigoriev A.I. On the correct form of recording the law of preserving the amount of matter at the moving interface of two liquid media // ZHTF, 2004, vol.74, no. 11, pp. 22–27.
- Belotserkovsky A.V., Divinsky L.I. Active-Passive Radar of Thunderstorm and Menacing Foci in the Clouds. St. Petersburg: Hydrometeoizdat, 1992. 286 p.
- Kachurin L.G. Physical Basis of Impact on Atmospheric Processes. Leningrad: Hydrometeoizdat, 1990. 463 p.
- Gorelik AG, Kozlov A.I., Sterlyadkin V.V. Scattering of radio waves on non-spherical and oscillating rain drops // Sci. Vest. MSTU GA, 2012, no. 176, pp. 25–30.
- Zhukov V.Yu., Shchukin G.G. Recognition of weather hazards in modern meteorological radar // Problems of Prisoner of War Geophysics and Monitoring of the State of the Natural Environment: Materials All Russian Sci. Conf. SPb.: 2020. pp. 40–50.
- Grigoriev A.I., Shiryaeva S.O. Etudy about a Thunderstorm. Lights of St. Elm, the Glow of Tornadoes, Different Lightning. Moscow;Berlin: Durect-MEDIA. 223 p.
- Sterlyadkin V.V. Full-scale measurements of precipitation drop oscillations // Izv. USSR AS. FAO, 1988, vol.24, no.6, pp. 613–621.
- Beard K.V., Tokay А. A field study of small raindrop oscillations // Geophys. Res. Lett., 1991, vol.18, no.12, pp. 2257–2260.
- Naife A.H. Methods of Perturbation. Moscow: Mir, 1976. 455 p.
- Lazaryants A.E., Shiryaeva SO, Grigoriev A.I. Scalarization of Vector Boundary Problems. Moscow: Russines, 2020. 140 p.
- Arfken G. Mathematical Methods in Physics. Moscow: Atomizdat, 1970. 712 p.
- Varshalovich D.A., Moskalev A.N., Kherson V.K. Quantum Theory of Angular Momentum. Leningrad: Nauka, 1975. 436 p.
- Abramowitz M, Steegan I. Handbook of Mathematical Functions. Moscow: Nauka, 1979. 830 p.
- Rayleigh. On the equilibrium of liquid conducting masses charged with electricity // Phil. Mag., 1882, vol.14, pp.184–186.
- Grigor’ev A.I. On some regularities of the implementation of instability of a highly charged viscous drop // ZHTF, 2001, vol.71, no.10, pp.1–7.
- Mazin I.P., Hrgian AH, Imyanitov I.M. Clouds and Cloud Atmosphere. Leningrad: Hydrometeoizdat, 1989. 647 p.
- Grigor’ev A.I., Kolbneva N.Yu., Shiryaeva S.O. On the effect of the liquid viscosity relaxation effect on the electromagnetic radiation intensity of an oscillating charged drop // Colloid J., 2023, vol.85, no. 4, pp. 483–501.
Arquivos suplementares
