SUBROUTINE CALC C COMMON /b_1/TEMPD(5,5),A(5,5),AX(5,5),B(5,5),BX(5,5), 1 C(5,5),CX(5,5),D(5,5),DX(5,5),ND,NTHE(6),FRAC(10,5) COMMON /b_2/AMPSSP(10),AMPSCR(10),OHM(10),RT(10), 1 FK(10,10),RK(10,10) COMMON /b_3/AMPS(10),VELZ(10),P1(10),P2(10),P3(10),P4(10), 1 X(10),VOL(10,800),IND(10),DT COMMON /b_4/NC,NT,AREA(10),AREARS(10),TSP(10),TQ(10),IQ(10), 1 TEMPZP,TEMP(10,800) COMMON /b_5/TOTVOL(10),ARCON(10),PROVEL(6,5),ALPHA(10) COMMON /b_6/A1(6),B1(6),C1(6),D1(6),A2(6),B2(6),C2(6),D2(6) COMMON /b_7/A3(6),B3(6),C3(6),D3(6),A4(6),B4(6),C4(6),D4(6) COMMON /b_8/A5(6),B5(6),C5(6),D5(6),A6(6),B6(6),C6(6),D6(6) COMMON /b_9/PROV1(6),PROV2(6),PROV3(6),PROV4(6),PROV5(6),PROV6(6) COMMON /b_10/FIELD1(6),FIELD2(10),BETA(10) C COMMON/CONN/IT(10),THER(10),DUMB(10,800),XPT(10) COMMON/PP/ A7(6),B7(6),C7(6),D7(6),VELP(6),CURR1(6) COMMON/TMSCAL/DTOLD,TMPMRK,DTNU,ASTOP COMMON/SAVEP/IC3,IC4,IC5,IC6 DIMENSION ENTH(10),ENMAG(10),ENOHM(10) DIMENSION Y(10),SUMV(10) DIMENSION AMPS1(10),AMPS2(10),AMPSM(10),RTM(10) DIMENSION EST(10),TAMP(800),AVOL(10),SUMVN(10) C C TIME =0.0 DO 4786 KC=1,NC,1 AVOL(KC)=3.1416*((P2(KC)+P3(KC))**2-P2(KC)**2)*P4(KC) 4786 XPT(KC)=0.0 C IC3=0 IC4=0 IC5=0 IC6=0 IC478=0 C TEMPZ=4.2 DT=DTOLD WRITE(6,101)NC WRITE(7,101)NC C 2 DO 8 KC=1,NC 4 DO 6 I =1,800 TEMP(KC,I)=TEMPZ VOL(KC,I)=0.0 TOTVOL(KC)=0.0 6 CONTINUE AMPS(KC)=AMPS SP(KC) AREA R=1.0/AREA(KC) AREARS(KC)=AREA R*AREA R X(KC)=0.0 RT(KC)=0.0 ENOHM(KC)=0.0 ENTH(KC)=0.0 IQ(KC)=0 IND(KC)=0 8 CONTINUE C RSUM=0.0 VOLSUM=0.0 VTSUM=0.0 TENTH=0.0 TENOHM=0.0 TEMAG=0.0 TOTENG=0.0 IF(NC.GT.1)GO TO 9898 WRITE(1,*)TIME,TEMPZ WRITE(2,*)TIME,AMPS SP(1) WRITE(3,*)TIME,TIME WRITE(4,*)TIME,TIME 9898 DO 200 IFI=1,NC,1 IF(FRAC(IFI,1).NE.0.0) GO TO 200 IND(IFI)=1 RHO1=RHO(TEMPZ,IFI) CP1=CP(TEMPZ,IFI) VOL(IFI,1)=3.142*((P2(IFI)+P3(IFI))**2-P2(IFI)**2)*P4(IFI) RT(IFI)=3.142*P2(IFI)*2.0*RHO1/(P4(IFI)*P3(IFI)) IQ(IFI)=1 200 CONTINUE WRITE(6,103)TIME C 10 DO 20 KC=1,NC ENMAG(KC)=0.0 14 DO 16 JC=1,NC ENMAG(KC)=ENMAG(KC)+0.5*FK(KC,JC)*AMPS(KC)*AMPS(JC) 16 CONTINUE IF(TSP(KC).LE.1.0E-06) AMPS SP(KC)=0.0 VOLTS=-(AMPSSP(KC)-AMPS(KC))*OHM(KC) WRITE(6,104)AMPS(KC),ENMAG(KC),ENOHM(KC),VOLTS,TEMPZ C VOL_START=0.0 IF(KC.EQ.1)WRITE(7,114) WRITE(7,116)KC,TIME,AMPS(KC),VOL_START,ENOHM(KC),VOLTS,TEMPZ, *VOL_START 114 FORMAT(1H ,'No TIME AMPS V INT E(INT) E(EXT) V(EXT) TEMP *EMP NORM VOL') 115 FORMAT(1X,I2,1X,F8.4,1X,1G9.2,6X,1G9.2,1X,1G9.2,3(1X,F6.1)) C VTSUM=VTSUM+VOLTS TENMAG=TENMAG+ENMAG(KC) TOTENG=TOTENG+ENMAG(KC) Y(KC)=AMPS(KC) AMPS1(KC)=AMPS(KC) SUMV(KC)=0.0 20 CONTINUE C IF(NC.EQ.1)GO TO 22 C WRITE(6,105)TOTENG,TENTH,TENMAG,TENOHM,VTSUM,VOLSUM 22 CONTINUE C C C 24 DO 90 I=1,NT CALL DEKM(Y,EST,TIME,DT,NC,AMPS1) IF(DT.EQ.0.0) RETURN C C WRITE(6,103)TIME REMP=FLOAT(I)/20. IREMP=INT(REMP+.01) PTZ=ABS(REMP-FLOAT(IREMP)) IF(PTZ.LT..001)WRITE(6,110) C 26 DO 28 KC=1,NC AMPS(KC)=Y(KC) AMPS2(KC)=AMPS(KC) AMPSM(KC)=0.50*(AMPS1(KC)+AMPS2(KC)) 28 CONTINUE C RSUM=0.0 VOLSUM=0.0 VTSUM=0.0 TENTH=0.0 TENMAG=0.0 TENOHM=0.0 TEMPMX=0.0 TOTENG=0.0 C C 30 DO 80 KC=1,NC ENMAG(KC)=0.0 32 DO 34 JC=1,NC ENMAG(KC)=ENMAG(KC)+0.5*FK(KC,JC)*AMPS(KC)*AMPS(JC) 34 CONTINUE C C APB=ABS(AMPS(KC)/AMPS1(KC)) IF(APB.NE.1.0) GO TO 4783 SS=-1.0/DT GO TO 4784 4783 SS=ALOG(APB) SS=-SS/DT 4784 CONTINUE ER=(AMPS1(KC)**2-AMPS2(KC)**2)*RT(KC)/(2.0*SS) ENTH(KC)=ENTH(KC)+ER EX=(AMPS1(KC)**2-AMPS2(KC)**2)*OHM(KC)/(2.0*SS) IF(TIME.LE.TSP(KC))EX=0.0 ENOHM(KC)=ENOHM(KC)+EX VOLTS=-(AMPSSP(KC)-AMPS(KC))*OHM(KC) C ICON=0 IF(IT(KC).NE.0.AND.TEMP(IT(KC),1).GT.THER(KC)) ICON=1 IF(ICON.EQ.1) GO TO 36 IF(IQ(KC).GE.1) GO TO 36 IF(TIME.GE.TQ(KC)) GO TO 36 IF(ABS(AMPS(KC)).LT.AMPSCR(KC)) GO TO 76 36 IQ(KC)=1 RT(KC)=0.0 C M=IND(KC) IF(IND(KC).GE.1) GO TO 40 IF(ICON.EQ.1) GO TO 469 CALL NORVOL(KC,I) GO TO 470 469 CALL NEWVOL(KC,I) 470 CONTINUE 40 CONTINUE M=I C C DO 70 J=1,M,1 C C TEMP1=TEMP(KC,J) LL=2 IF(TEMP1.LE.150.0) LL=5 IF(TEMP1.LE. 50.0) LL=10 DDT=DT/LL 50 DO 60 L=1,LL RHO1=RHO(TEMP1,KC) CP1 = CP(TEMP1,KC) F=1.0 IF(J.EQ.M) F=0.5 TEMP2=TEMP1+AMPS(KC)*AMPS(KC)*AREARS(KC)*F*DDT*RHO1/CP1 TEMP1=TEMP2 60 CONTINUE TEMP(KC,J)=TEMP2 C C RT(KC)=RT(KC)+RHO(TEMP(KC,J),KC)*VOL(KC,J)*AREARS(KC) 70 CONTINUE SUMV(KC)=SUMV(KC)+VOL(KC,I) SUMVN(KC)=(SUMV(KC)/AVOL(KC))*100. IF(ICON.EQ.1.AND.IC478.EQ.0) WRITE(6,478) 478 FORMAT(X,' THE NORMAL REGION OF THIS COIL IS NOW A PLANE FRONT 1 STARTING AT THE INSIDE FACE') IF(ICON.EQ.1)IC478=1 VINT=AMPS(KC)*RT(KC) C WRITE(6,107)TIME, * AMPS(KC),RT(KC),VINT,ENTH(KC),ENMAG(KC), 1 ENOHM(KC),VOLTS,TEMP(KC,1),VOL(KC,I),SUMVN(KC) WRITE(7,116)KC,TIME,AMPS(KC),VINT,ENTH(KC),VOLTS, *TEMP(KC,1),SUMVN(KC) 116 FORMAT(1X,I2,1X,F8.4,1X,F8.1,1X,F8.2,2(1X,1PE9.2),1X,0PF6.1, *1X,1PE9.2) c ********* screen monitering if(m.eq. 50)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq. 100)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.150)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.200)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.250)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.300)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.350)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq. 400)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.450)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.500)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.550)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.600)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.650)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq. 700)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) if(m.eq.750)PRINT 1050,m,kc,time,amps(kc),temp(kc,1) 1050 format(1h ,'step= ',i3,' coil ',i2,' time= ',f8.4, 1' cur= ',f8.1,' temp= ',0pf6.1) c ******* IF(KC.GT.1)GO TO 999 WRITE(1,*)TIME,TEMP(1,1) WRITE(2,*)TIME,AMPS(1) WRITE(3,*)TIME,VINT WRITE(4,*)TIME,SUMVN(1) 999 CONTINUE C GO TO 78 C 76 RT(KC)=0.0 WRITE(6,104)AMPS(KC),ENMAG(KC),ENOHM(KC),VOLTS,TEMPZ C C 78 TOTENG=TOTENG+ENTH(KC)+ENMAG(KC)+ENOHM(KC) RSUM=RSUM+RT(KC) VOLSUM=VOLSUM+SUMV(KC) VTSUM=VTSUM+VOLTS TENTH=TENTH+ENTH(KC) TENMAG=TENMAG+ENMAG(KC) TENOHM=TENOHM+ENOHM(KC) TEMPMX=AMAX1(TEMPMX,TEMP(KC,1)) AMPS1(KC)=AMPS2(KC) IF(TSP(KC).LE.TIME ) AMPS SP(KC)=0.0 80 CONTINUE C IF(NC.GT.1) * WRITE(6,105)TOTENG,TENTH,TENMAG,TENOHM,VTSUM,VOLSUM C IF(TEMPMX.LT.2000.) GO TO 84 WRITE(6,102) GO TO 95 84 CONTINUE 86 JC=0 DO 88 KC=1,NC IF(DT.GE.DTNU)GO TO 91 IF(TEMP(KC,1).GT.TMPMRK) DT=DTNU 91 IF(AMPS(KC).GE.1.0) JC=1 88 CONTINUE TAMP(I)=AMPS(1) IF(I.LT.10)GO TO 89 IF(TAMP(I).GT.(0.15*AMPS SP(1)))GO TO 89 IF((TAMP(I-1)-TAMP(I)).LE.ASTOP)DT=DT*2. 89 CONTINUE IF(JC.EQ.0)GO TO 95 90 CONTINUE 95 RETURN C C C 101 FORMAT(/,20X,41HCALCULATION OF CURRENT DECAY IN SYSTEM OF, 1 I3,6H COILS) 102 FORMAT('0 TEMPERATURE TOO HIGH') 103 FORMAT('0'F8.4,'S. AMPS'3X,'COIL R'6X,'V INT'7X, & 'E INTERNAL EMAG' 1 8X,'E EXT.'6X,'V EXT'7X,' TEMP DELTA VOL. NORM VOL'4X) 110 FORMAT('0'8X,'S. AMPS'3X,'COIL R'6X,'V INT'7X, & 'E INTERNAL EMAG' 1 8X,'E EXT.'6X,'V EXT'7X,' TEMP DELTA VOL. NORM VOL'4X) 104 FORMAT(9X,F9.3,1X,3(12X),3(G11.5,1X),F7.1,1X) 105 FORMAT(' ETOT='F12.0,4X,' TOTALS:'1X,12X,4(G11.5,1X),20X,G11.5) C 107 FORMAT(1X,F17.3,1X,6(G11.5,1X),F7.1,1X,2(G11.5,1X)) 107 FORMAT(1X,F8.4,G9.3,1X,6(G11.5,1X),F7.1,1X,2(G11.5,1X)) C C C C C TEMPZ = INITIAL TEMPERATURE OF SYSTEM C AMPS(KC) = CURRENT AT ANY TIME IN COIL KC C ( Y IS ALSO USED AS A SYNONYM FOR AMPS ) C VOLTS = VOLTAGE AT ANY TIME ACROSS A COIL C RT(KC) = INTERNAL RESISTANCE AT ANY TIME OF COIL KC C X(KC) = EXTENT OF PROPAGATION IN COIL KC,AS DEFINED IN C SUBROUTINE NORVOL C ENTH(KC) = TOTAL THERMAL ENERGY AT ANY TIME IN COIL KC C ENOHM(KC) = TOTAL THERMAL ENERGY AT ANY TIME IN EXTERNAL C RESISTANCE ACROSS COIL KC C ENMAG(KC) = TOTAL MAGNETIC ENERGY AT ANY TIME IN COIL KC C TENTH = TOTAL THERMAL ENERGY AT ANY TIME IN ALL COILS C TENOHM = TOTAL THRMAL ENERGY AT ANY TIME IN ALL RESISTANCES C TENMAG = TOTAL MAGNETIC ENERGY AT ANY TIME IN ALL COILS C TOTENG = TOTAL ENERGY AT ANY TIME IN WHOLE SYSTEM C IQ(KC) = INDICATION THAT COIL KC HAS QUENCHED, IF GREATER C THAN OR EQUAL TO UNITY C IND(KC) = INDICATION THAT PROPAGATION HAS OCCURED THROUGH WHO C VOLUME OF COIL, IF GREATER THAN OR EQUAL TO UNITY C AMPS1(KC) = CURRENT IN COIL KC BEFORE EACH TIME INTERVAL C AMPS2(KC) = CURRENT IN COIL KC AFTER EACH TIME INTERVAL C AMPSM(KC) = MEAN CURRENT IN COIL KC DURING EACH TIME INTERVAL C VOL(KC,I) = VOLUME OF NORMAL CONDUCTOR IN REGION I OF COIL KC C ( REGION I IS THAT VOLUME GENERATED BY PROPAGATION C DURING TIME INTERVAL I ) C SUMV(KC) = TOTAL VOLUME OF NORMAL CONDUCTOR AT ANY TIME IN COI C VTSUM = TOTAL VOLTAGE ACROSS ALL COILS AT ANY TIME C RSUM = TOTAL RESISTANCE OF ALL COILS AT ANY TIME C VOLSUM = TOTAL VOLUME OF NORMAL REGIONS IN ALL COILS AT ANY C TEMP(KC,I)= TEMPERATURE AT ANY TIME OF REGION I IN COIL KC C TEMPMX = MAXIMUM TEMPERATURE IN ALL COILS AT ANY TIME C C END C C FILE MN.CQUENCH.NORVOL(NO) C C PART OF THE CULHAM VERSION OF QUENCH. SEE CQMAIN. C SUBROUTINE NORVOL(KC,I) COMMON /b_3/AMPS(10),VELZ(10),P1(10),P2(10),P3(10),P4(10) 1, X(10),VOL(10,800),IND(10),DT COMMON /b_4/NC,NT,AREA(10),AREARS(10),TSP(10),TQ(10),IQ(10) 1,TEMPZ,TEMP(10,800) COMMON /b_5/TOTVOL(10),ARCON(10),PROVEL(6,5),ALPHA(10) COMMON /b_10/FIELD1(6),FIELD2(10),BETA(10) C COMMON/SAVEP/IC3,IC4,IC5,IC6 C DIMENSION A7(6),B7(6),C7(6),D7(6),VELP(6),CURR1(6) C C THIS ROUTINE CALCULATES THE INCREASE IN THE NORMAL VOLUME IN COIL C VOLUME CALCULATED AS FOR WILSON RHEL/M 151 PAGE 2 666 DELTAX=VELZ(KC)*ABS(AMPS(KC))*DT VELZ(KC)=VELZ(KC) X(KC)=X(KC)+DELTAX C P2 IS THE INTERNAL RADIUS OF THE COIL C P3 IS THE RADIAL DEPTH OF THE WINDING C P4 IS THE AXIAL LENGTH OF THE WINDING C C ZMAX IS USUALLY EQUAL TO THE MEAN CIRCUMFERENCE OF THE WINDING C XMAX AND YMAX ARE USUALLY SET TO HALF THE WIDTH AND P2 C OF THE COIL RESPECTIVELY C ALPHA=RATIO OF PERPENDICULAR PROPOGATION VELOCITY TO VELOCITY C ALONG THE WIRE. TYPICAL VALUE-ALPHA=0.05 C BETA DETERMINES THE POSITION OF THE ORIGIN OF THE NORMAL REGION C IN THE CROSS SECTION OF THE COIL C BETA=1.0,PROPOGATION STARTS IN THE CENTRE OF THE COIL C BETA=0.5,PROPOGATION STARTS AT THE CENTRE OF THE INSIDE OR C OUTSIDE FACE OF THE COIL C BETA =0.25,ORIGIN AT THE CORNER OF THE COIL C TOTST=TOTVOL(KC) PI=3.14159 IF((BETA(KC).GT.0.99).AND.(BETA(KC).LT.1.01))GO TO 90 IF((BETA(KC).GT.0.49).AND.(BETA(KC).LT.0.51))GO TO 91 IF((BETA(KC).GT.0.24).AND.(BETA(KC).LT.0.26))GO TO 92 WRITE(6,8)KC GO TO1 90 XMAX=0.5*P4(KC) YMAX=0.5*P3(KC) ZMAX=PI*(P2(KC)+0.5*P3(KC)) GO TO 93 91 XMAX=0.5*P4(KC) YMAX=P3(KC) ZMAX = PI*(P2(KC)+0.5*P3(KC)) GO TO 93 92 XMAX=P4(KC) YMAX=P3(KC) ZMAX = PI*(P2(KC)+0.5*P3(KC)) 93 DIAG=SQRT(XMAX**2.0+YMAX**2.0) C DIAG IS THE DIAGONAL OF THE CROSS SECTION OF THE COIL C DEFINED BY XMAX AND YMAX AND IS USED TO DETERMINE THE POINT C AT WHICH THE WHOLE CROSS SECTION IS JUST ENCLOSED BY THE C AT WHICH THE WHOLE CROSS SECTION IS JUST ENCLOSED BY THE C NORMAL CONE.IT IS ALSO USED AFTER STATEMENT 15 TO DETERMINE C THE TIME AT WHICH THE WHOLE COIL IS NORMAL RR=X(KC)*ALPHA(KC) RRR=0.0 C DOES THE BASE OF THE CONE OF NORMAL SUPERCONDUCTER EXCEED THE C CROSS SECTION PREVIOUSLY DEFINED C SET RADIUS OF CONE TO DIAGONAL PREVIOUSLY DEFINED IF((RR.LE.XMAX).AND.(RR.LE.YMAX))GO TO 10 IF(RR.LE.DIAG)GO TO 9 RRR=RR RR=DIAG 9 IF(RR.LE.XMAX)GO TO 20 C CALCULATE THE VOLUME OF ONE QUADRANT OF THE DEFINED CROSS SECTION C WHERE THE NORMAL REGION EXCEEDS XMAX C THE NORMAL REGION IS BOUNDED BY A LINE OF LENGTH XMAX FROM THE C CENTRE OF THE CONE AND PERPENDICULAR TO THE Y DIMENSION OF THE C COIL AND A LINE FROM THE CONE CENTRE JOINING THE JUNCTION OF THE C CIRCUMFERENCEOF THE CONE BASE AND A LINE PERPENDICULAR TO THE C FIRST BOUNDRY TO WHICH IT IS JOINED AT POINT XMAX FROM THE CONE C BASE CENTRE ROOTXM=SQRT(RR**2.0-XMAX**2.0) VOL2=XMAX*RR*ROOTXM/(3.0*ALPHA(KC))-XMAX**3.0/(12.0*ALPHA(KC)) 1*ALOG((RR+ROOTXM)/(RR-ROOTXM)) 20 IF(RR.LE.YMAX)GO TO 21 C REPEAT WHERE YMAX IS EXCEEDED ROOTYM=SQRT(RR**2.0-YMAX**2.0) VOL3=YMAX*RR*ROOTYM/(3.0*ALPHA(KC))-YMAX**3.0/(12.0*ALPHA(KC)) 1*ALOG((RR+ROOTYM)/(RR-ROOTYM)) 21 IF((RR.GT.XMAX).AND.(RR.LE.YMAX))GO TO 11 IF((RR.LE.XMAX).AND.(RR.GT.YMAX))GO TO 12 IF((RR.GT.XMAX).AND.(RR.GT.YMAX))GO TO 13 WRITE(6,7)KC GO TO 1 11 VOL1=X(KC)/3.0*RR**2.0*(PI- 2.0*ATAN(ROOTXM/XMAX)) VOL1=VOL1+4.0*VOL2 C IF(IC3.EQ.0)WRITE(6,3)KC C WRITE(6, 3)KC IC3=1 C GO TO 14 12 VOL1=X(KC)/3.0*RR**2.0*(PI-2.0*ATAN(ROOTYM/YMAX)) VOL1=VOL1+4.0*VOL3 C IF(IC5.EQ.0)WRITE(6,5)KC C WRITE(6, 5)KC IC5=1 C GO TO 14 13 IF(RRR.GT.DIAG)GO TO 16 VOL1=1.0/3.0*RR**2.0*X(KC)*(PI-2.0*ATAN(ROOTXM/XMAX)-2.0* 1ATAN(ROOTYM/YMAX)) VOL1=VOL1+4.0*VOL2+4.0*VOL3 C IF(IC6.EQ.0)WRITE(6,6)KC C WRITE(6, 6)KC IC6=1 C GO TO 14 16 CONEL=DIAG/ALPHA(KC) RECL=X(KC)-CONEL RECVOL=XMAX*YMAX*RECL*4.0 VOL1=4.0*VOL2+4.0*VOL3+RECVOL C IF(IC6.EQ.0)WRITE(6,6)KC C WRITE(6, 6)KC IC6=1 C 14 VOL1=2.0* BETA(KC)*VOL1 VOL14=VOL1 GO TO 15 10 VOL1=2.0*BETA(KC)*PI/3.0*RR**2.0*X(KC) DTOT=VOL1 15 IF(X(KC).LE.ZMAX)GO TO 37 XD=X(KC)-ZMAX RREX=XD*ALPHA(KC) IF(RREX.GE.DIAG)GO TO 38 CALL OVLAP(KC,I,XD,RREX,VOLEX1,XMAX,YMAX) VOLEX1=VOLEX1*2.0*BETA(KC) VOL1=VOL1-VOLEX1 C IF(IC4.EQ.0)WRITE(6,4)KC C WRITE(6, 4)KC IC4=1 C 37 VOL(KC,I)=VOL1-TOTVOL(KC) 999 FORMAT(2X,'VOL14='F12.3,5X,'VOLEX1='F12.3, 15X,'VOL1='F12.3,5X,'VOL='F12.3) TOTVOL(KC)=TOTVOL(KC)+VOL(KC,I) GO TO 1 38 DO 2 K=I,NT 2 VOL(KC,K)=0.0 IND(KC)=I-1 WRITE(6, 7)KC C WRITE(7,7)KC C 1 RETURN 4 FORMAT(47HNORMAL REGION OVERLAPS ON OPPOSITE SIDE OF COIL,I3) 3 FORMAT(30X,38HPROPOGATION REACHES X BOUNDARY IN COIL,I3) 5 FORMAT(30X,38HPROPOGATION REACHES Y BOUNDARY IN COIL,I3) 6 FORMAT(30X,44HPROPOGATION REACHES X AND Y BOUNDARY IN COIL,I3) 7 FORMAT(30X,7HCOIL NO,I3,17HCOMPLETELY NORMAL) 8 FORMAT(33HVALUE OF BETA IN DATA FOR COIL NO,I3,1X,7HINVALID) END C C FILE MN.CQUENCH.OVLAP(NO) C C PART OF THE CULHAM VERSION OF QUENCH. SEE CQMAIN. C SUBROUTINE OVLAP(KC,I,XD,RREX,VOLEX1,XMAX,YMAX) COMMON /b_3/AMPS(10),VELZ(10),P1(10),P2(10),P3(10),P4(10) 1, X(10),VOL(10,800),IND(10),DT COMMON /b_4/NC,NT,AREA(10),AREARS(10),TSP(10),TQ(10),IQ(10) 1,TEMPZ,TEMP(10,800) COMMON /b_5/TOTVOL(10),ARCON(10),PROVEL(6,5),ALPHA(10) C C THIS ROUTINE CALCULATES THE INCREASE IN THE NORMAL VOLUME IN COIL C DOES THE BASE OF THE CONE OF NORMAL SUPERCONDUCTER EXCEED THE CROS C SET RADIUS OF CONE TO HALF DIAGONAL PREVIOUSLY DEFINED PI=3.14159 IF((RREX.LE.XMAX).AND.(RREX.LE.YMAX))GO TO 35 IF(RREX.LE.XMAX)GO TO 30 C CALCULATE THE VOLUME OF ONE QUADRANT OF THE DEFINED CROSS SECTION C WHERE THE NORMAL REGION EXCEEDS XMAX C THE NORMAL REGION IS BOUNDED A LINE OF LENGTH XMAX FROM THE CENTRE C OF THE CONE AND PERPENDICULAR TO THE Y DIMENSION OF THE COIL,AND C A LINE FROM THE CONE CENTRE JOINING THE JUNCTION OF THE CIRCUMFERE C OF THE CONE BASE AND A LINE PERPENDICULAR TO THE FIRST BOUNDRY TO C C C WHICH IT IS JOINED AT POINT XMAX FROM THE CONE BASE CENTRE ROEXXM=SQRT(RREX**2.0-XMAX**2.0) VOLEX2=XMAX*RREX*ROEXXM/(3.0*ALPHA(KC))-XMAX**3.0/(12.0*ALPHA(KC)) 1*ALOG((RREX+ROEXXM)/(RREX-ROEXXM)) 30 IF(RREX.LE.YMAX)GO TO 31 C REPEAT WHERE YMAX IS EXCEEDED ROEXYM=SQRT(RREX**2.0-YMAX**2.0) VOLEX3=YMAX*RREX*ROEXYM/(3.0*ALPHA(KC))-YMAX**3.0/(12.0*ALPHA(KC)) 1*ALOG((RREX+ROEXYM)/(RREX-ROEXYM)) 31 IF((RREX.GT.XMAX).AND.(RREX.LE.YMAX))GO TO 32 IF((RREX.LE.XMAX).AND.(RREX.GT.YMAX))GO TO 33 IF((RREX.GT.XMAX).AND.(RREX.GT.YMAX))GO TO 34 GO TO 36 32 VOLEX1=XD/3.0*RREX**2.0*(PI- 2.0*ATAN(ROEXXM/XMAX)) VOLEX1=VOLEX1+4.0*VOLEX2 VOLA=VOLEX1 GO TO 36 33 VOLEX1=XD/3.0*RREX**2.0*(PI-2.0*ATAN(ROEXYM/YMAX)) VOLEX1=VOLEX1+4.0*VOLEX3 VOLB=VOLEX1 GO TO 36 34 VOLEX1=1.0/3.0*RREX**2.0*XD*(PI-2.0*ATAN(ROEXXM/XMAX) 1-2.0*ATAN(ROEXYM/YMAX)) VOLEX5=VOLEX1 VOLEX1=VOLEX1+4.0*VOLEX2+4.0*VOLEX3 VOLC=VOLEX1 GO TO 36 35 VOLEX1=PI/3.0*RREX**2.0*XD VVV=VOLEX1 36 CONTINUE 999 FORMAT(2X,'VOLA='F12.3,5X,'VOLB='F12.3,5X,'VOLC='F12.3, 15X,'VVV='F12.3) RETURN END C ******************************************************************* C FL=KE.QUENCH.NEWVOL C FL TO CALC PROPAGATION OF QUENCH FROM A COMPLETE SURFACE OF A COI C DESIGNED TO ENABLE USE OF QUENCHH WITH CONDUCTING RINGS THAT C TRIGGER THE SPREAD OF A NORMAL REGION. VERSION 1 IN THIS PROGR C THE INTERACTION OF THE CONVENTIONAL CONICAL NORMAL REGION AND C THE LAYER REGIONS IS IGNORED AND BOTH ARE ASSUMED TO EXIST INDEPE C THOUGH THE SPRAED OF THE CONE IS STOPPED WHEN THE LAYER SPREAD C BEGINS C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ SUBROUTINE NEWVOL(KC,I) COMMON /b_3/AMPS(10),VELZ(10),P1(10),P2(10),P3(10),P4(10) 1, X(10),VOL(10,800),IND(10),DT COMMON /b_4/NC,NT,AREA(10),AREARS(10),TSP(10),TQ(10),IQ(10), 1 TEMPZP,TEMP(10,800) COMMON /b_5/TOTVOL(10),ARCON(10),PROVEL(6,5),ALPHA(10) COMMON/CONN/IT(10),THER(10),DELYY(10,800),XX(10) DIMENSION A7(6),B7(6),C7(6),D7(6),VELP(6),CURR1(6) PI=3.14159265 DELTAX=VELZ(KC)*ABS(AMPS(KC))*DT IF(XX(KC).GE.P3(KC)) RETURN DELYY(KC,I)=DELTAX*ALPHA(KC) XX(KC)=XX(KC)+DELYY(KC,I) IF(XX(KC).GE.P3(KC))GO TO 50 VOL(KC,I)=P4(KC)*PI*(P2(KC)+0.5*P3(KC))*2.0*DELYY(KC,I) TOTVOL(KC)=TOTVOL(KC)+VOL(KC,I) RETURN 50 CONTINUE XX(KC)=XX(KC)-DELYY(KC,I) DELYY(KC,I)=P3(KC)-XX(KC) XX(KC)=P3(KC) VOL(KC,I)=P4(KC)*PI*(P2(KC)+0.5*P3(KC))*2.0*DELYY(KC,I) TOTVOL(KC)=TOTVOL(KC)+VOL(KC,I) I1=I+1 DO 2 K=I1,NT,1 2 VOL(KC,K)=0.0 IND(KC)=I-1 WRITE(6,78) KC 78 FORMAT(X,'COIL NO ',I5,' IS ALL NORMAL PLANE WAVE PROPGTN') RETURN END