Математический форум 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/ |