Математическое моделирование распространения радиоактивных веществ в воздушной среде в районах объектов энергетики
На правах рукописи
Зубов Владимир Николаевич
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ РАДИОАКТИВНЫХ ВЕЩЕСТВ В ВОЗДУШНОЙ СРЕДЕ В РАЙОНАХ ОБЪЕКТОВ ЭНЕРГЕТИКИ
05.13.18 – Математические моделирование,
численные методы и комплексы программ
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата физико-математических наук
Ростов-на-Дону
2009
Работа выполнена на кафедре высокопроизводительных вычислений и информационно-коммуникационных технологий факультета математики, механики и компьютерных наук ФГОУ ВПО «Южный федеральный университет»
Научный руководитель: доктор физико-математических наук,
профессор Крукиер Лев Абрамович
Официальные оппоненты: доктор физико-математических наук,
профессор Сухинов Александр Иванович
доктор физико-математических наук,
доцент Еремеев Виктор Анатольевич
Ведущая организация: Институт математического моделирования
РАН г. Москва
Защита диссертации состоится 26 марта 2009 г в 14 ч. 20 мин. на заседании диссертационного совета Д212.208.22 при ФГОУ ВПО «Южный федеральный университет» по адресу: Таганрог, пер. Некрасовский, 44,
ауд. Д-406.
С диссертацией можно ознакомиться в Зональной научной библиотеке ФГОУ ВПО «Южный федеральный университет» по адресу: Ростов-на-Дону, ул. Пушкинская, 148.
Автореферат разослан “___” февраля 2009 г.
Ученый секретарь
диссертационного совета Д212.208.22,
доктор технических наук, профессор Целых А.Н.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Математическое моделирование распространения веществ в движущейся воздушной среде является одной из наиболее значимых задач при решении проблемы определения района рассеивания и предупреждения загрязнения воздуха, воды и почвы загрязняющимися веществами от одного или нескольких источников. Объектами, загрязняющими природную среду, могут являться как объекты энергетики (электростанции), так и промышленные предприятия (заводы и фабрики). Пространственно-временные картины гипотетических выбросов одного или нескольких веществ из различных источников при разных погодных условиях могут быть использованы экспертами для определения потенциально опасных зон, зараженных соответствующими веществами, подготовки необходимых и своевременных мероприятий, адекватных складывающейся ситуации.
Так, в случае действительного выброса или утечки радионуклидов угроза заражения может потребовать незамедлительных действий по защите биосферы в целом и обеспечения безопасности здоровья людей в частности. Кроме того, модельный прогноз распространения концентраций загрязняющих веществ имеет огромное значение на этапе устранения последствий чрезвычайных ситуаций, которые могут иметь продолжительное воздействие на окружающую среду после относительно кратковременного выброса. В случае возникновения долгосрочных и продолжительных выбросов с потенциально длительным воздействием этапы планирования, разработки ответных действий и восстановления экосистемы переплетаются и оказывают значительное воздействие друг на друга.
Важное практическое требование, предъявляемое к большинству современных моделей распространения радионуклидов, связано с необходимости оценки накопленной дозы облучения на основе пространственно-временного распределения активности. На территориях, заселенных людьми, необходимо своевременно принимать решения об оказании помощи населению.
Разработке и исследованию математических моделей транспортно-диффузионного распространения радиоактивного загрязнения в атмосфере посвящены многочисленные научные труды как российских, так и зарубежных ученых: Алояна А.Е., Бакланова А.А., Борисевича М.Д., Брандта Д., Головизнина В.М., Златева З., Пененко В.В., Печингера У., Семенчина Е.А., Соренсона Д., Сухинова А.И., Тишкина В.Ф., Хирота М., и многих других.
Цель и задачи исследования. Целью работы является разработка математической модели распространения радиоактивного вещества в движущейся среде под воздействием процессов конвекции, диффузии и радиоактивного распада для района Волгодонской АЭС, создание экономичных вычислительных алгоритмов реализации этой модели и их реализация на многопроцессорных вычислительных комплексах.
В соответствии с выбранной целью необходимо решить следующий ряд задач:
- разработать математическую модель распространения радиоактивного вещества в движущейся среде под воздействием процессов конвекции, диффузии и радиоактивного распада для района со слабохолмистым рельефом.
- разработать и исследовать попеременно-треугольных кососимметричных конечно-разностных схемы решения динамических задач конвективно-диффузионного переноса в условиях преобладания конвективных процессов над диффузионными
- создать программный модуль для моделирования распространения радиоактивных примесей в воздушной среде на основе рассматриваемых алгоритмов и реализовать его на высокопроизводительных вычислительных системах.
- провести вычислительные эксперименты по моделированию распространения радиоактивного вещества
на основе созданной математической модели с учетом метеорологических данных, рельефа и свойств подстилающей поверхности района наблюдений.
Методы исследования. Основу методологии теоретического исследования работы составляют фундаментальные положения и общие принципы теории операторно-разностных схем, а так же теорий итерационных методов и матричных вычислений.
Научная новизна. Предложен класс кососимметричных попеременно-треугольных разностных схем численного решения динамической задачи конвективно-диффузионного переноса с преобладающей конвекцией.
Теоретически обоснована и численно проверена устойчивость разработанных разностных схем.
Достоверность. Представленные в диссертационной работе теоремы и утверждения имеют строгое математическое доказательство. Результаты проведенных вычислительных экспериментов хорошо согласуются с полученными теоретическими результатами.
Практическая значимость. Разработанные кососимметричные попеременно-треугольные конечно-разностные схемы обладают высокой экономичностью при решении нестационарного уравнения конвекции-диффузии, особенно в условиях преобладания конвективных процессов над диффузионными.
Реализация используемой математической модели как на квадратной, так и на треугольной сетке позволяет наиболее корректно аппроксимировать различные области с нерегулярным рельефом местности.
Реализованный на многопроцессорной вычислительной системе программный модуль позволяет выполнять расчёты распространения радиоактивных примесей в атмосфере для поддержки принятия решений в случае возникновения нештатных ситуаций на объектах энергетики.
Апробация работы
Основные результаты диссертационной работы докладывались на III и IV Всероссийской конференции “Актуальные проблемы прикладной математики и механики” (п. Абрау-Дюрсо, 2006г., 2008г.), на XI и XII Всероссийских школах-семинарах молодых ученых “Современные проблемы математического моделирования” (п. Абрау-Дюрсо, 2005г., 2007г.),
на Всероссийской молодёжной школе-конференции “Численные методы решения задач математической физики” (г. Казань, 2006г.),
на Международной конференции “Тихонов и современная математика” (г. Москва, 2006г.), на II Международной конференции “Matrix methods and operator equations” (г. Москва, 2007г.), на международном семинаре NASCom08 (г. Ростов на Дону, 2008г.)
Публикации. По теме диссертации опубликовано 14 печатных работ [1]-[14], в том числе 10 в соавторстве. Из них 1 статья в российском реферируемом журнале, 1 свидетельство о регистрации программы на ЭВМ, 1 монография, 7 статей в сборниках трудов, 4 в тезисах докладов всероссийских и международных конференций.
Структура и объем диссертации.
Диссертационная работа состоит из введения, трех глав, заключения и списка литературы. Общий объем диссертации составляет 120 страниц, в том числе 37 рисунков и 6 таблиц. Список литературы состоит из 127 наименований.
Автор работы выражает глубокую признательность своему научному руководителю д.ф.-м.н., проф. Крукиеру Л.А., научным сотрудникам
д.ф.-м.н. Муратовой Г.В и к.ф.-м.н. Чикиной Л.Г., благодарит научный коллектив ЛВЭ ЮГИНФО ЮФУ за помощь и ценные советы.
СОДЕРЖАНИЕ РАБОТЫ
Во введении раскрывается актуальность темы работы, дана постановка основных задач и изложены цели диссертационного исследования, представлено краткое содержание работы и сформулированы основные результаты, представленные к защите.
Первая глава посвящена построению модели распространения и осаждения радионуклидов в атмосфере, включающей диагностическую модель построения поля ветра на основе наблюдения. Глава состоит из трех разделов.
В первом разделе главы рассмотрено общее описание моделей распространения загрязнения в атмосфере. Выделены и охарактеризованы основные компоненты. Для каждого из описываемых физических процессов, участвующих в распространении радионуклидов, приведен необходимый набор метеорологических данных[1].
Обозначены режимы использования модели: в подготовительном планировании природоохранных мероприятий, при возникновении чрезвычайной ситуации (ЧС) в районе объекта энергетики, а так же на этапе устранения последствий ЧС. Изложены требования, предъявляемые к модели на различных этапах эксплуатации[2].
Второй раздел посвящен классификации моделей, проводимой на основании различных факторов. Рассмотрены два типа моделей – лагранжева и эйлерова. Приведены основные способы описания физических процессов моделями. Представлена пространственная классификация моделей (Рис. 1). Для каждого из классов моделей перечислены наиболее значимые физические явления и процессы, оказывающие существенное влияние на распространение загрязнения[3].
Приведены некоторые атмосферные траспортно-диффузионные модели разных масштабов, используемые в построении оперативных прогнозов различными метеорологическими центрами и исследовательскими институтами.
Рис. 1. Пространственные масштабы и физические процессы, имеющие к ним набольшее отношение
В третьем разделе главы представлено описание и структура отдельных компонент разрабатываемой мезо-
масштабной эйлеровой атмосферной траспортно-диффузионной модели распространения газообразной радиоактивной примеси. Перечислены используемые метеорологические данные:
- измерение скорости и направления ветра с помощью анемометра;
- сведения об осадках и классе стабильности атмосферы;
- карта местности со слабохолмистым рельефом.
Описана модель построения стационарного солиноидального поля скоростей ветра на основе данных наблюдений. Модель представляет собой модификацию модели URBAN[4], разработанной в ИММ РАН. В настоящей работе модель была модифицирована, в неё внесены изменения в соответствии с классом решаемых задач, связанные с включением в расчёт вертикальной компоненты скорости, необходимой для описания поля ветра над слабохолмистым рельефом.
Процесс построения стационарного поля ветра является итерационным, начинается с задания начального распределения поля ветра, аналогично модели URBAN, после чего циклически повторяются три процедуры, в которые внесены необходимые изменения:
1. Построение поля дивергенций по формуле(1):
, (1)
где 
– 
-ое приближение трехмерного поля ветра, 
.
Изменения в данной процедуре связаны с добавлением компоненты 
в формулу (1) построения поля дивергенций.
2. Минимизация дивергенции по формуле:
, (2)
где 
.
Модификация процедуры URBAN состоит в пересчёте величины 
по формуле (2).
3. Итерационное сглаживание поля ветра (порядка 2-3 итераций):
, (3)
где 
.
Изменения в третьей компоненте связаны с добавлением слагаемого с множителем 
, отвечающего вертикальному сглаживанию.
Итерационный процесс останавливается, как только абсолютная величина дивергенции, найденная по формуле (1), становится меньше 
или число итераций превышает установленный верхний порог.
Основу эйлеровой модели переноса радиоактивного загрязнения, разработанной в данной работе, составляет дифференциальная форма закона сохранения количества вещества в движущейся среде с постоянной плотностью.
Будем считать, что воздушная среда несжимаемая, то есть
, (4)
Процессы конвективного распространения, диффузии, радиоактивного распада и влажного осаждения в области с границей описываются уравнением:
, (5)
где 
– объемная активность радионуклидов (
), 
– коэффициент радиоактивного распада (
, где 
– период полураспада (с)), 
– диффузионные коэффициенты, включающие в себя турбулентную диффузию, 
, 
– источник радионуклидов, 
– коэффициент влажного осаждения.
С учетом (4) уравнение конвекции-диффузии (5) запишем в симметричной форме:
, (6)
К (6) добавим начальные и граничные условия в области проведения расчётов 
:
(7)
, 
, (8)
где 
– нормаль к границе.
Выбор параметров 
и 
в проводимых вычислительных экспериментах осуществляется, исходя из сведений о естественном радиационном фоне, других условий физической задачи, включая направление ветра у границы.
Для определения параметра 
влажного осаждения используется методика, применяемая в мезо-макромасштабной модели DERMA[5], разработанной в Датском метеорологическом институте:
, (9)
где 
– коэффициент осаждения, соответствующий вымыванию загрязняющих веществ дождем (washout);
– коэффициент осаждения, соответствующий вымыванию загрязняющих веществ дождем или снегом (snow scavenging);
– коэффициент осаждения, соответствующий растворению радионуклидов в тумане и облачности с последующим осаждением (rainout).
Отметим, что сухое мелкозернистое осаждение газообразной примеси не принимается во внимание ввиду незначительности данного процесса в мезо-
масштабе (Рис. 1).
Таким образом, описанная эйлерова модель (6)-(9) используется для моделирования распространения мелкозернистых и газообразных инертных (без учета химического взаимодействия веществ) радионуклидов на мезо-макромасштабном уровне.
Во второй главе диссертационной работы представлены основные теоретические результаты, сопровождаемые вычислительными экспериментами.
Глава состоит из трех разделов.
В первом разделе приведены общие сведения, определения и теоремы теории линейных операторов и разностных схем, используемые в последующих разделах второй главы.
Во втором разделе представлено описание базового алгоритма решения нестационарных задач, названного кососимметричными попеременно-треугольными конечно-разностными схемами (КПТС). В канонической форме алгоритм имеет следующий вид:
, 
, 
, (10)
где 
– шаг по времени. Операторы 
и 
, действующие в конечномерном гильбертовом пространстве 
, представлены соответсвующими матрицами в стандартном базисе. 
– операторно-разностный аналог пространственного дифференциального оператора (6). Оператор 
имеет вид:
, (11)
где 
– оператор с диагональной матрицей, 
– константа, 
и 
нижне и верхнетреугольные части матрицы 
, где 
, 
или 
Получены следующие достаточные условия устойчивости кососимметричных разностных схем (10)-(11).
Теорема 1. Пусть оператор 
диссипативный, оператор 
имеет вид (11). Если 
и выполнено операторное неравенство 
, то схема (10) устойчива в 
, 
.
Следствие 1. Если оператор 
диссипативен, оператор 
имеет вид (11), выполнено неравенство 
и операторное неравенство 
, где 
, то схема (10) устойчива в 
, 
.
Для частного случая BS=E, т.е.
(12)
доказано следствие.
Следствие 2. Пусть оператор 
диссипативный, оператор 
имеет вид (12), выполнено условие 
где 
, тогда схема (10) устойчива в 
, 
.
В третьем разделе второй главы диссертационной работы представлены результаты численного исследования устойчивости кососимметричной попеременно-треугольной разностной схемы (10) с 
, а так же сравнения данной схемы с кососимметричной треугольной схемой[6] 
(КТС) и явной схемой 
.
Для сравнения схем решалась следующая модельная задача распространения радиоактивного загрязнения в воздухе: краевая задача (6)-(8) на равномерной прямоугольной пространственной сетке размера 
в параллелепипеде 
(м3), на временном интервале 
(с), с краевыми условиями первого рода, с тремя диффузионными коэффициентами: 
. Коэффициент 
соответствует сильному конвективному преобладанию, а 
– слабому диффузионному преобладанию.
В качестве начальных и краевых условий выбираются функции, удовлетворяющие следующему аналитическому решению задачи:
, (13)
Результаты вычислительного эксперимента для ветрового поля
. (14)
представлены в таблице 1.
В результате проведенных вычислительных экспериментов для разных полей скоростей были сделаны следующие выводы:
- Точность приближенного решения КПТС значительно выше, чем у КТС, сравнима с явной схемой
- В случае конвективного преобладания:
- КПТС решает задачу с шагом по времени в 2 раза больше, чем у ТКС, и на порядок больше, чем у явной схемы
- Точностью приближенного решения с КПТС сравнима с явной схемой, выше, чем КТС, тогда как время расчёта меньше, чем у КТС и явной схемы
Таблица 1. Численное сравнение точности приближенного решения и схемной устойчивости для поля ветра (14)
| Расчётная схема |  (м2/с) | ||||
| 10 | 100 | 1000 | |||
| Явная схема | Время расчёта (с) | 42.9 | 6.3 | 14.5 | |
| Отн. погрешность прибл.решения (%) | 1.8 | 1.5 | 0.89 | ||
| Макс. шаг, сохраняющий устойчивость схемы (с) | 0.5 | 3.4 | 1.5 | ||
| КПТМ | Время расчёта (с) | 2.4 | 3.5 | 34.41 | |
| Погрешность прибл.решения (%) | 1.6 | 0.95 | 0.85 | ||
| Макс. шаг, сохраняющий устойчивость схемы (с) | 19.9 | 13.5 | 1.4 | ||
| Теоретический шаг, сохраняющий уст-ть (из дост. условия устойчивости) (с) | 18.9 | 10.5 | 1.4 | ||
| КТМ | Время расчёта (с) | 3.5 | 4.4 | 23.3 | |
| Погрешность прибл.решения (%) | 4.2 | 2.6 | 0.97 | ||
| Макс. шаг, сохраняющий устойчивость схемы (с) | 9.9 | 7.8 | 1.5 | ||
| Теоретический шаг, сохраняющий уст-ть (из дост. условия устойчивости) (с) | 9.8 | 7.7 | 1.5 | ||
Третья глава посвящена практическому применению разработанных методов и содержит результаты численных экспериментов на основе используемых моделей. В вычислительных экспериментах моделируется краткосрочное распространение инертного газообразного радионуклида для района Волгодонской АЭС со слабохолмистым рельефом подстилающей поверхности. Криптон является одной из основных компонент газообразного радиоактивного выброса данного энергоблока. Для регуляризации решения, получаемого с помощью используемых немонотонных центрально-разностных схем, в аппроксимацию транспортно-диффузионного уравнения вносится стандартный член искусственной вязкости, гарантирующий первый порядок точности по пространству аппроксимации краевой задачи перед “фронтом” решения.
Глава состоит из четырех разделов.
В первом разделе представлено описание входных параметров и начальных данных проведенных вычислительных экспериментов. В работе проделаны две серии расчетов распространения радиоактивного загрязняющего изотопа инертного газа 85Kr в воздухе в районе расположения Волгодонской AЭС в различных погодных условиях.
Во втором разделе описан параллельный алгоритм реализации транспортно-диффузионной модели. Выполнение расчётов в обоих последующих вычислительных экспериментах осуществлено в многопроцессорной вычислительной среде с распределенной памятью. Алгоритмической основой распределения вычислений на несколько вычислительных узлов является применение геометрического метода декомпозиции области. На рисунке 2 условно обозначен механизм взаимодействия соседствующих подобластей D(I) и D(II) с помощью граничных узлов. В используемом способе распределения вычислений каждая подобласть обрабатывается на отдельном узле вычислительной системы. При этом нумерация подобластей осуществляется по методике красно-черного разбиения, что позволяет эффективно реализовать предложенные в данной работе кососимметричные конечно-разностные схемы.
Рис 2. Механизм взаимодействия соседних подобластей и обмена данными между узлами вычислительной системы
В результате параллельной реализации математической модели на вычислительном кластере “IBM Cluster 1350”, функционирующем в Южном федеральном университете, удалось в значительной мере ускорить выполнение производимых вычислений и добиться многократного увеличения производительности реализации задачи за счет относительно большого размера собственной памяти процессора узла ВС и невысоких требований к выделенной памяти самой задачи (Таблица 2). В таблице представлено время расчётов каждой из двух основных компонент модели распространения загрязнения на одном и нескольких узлах имеющегося в наличии вычислительного кластера IBMX.
Данные результатов расчётов, приведенные в таблице 2, характеризуют задачу как недостаточно полно задействующую все вычислительные возможности высокопроизводительного кластера, что отражается более чем в двукратном увеличении производительности системы при увеличении числа узлов вдвое.
Таблица 2. Время расчёта двух основных компонент транспортно-диффузионной модели
| Число узлов ВС | Время расчёта (сек) | |
| Диагностическая модель | Транспортно-диффузионная компонента | |
| 1 | 45 | 2004 |
| 2 | 16,7 | 423 |
| 4 | 6,8 | 208 |
| 8 | 3,3 | 102 |
В третьем разделе представлены результаты расчётов с помощью диагностической модели метеорологической обстановки в рассматриваемой области размером 2060.5 км3 (Рис. 3). В условиях С-В ветра 4-6 м/с на высоте 140 м над уровнем водохранилища, построена стационарная трехмерная картина движения воздушных масс. В рассматриваемой системе координаты 
ветер в направлен вдоль положительного направления оси 
.
Рис 3. Область с рельефом подстилающей поверхности.
На рисунке 4 приведена зависимость скорости ветра (компонента 
) от высоты в вертикальном сечении, проходящем через трубу энергоблока АЭС.
Рис 4. Скорость ветра (компонента 
) в вертикальном сечении, проходящем через трубу энергоблока АЭС 
, (м/с).
На основе построенного стационарного поля ветра проведены два вычислительных эксперимента в различных метеорологических условиях – с осадками в виде дождя и без осадков.
В качестве начальных данных, общих для обоих экспериментов, рассматриваются следующие:
- класс атмосферной стабильности – “E”[7] (высота пограничного слоя – 300 м);
- источник выброс – вершина трубы энергоблока АЭС
с координатами ( 4000 , 3000 , 150 ) м; - активность источника – 208333 Бк (соответствует штатному режиму работы энергоблока);
- время выброса – 1 час;
- начальная активность
в области – 0 Бк/м3; - общее время вычислительных экспериментов – 2 часа.
В первом эксперименте моделируется перемещение загрязняющего облака в расчётной области на протяжении всего расчётного времени (Рис. 5 а) ) без учета влажного осаждения.
Во втором эксперименте в начальный момент времени на высоте верхней границы расчётной области в “вдувает” дождевое облако размером 25 8 
с изрезанной границей и скругленными краями. Облако движется со скоростью ветра у верхней границы расчётной области. При этом выпадают осадки в виде дождя с постоянной скоростью 5 мм/час, вызывая влажное осаждение радиоактивных веществ на подстилающую поверхность (Рис. 5 б) ).
а) б)
Рис 5. Объемная активность 85Kr (Бк/м3) в различные моменты времени: а) – эксперимент без осадков, б) – с дождевым облаком.
Картины активности, представленные на рис. 5, позволяют произвести численное сравнение присутствия в области влажных осадков в виде дождя, оценить влияние влажного осаждения на распространение радиоактивного облака как во время выброса, так и после него.
Во четвертом разделе представлены результаты краткосрочных расчётов распространения радиоактивного загрязнения в меняющихся погодных условиях.
Рис 6. Область с рельефом подстилающей поверхности.
Исходя из предположения о непрерывности изменения направлений движения воздушных масс в области под влиянием различных факторов и невозможности в среднесрочной перспективе пренебрежения динамическими процессами, протекающими в атмосфере, оказывающими непосредственное влияние на поле ветра, в данном разделе представлены результаты расчётов продолжительного транспортно-диффузионного распространения радиоактивной примеси 
в условиях периодичного обновления поля ветра.
В качестве рассматривается так же мезо- масштабная расчётная область 40400.5 км3 с изрезанной нижней гранью (Рис. 6).
Расчёты производились как на квадратной, так и на треугольной сетках со следующими входными данными задачи:
- класс атмосферной стабильности – “E”[8] (высота пограничного слоя – 300 м);
- источник выброс – вершина трубы энергоблока АЭС
с координатами ( 0 , 0 , 150 ) м; - активность источника – 208333 Бк (соответствует штатному режиму работы энергоблока);
- время выброса – 3 час;
- начальная активность
в области – 0 Бк/м3; - общее время вычислительных экспериментов – 4 часа.
- начальная скорость ветра – западный 10-12 м/с на высоте 150 м;
- скорость изменения направления ветра – (–/180) рад/мин;
- скорость обновления поля ветра – 1 мин.
Рис 7. Распределение активности (Бк) 
в области , квадратная сетка, вид сверху.
Результаты вычислительных экспериментов представлены на рисунках 7-10.
Рис 8. Распределение активности (Бк) 
в области , квадратная сетка, вид сбоку.
Рис 9. Распределение активности (Бк) 
в области , треугольная сетка, вид сверху.
Рис 10. Распределение активности (Бк) 
в области , треугольная сетка, вид сверху.
Результаты расчётов позволяют сформулировать следующие выводы:
- дождевое облако, находящееся на пути следования радиоактивного загрязнения, обладает способностью значительного уменьшения объемной активности радионуклидов в результате влажного осаждения на поверхность Земли;
- осадки оказывают наиболее значимое влияние на картину активности радионуклидов после завершения выброса, во время свободного перемещения загрязняющего облака по направлению к границе расчётной области.
- результаты вычислительных экспериментов, полученные на треугольной и квадратной сетках, показывают качественное совпадение;
- применение треугольной сетки в аппроксимации транспортно-диффузионного уравнения во многих случаях позволяет более точно аппроксимировать расчётную область с изрезанным рельефом, однако, не дает существенного увеличения точности решения задачи по сравнению с квадратной сеткой.
В заключении приведены основные результаты, полученные в диссертационной работе.
Положения, выносимые на защиту
- Разработана математическая модель распространения радиоактивного загрязнения в воздухе для районов со слабохолмистым рельефом местности.
- Предложен, теоретически и численно исследован новый класс кососимметричных попеременно-треугольных разностных схем аппроксимации динамической задачи конвективно-диффузионного переноса в несжимаемой среде.
- Программно реализована математическая модель транспортно-диффузионного распространения радиоактивного вещества в воздушной среде на многопроцессорной вычислительной системе с распределенной памятью с применением треугольной и прямоугольных сеток.
- Проведены расчёты модельного рассеивания радиоактивного облака в районе Волгодонской АЭС в различных погодных условиях.
Список работ, опубликованных по теме диссертации
I. Публикации в центральных изданиях, включенных в перечень периодических изданий ВАК РФ
- Зубов В. Н. Моделирование процессов радиоактивного загрязнения воздушной среды в районе Волгодонской АЭС / Зубов В. Н., Крукиер Л. А., Муратова Г. В., Субботина Т. Н. // Математическое моделирование. 2008. Т. 20. № 7. С.85-92
II. Монографии
- Зубов В. Н. Математическое моделирование процессов конвективно-диффузионного переноса в задачах экологии / Крукиер Л. А., Муратова Г. В., Андреева Е. М., Зубов В. Н., Субботина Т. Н., Никитенко О. Б. – Ростов-на-Дону: Издательство Южного федерального университета, 2008. – 304 c.
III Свидетельства о регистрации программ на ЭВМ
- Зубов В. Н. Расчёт движения атмосферы в районе Волгодонской АЭС / Зубов В. Н., Субботина Т. Н. РОСПАТЕНТ. Свидетельство №2008611407 от 20.03.2008.
IV. Публикации в других изданиях
- Зубов В. Н. Численная реализация математической модели процессов радиоактивного загрязнения воздушной среды / Субботина Т. Н., Зубов В. Н. // Сборник трудов XI Всероссийской школы-семинара “Современные проблемы математического моделирования”, Ростов-на-Дону: РГУ, 2005.
- C. 377-386 - Зубов В. Н. Модульный комплекс моделирования распространения радиоактивных веществ в несжимаемой среде / Зубов В. Н. // Сборник трудов научно-методической конференции “Современные информационные технологии в образовании: Южный Федеральный округ”, Ростов-на-Дону: ЦВВР, 2006. - С. 113-114
- Zubov V. N. Modelling of Radioactive Pollution Transport from a Point Source Air / Zubov V. N., Subbotina T. N. // In Proceeding of the International Conference “Tikhonov and Contemporary Mathematics”, Moscow, 2006, Vol.2. - P.207-209.
- Зубов В. Н. Решение трехмерного нестационарного уравнения конвекции-диффузии двухслойными треугольными кососимметричными разностными схемами / Зубов В. Н., Субботина Т. Н. // Тезисы докладов III Всероссийской конференции «Актуальные проблемы прикладной математики и механики», Екатеринбург: УрО РАН, 2006 г., - с. 41-42.
- Зубов В. Н. Моделирование переноса радиоактивного вещества в атмосфере / Зубов В. Н., Муратова Г. В., Субботина Т. Н. // Тезисы докладов III Всероссийской конференции «Актуальные проблемы прикладной математики и механики», Екатеринбург: УрО РАН, 2006 г., - с. 40-41.
- Зубов В. Н. Применение двухслойных кососимметричных разностных схем для решения нестационарной задачи конвекции-диффузии / Зубов В. Н. // Труды Математического центра им. Н.И. Лобачевского. Т.31. Численные методы решения задач математической физики. Казань: Изд-во Казанского мат. общества, 2006
- Zubov V. N. The radioactive pollution transport modeling in air nearby Volgodonsk nuclear station / Zubov V. N., Krukier L. A., Muratova G. V., Subbotina T. N. // Proceedings of II international conference on matrix methods and operator equations, Moscow, Institute of Numerical Mathematics Russian Academy of Science, 2007. – p. 94-95.
- Зубов В. Н. Решение стационарной задачи конвекции-диффузии с преобладающей конвекцией на треугольной и квадратной сетках / Крукиер Л. А., Зубов В. Н., Субботина Т. Н. // Сборник трудов XII Всероссийской школы-семинара «Современные проблемы математического моделирования», Ростов-на-Дону: РГУ, 2007. – C. 166-175
- Зубов В. Н. Математическое моделирование и вычислительный эксперимент для задач экологии воздушной среды / Зубов В. Н. // Сборник трудов научно-методической конференции «Современные информационные технологии в образовании: Южный Федеральный округ», Ростов-на-Дону: ЦВВР, 2007. – С. 110-111
- Зубов В. Н. Моделирование радиоактивного загрязнения атмосферы с помощью кососимметричных попеременно-треугольных разностных схем / Зубов В. Н. // Тезисы докладов III Всероссийской конференции «Актуальные проблемы прикладной математики и механики», Екатеринбург: УрО РАН, 2008. –С. 24-25.
- Зубов В. Н. Моделирование процессов распространения примеси в воздушной среде в районе объектов энергетики / Муратова Г. В., Крукиер Л. А., Дацюк В. Н., Дацюк О. В., Чикин А. Л., Зубов В. Н. // Сборник трудов XV Всероссийской научно-методической конференции “Телематика-2008”. Т.1. С.-П:, 2008. – с.80-82
Лично автором в работах [2,4,5,6] разработана модель переноса загрязнения в атмосфере в условиях преобладания конвективных процессов. В работах [3,7,8] автор разработал и применил кососимметричную попеременно-треугольную разностную схему в решении краевой задачи конвекции-диффузии. В работах [9,10] автором произведена адаптация модели транспортно-диффузионного переноса к району Волгодонской АЭС со слабохолмистым рельефом в холодных погодных условиях. В [11] автор построил аппроксимацию краевой задачи конвекции-диффузии со смешанными краевыми условиями кососимметричной попеременно-треугольной схемой на треугольной сетке. В [1] автором теоретически обоснована устойчивость предложенных кососимметричных схем, сформулированы и доказаны достаточные условия устойчивости. В [3,12,13] произведена модификация вычислительного алгоритма модели URBAN, адаптация его к району со слабохолмистым рельефом. В [14] предложен алгоритм реализации разработанной ранее модели переноса на многопроцессорных вычислительных системах с распределенной памятью на основе методики красно-черного разбиение расчётной области.
[1] Cochairperson W., Cochairperson N., Hicks B., Payton D. Federal Research and Development Needs and Priorities for Atmospheric Transport and Duffusion Modeling. Washington DC, 2004 – 207 p.
[2] Serafin R.J., Barron E.J., Clifford S.F., Duncan L.M., LeMone M.A. Tracking and Prediction the Atmospheric Dispersion of Hazardous Releases. Washington D.C.: National Academies Press, 2003. – 101 p.
[3] Borysiewicz M.J., Borysiewicz M.A., Garanty I., Kozubal A., Potempski S., Rowinski P.M., Zheleznyak. M. Models and techniquies for health and environmental hazard assessment and management, Warsaw, 2006. – 1360p.
[4] Тишкин В.Ф., Клочкова Л.В., Сузан Д.В.. Применение нового метода зануления дивергенции ветрового поля для моделирования процессов распространения летучих газов при огибании крутых препятствий. // Сборник трудов Всероссийской конференции "Математическое моделирование и проблемы экологической безопасности", 2000. – с.208215.
[5] Sorenson J.H., Baklanov A., Hoe S. The Danish emergency response model of the atmosphere (DERMA) // Journal of Environmental Radioactivity, 96, 2007. – pp. 122-129.
[6] Крукиер Л.А., Муратова Г.В., Субботина Т.Н. Эффективные разностные схемы решения нестационарного уравнения конвекции-диффузии // Мат. моделирование. – 2005 – Т.17, №12 – С. 80-86.
[7] Pasquill F. The Estimation of the dispersion of windborne material // Meteorological Magazine Vol. 90, 1961. – pp. 33-49
[8] Pasquill F. The Estimation of the dispersion of windborne material // Meteorological Magazine Vol. 90, 1961. – pp. 33-49
