Двумерное уравнение теплопроводности. Фролова Ксения. 6 курс

Материал из Department of Theoretical and Applied Mechanics
Версия от 14:50, 19 января 2016; Фролова Ксения (обсуждение | вклад) (Реализация)

(разн.) ← Предыдущая | Текущая версия (разн.) | Следующая → (разн.)
Перейти к: навигация, поиск

Постановка задачи

Необходимо решить задачу Коши для двумерного уравнения теплопроводности (дифференциальное уравнение в частных производных второго порядка, которое описывает распределение температуры в заданной области пространства и его изменение во времени.) с использованием средств параллельного программирования на основе MPI. Задача решается для однородного уравнения теплопроводности (система теплоизолирована) в области [0..L]x[0..L]:
[math]\frac{\partial U}{\partial t} - a^2(\frac{\partial^2 U}{\partial x^2}+\frac{\partial^2 U}{\partial y^2}) = 0[/math]
[math]U(t=0) = U_0[/math]
при граничных условиях:
[math] U(x,y,t) = \begin{cases}T_0, x=0\\ T_1, x!=0 \end{cases}[/math]
Используемые величины параметров:
[math]L=1, T_0=100, T_1=0[/math]

Используемый метод

В вычислительных системах с распределенной памятью процессоры работают независимо друг от друга. Для организации параллельных вычислений в таких условиях необходимо иметь возможность распределять вычислительную нагрузку и организовать информационное взаимодействие (передачу данных) между процессорами. Параллельное программирование служит для создания программ, эффективно использующих вычислительные ресурсы за счет одновременного исполнения кода на нескольких вычислительных узлах. Для создания параллельных приложений используются параллельные языки программирования и специализированные системы поддержки параллельного программирования, такие как MPI и OpenMP. Итак, MPI - это библиотека передачи сообщений, собрание функций на C/C++ (или подпрограмм в Фортране), облегчающих коммуникацию (обмен данными и синхронизацию задач) между процессами параллельной программы с распределенной памятью. Акроним MPI установлен для Message Passing Interface (интерфейс передачи сообщений). Под параллельной программой в рамках MPI понимается множество одновременно выполняемых процессов. Все процессы порождаются один раз, образуя параллельную часть программы. Каждый процесс работает в своем адресном пространстве, никаких общих переменных или данных в MPI нет. Процессы могут выполняться на разных процессорах, но на одном процессоре могут располагаться и несколько процессов (в этом случае их исполнение осуществляется в режиме разделения времени).

Реализация

При решении поставленной задачи будем использовать замену частных производных в дифференциальных уравнениях их разностными аналогами. Сеточный метод, основанный на замене в дифференциальном уравнении производных конечными разностями, называют методом конечных разностей, а сеточную схему такого метода - конечно-разностной.
По аналогии с одномерной задачей для уравнения теплопроводности вводим явную конечно-разностную схему. Область [0..L]x[0..L] разбивается на подобласти согласно количеству процессов в выполняемой параллельной программе. На каждом полученном таким способом интервале процесс интегрирования осуществляется отдельным процессом, при этом в связи с использованием явной схемы соседние процессы должны обмениваться крайними значениями, полученными на предыдущем шаге, для выполнения следующего шага.
Программа для решения двумерного уравнения теплопроводности: программа

Результаты

Найдено решение однородного уравнения теплопроводности в двумерной постановке для следующей сетки узлов: 300х300.

Количество процессов [-] Время рассчета [сек]
1 40.2082
3 13.7626
5 8.38831
7 6.56195
15 3.08675
35 3.90614

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