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

UptoLike

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

Рубрика:

real*8 Axi(0:Nx), Bxi(0:Nx), Cxi(0:Nx), Dxi(0:Nx),

1 TY1(1:Nx–1)

real*8 a0, b0, d0, an, cn, dn

!г.у. для левой границы ao To = bo T1+do

Axi(0) = a0

Bxi(0) = b0

Dxi(0) = d0

!г.у. для правой границы an Tn = cn Tn–1+dn

Axi(Nx) = an

Cxi(Nx) = cn

Dxi(Nx) = dn

do i = 1, Nx–1

Axi(i) = 2. * (Cps * Ros / tau + Lams / hx**2)

Bxi(i) = Lams / hx**2

Cxi(i) = Lams / hx**2

Dxi(i) = 2. * Cps * Ros * TY1(i) / tau

end do

End

!прогоночные коэффициенты по оси OY

Subroutine ABCD_Y(Ny, Ayj, Byj, Cyj, Dyj, Lams, Ros,

1 Cps, hy, tau, TX1, a0, b0, d0,

2 an, cn, dn)

integer Ny

real*8 Lams, Ros, Cps, hy, tau

real*8 Ayj(0:Ny), Byj(0:Ny), Cyj(0:Ny), Dyj(0:Ny),

1 TX1(1:Ny–1)

real*8 a0, b0, d0, an, cn, dn