CC BY
46
Вестник КРАУНЦ. Физ.-мат. науки. 2015. № 1(10). C. 12-17. ISSN 2079-6641
МAТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 517.958 [550.3 + 551.5]
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗАКОНА ИЗМЕНЕНИЯ ЗАРЯДА ОБЛАЧНЫХ КАПЕЛЬ ВО ФРАКТАЛЬНОЙ СРЕДЕ
Т.С. Кумыков1, Р.И. Паровик2' 3
1 Институт прикладной математики и автоматизации, 360000, Республика Кабардино-Балкария, г. Нальчик, ул. Шортанова, 89а
2 Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, п. Паратунка, ул. Мирная, 7
3 Камчатский государственный университет имени Витуса Беринга, 683032, г. Петропавловск-Камчатский, ул. Пограничная, 4
E-mail: macist20@mail.ru; romanparovik@gmail.com
В работе предложена новая математическая модель изменения заряда облачных капель в грозовых облаках. Модель учитывает фрактальные свойства грозовых облаков, а ее решение было получено с помощью аппарата дробного исчисления.
Ключевые слова: фрактальная размерность, математическая модель, оператор Римана-Лиувилля, оператор Капуто
(с) Кумыков Т.С., Паровик Р.И., 2015
MATHEMATICAL MODELING
MSC 37C70
MATHEMATICAL MODELING OF CHANGES IN THE CHARGE CLOUD DROPLETS IN A FRACTAL ENVIRONMENT
T.S. Kumykov1, R.I. Parovik2' 3
1 Institute of Applied Mathematics and Automation, 360000, Nalchik, Shortanova st., 89a, Russian
2 Institute of Cosmophysical Researches and Radio Wave Propagation Far-Eastern Branch, Russian Academy of Sciences, 684034, Kamchatskiy Kray, Paratunka, Mirnaya st., 7, Russia
3 Vitus Bering Kamchatka State University, 683031, Petropavlovsk-Kamchatsky, Pogranichnaya st., 4, Russia
E-mail: macist20@mail.ru; romanparovik@gmail.com
The paper proposed a new mathematical model of the variation of the charge cloud drops in storm clouds. The model takes into account the fractal properties of storm clouds, and the solution was obtained using the apparatus of fractional calculus.
Key words: fractal dimension, the mathematical model, operator Riemann-Liouville, operator Caputo
© Kumykov T.S., Parovik R.I., 2015 .
Введение
В течение последнего десятилетия многими геофизиками интенсивно изучается фрактальность природных структур и его влияние на различные геофизические процессы. К таким природным явлениям относится, и облако в котором вопрос образования и разделения электрических зарядов является актуальной. Изучению закономерностей разделения электрических зарядов в облаках посвящены многочисленные исследования, основные результаты которых обобщены в классических работах [1]-[9], где предложено большое количество объяснений без учета фрактальности среды. Результаты исследований в этой области показывают, что одной из важных предпосылок разделения электрических зарядов в облаках, является наличие ледяной фазы (кристалликов льда, крупы и градин) и переохлажденных водяных капель [10].
Известно, что облака с мощными конвективными токами имеют фрактальную структуру и облако является фрактальной средой [11]. Поэтому можно утверждать, что процессы, протекающие в такой среде хорошо описываются с помощью аппарата дробного исчисления.
Постановка и решение задачи
Согласно классическим представлениям, например, в работе [12], следуя теории Френкеля [13], для облачных капель, находящихся в слабо ионизированной воздушной среде, был получен средний заряд , который создается одной облачной каплей радиуса г, в виде
= 4пе°п£ а, (1)
где £° - электрическая постоянная; а - радиус пузырька; £ - электрокинетический потенциал; п - количество пузырьков радиуса а, образующихся в облачной капле радиуса г. Поэтому опираясь на теорию Френкеля, полный заряд капли можем записать в следующей форме:
q (х, г ) = 4пе°£ Я (х, г), (2)
где Я(х,г) - радиус капли. Закон изменения заряда капли, с учетом уравнения (2) можем записать в виде:
дq (хг) „„ г Я (x, г) ,оч
В уравнении (3) величина ] (х, г) = дq^г) представляет собой плотность потока
д г
Я (х, г)
заряда, который зависит от скорости изменения радиуса капли —-— и совпадает с
дг
диффузионным потоком у поверхности капли, при условии ее роста за счет диффузии из окружающей среды [14].
В случае, когда процесс протекает во фрактальной среде, то вместо модели (3)необходимо рассмотреть обобщенную модель изменения заряда капли, который учитывает фрактальные свойства среды. Изменение заряда капли зависит от изменения ее размера (3). Поэтому рассмотрим сначала как изменяется размер капли в среде с фрактальной структурой.
Согласно работе [15], диффузионный поток можно определить следующим образом:
q (х, г) = -Ш^и (х, г), 0 < а < 1, (4)
где к - коэффициент диффузии вещества; и (х,г) - концентрация вещества; - оператор дробного дифференцирования в смысле Римана-Лиувилля дробного порядка а с началом в точке а, который определяется следующим образом [16]:
х
а _ 1 д ги (д, г) йд (д, г) = Г (1 - а) д*У (х - д )а •
а
В случае, когда а = 1 оператор Римана-Лиувилля переходит в класическую производную первого порядка.
Отметим, что замена оператора д/дг на дробный оператор в дифференциальных уравнениях неявно включает дополнительные факторы взаимодействия физической системы. Дробный параметр а зависит от фракталльной размерности среды, поэтому можно считать, что уравнение (4) описывает фрактальный процесс [15].
Пусть а = 0, тогда из формул (3) и (4) получим:
к
3 (х, г ) = - ^ ^ (х, г) • (5)
к
Обозначим через X =---т и подставим в формулу (5) соотношение для плотности
4п£о<;
потока 3 (х,г) с учетом (3) приходим к уравненнию:
д- X (х, г ) = 0. (6)
Формула (6) представляет собой уравненние с дробной производной порядка а для которого справедливы следующие нелокальные условия [11]:
Я (х, 0) = г1 (х), х е [0, Ь], (7)
lim (x, t) = Г2 (t), t E [0, T],
x-» 0
(8)
Решение задачи (6)-(8) имеет вид [17]:
л
R (x, t ) = J
ri(s)
-1
x — s
,1,0 '1, а
-X t (x - s)'
-jds + X J
Г2 (n) 1,0 f -X (t - n)
x
1,
d n.
(9)
где e ав(z) = i
- функция типа Райта.
п=о Г ( ап + V) Г (8 - вп) Подставляя (9) в формулу (2), получаем соотношение для заряда капли с учетом фрактальной геометрии среды.
q (x, t) = 4ne0Z
/x-)«К fe^)ds + X /e1,a dn
С учетом (10) поток заряженных капель примет вид:
j (x, t) = 4n£0Cdt
r1(s) 1,0 f -Xt x-se1, а V(x-s)a
jds + X I
Г2 (n) 1,0 ^ -X (t - n)
x
1,
dn
(10)
. (11)
x
n
z
x
Упростим формулу (11), принимая во внимание следующую известную формулу дифференцирования интеграла с переменными пределами интегрирования:
х2(г) д
I' (г) = | (х, г) ¿х+/ (х2 (г), г) х2 (г) - / (х1 (г), г) х1 (г),
получим
j (x, t) = 4neoZ
xi(t)
ri (s) 1,0 f -Xt x-s dЛ« V(x-s)a
ds+X*/
i
d i Г2 (n) 1,0 -X (t - n)
x 1,n x«
/I „Лг1 (s) д 1,0 / -At \ Г2 (t) 1,0
4m4T-j w, M ÖVds+X ~Г" e1-(0)+
dn
(12)
+X/ Г2.М l ец( dn.
x dt 1,a
x
С учетом того, что выполняются свойства [17]:
*1£ (0) = 1, ^-1<(Яга) = г5/ (Яга) окончательно получим:
j (x, t) = 4л£оС
л-
X Г2 (t) [ Г1 (s)
е
0,0 1,а
x У (x - s) t 1,а V (x - s)
-X t
ds I +
(13)
+4пе0^Х
r2 (n) o,o f -X (t - n)
xt V xa
dn .
Формула (13) представляет собой закон изменения заряда облачных капель с учетом фрактальных свойств среды записанный с помощи функции типа Райта.
В работе [11], было получено уравнение типа (4) с оператором дробного производного по Капуто:
q (х, г) = уд°0и (х, т), 0 < а < 1, (14)
а-1 du (x Т)
0tu (x, v = D0t dТ
где у > 0, д^и (x, т) = D(
- регуляризованная дробная производная порядка
о от функции и (х, т) с началом в точке 0 и т .С учетом формулы (14) и закона изменения размера капли, считая функцию Я зависящей только от времени аналогично формуле (3) перепишем в виде:
за* (t) - kR (t )=o,
(15)
1
где к = -. Уравнение(15)- обыкновенное дифференциальное уравнение дробного по-7
рядка, Для уравнения (15) справедливо начальное условие вида:
R (x, 0) = R0.
(16)
x
t
Следуя работе [17], решение задачи может быть записано в терминах функции типа Миттаг-леффлера следующим образом:
R (t)= RoEa,i (ktа), (17)
zk
где Ea в (z) = L —;-¡5Т - функция типа Миттаг-Леффлера. Подставляя решение
k=o Г (ak + в) (17) в формулы (2) и (3), получим:
q (t) = 4neoCRoEa,i (kta), j (t) = 4n£oCRota-1Ea,a (kta). (18)
Формула (18) представляют собой закон изменения заряда капли в обобщенной теории Френкеля в облачной среде через функцию Миттаг-Леффлера.
Заключение
Рассматривая облака, которые, как известно, имеют разную структуру, и имеют свою классификацию по происхождению и морфологическим признакам, к которой можно добавить данные об их фрактальной структуре, возможно в дальнейшем -формирование более общей картины состояния физики облаков.
В данной работе предложена математическая модель закона изменения заряда капель во фрактальной облачной среде, обобщающая теорию Френкеля. Получено решение этой модели с учетом функции типа Райта и функции Миттаг-Леффлера.
Библиографический список
1. Качурин Л. Г., Морачевский В.Г. Кинетика фазовых переходов воды в атмосфере. Л.: Гидрометео-издат, 1965. 114 с.
2. Мейсон Б. Дж. Физика облаков. Л.: Гидрометеоиздат, 1961. 542 с.
3. Мучник В. М. Физика грозы. Л: Гидрометеоиздат, 1974. 252-257 c.
4. Чалмерс Дж. А. Атмосферное электричество. Л.: Гидрометеоиздат, 1974. 420 с.
5. Юман М. Молния. М.: Мир, 1972. 328 с.
6. Ribeira J.C. On the thermo-dielectric effect // Ann. Acad. Brasil. Sci. 1950. vol. 22. №3. P. 547-556.
7. Workman E. J., Reynold S.E. Electrical phenomena occurring during the freezing of delute aqueous solution and their possible relationship to thunderstorm electricity // Phys.Rev. 1956. vol. 94. №4. P. 1073-1075.
8. Workman E. J. The possible role of ammonia in thunderstorm electrification. In: Proc. Intern. Conf. Cloud Phys.. Toronto. 1968. P. 653-656.
9. Имянитов И. М. Электрическая структура конвективных облаков и ее связь с движением воздуха в облаках. Исследование облаков, осадков и грозового электричества // М., Гидрометеоиздат, 1961, с. 225-238.
10. Аджиев А.Х. Куповых Г.В. Атмосферно-электрические явления на Северном Кавказе. Таганрог: 2004. 122 c.
11. Нахушев А.М. Дробное исчисление и его применение. М: Физматлит, 2003. 272 с.
12. Кумыков Т.С. Жекамухов М.К., Каров Б.Г. Электризация и пространственное разделение зарядов при выделении пузырьков воздуха в процессе коагуляционного роста градин в облаке II. Генерирование грозового электричества за счет выделения заряженных пузырьков при намерзании переохлажденных облачных капель на поверхности градин // Метеорология и гидрология. 2008. №12. С. 15-24.
13. Френкель Я.И. Теория основных явлений атмосферного электричества // Сборник избранных трудов. М.: Наука. 1958. Т.2. С. 538-567.
14. Лифшиц Е.М., Питаевский Л.П. Физическая кинетика: М.: Наука, 1979. 527 с.
15. Шогенов В.Х., Шхануков-Лафишев М.Х., Бештоев Х.М. Дробные производные: интерпретация и некоторые применения в физике. Сообщения объединенного института ядерных исследований. Дубна. 1997. 20 с.
16. Самко С.Г., Килбас А. А., Маричев О. И. Дробные интегралы и производные и некоторые их приложения. Минск: Наука и техника, 1987. 688 с.
17. Псху А.В. Краевые задачи для дифференциальных уравнений с частными производными дробного и континуального порядка. Нальчик: КБНЦ РАН, 2005. С. 185.
Поступила в редакцию / Original article submitted: 17.05.2015
