Численное решение двумерных нестационарных уравнений теплопроводности. Маслов Е.А. - 19 стр.

UptoLike

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

Рубрика:

Do i = 1, Nx–1 !направление 0Y

!г.у.(III рода)на нижней г. (AE)ao To=bo T1+do

a0 = 1. + Alfae * hy / Lams

b0 = 1.

d0 = Te * Alfae * hy / Lams

!г.у.(III рода)на г.(BC)и(AD) an Tn=cn Tn–1+dn

if (i < = Nxlg) then !граница (BC)

an = 1. + Alfag * hy / Lams

cn = 1.

dn = Tg * Alfag * hy / Lams

else !граница (CD)

an = 1. + Alfae * hy / Lams

cn = 1.

dn = Te * Alfae * hy / Lams

end if

TX1(1:Ny–1) = TXY(i,1:Ny–1) !Ti,1:Ny–1

call ABCD_Y(Ny, Ayj, Byj, Cyj, Dyj, Lams, Ros, Cps,

1 hy, tau, TX1, a0, b0, d0, an, cn, dn)

call TDMA(Ny, Ayj, Byj, Cyj, Dyj, TY)

TXY(i, 0:Ny) = TY(0:Ny)

end Do

TXY(0, 0:Ny) = TXY(1, 0:Ny)

TXY(Nx, 0:Ny) = TXY(Nx–1, 0:Ny)

write(*, '(A, f5.2, A)') ' Time= ', L*tau, '.c'

end DO