Главная Журналы поверхностях вырезов дополнительные источниковые члены прибавляются к заданным в тех контрольных объемах, которые получают тепловой поток из вырезов. Детали такого представления граничных условий подробно описаны в п. 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.,
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 |