Математический форум Math Help Planet
http://mathhelpplanet.com/

Разностная схема для плоскости и сферы
http://mathhelpplanet.com/viewtopic.php?f=64&t=55507
Страница 1 из 1

Автор:  mihailothepooh [ 24 авг 2017, 22:19 ]
Заголовок сообщения:  Разностная схема для плоскости и сферы

Имеется уравнение переноса в приближении метода сферических гармоник в одномерной геометрии в плоской симметрии:
[math]\left\{\!\begin{aligned}
& \frac{ d f_{1} (x) }{ dx } = G_{0} f_{0} (x) + Q(x) \\
& 2\frac{ d f_{2} (x) }{ dx } + \frac{ d f_{0} (x) }{ dx } = 3G_{1} f_{1} (x) \\
& 3\frac{ d f_{3} (x) }{ dx } + 2\frac{ d f_{1} (x) }{ dx } = 5G_{2} f_{2} (x) \\
& 3\frac{ d f_{2} (x) }{ dx } = 7G_{3} f_{3} (x)
\end{aligned}\right.[/math]


и сферической симметрии:

[math]\left\{\!\begin{aligned}
& \left( \frac{ d }{ dr } + \frac{ 2 }{ r } \right) f_{1} (r) = G_{0} f_{0} (r) + Q(r) \\
& 2\left( \frac{ d }{ dr } + \frac{ 3 }{ r } \right) f_{2} (r) + \frac{ d }{ dr }f_{0} (r) = 3G_{1} f_{1} (r) \\
& 3\left( \frac{ d }{ dr } + \frac{ 4 }{ r } \right) f_{3} (r) + 2\left( \frac{ d }{ dr } - \frac{ 1 }{ r } \right) f_{1} (r) = 5G_{2} f_{2} (r) \\
& 3\left( \frac{ d }{ dr } - \frac{ 2 }{ r } \right) f_{2} (r) = 7G_{3} f_{3} (r)
\end{aligned}\right.[/math]


Можно ли построить разностную схему, которая бы описывала и плоский и сферический случай одновременно?
Имеется ввиду, что в схему входят такие характеристики ячеек как площади граней, объёмы и т.д. и если мы подставляем величины для плоских ячеек, то получаем схему для плоскости, а если для сферических ячеек, то схема начинает описывать сферический случай

Автор:  Xmas [ 25 авг 2017, 04:26 ]
Заголовок сообщения:  Re: Разностная схема для плоскости и сферы

Едва ли возможно совместить в одной схеме плоский и сферически-симметричный случаи. Даже если вы привлекаете для каждого узла сетки "площади граней" и "объёмы", вам ещё нужно сообщить в расчётную схему, каким образом пересчитывать потоки от одной ячейки к другой. Там все эти [math]1\slash r[/math] и возникнут. В уравнениях они просто сразу в явном виде представлены.

Ещё довод.
Если взять, например, уравнение теплопроводности, и рассчитывать по нему теплоизоляцию, то для плоских и цилиндрических поверхностей всегда найдётся толщина изоляции, при которой потери тепла сколь угодно малы. А для горячего шара - не найдётся. Даже бесконечная теплоизоляция, причём уже прогретая полностью (т.е., в стационарном состоянии), не помогает, шар будет вечно терять тепло (с потоком [math]Q=4\pi r \lambda \Delta T[/math]).

Опять же, принципиально разное поведение, которое обеспечивается за счёт членов с [math]1\slash r[/math], независимо от размеров ячеек.

Автор:  mihailothepooh [ 25 авг 2017, 06:19 ]
Заголовок сообщения:  Re: Разностная схема для плоскости и сферы

Спасибо за ответ, но для уравнения теплопроводности такую схему построить можно, если записать уравнение в следующем виде

[math]\left\{\!\begin{aligned}
& \frac{ 1 }{ r^{2} } \left( \frac{ d }{ dr }r^{2} f_{1}(r) \right) = G_{0}f_{0}(r) + Q(r) \\
& \frac{ d }{ dr }f_{0}(r) = 3G_{1}f_{1}(r)
\end{aligned}\right.[/math]


Затем первое уравнение умножить на [math]4\pi r^{2}[/math] и проинтегрировать по i-той ячейке от [math]r_{i - 1\slash2 }[/math] до [math]r_{i + 1\slash2 }[/math]. Второе уравнение просто записать на грани ячейки , например [math]r_{i + 1\slash2 }[/math]. Получим

[math]\left\{\!\begin{aligned}
& 4 \pi r_{i + 1\slash2 }^{2}f_{1}^{i + 1\slash2 } - 4 \pi r_{i - 1\slash2 }^{2}f_{1}^{i - 1\slash2 } = \left( \frac{ 4 }{ 3 } \pi r_{i + 1\slash2 }^{3} - \frac{ 4 }{ 3 } \pi r_{i - 1\slash2 }^{3} \right) \left( G_{0}^{i}f_{0}^{i} + Q^{i} \right) \\
& \frac{ f_{0}^{i + 1} - f_{0}^{i} }{ l_{i + 1, i} }= 3G_{1}^{i + 1\slash2 } f_{1}^{i + 1\slash2 }
\end{aligned}\right.[/math]


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

Страница 1 из 1 Часовой пояс: UTC + 3 часа [ Летнее время ]
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
http://www.phpbb.com/