Главная  Журналы 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 [ 45 ] 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99

поверхностях вырезов дополнительные источниковые члены прибавляются к заданным в тех контрольных объемах, которые получают тепловой поток из вырезов. Детали такого представления граничных условий подробно описаны в п. 7.7.2. Заметим, что при определении источниковых членов в приграничных контрольных объемах идентификаторы массивов SC(I,J) и SP(I, J) присутствуют и в правых частях выражений на ФОРТРАНе, что необходимо для учета дополнительных источниковых членов в этих контрольных объемах. Если опустить SC(I,J) или SP(I,J) в правых частях соответствующих выражений, то обычные источниковые члены [заданные соотношением (8.76)] для этих контрольных объемов не будут учтены.

Наконец, короткий теплоизолированный участок на верхней границе реализуется установкой КВС = 2.

8.4.3. Дополнительные имена на ФОРТРАНе

ABYV - отношение площади грани к объему КО;

COND - теплопроводность к [см. (8.7а)];

CONS - ABYV/RES;

DX, DY - размеры контрольных объемов по осям х и у;

GBYD - тепловая (термическая) проводимость;

НЕ - коэффициент теплоотдачи [см. (8.7а)];

QW - постоянная плотность теплового потока д,, через нижнюю границу [см. (8.7а)];

RES - тепловое (термическое) сопротивление;

Т (I, J) - температура Т;

TINF - температура окружающей среды [см. (8.7)];

TW1, TW2- температуры границ Г„, и Т2 [м- (8-7)];

YM - координата у для середины грани контрольного объема.

8.4.4. Листинг подпрограммы ADAPT

сссссссссссссссссссссссссссссссссссссссссссссссссссссссссс

SUBROUTINE ADAPT

с----

с- EXAMPLE 4 - CONDUCTION IN А COMPLEX CYLINDRICAL GEOMETRY

С------

$INCLUDE:COMMON

DIMENSION T(NI,NJ) EQUIVALENCE (F(1,1,1),T(1,1))



ENTRY GRID

HEADER=CONDUCTION IN A COMPLEX CYLItiDRICAL GEOMETRY

PRINTF=PRINT4

PLOTF=PL0T4

M0DE=2

R(l)=0.5

Ll = 14

Ml=li

YV(2)=0.

DY=1.5/FLOAT{Ml-2) DO 10 J--3,M1

YV{J)=YV(J-1)rDY 10 CONTINUE XU(2)=--0. DX=1./6. DO 20 1=3,8

XU(I)-XU(I-l)+DX 20 CONTINUE

CONSTRUCT CONTROL-VOLUME FACES TO MATCH CIRCULAR BOUNDARY DO 30 1=9,13 J=19-I

YM--0.5* (YV(J) +YV(J-1) )

XU(I)=2.-SQRT(1.-(YM-1.5)* * 2.) 30 CONTINUE

XU(Ll)=2. RETURN

ENTRY BEGIN

TITLE(1)= TEMPERATURE

CALL INTA4(KSOLVE(1),1,KPRINT(1),1,KPLOT(1),1,LAST,3) CALL DATA6(COND,2.,TWl,200.,TW2,100.,QW,1.,HE,5.,

TINF, 20 .

DO 100

= 1,M1

DO 100 I

= 1,L1

T (I, J)

=TW2

CONTINUE

DO 110

=2,M2

T(l,J)

=TW1

CONTINUE

DO 120

= 2, L2

T(I,1)

=TW1

CONTINUE RETURN

С *-

* ** * *

ENTRY OUTPUT

DO 200 lUNIT-IUl,IU2

IFdTER.EQ. 0) WRITE (lUNIT, 210) 210 FORMAT(2X, ITER,2X, T(4,5) ,4X, 1(6,10) , 1 3X,T(10,4),3X,T(12,4))

WRITE(lUNIT,22 0) ITER,T(4,5),T(6,10),T(10,4),T(12,1! 220 FORMAT(3X,12,1P4E10.2) 200 CONTINUE



IF(ITER.EQ.LAST) THEN CALL PRINT

COME HERE TO FILL IBLOCK(I,J) BEFORE CALLING PLOT DO 230 J=2,M2 DO 230 1=2,L2

IF(X(I).LT.0.5.AND.Y(J).GT.1.0) IBLOCK(I,J)=1 IF(X(I).LT.0.5.AND.Y(J).LT.0.5) IBLOCK(I,J)=1 IF(X(I).GT.l.O.AND.Y(J).GT.0.5) IBLOCK(I,J)=1 230 CONTINUE

DO 240 J=5,9 DO 240 1=8,17-J IBLOCK(I,J)=0 240 CONTINUE

CALL PLOT

ENDIF RETURN

ENTRY PHI DO 300 J=2,M2 DO 300 1=2,L2 GAM(I,J)=COND

IF(X(I).LT.0.5.AND.Y(J).GT.1.0) GAM(I,J)=0.

IF(X(I).LT.0.5.AND.Y(J).LT.0.5) GAM(I,J)=0.

IF(X(I).GT.l.O.AND.Y(J).GT.0.5) GAM(I,J)=BIG

SC(I,J)=50.

SP(I,J)=-4. 300 CONTINUE

DO 310 J=5,9 DO 310 1=8,17-J

GAM(I,J)=COND 310 CONTINUE

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 320 J=8,M2

SC(5,J)=SC(5,J)+QW/XCV(5) 320 CONTINUE

DO 330 J=2,4

ABYV=1./XCV(5)

RES=1./НЕ+0.5*XCV(5)/GAM(5,J) CONS=ABYV/RES SC(5,J)=SC(5,J)+CONS*TINF SP(5,J)=SP(5, J)-CONS 330 CONTINUE

DO 340 1=2,4

GBYD=GAM(I,7)/(0.5*YCV(7) )

ABYV=RV(8)/YCVR(7)

SC(I,7)=SC(I,7)+ABYV*GBYD*TW1

SP(I,7)=SP(I,7)-ABYV*GBYD

ABYV=RV(5)/YCVR(5)

RES=1./HE+0.5*YCV(5)/GAM(I,5)

CONS=ABYV/RES

SC(I,5)=SC(I,5)+C0NS*TINF





0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 [ 45 ] 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99