Весовые алгоритмы метода мажорантной частоты для статистического моделирования решения пространственно-однородных нелинейных кинетических уравнений больцмановского типа тема автореферата и диссертации по математике, 01.01.07 ВАК РФ

На правах рукописи

Блощицьша Ольга Витальевна

ВЕСОВЫЕ АЛГОРИТМЫ МЕТОДА МАЖОРАНТНОЙ ЧАСТОТЫ ДЛЯ СТАТИСТИЧЕСКОГО МОДЕЛИРОВАНИЯ РЕШЕНИЯ ПРОСТРАНСТВЕННО-ОДНОРОДНЫХ НЕЛИНЕЙНЫХ КИНЕТИЧЕСКИХ УРАВНЕНИЙ БОЛЬЦМАНОВСКОГО ТИПА

01.01.07 - вычислительная математика

Автороферат

диссертации на соискание ученой степени кандидата физико-математических наук

3 О МАЙ 2013

Новосибирск - 2013

005060261

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте вычислительной математики и математической геофизики Сибирского отделения Российской академии наук (ИВМиМГ СО РАН).

Научный руководитель:

доктор физико-математических наук, Рогазинский Сергей Валентинович

Официальные оппоненты:

Иванов Михаил Самуилович, доктор физико-математических наук,

Федеральное государственное бюджетное учреждение науки Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук, заведующий лабораторией

Плотников Михаил Юрьевич, кандидат физико-математических наук,

Федеральное государственное бюджетное учреждение науки Институт теплофизики им. С.С. Кутателадзе Сибирского отделения Российской академии наук, старший научный сотрудник

Ведущая организация:

Омский филиал Федерального государственного бюджетного учреждения науки Института математики им. С.Л. Соболева Сибирского отделения Российской академии наук

Защита состоится 19 июия 2013 года в 15 часов на заседании диссертационного совета Д 003.061.01 на базе Федерального государственного бюджетного учреждения науки Института вычислительной математики и математической геофизики Сибирского отделения Российской академии наук по адресу: 630090, Новосибирск, пр. академика Лаврентьева, 6.

С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Института вычислительной математики и математической геофизики Сибирского отделения Российской академии наук.

Автореферат разослан мая 2013 года.

Ученый секретарь

диссертационного совета Д 003.061.01 д.ф.-м.н.

^ С.В. Рогазинский

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

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

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

Нелинейное кинетическое уравнения Больцмана и нелинейное кинетическое уравнение Смолуховского называют нелинейными кинетическими уравнениями больцмановского типа из-за того, что они оба имеют одинаковый тип нелинейности. Их сложная нелинейная структура делает, в большинстве случаев, невозможным их аналическое решение. Поэтому единственным способом нахождения решения задач для уравнений этого типа являются численные методы. В силу того, что физическая интерпретация нелинейных кинетических уравнений больцмановского типа носит вероятностный характер, эффективным методом численного решения задач для уравнений такого типа, является численное статистическое моделирование или метод Монте-Карло.

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

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

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

Основные цели работы.

• Разработка весовой модификации метода мажорантной частоты для численного решения уравнения Больцмапа и уравнения Смолуховского.

• Обоснование применимости метода путем доказательства конечности дисперии стандартных весовых оценок.

• Исследование эффективности метода с помощью проведения численных экспериментов.

Методы исследования базируются на теории кинетических уравнений, теории интегральных уравнений второго рода и теории весовых методов статистического моделирования.

Научная новизна и практическая значимость работы.

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

Личный вклад соискателя заключается в обосновании применимости метода путем доказательства теорем о конечности дисперсии весовых оценок метода Монте-Карло, создании программных алгоритмов и подтвержде-

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

Апробация работы. Результаты, изложенные в диссертации, докладывались и обсуждались на семинаре Отдела статистического моделирования в физике ИВМиМГ СО РАН, а также на следующих конференциях:

— Международная научная студенческая конференция "Студент и научно-технический прогресс" (2007, 2008 гг.)

- 6th St. Peterburg Workshop on Simulaton, (г. Санкт-Петербург, 2009 г.)

- Конференция молодых ученых ИВМиМГ СО РАН (2009, 2010 гг.)

— Всероссийская конференция по вычислительной математике КВМ-2011 (г. Новосибирск, 2011 г.)

Публикации. По тематике диссертации автором опубликовано 7 работ, среди которых 2 работы в изданиях из списка ВАК. Список опубликованных работ помещен в конце автореферата.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка использованной литературы. Объем работы -87 страниц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Во Введении обоснована актуальность, цель и задачи исследований, изложено краткое содержание диссертации по главам и параграфам.

Первая глава носит справочный характер. В пей приведены основные сведения из теории метода Монте-Карло для решения интегральных уравнений, которые будут использоваться в дальнейшем.

В параграфе 1.1 рассматривается интегральное уравнение второго рода и дается краткое обоснование представления функционалов от его решения рядом Неймана.

В параграфе 1.2 вводится понятие цепи Маркова и рассматриваются весовые оценки метод а Монте-Карло для оценки функционалов от решения интегрального уравнения второго рода и приводится обоснование их несмещенности.

В параграфе 1.3 приводятся условия конечности дисперсии весовых оценок.

Во Второй главе представлен весовой метод Монте-Карло для решения нелинейного уравнения Больцмана и приведены результаты численных расчетов. Глава состоит из семи разделов.

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

/ОМ) | =/оМ, ¿е(о,Т].

/=о

Под численным решением задачи Коши для уравнения Больцмана понимается нахождение линейных функционалов от функции /(у, вида:

щ = I /1(у)/(у,<)сгу.

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

В разделе 2.3 описывается метод прямого статистического моделирования.

В разделе 2.4 формулируется метод мажорантной частоты и приводится алгоритм прямого статистического моделирования для приближенного решения уравнения Больцмана с использованием принципа мажорантной частоты.

Раздел 2.5 посвящен описанию и обоснованию весовой модификации метода мажорантной частоты для численного решения уравнения Больцмана. В подразделе 2.5.1 введена модификация фазового пространства, позволяющая вывести интегральное уравнение, на основе которого можно строить весовые модификации прямого статистического уравнения. В число координат фазового пространства вводится номер взаимодествующей пары. В подразделе 2.5.2 в расширенном фазовом пространстве строится базовое интегральное уравнение. В подразделе 2.5.3 в число координат фазового пространства помимо номера взаимодействующей пары вводитстея также тип столкновения, что позволяет сформулировать новое интегральное уравнение, на основе которое можно стандартным способом строить весовые модификации прямого статистического моделирования. В разделе вводится цепь Маркова и случайные веса, на основе которых строятся стандартные весовые оценки метода

Монте-Карло для вычисления функционала от решения задачи Коши для уравнения Больцмана.

В разделе 2.6 формулируется теорема о конечности дисперсии стандартных весовых оценок метода Монте-Карло и приводится ее док;1зательство.

Теорема 1. Пусть для некоторого а > 0 выполняются условия

<Си/ f0(V)(l + V2fdV = Со < оо,

V _

тогда оценка по столкновениям £ = QnHn и оценка по поглощениям 1] =

п=О

QUHU, где v - номер последнего состояния цепи, имеют конечную дисперсию, а веса Qn вычисляются но формулам, приведенным в разделе 2.5.

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

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

Третья глава посвящена весовому методу Монте-Карло для решения нелинейного уравнения коагуляции. Глава состоит из семи разделов.

В разделе 3.1 формулируется задача Коши для нелинейного уравнения Смолуховского в случае чистого парного слипания частиц.

d 1

■Щп(1'Ь) = § Y, k(i,j)n(i,t)-n(l,t)^2k(l,i)n(i,t),

i+j=l 1=1

п{1, t — 0) = щ{1).

Здесь n(l, t) - числовая плотность частиц размера I в момент времени t, причем размер частицы I принимает натуральные значения; k(i,j) - коэффициенты коагуляции, которые являются заданными величинами. Под численным решением задачи Коши для уравнения Смолуховского понимается нахождение линейных функционалов от функции п{1. t).

В разделе 3.2 рассматривается математическая модель физического процесса однородной коагуляции.

Раздел 3.3 посвящен методу прямого статистического моделирования.

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

Раздел 3.5 посвящен описанию и обоснованию весовой модификации метода мажорантной частоты для численного решения уравнения коагуляции. В подразделе 3.5.1 введена модификация фазового пространства, позволяющая вывести интегральное уравнение, на основе которого можно строить весовые модификации прямого статистического уравнения. В число координат фазового пространства вводится номер взаимодествующей пары. В подразделе 3.5.2 в расширенном фазовом пространстве строится базовое интегральное уравнение. В подразделе 3.5.3 в число координат фазового пространства помимо номера взаимодействующей пары вводитстся также тип столкновения, что позволяет сформулировать новое интегральное уравнение, на основе которое можно стандартным способом строить весовые модификации прямого статистического моделирования. В разделе вводится цепь Маркова и случайные веса, на основе которых строятся стандартные весовые оценки метода Монте-Карло для вычисления функционала от решения задачи Коши для уравнения коагуляции.

В разделе 3.6 приводится доказательство конечности дисперсии стандартных весовых оценок метода Монте-Карло.

Теорема 2. Пусть для некоторого а > 0 выполняются условия

^С и I п0{Ь){ 1 + Ь)айЬ = С0 < оо, тогда оценка по столкновениям £ = £ ЯпНп « оценка по поглощениям

п—О

Т] = где V - номер последнего состояния цепи, имеют конечную дис-

N

Персию. Здесь Ь = /¿, а веса <5„ вычисляются по формулам, приведенным ¿=1

в разделе 3.5.

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

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

в отличие от оценки, полученной методом прямого моделирования, для которой наблюдается занижение, начиная с некоторого момента времени, что объясняется "исчезновением" мономеров.

В Четвертой главе описаны ценностные алгоритмы моделирования для приближенного решения уравнения Больцмана и уравнения коагуляции и приведены результаты численных экспериментов. Глава состоит из трех разделов.

В разделе 4.1 определяется функция ценности и метод выборки по важности. Приведены основные теоремы о несмещенности оценок, построенных методом выборки по важности. Также вводятся вспомогательные функции ценности.

В разделе 4.2 формулируется ценностная весовая модификация метода мажоратной частоты для численного решения уравнения коагуляции. Вводится приближение к временной функции ценности и ценностное моделирования номера пары. Приводится алгоритм ценностного моделирования цепи Маркова с использованием весовой модификации метода мажорантной частоты. В подразделе 4.2.1 приведены результаты численных расчетов для уравнения коагуляции.

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

Раздел 4.3 посвящен ценностной весовой модификации метода мажоратной частоты для численного решения уравнения Больцмана. Введено приближение к временной функции ценности. Приводится алгоритм частично-ценностного моделирования цепи Маркова с использованием метода мажорантной частоты. В подразделе 4.3.1 приведены результаты численных расчетов для уравнения Больцмана.

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

Заключение содержит перечень основных результатов диссертационной работы.

ЗАКЛЮЧЕНИЕ

• Построен весовой алгоритм метода Монте-Карло для приближенного решения нелинейных кинетических уравнений, который с одной сто-

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

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

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

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

СПИСОК РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в изданиях списка ВАК

1. S.V. Rogazinskii, O.V. Bloshchitsyna. Weighted modification of the majorant frequency method for solution of nonlinear kinetic equations // Russ.J.Numer.Anal.Math.Modelling, 2011, Vol. 26, No. 5, pp. 501-514.

2. O.B. Блощицыиа, C.B. Рогазинский. Весовая модификация метода мажорантной частоты для решения нелинейных кинетических уравнений // Журнал вычисл. матем. и матем. геофизики. 2012, Т.52, № 9, С. 1724-1734.

Прочие публикации

3. Bloshchitsyna O.V., Rogasinsky S.V. Weighted modification of the majorant frequency method for solving the Boltzrriarm equation. 6th St. Peteiburg Workshop on SiinuJaton: Proceedings, St. Petersburg, 2009, Vol. 1, C. 14.3-148.

4. Блощицыпа O.B. Весовая модификация метода мажорантной частоты для решения нелинейного уравнения Больцмана. Труды Конференции молодых ученых, Новосибирск: Издательство ИВМ и МГ СО РАН, 2009, С.3-9.

5. Блощицыпа О.В. Весовой вариант метода мажорантной частоты для приближенного решения нелинейного уравнения коагуляции. Труды конференции молодых ученых, Новосибирск: Издательство ИВМ и МГ СО РАН, 2010, С.3-10.

6. Блощицыпа О.В. Весовая модификация метода мажорантной частоты для решения уравнения Больцмана. Материалы XLV Международной научной студенческой конференции "Студент и научно-технический прогресс": Математика. Новосибирский государственный университет, Новосибирск, 2007, С.201.

7. Блощицыпа О.В. Весовая модификация алгоритма моделирования решения уравнения Больцмана с использованием схемы мажорантной частоты. Материалы XLVI Международной научной студенческой конференции "Студент и научно-технический прогресс": Математика. Новосибирский государственный университет, Новосибирск, 2008, С.311.

Подписано в печать 13.05.2013 Формат 60x84 1\16 Усл. печ. л. 1 Объем 16 стр. Тираж 100 экз. Заказ №84 Отпечатано Омега Принт 630090, г. Новосибирск, пр. Ак.Лаврентьева,6 email: omegap@yandex.ru

Федеральное государственное бюджетное учреждение науки Институт вычислительной математики и математической геофизики

Сибирского Отделения РАН

04201360318

На правах рукописи

Блощицына Ольга Витальевна

Весовые алгоритмы метода мажорантной частоты для статистического моделирования решения пространственно-однородных нелинейных кинетических уравнений больцмановского типа

01.01.07 - Вычислительная математика

Диссертация на соискание ученой степени кандидата физико-математических наук

/

Научный руководитель: д.ф.-м.н. C.B. Рогазинский

Новосибирск - 2013

ОГЛАВЛЕНИЕ

Введение......................................................................................................................................................................................................4

Глава 1. Решение интегральных уравнений методом Монте-Карло

1.1. Интегральные уравнения ..............................................................................................................................11

1.2. Цепи Маркова. Весовые оценки ..........................................................................................................13

1.3. Дисперсии оценок......................................................................................................................................................15

Глава 2. Весовые алгоритмы метода Монте-Карло для решения

нелинейного кинетического уравнения Больцмана ... 16

2.1. Уравнение Больцмана........................................................................................................................................16

2.2. Математическая модель многочастичной системы для решения уравнения Больцмана........................................................................................................................................18

2.3. Метод прямого статистического моделирования ....................................................21

2.4. Метод мажорантной частоты ................................................................................................................26

2.5. Весовая модификация прямого статистического моделирования 30

2.5.1. Модификация фазового пространства......................................................................30

2.5.2. Построение базового интегрального уравнения

в модифицированном фазовом пространстве..............................................31

2.5.3. Весовая модификация метода мажорантной частоты ....................35

2.6. Доказательство конечности дисперсии весовых оценок ............................38

2.7. Результаты численных расчетов ......................................................................................................43

Глава 3. Весовой алгоритм метода Монте-Карло для решения

нелинейного уравнения коагуляции..........................................................46

3.1. Уравнение Смолуховского ..........................................................................................................................46

3.2. Математическая модель системы с переменным количеством частиц для решения уравнения Смолуховского ..............................................................................47

3.3. Метод прямого статистического моделирования ..................................................52

3.4. Метод мажорантной частоты ................................................................................................................55

3.5. Весовая модификация прямого статистического моделирования 58

3.5.1 Модификация фазового пространства......................................................................58

3.5.2. Построение базового интегрального уравнения

в модифицированном фазовом пространстве..............................................59

3.5.3. Весовая модификация алгоритма мажорантной частоты .... 62

3.6. Доказательство конечности дисперсии весовых оценок............................64

3.7. Результаты численных расчетов......................................................................................................69

Глава 4. Алгоритмы моделирования по ценности ......................................72

4.1. Ценностное моделирование для кинетических уравнений ....................72

4.2. Ценностная весовая модификация метода мажорантной частоты для численного решения уравнения коагуляции ..................................................74

4.2.1 Результаты численных расчетов ..........................................................................................77

4.3. Ценностная весовая модификация метода мажорантной частоты для численного решения уравнения Больцмана..................................................80

4.3.1. Результаты численных расчетов........................................................................................82

Заключение........................................................................................................................................................................................84

Список литературы............................................................................................................................................................85

ВВЕДЕНИЕ

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

В кинетической теории газ описывается с помощью функции распределения, которая содержит информацию как о распределении самих молекул внутри рассматриваемой системы, так и о распределении молекулярных скоростей [2]. Функция распределения в общем случае изменяется со временем, и ее изменение для случая разреженного газа описывается интегро-дифференциальным уравнением Больцмана [3], [4]. Уравнение Больцмана лежит в основе кинетической теории газов и находит широкое применение при изучении таких явлений как перенос электронов в твердых телах и плазме, перенос нейтронов в ядерных реакторах, перенос фотонов в сверхтекучих жидкостях, перенос излучения [5]. В случае задач, связанных с рассмотрением процессов коагуляции (слипания) в газообразных или жидких средах, возникает необходимость решения нелинейного кинетического уравнения Смолуховского [6].

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

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

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

Весовой вариант статистического моделирования для численного решения нелинейных кинетических уравнений Больцмана и Смолуховского с введением глобальных весов ранее был исследован в работах [8], [9]. Для уравнения Больцмана рассматривалась модель газа - максвелловские молекулы (случай, когда произведение относительной скорости на сечение -константа); для уравнения Смолуховского был рассмотрен случай, когда коэффициенты коагуляции постоянны.

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

Целью диссертационной работы является построение и обоснование весовых модификаций метода мажорантной частоты для приближенного решения нелинейных кинетических уравнений Больцмана и Смолуховского.

Основным результатом работы является построение весового алгоритма

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

Диссертация состоит из введения, четырех глав, заключения и списка использованной литературы.

Во введении обоснована актуальность, цель и задачи исследований, изложено краткое содержание диссертации по главам и параграфам.

Первая глава носит справочный характер. В ней приведены основные сведения из теории метода Монте-Карло для решения интегральных уравнений и которые будут использоваться в дальнейшем.

В параграфе 1.1 рассматривается интегральное уравнение второго рода и дается краткое обоснование представления функционалов от его решения рядом Неймана.

В параграфе 1.2 вводится понятие цепи Маркова и рассматриваются весовые оценки метода Монте-Карло для оценки функционалов от решения интегрального уравнения второго рода и приводится обоснование их несмещенности.

В параграфе 1.3 приводятся условия конечности дисперсии весовых оценок.

Во второй главе представлен весовой метод Монте-Карло для решения нелинейного уравнения Больцмана и приведены результаты численных расчетов. Глава состоит из семи параграфов.

В параграфе 2.1 определяется одночастичная функция распределения

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

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

В параграфе 2.3 описывается метод прямого статистического моделирования.

В параграфе 2.4 формулируется метод мажорантной частоты и приводится алгоритм прямого статистического моделирования для приближенного решения уравнения Больцмана с использованием принципа мажорантной частоты.

Параграф 2.5 посвящен описанию и обоснованию весовой модификации метода мажорантной частоты для численного решения уравнения Больцмана. В разделе 2.5.1 введена модификация фазового пространства, позволяющая вывести интегральное уравнение, на основе которого можно строить весовые модификации прямого статистического уравнения. В число координат фазового пространства вводится номер взаимодействующей пары. В разделе 2.5.2 в расширенном фазовом пространстве строится базовое интегральное уравнение. В разделе 2.5.3 в число координат фазового пространства, помимо номера взаимодействующей пары, вводится также тип столкновения, что позволяет сформулировать новое интегральное уравнение, на основе которого можно стандартным способом строить весовые модификации прямого статистического моделирования. В разделе вводится цепь Маркова и случайные веса, на основе которых строятся стандартные

весовые оценки метода Монте-Карло для вычисления функционала от решения задачи Коши для уравнения Больцмана.

В параграфе 2.6 приводится доказательство конечности дисперсии стандартных весовых оценок метода Монте-Карло.

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

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

Третья глава посвящена весовому методу Монте-Карло для решения нелинейного уравнения коагуляции. Глава состоит из семи параграфов.

В параграфе 3.1 формулируется задача Коши для нелинейного уравнения Смолуховского в случае чистого парного слипания частиц.

В параграфе 3.2 рассматривается математическая модель физического процесса однородной коагуляции.

Параграф 3.3 посвящен методу прямого статистического моделирования.

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

Параграф 3.5 посвящен описанию и обоснованию весовой модификации метода мажорантной частоты для численного решения уравнения коагуля-

ции. В разделе 3.5.1 введена модификация фазового пространства, позволяющая вывести интегральное уравнение, на основе которого можно строить весовые модификации прямого статистического уравнения. В число координат фазового пространства вводится номер взаимодействующей пары. В разделе 3.5.2 в расширенном фазовом пространстве строится базовое интегральное уравнение. В разделе 3.5.3 в число координат фазового пространства, помимо номера взаимодействующей пары, вводится также тип столкновения, что позволяет сформулировать новое интегральное уравнение, на основе которого можно стандартным способом строить весовые модификации прямого статистического моделирования. В разделе вводится цепь Маркова и случайные веса, на основе которых строятся стандартные весовые оценки метода Монте-Карло для вычисления функционала от решения задачи Коши для уравнения коагуляции.

В параграфе 3.6 приводится доказательство конечности дисперсии стандартных весовых оценок метода Монте-Карло.

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

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

В четвертой главе описаны ценностные алгорр1тмы моделирования для приближенного решения уравнения Больцмана и уравнения коагуляции и приведены результаты численных экспериментов. Глава состоит из трех

параграфов.

В параграфе 4.1 определяется функция ценности и метод выборки по важности. Приведены основные теоремы о несмещенности оценок, построенных методом выборки по важности. Также вводятся вспомогательные функции ценности.

В параграфе 4.2 формулируется ценностная весовая модификация метода мажоратной частоты для численного решения уравнения коагуляции. Вводится приближение к временной функции ценности и ценностное моделирования номера пары. Приводится алгоритм ценностного моделирования цепи Маркова с использованием весовой модификации метода мажорантной частоты. В разделе 4.2.1 приведены результаты численных расчетов для уравнения коагуляции.

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

Параграф 4.3 посвящен ценностной весовой модификации метода мажоратной частоты для численного решения уравнения Больцмана. Введено приближение к временной функции ценности. Приводится алгоритм частично-ценностного моделирования цепи Маркова с использованием метода мажорантной частоты. В разделе 4.3.1 приведены результаты численных расчетов для уравнения Больцмана.

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

Заключение содержит перечень основных результатов диссертационной работы.

Глава 1

Решение интегральных уравнений методом Монте-Карло

1.1. Интегральные уравнения

Рассмотрим интегральное уравнение второго рода

<р(х) = J к{х', x)ip{x')dx' + f{x), ip = Kp + f, (1.1)

где X — I-мерное евклидово пространство; f,ip € L\ L — банахов