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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

126
Readln(k0);
Writeln('Введите энергию активации химической реакции, E');
Readln(E);
Writeln('Введите начальную температуру, T0');
Readln(T0);
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);
{определяем расчетный шаг сетки по времени}
tau:=t_end/100.0;
{определяем поле температуры в начальный момент времени}
for i:= 1 to N do
T[i]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия, используя соотношения (2
4)}
alfa[1]:=2.0*tau*lamda/(2.0*tau*(lamda+kapa*h)+ro*c*sqr(h));
beta[1]:=(ro*c*sqr(h)*T[1]+2.0*tau*kapa*h*Te)/(2.0*tau
*(lamda+kapa*h)+ro*c*sqr(h));
{запоминаем поле температуры на предыдущем временном слое}
for i:=1 to N do
Tn[i]:=T[i];
{цикл с постусловием, позволяющий итерационно вычислять поле
температуры, вследствие наличия нелинейности в самом уравнении
теплопроводности (56)}
repeat
{запоминаем поле температуры на предыдущей итерации}
for i:=1 to N do
Ts[i]:=T[i];
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2
to N-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления системы
уравнений с трехдиагональной матрицей}
ai:=lamda/sqr(h);

Страницы