Моделирование пространственно однородного процесса коагуляции для больших систем

Явление коагуляции (или слияния) частиц есть одна из основных причин эволюции дисперсных систем, под которыми понимают механическую смесь среды (газообразной или жидкой) с частицами (твердыми или жидкими). Это явление наблюдается в различных физических ситуациях: в растворах — броуновская коагуляция, при образовании капель дождя — гравитационная коагуляция. Процесс коагуляции оказывает воздействие на рост кристаллов в растворах, на рост газовых пузырьков и пор в твердом теле. Серьезное влияние оказывает коагуляция продуктов горения топлива на тягу реактивных двигателей (см. [7]).

1. Физические аспекты процесса коагуляции.

Остановимся более подробно на описании физических аспектов изучаемого явления коагуляции. Подробный анализ механизмов, приводящих к движению и столкновению частиц в дисперсных системах, приведен в работе [7], выдержка из которой послужит кратким обзором возможных физических причин коагуляции в дисперсных системах.

Коагуляция аэрозольных частиц играет важную роль в протекании различных атмосферных процессов, включая эволюцию спектра тонкодисперсной аэрозольной компоненты и образование осадков. Поэтому при исследовании физики атмосферы необходимо иметь довольно четкое представление о физических механизмах, приводящих к коагуляции или препятствующих ей, правильно оценивать влияние различных факторов на взаимодействие частиц при отдельных актах коагуляции и уметь рассчитать с требуемой точностью эволюцию функции распределения по размерам в результате коагуляции.

Рассмотрим возможные механизмы, способствующие или препятствующие коагуляции частиц в атмосфере. Надо, однако, иметь в виду, что они зачастую действуют совместно.

Броуновская коагуляция. Мелкие аэрозольные частицы (для нижней тропосферы Земли при Я < 1 мкм) реагируют на случайные молекулярные флуктуации плотности и средней скорости молекул воздуха, поэтому пребывают все время в нерегулярном броуновском) движении. Броуновское блуждание приводит к их взаимному столкновению, поэтому является одним из основных, постоянно действующих механизмов, способствующих коагуляции мелких аэрозольных частиц.

Броуновская диффузия. Броуновское блуждание частиц приводит к так называемой броуновской диффузии и способствует также их осаждению на более крупные предметы, например на растения, деревья и т. д. Этот процесс является одной из основных причин сухого вымывания тонкодисперсной аэрозольной компоненты из атмосферы. Естественно, что броуновская диффузия возможна и на крупные атмосферные частицы, физически не находящиеся в состоянии броуновского блуждания, например, на крупные частицы атмосферной пыли, частицы облаков и осадков. Однако, вследствие седиментации крупных частиц в атмосфере последние обтекаются воздухом. Следовательно, мы имеем дело одновременно с конвективным переносом мелких аэрозольных частиц и броуновской диффузией, т. е. с конвективной броуновской диффузией. Этот процесс является одной из основных причин влажного вымывания тонкодисперсной аэрозольной компоненты из атмосферы.

Эффект «зацепления». Конвективный перенос мелких аэрозольных частичек в окрестности падающей крупной частицы сам по себе, без влияния других причин, практически не может привести к захвату. Это связано с тем, что переменная составляющая скорости воздуха на поверхности этой частицы исчезает. Следовательно, безынерционная неброуновская частичка должна, вообще говоря, крупную частицы обойти таким же образом, как ее обтекает воздух. Однако частички обладают конечным размером. И поскольку на расстоянии порядка радиуса частички от поверхности крупной частицы нормальная составляющая скорости воздуха конечна, при конвективном переносе возможно «зацепление» мелкой частички за крупную. Эффект зацепления играет основную роль при влажном вымывании каплями облаков и осадков мелких аэрозольных частиц в диапазоне размеров, где конвективная броуновская диффузия уже неэффективна, а влияние инерции на осаждение еще очень мало.

Влияние седиментации осаждающихся частичек. Падение частичек под влиянием силы тяжести (седиментация) может приводить как к уменьшению, так и к увеличению их осаждения на закрепленные тела. Это зависит, естественно, от направления действия силы тяжести. Если же и крупные частицы и частички находятся в состоянии падения, то сила тяжести всегда приводит к уменьшению осаждения, поскольку она направлена против конвективного переноса частичек воздухом, обтекающим крупную частицу. Влияние седиментации осаждающихся частичек сравнимо по порядку величины с влиянием эффекта зацепления.

Инерционное осаждение. Поскольку гидродинамическое поле крупной частицы неоднородно, на движение мелких частичек в этом поле всегда в той или иной мере будет оказывать влияние их инерция. Влияние инерции проявляется в уменьшении в некоторых областях течения кривизны траекторий частичек по сравнению с кривизной линий тока среды. Возможны два различных режима движения частичек под влиянием инерции в неоднородном поле среды: докритический, когда траектории частичек и линий тока среды в окрестности поверхности крупной частицы не совпадают, но их поведение аналогично, и сверхкритический, когда влияние инерции столь велико, что траектории частиц пересекают поверхность крупных частиц. В первом случае инерция частиц либо усиливает, либо уменьшает действие других механизмов коагуляции, например, эффекта зацепления. При сверхкритическом режиме мы имеем дело с новым механизмом коагуляции, действующим самостоятельно, — инерционным осаждением.

Все три рассмотренных выше процесса действуют совместно при гравитационной коагуляции падающих в атмосфере частиц сравнимых размеров. Гравитационная коагуляция является одним из основных микрофизических механизмов образования осадков.

Анализ характера протекания коагуляционных процессов в среде (в особенности, для частиц сравнимых размеров) сильно осложнен из-за влияния гидродинамического взаимодействия. При вращении или поступательном движении частиц в среде вокруг них образуются неоднородные гидродинамические поля. Взаимное искажение гидродинамических полей, в особенности при сближении частиц, и приводит к появлению гидродинамического взаимодействия.

Силы гидродинамического взаимодействия пары частиц при малых скоростях движения (малых числах Рейнольдса) препятствует сближению частиц, при больших числах Рейнольдса и определенных углах между линией центров и направлением движения способствует ему. Заметим, что силы гидродинамического взаимодействия, по крайней мере для малых чисел Рейнольдса, неограниченно возрастают при сближении частиц. Это явление. приводит к гидродинамическому «запрету» коагуляции.

Следовательно, в каждом конкретном случае необходим строгий физико-химический анализ эффектов, возникающих на малых расстояниях между поверхностями частиц и способствующих снятию гидродинамического запрета коагуляции.

Эффект «втягивания в гидродинамический след». При движении тела за ним всегда, кроме случая очень малых чисел Рейнольдса, образуется гидродинамический след. Режим течения в следе при разных числах Рейнольдса различен. Однако всегда среда в следе в той или иной мере увлекается падающим телом. Поэтому частица, попавшая в гидродинамический след движущегося в атмосфере тела, будет падать несколько быстрее, чем в невозмущенной атмосфере. Этот эффект должен приводить к гравитационной коагуляции частиц одинаковых масс, если они падают вдоль линии, близкой к линии центров. Даже в случае, когда в среде находится другая частица, с несколько меньшей массой, гравитационная коагуляция возможна. Если бы эффекта втягивания в гидродинамический след не существовало, в последнем случае частицы просто разошлись бы. Однако здесь становится нетривиальным вопрос о причинах появления таких частиц в гидродинамическом следе частиц с большей (или равной) массой, ведь на очень большом расстоянии взаимное влияние их гидродинамических полей на седиментацию исчезает. Поэтому роль эффекта втягивания в гидродинамический след в процессах коагуляции может проявляться в случае, когда легкие частицы образуются непосредственно в гидродинамическом следе тяжелой, либо в случае, когда в дисперсных системах существуют силы, забрасывающие легкие частицы в гидродинамический след тяжелой.

Электростатическая коагуляция имеет место, естественно, когда существуют частицы, несущие заряды противоположных знаков. Поскольку электризация частиц в атмосфере — довольно обычное явление, электростатическая коагуляция может оказывать определенное влияние на эволюцию аэрозольной компоненты атмосферы. Анализ электростатической коагуляции частиц в атмосфере сильно усложняется по следующим двум причинам. Первая — из-за седиментации частиц мы, как правило, всегда имеем дело одновременно и с гравитационной коагуляцией (или, в предельных случаях, с каким-нибудь из перечисленных выше первичных процессов, ее определяющих). Вторая — необходимость учета влияния индукционных (зеркальных) сил, возникающих всегда, когда перемещение частиц в зарядах возможно.

Индукционные силы, независимо от знаков зарядов частиц, всегда способствуют коагуляции. Под их влиянием даже одноименно заряженные частицы могут коагулировать. Однако в последнем случае всегда существуют две различных области: внешняя, где частицы отталкиваются, и внутренняя (зона доминирующего влияния индукционных сил), где они притягиваются. Поэтому, как и эффект втягивания в гидродинамический след, роль индукционных сил в процессах коагуляции одноименно заряженных частиц может проявляться либо в случае, когда частицы образуются непосредственно в зоне доминирующего влияния индукционных сил, либо в случае, когда в дисперсных системах существуют силы, забрасывающие их в эту зону. Влияние индукционных сил уменьшается, если по каким-нибудь причинам либо перемещение зарядов в частицах затруднено, либо число носителей элементарных зарядов в частицах сравнительно мало.

Для сравнительно плотно упакованной дисперсной системы электростатическая коагуляция может быть затруднена из-за влияния экранировки.

Влияние турбулентности на коагуляцию частиц в дисперсных системах в различных ситуациях может сказываться совершенно различным образом, однако, по-видимому, всегда турбулентность способствует коагуляции. Если внутренний (колмогоровский) масштаб турбулентности Я меньше или сравним с размерами частиц, то мы имеем дело с турбулентным блужданием частиц, которое аналогично броуновскому. Это блуждание будет приводить к взаимному столкновению частиц, т. е. к непосредственной турбулентной коагуляции. Необходимо, однако, отметить, что в атмосферных условиях Земли внутренний масштаб турбулентности обычно равен 0,1 см, поэтому непосредственная турбулентная коагуляция будет иметь место только для очень крупных частиц, например хлопьев снега.

Влияние турбулентности на эволюцию тонкодисперсной аэрозольной компоненты осадков в обычных атмосферных условиях может проходить по иным каналам. Размеры частиц в этом случае значительно меньше X. Из-за действия молекулярной вязкости среды в областях, размеры которых меньше Я, будет происходить ламинаризация течения. Существование турбулентных флуктуаций с размерами, меньшими Я, возможно, но мало вероятно. В этом случае имеет смысл говорить не о непосредственной турбулентной коагуляции частиц, а об усилении коагуляционных процессов в дисперсных системах, вызванных другими причинами. Такое усиление будет происходить, с одной стороны по причине интенсивного турбулентного перемешивания частиц на расстояниях, больших Я. При этом сближение частиц на масштабах, меньших Я, можно рассматривать, считая течение ламинарным. Характеристики турбулентности среды будут входить во внешние граничные условия. С другой стороны, турбулентность приводит к появлению дополнительной регулярной неоднородности гидродинамического поля на масштабах, меньших Я. Эта внешняя неоднородность может эффективно влиять на сближение частиц в ламинарных областях течения.

Заметим, что турбулентные пульсации могут способствовать втягиванию частиц в гидродинамический след или в зону действия индукционных сил в случае одноименно заряженных частиц. В плотно упакованной дисперсной системе они могут способствовать электростатической коагуляции также путем разрушения экранировки. Турбулентность атмосферы приводит к уменьшению размеров гидродинамического следа достаточно крупной частицы, т. е. к увеличению установившейся скорости ее падения. Этот эффект может существенно увеличить эффективность гравитационной коагуляции частиц осадков с облачными.

Интересные эффекты, способствующие коагуляции частиц, возникают в среде, температура которой неоднородна." .

1. Алексеева И. В., Кабанов A.C. Численная модель гомогенного зарождения капельной фазы с учетом коагуляции частиц // Коллоидный журнал, 1990, т. 52, № 2, с. 227−234.

2. Арнольд В. И. Обыкновенные дифференциальные уравнения. -М.: Наука, 1971.

3. Багдасарова И. Р., Галкин В. А. Моделирование периодических структур в распределении дефектов, возникающих в конструкционных материалах ЯЭУ, под действием стационарного источника // Известия вузов, Ядерная энергетика, 1999, № 1, с. 85−93.

4. Багдасарова И. Р., Галкин В. А. Моделирование процесса коагуляции в пространственно однородном случае // Математическое моделирование, 1999, №, с.

5. Бегалишвили H.A., Енукашвили И. М. Интегральные уравнения кинетической теории коагуляции в смешанных атмосферных облаках//Труды ЗакНИГМИ, 1971, вып. 36(42), с. 15−21.

6. Волощук В. М., Седунов Ю. С. Процессы коагуляции в дисперсных системах. JL: Гидрометеоиздат, 1975.

7. Галкин В. А. Итерационный метод решения одного класса эволюционных уравнений, связанных с физической кинетикой //ЖВМиМФ, 1981, т. 21, № 2, с. 385−399.

8. Галкин В. А. О решении уравнения коагуляции с ядром Ф = ху II Метеорология и гидрология, 1984, № 5, с. 33−39.

9. Галкин В. А. О существовании и единственности решения уравнения коагуляции // Дифференциальные уравнения, 1977, т. 13, № 8, с. 1460−1470.

10. Галкин В. А. Об одном свойстве процесса коагуляции атмосферного аэрозоля // Метеорология и гидрология, 1983, № 12, с. 11−19.

11. Галкин В. А. Об устойчивости и стабилизации решения уравнения коагуляции // Дифференциальные уравнения, 1978, т. 14, № 10, с. 1863−1874.

12. Галкин В. А. Теория функциональных решений квазилинейных систем законов сохранения и ее приложения // Труды семинара им. И. Г. Петровского, 1997, № 20, с. 81−120.

13. Галкин В. А., Тупчиев В. А. Об асимптотическом поведении решения уравнения коагуляции // Труды Института Экспериментальной Метеорологии, 1978, вып. 1972, с. 31−41.

14. Годунов С. К., Рябенький B.C. Разностные схемы. М.: Наука, 1973.

15. Головин A.M. К вопросу о решении уравнения коагуляции дождевых капель с учетом конденсации // ДАН СССР, 1963, т. 148, № 6, с. 1290−1293.

16. Головин A.M. О кинетическом уравнении коагулирующих облачных капель // Изв. АН СССР, Сер. геофизика, 1963, № 10, с. 1571−1580.

17. Головин A.M. О спектре коагулирующих облачных капель // Изв. АН СССР, Сер. геофизика, 1963, № 9, с. 1438−1447.

18. Головин A.M. Решение уравнения коагуляции облачных капель в восходящем потоке воздуха // Изв. АН СССР, Сер. геофизика, 1963, № 5, с. 783−791.

19. Егоров В. И. К теории стационарной коагуляции пространственно неоднородных систем // Сообщения АН ГССР, 1979, № 3, с. 601−604.

20. Егоров В. И. Теоретические аспекты кинетики коагуляции в дисперсных системах // Сообщения АН ГССР, 1979, № 2, с. 325−328.

21. Енукашвили И. М. К теории конвективного переноса облачных частиц // Труды ЗакНИГМИ, 1971, вып. 41(47), с. 3−18.

22. Енукашвили И. М. О влиянии коагуляции и пространственного перераспределения облачных частиц на спектр выпадающих осадков // Труды Института геофизики АН ГССР, т. XXVIII, 1972, с. 189−195.

23. Енукашвили И. М., Бегалишвили H.A. Интегральное уравнение кинетической теории коагуляции облачных частиц с учетом ихконденсационного роста // Труды ЗакНИГМИ, 1971, вып. 36(42), с. 3−10.

24. Енукашвили И. М., Бегалишвили H.A. О численном решении нелинейных кинетических уравнений коагуляции методом Монте-Карло для пространственно-неоднородной стационарной задачи // Труды ЗакНИГМИ, 1974, вып. 55(61), с. 30−41.

25. Енукашвили И. М., Бегалишвили H.A. Об одном алгоритме численного решения интегральных кинетических уравнений коагуляции // Труды ЗакНИГМИ, 1974, вып. 55(61), с. 3−11.

26. Енукашвили И. М., Джапаридзе Н. Д., Надибаидзе Г. А., Бегалишвили H.A. О численном решении методом Монте-Карло линеаризованных кинетических уравнений коагуляции с учетом дробления облачных капель // Труды ЗакНИГМИ, 1975, вып. 63(69), с. 124−134.

27. Енукашвили ИМ., Джапаридзе Н. Д., Надибаидзе Г. А., Бегалишвили H.A. К вопросу численного моделирования лэнгмюровского цепного процесса в конвективных облаках методом Монте-Карло // Сообщения АН ГССР, 1974, вып. 76, № 1, с. 73−76.

28. Енукашвили И. М., Джапаридзе Н. Д., Цицвашвили Ш. И., Бегалишвили H.A. О численном моделировании методом Монте-Карло активных воздействий на кинетику осадкообразования в конвективных облаках // Труды ЗакНИГМИ, 1974, вып. 55(61), с. 12−29.

29. Кружков С. Н. К методам построения обобщенных решений задачи Коши для квазилинейного уравнения первого порядка // Успехи мат. Наук, 1965, т. 20, № 6, с. 112−118.

30. Кружков С. Н. Квазилинейные уравнения первого порядка со многими независимыми переменными // Мат. сб., 1970, т. 81, № 2, с. 228−255.

31. Сафронов B.C. Частный случай решения уравнения коагуляции // ДАН СССР, 1962, т. 147, № 1, с. 64−67.

32. Сафронов B.C. Эволюция допланетарного облака и образования Земли и планет. М.: Наука, 1969.

33. Туницкий Н. О коагуляции полудисперсных систем // Журнал экспериментальной и теоретической физики, 1938, т. 8, вып. 4, с. 418−424.

34. Тупчиев В. А. Об асимпротических свойствах решения уравнения коагуляции // Труды ИЭМ, 1971, вып. 23, с. 17−27.

35. Тупчиев В. А. Об асимптотических свойствах решения уравнения коагуляции // Труды Института Экспериментальной Метеорологии, 1971, вып. 23, с. 17−27.

36. Филиппов А. Ф. Обыкновенные дифференциальные уравнения с разрывными правыми частями // Мат. сб., 1960, т. 51, № 4, с. 101−128.

37. Хартман Ф. Обыкновенные дифференциальные уравнения. -М.: Мир, 1970.

38. Чандрасекар С. Стохастические проблемы в физике и астрономии. М.: Издательство иностранной литературы, 1947, 170 с.

39. Эдварде P.E. Функциональный анализ. Теория и приложения. -М.: Мир, 1969.

40. Эйнштейн А., Смолуховский М. Сб. статей «Броуновское движение». М.: ОНТИ, 1936, 607.

41. Klett James D. A Class of Solutions to the Steady-State, Source-Enhanced, Kinetic Coagulation Equation // Journal of the Atmospheric Sciences, 1975, v. 32, № 2, p. 380−389.

42. Lai F.S., Friedlander S.K., Pich J., Hidy G.M. The Self-Preserving Particle Size Distribution for Brownian Coagulation in the Free-Molecule Regim // J. Colloid and Interface Sei., 1972, v. 39, № 2, p. 395−405.

43. Lang Reihard, Xanh Nguen Xuan. «Smoluchovski's Theory of Coagulation in Colloids Holds Rigorously in the Boltzmann-Grad-Limit», Wahrscheinlichkeitstheorie.

44. Langevin P. Sur la theorie du mouvement brownien // Comptes Rendus, Pans, 1906, vol. 146, № 10, p. 530−533.

45. Lee K.W., Chen H. Coagulation Rate of Polidisperse Particles // Aerosol Science and Technology, 1984, v. 3, № 3, p. 327−334.

46. McLeod I.B. On the scalar transport equation I I Proc. London Math. Soc, 14, 1964, p. 445−458.

47. Melzak Z.A. A scalar transport education // Trans. Amer. Math. Soc., 1957, v. 85, p. 547−566.

48. Scott W.I. Analytical studies of cloud droplet coalescence // Journal of Atmospheric Sciences, 1968, v. 25, p. 54−65.