CC BY
44
УДК 551 .576: 551.594
ИССЛЕДОВАНИЕ НА ОСНОВЕ ЧИСЛЕННОЙ МОДЕЛИ ВЗАИМОДЕЙСТВИЯ МИКРОФИЗИЧЕСКИХ И ЭЛЕКТРИЧЕСКИХ ПРОЦЕССОВ В КОНВЕКТИВНЫХ ОБЛАКАХ ПРИ ИХ ЭВОЛЮЦИИ
© 2004 г. Е.А. Корчагина, А.В. Шаповалов
The two-dimension numerical model is used to research the interaction of microphysical and electrical processes in convective clouds during its micro-structuring. The results of numerical simulations of microstructural and electrical parameters’ evolution are given. The formation of precipitation in clouds regarding as well as negleching electrical phenomena is considered.
В исследованиях физических процессов в облаках и их взаимодействия в последние годы возрастает роль математического моделирования, которое в комплексе с лабораторными и натурными физическими экспериментами позволяет анализировать некоторые важные аспекты физических процессов в конвективных облаках.
Анализ теоретических и экспериментальных работ по физике облаков показывает, что еще мало изучено влияние электрических зарядов на частицах и электрического поля на протекающие в облаке микрофизиче-ские процессы и обратная связь между ними [1-5].
Целью данной работы является теоретическое исследование проблемы взаимодействия термодинамических, микрофизических и электрических процессов в конвективных облаках при их эволюции.
В качестве инструмента исследований выступает двумерная математическая модель конвективного облака с детальной термодинамикой, микрофизикой и учетом электрических процессов. Использование двух пространственных координат связано с существенным усложнением микрофизического блока в предложенной постановке. В реализованной авторами модели коэффициент коагуляции частиц не является постоянной величиной, а зависит от напряженности электростатического поля облака, рассчитываемой в каждый момент времени. Объем расчетов на каждом временном шаге значительно возрастает по сравнению со случаем использования коэффициента коагуляции частиц, зависящего только от их размера, поэтому применение трехмерной модели становится затруднительным.
В работе приведены результаты расчетов эволюции термодинамических, микроструктурных и электрических параметров конвективного облака. Рассмотрено образование осадков без электрических явлений и с их учетом. С помощью модели проведена количественная оценка влияния электрических процессов на формирование осадков, исследована положительная обратная связь между ансамблем частиц и электрическим полем облака.
Численная модель конвективного облака
Ниже представлено краткое описание разработанной авторами модели конвективного облака с детальным описанием термодинамических, микрофизических и электрических процессов. В отличие от имеющихся в настоящее время моделей [6-11] рассчитывались плотности объемных зарядов в облаке, потенциал и напряженность электрического поля, создаваемого этими зарядами, учитывалось влияние электрического поля облака на процессы слияния облачных частиц.
В основе предлагаемой модели конвективного облака с учетом электрических процессов лежат известные уравнения гидротермодинамики влажной конвекции [6,7,12]:
^ дп'
& дх
dw дп'
— =-------+ ЛЪ;
(1)
(
dt
du
dz
^ + 0,61q'- Qs
,00
A
d0
dt
L3
dw
- +------= aw ;
dx dz
Lk 0 f SM
Cp T
0 f SM3
St
+ Lç 0f SM,
Cp T
St
St
+ л'0;
(2)
(3)
SMk SMC
St
St
^'q;
u-
dx
w
dz:
d
+ СрТ
dt d д
где dt “ дt' д' = АкА + ^к- ,
дх дх д7 д7 u(x,z), w(x,z) - компоненты вектора скорости воздушных потоков в облаке; 9(х^) - потенциальная температура; л(х^)=ср 0 (p(x,z)/1000)R/Cp - безразмерное давление; 0 - средняя потенциальная температура; R-газовая постоянная; q(x,z) - удельная влажность воздуха; р8(х^) - суммарное отношение смеси жидкой и твердой фаз в облаке; ст^) - параметр, учитывающий изменение плотности воздуха с высотой; Р(х^) и Т(х^) - соответственно давление и температура; Ср -теплоемкость воздуха при постоянном давлении; Ьк, ЬС, Ьз - удельная теплота конденсации, сублимации и замерзания; я’(х^), 0’(х^), q,(x,z) - отклонения безразмерного давления, потенциальной температуры и удельной влажности от их фоновых значений в окружающей атмосфере Пф(х^), 0ф(х^) и qф(x,z); ( 5М к
Г 8 М С
I S t
I S t
изменения удельной влажности за счет
диффузии пара на капли и кристаллы;
f S МЗ
I S t
З I - масса
капельной воды, замерзающей в единицу времени в единице объема воздуха; К(х^) - коэффициент турбулентной диффузии. Для границ пространственной области используются обозначения 0,Ьх и 0,Ь;,.
Система (1)-(3) решается при следующих начальных и граничных условиях: и(х^,0)= и°(х^,0);
5в
w(x,z,0)=w (х^,0);
0(х^,0)= 0°(х^,0); (4)
и=и0^); 0=0o(z); p=Po(z); q=qo(z) при х=0,Ьх. (5)
u=w=0, 0=00(0), р=р0(0), q=q0(0) при z=0
u=u(Lz), w=0, 0=0o(Lz), p=po(Lz), q=qo(Lz) при z= Lz.
Предлагаемая модель отличается микрофизиче-ским блоком, который описывает процессы нуклеа-ции, конденсации, коагуляции капель с каплями, сублимации, аккреции, замерзания капель, осаждения облачных частиц в поле силы тяжести, их перенос воздушными потоками, а также взаимодействие облачных частиц под влиянием электрического поля облака. В отличие от моделей с параметризованной микрофизикой [9-11] нами рассматривается непрерывный спектр частиц в облаке. Система уравнений для функций распределения по массам капель ^(х^шД), ледяных частиц ^(х^шД) и осколков замерзания капель ^(х^шД) имеет следующий вид [13]:
—L + и —L + (w - V, ) —L = I
дt дх д z I
dt
KД
dfi
dt
Kr
dt
AK
+ ífI і + ЛЇ. +1,;
df-
dt
df
ДР
dt
(в)
df2
^ + u+ (w - Y2)-^ = | — I + | —
dt
dx
+fdf2
1 dt УЗ
dz
+ ЛІ2 + I2 + IAВ ;
df2
dt
df2
dt
AK
d f3 df3 df3
—— + u —— + (w - Y2 ) —— =
d t d x d z
d t
d t
+ ЛТ3
AK
где Vi(m), V2(m) - установившиеся скорости падения
жидких и твердых частиц;
dfl d t
KД
f d t
Kr
d t
d t
f
dt
- изменения функции
АК V д 1 У ДР V ^ у З
распределения капель за счет микрофизических процессов конденсации, коагуляции капель, аккреции капель и кристаллов, дробления и замерзания соот-
ветственно;
f
dt
f
dt
AK
dt
- измене-
ния функции распределения кристаллов за счет суб-
лимации, аккреции и замерзания капель;
d f3
f
dt
—— I - изменения функции распределения
д1: У АК
f3(r,ш,t) за счет образования осколков при спонтанном замерзании облачных капель и аккреции.
Для системы уравнений (6) используются следующие начальные и граничные условия:
^ (х^т^^^х^т^^^х^т^)^; (7)
fl(x,z,ш,t)= ^(х^,тД)= f3(x,z,ш,t)=0 при x=0,Lx;
f1(x,z,m,t)= f2(x,z,m,t)= f3(x,z,m,t)=0 при z=Lz ; (8)
df df2 д f3 „ 0
—L = —L = —3 = 0 при z=0.
dz dz dz
Для описания коагуляционных процессов в облаке применялось интегродифференциальное уравнение в виде, представленном в [7].
Расчет взаимодействия капель и кристаллов проводился на основе следующих соотношений: df ^ ”
—11 =-f1(x,z,m,t)Jp2(m,m')f2(x,z,m',t)dm'; (9) д / А 0
df ^
—~ I = -f2(x,z,m,t) f р2(m',m)f1(x,z,m',t)dm' +
dt / А 0
+J f1 (x,z,m',t)P2 (m',m - m')f2 (x,z,m - m',t)dm'; (10)
0
P2 (m, m') = n(r(m) + r(m '))2 • | V1 (m) - V2 (m ')| • E(m, m'), где r(m) и r(m’) - радиусы сталкивающихся частиц; E(m,m’) - коэффициент слияния, который рассчитывается по аппроксимационным формулам.
Изменение функции f1(x,z,m,t) вследствие дробления капель рассчитывалось по формуле [6]:
д t
= -f1 (x,z,m,t)P(m) + J P(m)Q(m,m')f1 (x,z,m',t)dm'. (11)
m
В этом выражении P(m) есть вероятность распада в единицу времени капли массой m; Q(m,m’) - вероятность образования капли массой m при распаде капли массой m’.
Расчет изменений функций f1 (x,z,m,t) и f2 (x,z,m,t) за счет кристаллизации капель проводился с использованием функции вероятности замерзания в единицу времени капли массой m в точке (x, z), обозначенной R (x,z,m) [6,12]:
dfL/) ^(^ff] -fi(x,z,m,t)R(x,z,m); (12)
R ( x, z, m ) = A exp [b ( Тм (m) - Тв (x, z) ) ] , (13) где А, В - параметры; Тм ( m ) - медианная температура замерзания капель массой m; Тв ( x, z ) - температура воздуха в точке (x,z).
Появление новых капель и ледяных частиц в естественных условиях учитывается с помощью функций I1 для капель и I2 - для кристаллов [12]
Слагаемые, описывающие испарение, конденсацию и сублимацию, имеют вид [7]:
d I f (r t) dm 1 ,
-----1 i, (r,m,t)------
dm í 1 dt
df2 I d f f t) dm I
it |С=-dm1 f2(r,m,t)_d^J.
(14)
(1З)
Более подробно выражения для слагаемых, входящих в уравнения модели (1) - (15), представлены в работах [6,7,12,13].
Численное моделирование процессов электризации конвективных облаков связано с трудностями принципиального характера, одной из которых является корректное описание микрофизических процессов электризации облачных частиц. При моделирова-
Зі
нии облакообразования с учетом разделения электрических зарядов большинство исследователей не учитывают или упрощают детальное описание процессов разделения зарядов [10,11].
Для исследования закономерностей взаимодействия поля и частиц в облаке нужна детализация этих процессов, что и предпринято в настоящей работе.
Как известно, в конвективных облаках одновременно реализуются ряд механизмов электризации облачных частиц. В разработанной модели учитываются процессы электризации облачных частиц на основе экспериментально исследованных закономерностей развития грозовой деятельности облаков и значений коэффициентов разделения зарядов, связанных с замерзанием капель воды, ростом крупы и градин и взаимодействием градин с кристалликами льда и переохлажденными каплями [3,14].
За счет микрофизических процессов замерзания капель и аккреции в облаке идет накопление отрицательного заряда на ледяных частицах. Одновременно формируется положительный заряд, состоящий из зарядов отдельных частиц - осколков замерзания капель.
Для замерзающих капель, диаметр которых больше 200 мкм, с достаточной точностью процесс электризации описывается выражением [14]
q(m)=alm, (16)
где т - масса замерзшей капли; а! - коэффициент пропорциональности, значение которого меняется в зависимости от содержания примесей в капле и температуры ее замерзания (а«3,5-10-1° Кл/г при Т= -8.. .-16 0С).
На крупных кристаллах, крупе и градинах накапливается электрический заряд за счет захвата переохлажденных капель. Заряд пропорционален массе, замерзшей на них воды. При этом коэффициент пропорциональности зависит от температуры растущей частицы, а также от концентрации и химического состава примесей в облачной воде и принимает значение от 10-10 до 10-8 Кл/г.
Образование осколков при замерзании капель учитывается следующим образом:
ЭТ3
dt
= | п(т,т')К(х,2,т',1)^(х,2,т',1)(1т', (17)
З т
где п(т,т') - число ледяных осколков массы т, образующихся при замерзании капель массы т'; Щх^т'Д) - вероятность замерзания капель массой т' в единицу времени.
Образование осколков в процессе аккреции ледяных кристаллов с облачными каплями описывается формулой
d f ^ w
I = J n(m,m')R2(x,z,m',t)dm';
dt У AK m
(18)
где
Я.2 (х,2,т'Д) = ^ (х,2,т',1) • | Р2 (т'^) (х^Д) -
0
число столкновений кристаллов массой £, с каплями массой т', приводящих к замерзанию последних и образованию осколков.
Число ледяных осколков п(т,т') определяется с помощью аппроксимаций экспериментальных зависимостей выбросов микрочастиц от размеров замер-
зающих капель. Для капель с г<75 мкм используются данные Лезема и Мейсона, приведенные в [3]. Для капель больших размеров применяются данные, полученные в работе [14].
Микроскопические осколки замерзания выносятся потоками в верхнюю часть облака, где образуется преимущественно положительный объемный заряд р+(х,2,1). Область сосредоточения отрицательно заряженных ледяных частиц образует зону преимущественно отрицательного объемного заряда р_(х,2,1).
Объемные заряды на временном шаге рассчитывались по формулам:
p- (x,z,t ) = a2 J mf2 (x,z,m,t )dm -p2 (x,z) :
(19)
р+ (х^Д) = а31 mfз (х,2,тД )с1т -р3 (х,2),
0
где а2 и а3 - коэффициенты разделения зарядов; р2(х^) и р3(х,7) - уменьшение объемных зарядов в результате тока проводимости атмосферы и разрядов.
Рассчитанные значения суммарных объемных зарядов рЭ(х,2,1) используются далее для определения потенциала И(х,2) создаваемого ими электрического поля. Для этого на каждом временном шаге (при фиксированном 1) решается уравнение Пуассона
рэ
d 2U
d 2U
дх2 &2 е0
при следующих граничных условиях
Ш
х=0,Ьх дг
—=0
dx
■ = 0
U = 0lz=0
(20)
(21)
E (>
,t Ь-(1
(22)
где 80 - диэлектрическая постоянная вакуума; рэ (х 2,1) = р+ (х z, 1) + р- (х z, 1) .
Напряженность электрического поля Е(х,2,1) определяется как градиент потенциала по формуле:
д и д И'
-----+ п2---
д х д 2
Значения напряженности электрического поля учитываются в модели при определении коэффициентов коагуляции облачных частиц согласно теоретическим и экспериментальным зависимостям Л.М.Левина и Н.В.Красногорской [4,5].
Влияние электрического поля на эффективность столкновения заряженных капель зависит от размеров частиц, знака и величины зарядов. Кроме эффекта поляризации капель, электрическое поле вызывает изменение скорости движения заряженных частиц, если их размеры достаточно малы. Чем меньше относительные скорости движения капель, тем больше импульс электрических сил. Поэтому влияние электрического поля тем значительнее, чем ближе размеры взаимодействующих частиц [4].
Для решения системы уравнений модели на временном интервале (1п, 1п+1) используются методы расщепления, прогонки, верхней блочной релаксации [15].
Алгоритм расчета уравнений переноса субстанций основан на прогонке с чередованием направлений (продольно-поперечная схема). Эллиптическое неоднородное уравнение для отклонения давления п' [7], полученное применением оператора дивергенции к уравнениям движения (1), с использованием уравне-
W
ния неразрывности (2), решается методом верхней блочной релаксации [15]. Этот же метод применяется при решении уравнения Пуассона для потенциала электрического поля (20).
Интегрирование системы уравнений модели проводилось по пространственно-временной области Б х Т : {0 < х < Ьх,0 < 2 < Ь2, т> 0} .
Размеры пространственной области составляли 50 км по горизонтали и 18 км по вертикали. Шаг сетки по оси х составлял 200 м, по оси 2 - 100 м, по времени - 10 с.
По массам облачных частиц использовалась сетка с переменным шагом, при этом градации задавались следующим образом [7]: щ=2,540-10 г; тп+1=1,414-тп; п е [1 : 60 ]. Для ледяных частиц использованы 75 категорий масс, которые перекрывают все размеры, включая град ^тах»5 см).
Исследование коагуляционного роста облачных частиц в электростатическом поле различной напряженности
Формирование осадков в облаке происходит при наличии в облаке достаточно крупных частиц (г > 50мкм), называемых зародышами осадков. Зародыши крупных капель не могут образоваться только в результате конденсации [2,3], поэтому для объяснения их роста привлекают различные механизмы, в том числе электрическую коагуляцию [3-5].
Для оценки влияния электрического поля на коагуляцию облачных частиц в данной работе были проведены численные эксперименты по моделированию образования осадков в пространственно однородной системе при наличии электрических полей различной интенсивности. Коэффициент коагуляции аппроксимировался в зависимости от напряженности электрического поля Е2 [4,5]. Некоторые результаты расчетов отображены на рис. 1, на котором представлена зависимость хода процесса коагуляции капель при различных значениях напряженности электрического поля в облаке. По оси ординат отложена водность по градациям радиусов, г/м3.
С увеличением напряженно -сти поля быстрее образуются зародыши осадков и спектр смещается в сторону крупных капель.
Из первоначального спектра облачных капель (кривая 1) образуются дождевые капли при отсутствии поля за 35 мин, а при напряженности Е2=1000 В/см за 23 мин, при этом мелкокапельная часть спектра коагулирует практически полностью. Таким образом, процесс осадкообразования сокращается в присутствии электрических полей приблизительно на одну треть общего
времени образования осадков. При больших значениях Е2 влияние поля может быть еще существеннее, если следовать зависимости коэффициента коагуляции от напряженности электрического поля [5].
Полученные в работе результаты, относящиеся к взаимодействию частиц, размером ниже известного в физике облаков предела в 19 мкм, согласуются с результатами, полученными в экспериментальных исследованиях [16]. В частности, в этих экспериментах наблюдалось заметное увеличение числа капель радиусами 25-45 мкм при высоковольтных разрядах, тогда как в отсутствие разрядов размер капель составлял 7-15 мкм.
Результаты расчетов полей термодинамических микроструктурных и электрических параметров конвективного облака
На основе уравнений разработанной модели (1)-(22) были выполнены расчеты по моделированию конвективных облаков при различных реальных стратификациях атмосферы (использовались данные зондирования атмосферы в аэропорту «Минеральные Воды» и на научно-исследовательском полигоне ВГИ «Кызбурун-2» в весенне-летний период).
Облако инициировалось заданием теплового импульса в начальный момент времени у поверхности земли с перегревом ДТ=1,0-1,5 °С.
Некоторые результаты моделирования термодинамических и микроструктурных параметров конвективных облаков представлены в [12], и они не противоречат данным полевых наблюдений.
Водность, г/м3
|°д(г)
Рис. 1. Изменение спектра облачных капель при коагуляции за 20 минут 1 - начальное распределение; 2 - в отсутствии электрического поля (Е2=0);
3 - в электрическом поле с напряженностью Е2=100 В/см; 4 - Е2=300 В/см;
5 - Е2=1000 В/см
Остановимся на других результатах исследований на основе модели облака.
Многие авторы высказывают мысль, что важную роль в облаках играет взаимодействие различных физических процессов [1-3, 6, 7]. В [12] исследовались некоторые виды взаимодействий процессов в облаках и оценена их роль, а в настоящей работе сделана попытка анализа одного из наиболее сложных видов взаимодействия: электрических и микрофизических с наличием положительной обратной связи между ними.
Формирование микроструктуры на начальной стадии развития облака определялось микрофизическими процессами конденсационного роста отдельных облачных капелек. По результатам расчетов облака состояли на начальной стадии (примерно до 20-й мин) из капель воды, позднее - из смеси капель и кристаллов. При отрицательных температурах воздуха наблюдалось смешанное фазовое состояние, при этом капли были переохлажденными.
Превращение мощных кучевых облаков в кучеводождевые сопровождалось оледенением их вершин, формированием частиц осадков, а также появлением значительных объемных электрических зарядов, которые вызывают молниевые разряды.
Накопление электрического заряда в облаках происходило вследствие замерзания капель и аккреции капель с ледяными частицами. Из-за разности скоростей падения осколков замерзания и более крупных частиц, крупы или градин, происходило пространственное разделение зарядов. В предвершинной части облака формировался положительный объемный заряд, ниже - отрицательный.
Электрическое строение облака на 20-й мин развития для одного из вариантов исходных данных (радиозондирования атмосферы) представлено на рис. 2. Интегральный положительный заряд в единице объема достигает значений 6-10-3 Кл км-3, отрицательный —1,1-10-2 Кл км-3, что согласуется с существующими представлениями [1-3].
10
5-
Z, км
__0.001
0.001 0.001 . ...-----------г
/:-ojöfl #=~
-0.003 г
;.0.003-
-0.005-
--0.003
Х,км
параметров для одного из облаков представлены на рис. 3.
Z, і^м
30 Х,км
—I-
20
Рис. 2. Распределение отрицательного (а), положительного (б) и суммарного (в) объемных зарядов в облаке на 20-й минуте развития, Кл/км3
На каждом временном шаге численно решалось уравнение для потенциала электрического поля, рассчитывались вертикальная Б2 и горизонтальная Е составляющие напряженности электрического поля. Некоторые результаты моделирования электрических
Рис. 3. Изолинии потенциала, В (а); горизонтальной (б) и вертикальной (в) составляющих электрического поля, В/см, на 20-й минуте эволюции облака
На 20-й мин развития облака (35 минута от начала расчета) потенциал электрического поля составляет величину 3,1-106 В (рис. 3а). Компоненты вектора напряженности поля имеют значения Бх «5 В/см (рис. 36), а Б2 «8 В/см (рис. 3е). С течением времени заряд в облаке и соответственно потенциал поля увеличиваются. По результатам расчетов облака на 40-й мин развития потенциал достигал значений 1,1-108 В, а напряженности - Бх « 100 В/см, Бг«240 В/см.
Максимальные значения потенциала электрического поля, полученные в расчетах, имели значения 300-500 МВ. Сравнение результатов, полученных в данной работе, с приведенными в [11], показало, что порядок величины потенциала совпадает.
Вследствие взаимодействия облачных частиц радиусами до 20 мкм, которое становится возможным при наличии зарядов на частицах и значительных электрических полей в облаке (Е > 50 В/см), в выполненных численных экспериментах наблюдалось ускорение процесса роста частиц осадков и увеличение электрического поля облака. На рис. 4 приведены поля вертикальной составляющей электрического поля для двух моментов времени (20 и 40 мин). Составляющие вектора напряженности электрического поля Б( и Бг при учете электрической коагуляции растут быстрее (рис. 4), чем без ее учета (рис. 3): на 40-й мин развития облака Б2 без учета электрической коагуляции составляла 240 В/см, а при учете - 650 В/см.
Сравнение времени образования осадков без учета и с учетом электрической коагуляции показало, что во втором случае осадки начинаются на 7-10 минут раньше.
Появляются изменения и в спектрах ледяных частиц - наблюдаются смещения максимумов в сторону крупных частиц.
Результаты расчетов напряженности электрического поля, полученные в данной работе, находятся в согласии с экспериментальными измерениями, представленными в [17], согласно которым в облаках имеются достаточно большие области со значениями Б порядка 1000 В/см. В этой же работе указывается на наличие нескольких областей с различными направлениями поля.
Рис. 4. Поле вертикальной составляющей напряженности электрического поля, В/см при учете процесса электрической коагуляции на 20-й (а) и 40-й (б) минутах эволюции облака
Было выполнено большое количество численных экспериментов, при различном термодинамическом состоянии атмосферы. Каждый эксперимент проводился без учета и с учетом электрической коагуляции.
В результате проведенных исследований подтверждена закономерность, что между микрофизическими, термодинамическими и электрическими процессами существует обратная связь (взаимодействие), т.е. взаимное влияние их друг на друга. В проведенных численных экспериментах получалось, что электрическое поле ускоряет рост частиц, с другой стороны при этом генерировалось большее количество электрического заряда, которое увеличивало само поле. Общая картина влияния электрического поля проявлялась в некотором ускорении прохождения облаком всех стадий развития, кроме начальной.
Применение двумерной модели, по мнению авторов, не нарушает общности полученных результатов, поскольку речь идет о взаимодействии микрофизики и электрического поля облака.
Заключение
Основные результаты работы:
1. Уточнены вопросы формирования микроструктуры конвективных облаков с учетом электризации облака и электрической коагуляции частиц. На основе численных экспериментов получило теоретическое подтверждение предположение о возможности быстрого образования зародышей осадков вследствие взаимодействия капель радиусом до 20 мкм при наличии электрических полей в облаках.
2. С помощью математического моделирования впервые проанализирован один из важных аспектов микрофизики конвективных облаков - взаимодействие микрофизических и электрических процессов.
Результаты численных экспериментов подтверждают существование положительной обратной связи между ростом частиц осадков в облаке и увеличением напряженности электростатического поля, которая заключается во взаимном влиянии их друг на друга. Анализ результатов моделирования позволил установить, что электрическое поле ускоряет рост частиц, при этом генерируется большее количество электрического заряда, которое увеличивает само поле. Общая картина взаимного влияния проявляется в ускорении прохождения облаком всех стадий развития, кроме начальной.
3. Из полученных в ходе расчетов данных следует, что взаимодействие микрофизических и электрических процессов в мощных конвективных облаках играет существенную роль в процессе осадкообразования. В частности, в грозо-градовых облаках период образования осадков сокращается на одну треть и более.
Таким образом, результаты исследований показали, что влиянием электрического поля и электрической коагуляции в мощных конвективных облаках пренебрегать нельзя (при разработке их моделей и концепций активного воздействия на них).
Следует отметить, что заметное влияние электрические поля оказывают на формирование микроструктуры, когда их напряженность достигает значений свыше 50 В/см.
Литература
1. Имянитов И.М. и др. Электричество облаков. Л., 1971.
2. Мазин И.П., Шметер С.М. Облака. Строение и физика образования. Л., 1983.
3. Мейсон Б.Дж. Физика облаков. Л., 1961.
4. Левин Л.М. // Тр. Эльбрусской высокогорной экспедиции. 1961. Т. 2. С. 5-42.
5. Красногорская Н.В. // Изв. АН СССР. ФАО. 1965. Т. 1. С. 339-345.
6. Ашабоков Б.А., Калажоков ХХ. Численное моделирование градовых облаков. М., 1992.
7. Коган Е.Л. и др. Численное моделирование облаков. М., 1984.
8. Хворостьянов В.И. // Изв. РАН. ФАО. 1994. Т. 30. № 4. С. 543-557.
9. Chen J.P., Lamb D. // J. Atmos. Sci., 1994, Vol. 51. P. 2613-2630.
10. Rawlins F. // Quart.Jour. of the Royal Met. Society. 1982. Vol. 108. P. 779-801.
11. Helsdon J. H., Jr., Farley R. D. // J. Geoph. Res. 1987. Vol. 92. P. 5661-5676.
12. Корчагина Е.А. и др. // Информационные системы и технологии: Межвед. сб. Вып. 1. Нальчик, 2000. С. 10-17.
13. Шаповалов А.В. // Мат. моделирование. 2003. Т. 15. № 4. C. 65-76.
14. Аджиев А.Х., Тамазов С.Т. // Метеорология и гидрология. 1987. № 7. С. 57-62.
15. Марчук Г.И. Методы вычислительной математики. М., 1977.
16. Першина Т.А. и др. // Тр. НИЦ ДЗА. Вып.1 (546). 1997. С. 62-67.
17. МашуковХ.М. и др. // Всерос. конф. по физике облаков и активным воздействиям на гидромет. процессы. Нальчик, 2001. С. 46-47.
ГУ«Высокогорный геофизический институт» Росгидромета, г. Нальчик___________________4 марта 2004 г.
