Рус Eng Cn Перевести страницу на:  
Please select your language to translate the article


You can just close the window to don't translate
Библиотека
ваш профиль

Вернуться к содержанию

Программные системы и вычислительные методы
Правильная ссылка на статью:

Апробация модели капельной фазы в эксперименте по моделированию образования кислотосодержащего облака

Пекунов Владимир Викторович

доктор технических наук

Инженер-программист, ОАО "Информатика"

153000, Россия, Ивановская область, г. Иваново, ул. Ташкентская, 90

Pekunov Vladimir Viktorovich

Doctor of Technical Science

Software Engineer, JSC "Informatika"

153000, Russia, Ivanovskaya oblast', g. Ivanovo, ul. Tashkentskaya, 90

pekunov@mail.ru
Другие публикации этого автора
 

 

DOI:

10.7256/2454-0714.2021.1.35104

Дата направления статьи в редакцию:

22-02-2021


Дата публикации:

10-05-2021


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


Ключевые слова:

численное моделирование, облако, кислотные дожди, водность, распределение капель, математическая модель, апробация, интерполяционные соотношения, конденсация, столкновения

Abstract: The problem of numerical modeling of the formation (as a result of condensation growth and droplet collisions) and development of primary acidulous cloud considers various factors: the presence of temperature gradients, turbulence, direct solar radiation heating the air and walls of buildings, diffuse solar radiation (which describes radiation cooling), transfer of gaseous pollutants and their absorption by droplets. The author earlier formulated the corresponding complex mathematical model that takes into account the aforementioned factors. This article sets the task of testing the droplet component of this model through numerical modeling of the processes in the emerging cloud, with subsequent comparison of the results with theoretical and empirical correlations. The author obtained the new results of numerical modeling of acidulous cloud in the air over a vast urban area with high-density development  on the basis of the comprehensive mathematical model that takes into account the above listed factors and relies on the interpolation-sectional submodel of droplet phase. The author models the dynamics and kinetics of such cloud that absorbs gaseous sulphur dioxide; and obtains results on the intensity of absorption of this pollutant in the forming cloud. The comparison of these results with the known data (Hrgian-Mazin droplet distribution and interpolation ratio for the water level of the cloud) demonstrated quite a coincidence of droplet distribution and water level of the cloud. The conclusion is made on sufficient adequacy of application of the ecological model that includes a special submodel of droplet phase.  


Keywords:

numerical simulation, cloud, acid rain, water content, droplet’s size distribution, mathematical model, approbation, interpolation relations, condensation, coalescence

Введение

Изучение закономерностей возникновения кислотных осадков является важной экологической задачей, точное решение которой почти невозможно без привлечения математического моделирования. Данная задача сводится к двум основным подзадачам: моделированию образования и развития первичного кислотосодержащего облака и моделированию его эволюции вплоть до выпадения осадков. Для решения данных подзадач автором ранее была предложена комплексная многофазная математическая модель, включающая подмодель капельной фазы с учетом конденсации, испарения, коалесценции капель, поглощения и высвобождения ими газообразных загрязнителей [1, 10]. Определим место данной подмодели среди известных подходов, выделив их достоинства и недостатки:

1. Дискретные модели [2], к которым относятся варианты лагранжева подхода с секционированием, а также метод крупных частиц (упоминается, например, в работе [3]). Это достаточно простой подход, но для получения адекватной точности требуется моделирование большого количества капель, что дает очень высокую трудоемкость расчета.

2. Эйлерово-эйлеровы модели (модели сплошных сред). Выделим пять основных подходов: а) применение фиксированных распределений капель (дает недостаточную точность [2]); б) подход А.Е.Алояна [4] с решением кинетического уравнения (имеет высокую трудоемкость); в) атмосферные модели (весьма приближенные, имеющие много допущений о характерных распределениях капель, см., например, [5]); г) метод Хилла (метод моментов [6, 7]), в котором рассматриваются лишь интегральные характеристики распределений (сами распределения явно не рассчитываются), при этом есть серьезное допущение – капли не имеют собственной скорости [8]; д) композитные методы, использующие элементы метода моментов для расчета локальных распределений в каждой расчетной ячейке. Здесь можно выделить метод аппроксимации распределений суммой δ-функций [7, 9], а также предложенный автором ранее метод на базе поиска кусочно-сплайновых распределений [10], используемый в данной работе. Эти методы отличаются сравнительно невысокой трудоемкостью и позволяют аппроксимировать произвольные, в том числе разрывные, распределения. Используемый нами метод [1, 10] позволяет ввести собственные скорости капель, учитывает конденсацию, испарение и столкновения капель, и является достаточно точным [1, 10] и универсальным.

Используемая в данной работе математическая модель ранее была апробирована на задаче о поглощении газообразных загрязнителей туманом [1], причем результаты представлялись качественно вполне правдоподобными. Однако нельзя не признать необходимость дополнительной апробации со сравнением результатов моделирования с некоторыми эмпирическими данными, например, о виде распределений капель в формирующихся облаках [11] и их водности [12], при этом интерес представляет протекание этих процессов именно в присутствии газообразных загрязнителей, например, диоксида серы, источником которого являются автомагистрали или предприятия.

Таким образом, целью данной работы является дополнительная апробация капельного компонента комплексной математической модели [1, 10]. Для достижения данной цели поставим следующие задачи: а) провести численный эксперимент по моделированию образования первичного кислотосодержащего облака в присутствии газообразного диоксида серы, б) сравнить полученные результаты с известными эмпирическими данными о распределении капель и водности облаков, в) сделать выводы о характере поглощения газообразного диоксида серы облаком

1. Постановка эксперимента

С применением математической модели [1, 10] был поставлен численный эксперимент в трехмерной расчетной области высотой 272 м над участком города Ганновер, Германия. Был задан вертикальный градиент температур от 27,7 градусов Цельсия на нижней границе до 22,7 градусов на верхней. В воздухе присутствовал газообразный диоксид серы с концентрацией 8,8×10-9 моль/см3 (200 млн‑1). Концентрация водяного пара была определена в размере 1,33×10-6 моль/см3, что соответствовало 100-процентной влажности на половине высоты области. Изначально в воздухе присутствовали ядра конденсации диаметром 1 микрон. Капельная фаза была поделена на семь секций: [10-6; 10-5), [10-5; 2×10-5), [2×10-5; 6×10-5), [6×10-5; 10-4), [10‑4; 5×10‑4), [5×10-4; 10-3), [10-3; 5×10-3) м.

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

2. Апробация

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

На рис. 1 показана кривая экспериментально полученного распределения капель [ед/см3/микрон] по радиусам r [микрон] (на высоте 192 м) в сравнении с аналитической кривой nаналит(r), формулы для которой приведены в работе [11]:

где - плотность воды [кг/м3], - плотность i-го компонента капельной фазы [кг/м3], n(r) – экспериментально полученная функция распределения капель по радиусам [ед/см3/микрон], rнач и rкон – минимальный и максимальный из рассматриваемых радиусов капель [микрон].

Из рис. 1 достаточно хорошо видно, что полученное распределение капель вполне соответствует ожидаемому, что подтверждает адекватность рассматриваемой математической подмодели капель.

Рис. 1. Распределение капель (на высоте 192 м) по радиусам (диаметры указаны в микронах). Сплошная линия – аналитическая кривая, пунктирная – экспериментальная кривая

На рис. 2 показана кривая экспериментально полученной водности облака [г/м3] (в зависимости от высоты h [м] над основанием облака) в сравнении с аналитической кривой wаналит(), формулы для которой приведены в работе [12]:

где H – высота облака, которая в наших экспериментах была подобрана из условия наилучшего соответствия экспериментальной и аналитической кривых водности и составила H = 370 м, а величина (отношение высоты точки с максимальной водностью к общей высоте облака) для большинства облаков [12].

Из рис. 2 видно, что характер изменения водности с высотой вполне соответствует ожидаемому, что дополнительно подтверждает адекватность примененной подмодели капельной фазы.

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

причем, при использовании данной зависимости для экстраполяции, вполне вероятен выбор значения m = 0. Хорошие результаты были получены при p = -1/8, q = 0,3.

Рис. 2. Водность облака (г/м3, абсциссы) в зависимости от высоты над основанием облака (м, ординаты). Сплошная линия – аналитическая кривая, пунктирная – экспериментальная кривая

Выводы

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

Библиография
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.

Результаты процедуры рецензирования статьи

В связи с политикой двойного слепого рецензирования личность рецензента не раскрывается.
Со списком рецензентов издательства можно ознакомиться здесь.

Математический инструментарий и методы математического моделирования находят применение в различных сферах научных исследований. Рецензируемая статья посвящена математическому моделированию в экологии и метеорологии – в ней исследуются закономерности возникновения кислотных осадков, решаются две взаимосвязанные задачи: во-первых, моделирования образования и развития первичного кислотосодержащего облака и, во-вторых, моделирования его эволюции и выпадения осадков.
Методология исследования базируется на использовании авторской математической модели, ранее апробированной на задаче о поглощении газообразных загрязнителей туманом, и дополнении ее эмпирическими данными о виде распределений капель в формирующихся облаках, их водности. При этом акцент сделан на протекании этих процессов в присутствии газообразного загрязнителя – диоксида серы, источником которого являются автомагистрали или предприятия.
Научная новизна представленного исследования, по мнению рецензента, заключается в апробации капельного компонента комплексной математической модели, в полученных результатах численного эксперимента по моделированию образования первичного кислотосодержащего облака в присутствии газообразного диоксида серы, их сравнении с известными эмпирическими данными о распределении капель и водности облаков, а также выводах о характере поглощения облаком газообразного диоксида серы.
Структура статьи включает Введение, два раздела: Постановка эксперимента, Апробация, Выводы и Библиографию.
Во введении показано место рассматриваемых моделей среди известных подходов, выделены их достоинства и недостатки, сформулирована цель исследования, приведена краткая история разработок. При описании постановки эксперимента указаны его условия, обоснован подбор примененных параметров. В разделе «Апробация» приведены использованные при моделировании формулы, показаны полученные графики распределения капель и водности облака.
В выводах отмечается, что менее чем за минуту роста облака было поглощено около 3% газообразного диоксида серы.
Представленный библиографический список статьи включает 12 источников: публикации отечественных и зарубежных ученых в периодических научных журналах и монографии, изданные за период 1984-2019 гг. На приведенные в списке литературы источники имеются адресные ссылки по тексту статьи, что свидетельствует о наличии в публикации апелляции к оппонентам. Статья может представлять интерес для читателей, интересующихся вопросами применения методов математического моделирования и компьютерного экспериментирования в экологии и метеорологии.
Однако, требуется доработка оформления формул и рисунков: во-первых, приведенные в работе формулы не пронумерованы; во-вторых, наименования рисунков представляются не лаконичными, поскольку включают элементы, которые обычно размещают в легенде графиков и подписях осей.
Отмеченные недоработки снижают ценность представленных материалов, но не настолько, чтобы изменить общую положительную оценку проведенной авторами работы – рецензируемый материал подготовлен на актуальную тему, обладает элементами научной новизны и практической значимости, ориентирован на применение методов математического моделирования для изучения образования кислотосодержащего облака. Считаю, что статья может быть опубликована после незначительной доработки.