Разностные методы решения задач теплопроводности. Кузнецов Г.В - 138 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

138
{определяем соответствующий шаг по времени}
tau:=(2.0*a1*a2*Qfr*w*0.5*(ro1+ro2)*sqr(h)-sqr(h)*(lamda1*a2+lamda2
*a1)*(Tfr-Tn[k]))/(2.0*a1*a2*(lamda1*(Tfr-T[k-1])-lamda2
*(T[k+1]-Tfr)));
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=0.0;
beta[1]:=Tc;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2 to k-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления системы
уравнений с трехдиагональной матрицей}
ai:=a1/sqr(h);
bi:=2.0*a1/sqr(h)+1.0/tau;
ci:=a1/sqr(h);
fi:=-Tn[i]/tau;
{alfa[i], beta[i] – прогоночные коэффициенты}
alfa[i]:=ai/(bi-ci*alfa[i-1]);
beta[i]:=(ci*beta[i-1]-fi)/(bi-ci*alfa[i-1]);
end;
{определяем значение температуры на границе фазового перехода}
T[k
]:=Tfr;
{используя соотношение (7) определяем неизвестное поле температуры
в промерзшей зоне грунта}
for i:= k-1 downto 1 do
T[i]:=alfa[i]*T[i+1]+beta[i];
{определяем начальные прогоночные коэффициенты для талой зоны
грунта на основе условия на границе раздела двух сред}
alfa[k]:=0.0;
beta[k]:=Tfr;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= k+1 to N-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления системы
уравнений с трехдиагональ
ной матрицей}
ai:=a2/sqr(h);
bi:=2.0*a2/sqr(h)+1.0/tau;
ci:=a2/sqr(h);

Страницы