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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

108
Writeln('Введите температуру внешней среды относительно границы х
= L, Te2');
Readln(Te2);
Writeln('Введите приведенную степень черноты на границе х = 0, eps1');
Readln(eps1);
Writeln('Введите приведенную степень черноты на границе х = L, eps2');
Readln(eps2);
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;
{определяем alfa начальный прогоночный коэффициент на основе левого
граничного условия, используя соотношение (48)}
alfa[1]:=2.0*tau*lamda/(2.0*tau*(lamda+kapa1*h)+ro*c*sqr(h));
{запоминаем поле температуры на предыдущем временном слое}
for i:=1 to N do
Tn[i]:=T[i];
{цикл с постусловием, позволяющий итерационно вычислять поле
температуры, всл
едствие наличия нелинейности в левом граничном
условии}
repeat
{определяем beta начальный прогоночный коэффициент на основе
левого граничного условия, используя соотношение (48), при этом
начинаем итерационный цикл по левому граничному условию}
d:=T[1];
beta[1]:=(ro*c*sqr(h)*Tn[1]+2.0*tau*kapa1*h*Te1+2.0*tau*eps1
*sigma*h*(sqr(sqr(Te1))-sqr(sqr(d))))/(2.0*tau*(lamda
+kapa1*h)+ro*c*sqr(h));

Страницы