CC BY
29
MS С 74Е20
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ИЗМЕНЕНИЯ ЗАРЯДА ПУЗЫРЬКОВ В КАПЛЯХ ОБЛАКОВ С УЧЕТОМ ФРАКТАЛЬНОСТИ СРЕДЫ
Т.С. Кумыков
Институт прикладной математики и автоматизации, ул. Шортанова, 89-А, Нальчик, 360004, Россия, e-mail: macist20@mail.ru
Аннотация. В работе предложена новая модель изменения заряда пузырьков в переохлажденных облачных каплях. Модель учитывает, фрактальные свойства облаков, а ее решение было получено с применением аппарата дробишх) исчисления.
Ключевые слова: фрактальная размерность, математическая модель, облачная капля, пузырьки.
Введение. Известно, что облака с мощными конвективными токами имеют фрактальную структуру |1|, то есть облако является фрактальной средой. Поэтому можно думать, что процессы, протекающие в такой среде, хорошо описываются с помощью аппарата дробного исчисления.
Несмотря на несомненные успехи в изучении процессов в облаках (Качурин Л.Г., Мейсоп Б.Дж., Мучник В.М., Чалмерс Дж.А., Юмап М,, Ribeira J.С., Workman E.J., Reynold S.E., Имянитов И.М., Френкель Я.И., и т.д.), многие из них до настоящего времени изучены на недостаточном уровне. Это относится и к процессам электризации облачных частиц (капель), к влиянию электрического ноля на микрофизические процессы, к влиянию микроструктуры облака на заряд и ноле, и ко многим другим процессам.
Таким образом, для физики облаков представляет большой интерес разработка математических моделей конвективных облаков с детальным учетом влияния фракталь-ности среды на различные геонроцессы в облаках, способствующих развитию общей картины физики облаков. Математическое моделирование позволяет детально изучать как отдельные физические процессы, так и их взаимодействие между собой, а применение аппарата дробного исчисления, позволяет неявно включать дополнительные факторы взаимодействия физической системы. Несомненным преимуществом моделирования является тот фактор, что оно позволяет изучать недоступные или малодоступные для экспериментального исследования процессы с учетом фракталыюсти среды.
В настоящей работе проведено моделирование изменения заряда пузырьков в облачных каплях с учетом фракталыюсти среды, которые играют важную роль в процессе электризации облачных частиц.
Постановка и решение задачи. В процессе увеличения пузырька за счет диффузионных механизмов, важным моментом в облачных каплях является случай утечки газа в общем балансе газа в капле. Утечка газа бурно протекает в тот момент, когда переохлажденная облачная капля радиуса r сталкивается с градиной [2, 3] в результате
которого, выходя на поверхность и разрушаясь, уносят с собой определенное количество заряда, определяемая формулой qa = 4п£0£а,
Известно [3], что средний заряд qr, который создается одной облачной каплей, радиуса г за счет пузырьков содержащихся в пей определяется в виде
qr = 4п£0£иа , (1)
где £0 - электрическая постоянная, а - радиус пузырька, £ - электрокинетический по-иа
В формуле (1) вводя переменную t, перепишем заряд одного пузырька в капле в следующей форме
qr (t) = 4п£о£а (t) , (2)
а изменение заряда пузырька образующегося в облачной капле примет вид
д д
-qr(t) = ^s^-a(t) , (3)
Учитывая фракталыюсть среды, (3) можно записать в виде
A50Xi) = 47re0e^«W , (4)
где А = const. > 0, д^а (t) = - регуляризовахшая дробная производная порядка
а от функции а (t) с началом и концом в точках 0 и t (производная по Капуто) [4], Dax - оператор интегро-дифференцирования в смысле Римана-Лиувилля дробного порядка а началом в точке а £ [A, B] , которая определяется следующим образом [5]:
Daaxu (t)
sign (x — a) [ и (t) dt
a < 0,
a+1 '
Г (—a) J |x — t|
a
и (t), a = 0, sign (x ~ a) ^[a-j + l Dax и (t) , OL > 0.
где символ sign z определяется равенствами signO = 0 , signz = z/ |z|, z = 0.
Замена d/dt на d0t или h a D0 в дифференциальных уравнениях, неявно включает дополнительные факторы взаимодействия физической системы. Поэтому можно утверждать, что уравнение (4) описывает фрактальный процесс, d p2
С учетом, что —a (t) = 2, 4ф-a (t). где а - размер п\'зырька. tp - пересыщение, р -
dt p0a
давление окружающей среды, p0 - атмосферное давление, а - поверхностное натяжение воды, перепишем (4)
\d%tqr (t) = 4тге0? • 2, 40 — a (i) , (5)
Poa
Обозначим 4пеоС ' 2, 4фр2/\р0а = В (в случае расчета полного заряда всех пузырьков образующихся в капле, B умножается на n - количество пузырьков радиуса а), тогда
(5) примет вид
Ö0tQr (t) = Ва (t) , (6)
К уравнению (6) добавим начальное условие
qr (0) = qo , (7)
да (t)
Решим задачу (6), (7). Для этого используя формулу d£ta(t) = DQt_1——— перепишем
(6) в виде
Dt-1 AU- (t) = Ва (t) , (8)
Применив обобщенную формулу Ньютона-Лейбница к уравнению (8) получим
|t - 01 -а+1-1 t-a
Daot-l+lqr (t) ~ г 1_J + 1_1) Й Dltlqr (t) = Daotqr (t) - r(1_Q^oV (t) =
t—a
= Daota (t) - —--q0 = Ba (t) , (9)
Г (1 — а)
Подействуем на обе части уравнения (9) оператором D—ta тогда получим следующее выражение
D—TDaotqr (t) = D
—о t
t-
+ В а (t)
Г(1 — а)
в которой поело применения обобщенной формулы Ньютона-Лейбница имеем
(10)
t 0 -1 t® - г м-u uDoo1(lr (t) = Dq"——--q0 + ZV [Ba (t)] , (11)
Г (1 + а — 1) Г (1 — а)
D -1 qr (t) = 0
t-
qr (t) = Dq" ——--qo + Düa [Ba (t)] , (12)
Г (1 — а)
Решение правой части уравнения (12) выглядит следующим образом
t t D-a t~a л I 1 f T~aqodT I 1 fp a{t)dT
Dot "F71-+ Dot [Ba (t)\ - —— / —--—-+ —— / В—--
Г(1 — аГ ' ut L WJ Г (а) 0 Г (1 — а) |t — т|1—а Г (а) J 11 — т|]
t t 90 ^ r~a (t - r)a~l dr + [ a (t) (t - r)a~l dr
Г (а) Г (1 — а) Г (а)
о 0
t
qo + f^T [ a (t) (t - rf-1 dr , (13)
Г (а) J о
где
т
(t - т)а-1 dT
s = -ь t'
т = ts dT = tds
t-as-a (t - st)a-1 ds = в (1 - a, a)
Г (a) Г (1 - a)
Г (1 - a + a)
Решение задачи (6), (7) в итоге принимает вид
qr (t) = qo + BD-taa (t)
Г (a) Г (1 - a) .
(14)
Уравнение (14) представляет собой обобщенное уравнение закона изменения заряда пузырьков создаваемое облачной каплей с учетом фракталыюети среды. Решение интегрального уравнения (14) имеет вид |6|
qr (t) = qo + Bao I (t - т)a-1 Ea,a (B (t - т)a) dT
ж
Вг (t — т)a
где с учетом Еаа (В (t. — т)а) = У / -¡^———— ~ функции тина Миттаг-Леффлера,
^ Г (ai + a)
i=о у '
окончательно (14) примет вид
(t - т г1 y,
о i=0
-1 ^ Bi (t - т)a
Г (ai + a)
Анализ результатов математического моделирования. Математическое моделирование проводилось с помощью математического пакета символьной математики Wolfram Mathematica 9. На рис. 1 приведены расчетные кривые изменения относитсль-
qr ( t)
ta
t
t
Рис. 1. Расчетные кривые определяющая область изменения дг (¿) в зависимости от параметров ¿и а, полученные согласно формуле (14): а=0,2 (кривая 1); а=0,4 (кривая 2); а= 0,6 (кривая 3); а=0,8 (кривая 4).
Видно, что при достаточно малых значениях а расчётные кривые перегруппировываются с бесконечно длинными «степенными хвостами». Можно также заметить, что «степенные хвосты» указывают на нелинейность фрактальных процессов в облачной среде. Полученные данные хорошо согласовываются с данными Ирибарпе и Мейсе-па |7| в котором па один пузырек радиусом более 0,1 мм образуются заряды порядка 10-14 - 10-15 Кл.
Заключение. Одним из примеров фрактальных объектов в природе являются облака. Облака па самом доле являются самоафиппыми фракталами, т.к. атмосфера стратифицирована. Это означает, что вертикальное направление неравномерно с горизонтальным, следовательно, облака не могут быть самоподобны но форме. Поэтому важным моментом является изучение процессов внутри облака с учетом влияния фракталыюети среды, другими словами фрактальные геонроцессы.
Облака имеют разную структуру, и имеют свою классификацию но происхождению и морфологическим признакам, к которым, добавляя данные об их фрактальной структуре и фрактальных процессах, в дальнейшем будет возможным формирование более общей картины состояния физики облаков.
В работе предложена математическая модель изменения заряда пузырьков образующиеся в облачных каплях с учетом фракталыюети среды. Получено решение этой модели с использованием аппарата дробного исчисления.
Дальнейшие исследования с помощью разработанной модели будут направлены па более глубокое изучение процессов электризации и пространственного разделения зарядов в грозовых облаках с учетом взаимодействия процессов внутри облака и фракталыюети среды.
Литература
1. Proceedings of the Sixth Trieste International Symposium on Fractals in Phvsies / Edited bv L.Pietronero, E. Tosatti / ICTP, Trieste, Italy, 198. P.644-649.
2. Кумыков T.C. Жекамухов M.K., Каров Б.Г. Электризация и пространственное разделение зарядов при выделении пузырьков воздуха в процессе коахуляционших) роста 1'радин в облаке. I. Кинетика процесса выделения пузырьков при повышении температуры переохлажденных облачных капель /7 Метеорология и Гидрология. 2008. №11. С.44-52.
3. Кумыков Т.С. Жекамухов М.К., Каров Б.Г. Электризация и пространственное разделение зарядов при выделении пузырьков воздуха в процессе коахуляционших) роста 1'радин в облаке II. Генерирование хрозовохх) электричества за счет выделения заряженных пузырьков при намерзании переохлажденных облачных капель на поверхности 1'радин /'/' Метеорология и Гидрология. 2008. №12. С.15-24.
4. Caputo М., Elastieita dc dissipazione, Zaniehelli / Bologna, Italy, (Links), 1969.
5. Нахушев A.M. Уравнения математической биологии / M: Высшая школа, 1995. 301 с.
6. Пеху А.В. Краевые задачи для дифференциальных уравнений с частными производными дробжнх) и континуальнохх) порядка / Нальчик: Издательство КБНЦ РАН, 2005. 185 е., С.22.
7. Iribarne J.V., Mason B..J. Electrication accompanying the bursting of bubbles in water and dilute aqueous solutions /7 Trans. Faraday Soc. 1967. 63. №537. P.143-151.
MATHEMATICAL MODELLING OF BUBBLE CHARGE CHANGE IN DROPLET CLOUDS WITH ACCOUNT OF ENVIRONMENT FRACTALITY
T.S. Kumykov
Institution of Applied Mathematics and Automation, Shortanova St., 89-A, Nalchik, 360004, Russia, e-mail: macist20@mail.rul
Abstract. The paper proposes a new model of the change charge bubbles in supercooled cloud droplets. The model takes into account, the fractal properties of clouds, and the solution was obtained with the application of the fractional calculus.
Key words: fractal dimension, mathematical model, cloud droplet, bubble.
