ML040500671
| ML040500671 | |
| Person / Time | |
|---|---|
| Site: | Catawba  |
| Issue date: | 02/23/2004 |
| From: | Eltawila F NRC/RES/DSARE |
| To: | Black S Division of Systems Safety and Analysis |
| Voglewede, J | |
| References | |
| Download: ML040500671 (3) | |
Text
C - 1 Attachment C Changes made to FRAPCON to predict the radioactive fission product release from the Duke MOX lead test assembly The following is a line by line listing of the changes made to FRAPCON in order to make the runs described in the above paper.
In the subroutine initial.f Add the variable ngasmod to the namelist /frpcon/
After the line igas=0 Add the line ngasmod=1 In the common file comms.h Add the variable ngasmod to the common /inpti/
In the subroutine frpcon.f Add the variable ngasmod and bu to the call to fgasre In the subroutine fgasre.f Add the variables ngasmod and bu to the subroutine statement Change the line if(ngasr.gt.10) go to 211 To if(ngasmod.eq.1) go to 211 Add the variables gasflg and bu to the call to ans54 After the call to ans54, add the lines, call ans54cs (it,brnup,dp,jpow,nt,im,rv,dv,nr,ProblemTime
+, qaxnorm,rc,rdotcs
+,ir1,qmpy,jst,releascs,den,dco,ngasr,crad,tfr,nr,pi,gasflg) if(gasflg.eq.1) then if(jpow.eq.1) then csrel=0.0 bpsum=0.0 endif csrel=csrel+bp*rdotcs bpsum=bpsum+bp if(jpow.eq.nt) then csrel=csrel/bpsum write(9,209) ProblemTime(it)/86400.0, csrel*100.0 209 format(2x,'Time(days), Cs Release(%)',2x,f10.1,2x,f10.4) endif endif In the subroutine ans54 Add the variables gasflg and bu to the subroutine statement In the dimension statements add the variables tempksub(50) and busub(50)
In the dimension statements change ansr(10) to ansr(50) f(10) to f(50) flxfac(10) to flxfac(50) pf(10,21) to pf(50,21)
C - 2 prdct(11,21,400) to prdct(50,21,400) flxfacb(11) to flxfacb(50) ansd(176000) to ansd(800000)
Change the line dimension decay(11), releas(11), half(11)
To dimension decay(5), releas(5), half(5), halflife(5) integer gasflg c Half-lives in seconds from Chart of Nuclides, 15th edition c Order is Kr-87, Kr-88, Xe-133, Xe-135, I-131 c
data (halflife(i),i=1,5)/4572.0,10224.0,452995.0,
& 32760.0, 692928.0/
After the line tempk=(tempf+459.65)/1.8 Add the line tempksub(i2)=tempk After the line bup = burnup*prdct(i2,jpow,nt)
Add the line bupsub(i2)=bup After the line 210 continue Add the line if(gasflg.eq.1) then Change the lines do 250 i=1,11 rtime = 1+(i-1)*deltim half(i) = 10.**rtime decay(i) = 1/half(i)
To do 250 i=1,5 half(i)=halflife(i) decay(i)=0.6931472/half(i)
After the line j1 = i+nreg*(jpow-1)+npow*(nt-1)*nreg Add the line ansdlock=ansd(j1)
Change the lines do 310 jx=1,11 xmu = decay(jx)/ansd(j1) taut = ansdrad*ProblemTime(nt)
To do 310 jx=1,5 ansdrad=ansdlock c
D is 2 times higher for Cs and 7 times higher for I c
NUREG/CR-2507 if(jx.gt.4) ansdrad=ansdlock*7.0 xmu = decay(jx)/ansdrad taut = ansdrad*dt(nt)
Change line frac=frac*wf To
C - 3 frac1=frac*wf Change line releas(jx) = releas(jx)+frac To releas(jx) = releas(jx)+frac1 Change line do 330 i=1,11 To do 330 i=1,5 After the line 330 continue Add the lines timday=ProblemTime(nt)/86400.0 if(iwrote.eq.0) write(7,1001) iwrote=1 1001 format(4x,'Time',7x,'Days',3x,'Rod Burnup',20x,'Gas Release',
& 'Fractions',/,4x,'Step',15x,'(MWd/MTU)',5x,
& 'Kr-87',6x,'Kr-88',5x,'Xe-133',5x,'Xe-135',6x,'I-131')
write(7,1002) nt-1, timday, bu, (releas(ii),ii=1,5) 1002 format(2x,i5,2x,2(f10.1,2x),8(e10.4,1x))
endif Add a new subroutine called ans54cs which is identical to ans54, except for the following changes In the subroutine statement, the name should be changed from ans54 to ans54cs After the line ansd(j1) = ansd(j1)*100.**(bup/28000.)
add the line ansd(j1) = ansd(j1)*2.0 Remove everything between the lines 210 continue And return