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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

167
Writeln(f,'Результат получен с шагом по координате x hx = ',hx:6:4);
Writeln(f,'Результат получен с шагом по координате y hy = ',hy:6:4);
Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4);
Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4);
close(f);
{определяем поле температуры в начальный момент времени}
for i:= 1 to Nx do
for j:= 1 to Ny do
T[i,j]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{запоминаем поле темпера
туры на предыдущем n временном слое}
for i:=1 to Nx do
for j:=1 to Ny do
Tn[i,j]:=T[i,j];
{решаем СЛАУ в направлении оси Ох для определения поля
температуры на промежуточном (n+1/2) временном слое}
for j:=1 to Ny do
begin
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=0.0;
beta[1]:=Th;
{цикл с постусловием, позволяющий итерационно вычислять значения
температуры в области решения, вследствие наличия нелинейности в
самом уравнении теплопроводности}
repeat
{запоминаем поле темпера
туры на предыдущей s итерации}
for k:=1 to Nx do
for l:=1 to Ny do
Ts[k,l]:=T[k,l];
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2 to Nx-1 do
begin

Страницы