2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




На страницу 1, 2  След.
 
 Трехдиагональная система для уравнения теплопроводности
Сообщение31.05.2013, 23:26 
Здравствуйте,
Пытаюсь решить уравнение теплопроводности, использую неявную разностную схему.
В учебнике Самарского написано, что все уравнения будут иметь один и тот же вид:
Цитата:
$\gamma y_{i+1}^{n+1} - (1+2\gamma) y_i^{n+1} + \gamma y_{i-1}^{n+1}=-F_i^n$, $i=1,2,\ldots, N-1$, $y_0^{n+1} = \mu_1(t_{n+1})$, $y_N^{n+1} = \mu_2(t_{n+1})$, (13)
где $\gamma = \tau/h^2$, $F_i^n=y_i^n +\tau \varphi_i^n$. Эту систему можно решить методом прогонки (см. п. 7 § 4 ч. I), так как условия устойчивости прогонки выполнены.

Изображение
Понятно, что трехдиагональная сисема получится
Не понятно, как получить правую часть $F$ этой системы. Ведь в нее входят пока еще не известные $y^n_i$
Помогите, пожалуйста, разобраться.

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 00:12 
Там крайне неудачные обозначения. Используйте тупо неявную схему: $\dfrac{y^{\text{нов}}_i-y^{\text{стар}}_i}{\tau}=\dfrac1{h^2}(y^{\text{нов}}_{i-1}-2y^{\text{нов}}_i+y^{\text{нов}}_{i+1})+\varphi_i$. Ну там с точностью до знака при фях -- лень вникать.

-- Сб июн 01, 2013 01:17:44 --

А, я не врубился:

shupenker в сообщении #731041 писал(а):
Ведь в нее входят пока еще не известные y_in

Что значит "неизвестные" -- ведь это же значения на старом временнОм слое и, соотв., давно (уже шаг прошёл) как известны.

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 06:12 
Тогда решить систему уравнений для всей квадратной области не получится? Нужно решать для каждого отсчета t отдельную систему? И не будет ли слишком быстро накапливаться погрешность в таком случае он нижнего к верхним слоям?

 
 
 
 Posted automatically
Сообщение01.06.2013, 07:54 
Аватара пользователя
 i  Тема перемещена из форума «Помогите решить / разобраться (М)» в форум «Карантин»
Причина переноса: формулы не оформлены $\TeX$ом

shupenker, наберите все формулы $\TeX$ом. Картинку, содержащую цитату из учебника, можете оставить. Инструкции по оформлению формул здесь или здесь (или в этом видеоролике).
После исправлений сообщите в теме Сообщение в карантине исправлено, и тогда тема будет возвращена.

 
 
 
 Posted automatically
Сообщение01.06.2013, 12:32 
 i  Тема перемещена из форума «Карантин» в форум «Помогите решить / разобраться (М)»

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 12:37 
Аватара пользователя
shupenker в сообщении #731141 писал(а):
Тогда решить систему уравнений для всей квадратной области не получится? Нужно решать для каждого отсчета t отдельную систему?
Почему решить систему уравнений для всей квадратной области не получится?

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 12:56 
shupenker в сообщении #731141 писал(а):
Тогда решить систему уравнений для всей квадратной области не получится? Нужно решать для каждого отсчета t отдельную систему?

Погуглите на "схему переменных направлений" или, что то же, "метод дробных шагов".

shupenker в сообщении #731141 писал(а):
И не будет ли слишком быстро накапливаться погрешность в таком случае он нижнего к верхним слоям?

Конечно, погрешность будет накапливаться, причём (в относительном масштабе) экспоненциально; это общее свойство для всех дифференциальных уравнений вообще. Но это -- лишь если нас интересует именно временнАя зависимость. Если же мы интересуемся лишь предельным значением, т.е. используем такого рода схему как итерационную для решения уравнения Пуассона, то никакого накопления погрешностей быть не может в принципе, и это -- тоже общее свойство всех итерационных схем вообще.

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 14:39 
Мне нужно сделать именно неявную схему для неоднородного уравнения теплопроводности.
В явной схеме можно даже без системы уравнений обойтисть, наверное, применяя шаблон р.с. снизу от начальных условий.
Но в случае неявной схемы, чтобы найти значения на $j+1$ временном слое необходимо вычислить значения на $j$-м. Следовательно получается по системе уравнений на каждый временнОй слой. Это нормальный подход или можно как-то лучше\проще сделать?

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 14:47 
Аватара пользователя
shupenker в сообщении #731263 писал(а):
Но в случае неявной схемы, чтобы найти значения на $j+1$ временном слое необходимо вычислить значения на $j$-м.

В случае явной схемы, чтобы найти значения на $j+1$ временном слое, тоже надо знать значения на $j$-м.

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 14:57 
TOTAL в сообщении #731269 писал(а):
В случае явной схемы, чтобы найти значения на временном слое, тоже надо знать значения на -м.


на $j$-м слое из начального условия $u(x,0)=...$можно найти же и никакой системы уравнений не будет

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 15:02 
Аватара пользователя
shupenker в сообщении #731275 писал(а):
на $j$-м слое из начального условия $u(x,0)=...$можно найти же и никакой системы уравнений не будет
Для неявной схемы тоже все можно найти из начального условия.

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 15:10 
TOTAL в сообщении #731277 писал(а):
Для неявной схемы тоже все можно найти из начального условия.


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

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 15:16 
Аватара пользователя
shupenker в сообщении #731280 писал(а):
Расскажите тогда, пожалуйста, как все сделать "по уму". В итоге мне нужно получить трехдиагональную систему и ее правую часть с числовыми коэффициентами, а то все общими фразами, не понятно :-)
Рассказываю: в приведенной выше цитате из Самарского трехдиагональная система уравнений записана. Чего ещё надо?

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 15:52 
TOTAL в сообщении #731282 писал(а):
Рассказываю: в приведенной выше цитате из Самарского трехдиагональная система уравнений записана. Чего ещё надо?


shupenker в сообщении #731041 писал(а):
Не понятно, как получить правую часть этой системы. Ведь в нее входят пока еще не известные $y^i_n$


Начальное условие задает чиловые значения правой части только для самого нижнего временнОго слоя.

shupenker в сообщении #731263 писал(а):
Но в случае неявной схемы, чтобы найти значения на временном слое $j+1$необходимо вычислить значения на $j$


Вот я и подумал, что для этого нужно решать систему для каждого слоя

 
 
 
 Re: Трехдиагональная система для уравнения теплопроводности
Сообщение01.06.2013, 15:59 
shupenker в сообщении #731263 писал(а):
Но в случае неявной схемы, чтобы найти значения на $j+1$ временном слое необходимо вычислить значения на $j$-м.

Зачем вычислять уже вычисленное?

 
 
 [ Сообщений: 22 ]  На страницу 1, 2  След.


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group