Содержание
Введение 4
Глава 1. Использование преобразования Фурье для решения задачи Дирихле для уравнения теплопроводности 8
1.1. Основная задача и использование преобразования Фурье 8
1.2. Решение задачи Дирихле для уравнения Гельмгольца с комплексным параметром............................. 11
1.3. Решение задачи Дирихле для уравнения теплопроводности ....................................................... 15
1.4. Решение многомерной задачи Дирихле для уравнения теплопроводности и уравнения Гельмгольца............. 21
1.5. Тестовые расчеты...................................... 26
Глава 2. Решение смешанных краевых задач для уравнения Гельмгольца с комплексным параметром методом Монте-Карло и применение преобразования Фурье при решении краевых задач третьего рода для уравнения теплопроводности 35
2.1. Решение смешанной краевой задачи для уравнения Гельмгольца с комплексным параметром........................... 35
2.2. Решение смешанной краевой задачи для уравнения теплопроводности ............................................... 43
2.3. Тестовые расчеты...................................... 46
2
Глава 3. Построение глобальных оценок решений
краевых задач 52
3.1. Построение глобальных оценок решения задачи Дирихле для уравнения Гельмгольца с комплексным параметром 52
3.2. Построение глобальной оценки решения смешанной краевой задачи для уравнения Гельмгольца с комплексным параметром............................................... 57
3.3. Построение глобальных оценок решения смешанной краевой задачи для уравнения теплопроводности с использованием дискретного преобразования Фурье.................. 60
3.4. Построение глобальной оценки решения смешанной краевой задачи для уравнения теплопроводности с использованием непрерывного преобразования Фурье................... 65
3.5. Тестовые расчеты.................................... 70
Глава 4. Вычисление производных от решений краевых
задач для уравнения теплопроводности 73
4.1. Вычисление производных по параметру ................ 73
4/2. Вычисление производной по пространственной переменной с использованием нецентральной функции Грина для оператора Лапласа........................................ 81
4.3. Вычисление производной по пространственной переменной с использованием нецентральной функции Грина для оператора Гельмгольца.................................... 85
4.4. Вычисление производной по времени от решения задачи для уравнения теплопроводности .......................... 91
4.5. Тестовые расчеты.................................... 92
Заключение 96
Литература 97
3
Введение
В последнее время методы Монте-Карло получили широкое развитие в задачах математической физики, в частности, в задаче теории переноса излучения [12, 31]. Преимуществами методов Монте-Карло являются, прежде всего, физическая наглядность и возможность решать задачи со сложной геометрией, а также возможность оценки функционалов и производных от решения задачи без оценки самого решения задачи.
Одним из подходов методов Монте-Карло является сведение исходной задачи к интегральному уравнению, используя формулу Грина. Затем, решение задачи оценивается методами статистического моделирования по локальным областям. Это легко реализуется через функцию Грина, которая может быть выписана явным образом для п ростейших об л астей.
Одной из разновидностей данных методов является алгоритм "блуждания по сферам”, который был впервые предложен Дж. Брауном [2] для приближенной оценки решения задачи Дирихле для уравнения Лапласа. Позднее данный алгоритм был расширен на решение краевых задач для уравнения Гельмгольца. Алгоритм ”блуждания по сферам” основан на моделировании точек последовательного выхода винеровского процесса из максимальных сфер, целиком лежащих в рассматриваемой области.
Методы Монте-Карло нашли широкое применение при решении краевых задач для уравнений в частных производных. Хорошо изучены следующие направления: построение оценок решений краевых задач для уравнения Лапласа, уравнения Гельмгольца с вещественным параметром [22]. Были разработаны так же алгоритмы оценки решений краевых задач для параболических уравнений. Краевые задачи для параболических уравнений решались, в основном, двумя
4
способами. Первый способ заключается в использовании алгоритма “блуждания по границе”. Второй способ основывается на применении преобразования Лапласа. Хотя, метод ”блуждания по границе” обладает безусловными преимуществами, но проигрывает методу ’ блуждания по сферам” при решении задач со сложной геометрией и не позволяет оценивать решение задачи с переменными коэффициентами. Решение задач с использованием преобразования Лапласа налагает дополнительные ограничения на обратное преобразование Лапласа, а именно, поведение решения преобразованной задачи должно убывать экспоненциально на бесконечности.
Все это делает привлекательным вопрос о применении преобразования Фурье при решении задач параболического типа. Во-первых, данный метод позволяет применить алгоритм ”блуждания по сферам" и, во-вторых, не налагает ограничений на поведение решения преобразованной задачи на бесконечности. Основная трудность использования преобразования Фурье связана с решением преобразованной задачи, которая является краевой задачей для уравнения Гельмгольца с комплексным параметром.
В настоящей работе приведены алгоритмы оценки решений краевых задач для уравнения Гельмгольца с комплексным параметром, а также построены оценки решений краевых задач для уравнения теплопроводности. Были также получены оценки производных по времени, пространственной переменной и параметру от решения задач для уравнения теплопроводности, а также производных по пространственной переменной и комплексному параметру для уравнения Гельмгольца. Причем, для оценки производных по пространственным переменным была использована нецентральная функция Грина краевой задачи для уравнения Гельмгольца в шаре.
Существует еще один метод применения статистического моделирования при решении краевых задач для эллиптических уравнений.
Это метод — ’'блуждания по решетке”, который может также интерпретироваться как дискретный вариант метода "блуждания по сферам”. В данной работе этот метод применяется для построения глобальных оценок решений задач для уравнения теплопроводности и уравнения Гельмгольца с комплексным параметром. Глобальные оценки решений краевых задач для уравнения теплопроводности были получены с использованием дискретного и непрерывного преобразований Фурье. Были получены также оценки трудоемкости и оптимальных параметров данных методов.
На защиту выводятся следующие результаты:
1. Получены оценки решений задачи Дирихле и смешанной краевой задачи для уравнения Гельмгольца с комплексным параметром и уравнения теплопроводности для трехмерного и многомерного случая.
2. Получены глобальные оценки решений краевых задач для уравнения теплопроводности и уравнения Гельмгольца с комплексным параметром с использованием метода "блуждания по решетке”. Причем, решение задачи для уравнения теплопроводности оценивалось с использованием:
(а) непрерывного преобразования Фурье
(б) дискретного преобразования Фурье.
Были получены оценки трудоемкостей и оптимальных параметров данных методов.
3. Получены оценки производных по времени, пространственной переменной и параметру для уравнения теплопроводности, а также оценки производных по параметру и пространственной переменной от решения краевой задачи для уравнения Гельмгольца с комплексным параметром.
б
Основные результаты работы были опубликованы в [13]-[16], [24] и [25].
Результаты, изложенные в диссертации, докладывались на семинаре отдела статистического моделирования в физике ИВМиМГ СО РАН, на конференциях молодых ученых ИВМиМГ СО РАН и на 3-ем международном семинаре по математическому моделированию в г. С.-Петербурге.
Пользуясь случаем, автор выражает искреннюю благодарность своему научному руководителю чл.-корр. РАН Г.А.Михайлову за постоянное внимание и руководство работой.
- Київ+380960830922