Document: NRC Regulatory Guide
Document ID: 2704425a-c58a-45c4-93ab-8761721c3e7a
Document Type: regulatory_guide
Title: Evaluation of Reactor Pressure Vessels with Charpy Upper-Shelf Energy Less Than 50 Ft-Lb
Source: NRC Regulatory Guide Division 1
Source URL: https://www.nrc.gov/docs/ML0037/ML003740038.pdf
Revision Date: 2023-06
Chapter: 
Section ID: RG-1.161
CFR Part: 
CFR Title: 

Content:
=1. 10 B-3 "TT = TMAX*TIMFJIO. C EQUATION FOR THE TEMPERATURE OF THE WATER TWATER(TIME) = TINT+CONST(I)+ CONST(2)*TT + CONST(3)*TT**2 + 1 CONST(4)*TT**3 + CONST(5)*TT**4 DO 150 CRACK= 1,5 K=KO I1O X = ZQ(CRACK)/TM TAU - K*TTfTH**2 DO 120M= 1, 5 W(M) = CONST(M) * 120 CONTINUE DO 130N- 1,4 ALNQ - AL(NQ) AN(N) = 2 * SIN(ALNQ)/(ALNQ + SIN(ALNQ)* COS(ALNQ)) CP(N) = COS(ALNQ (I-X)) Y(N,I) -I -EXP(-(ALNQ**2)*TAU) DO 130 M -2, 5 Y(NM) = TAU**(M-I) - (Y (NM-IyALNQ**2)*(M-1) 130 CONTINUE DO 140N= 1,4 , ALNQ - AL(NQ) SUM(N) - AN(N) * CP(N) * (S(I) * EXP(-(ALNQ**2*TAU)) + S(2) I * Y(N,IYALNQ**2 +2*S(3)* Y(N,2)ALNQ**2 + 3 *S(4) * Y(N,3) 2 /ALNQ**2 +4 *S(5)*Y(N,4)/ALNQ**2) 140 CONTINUE C EQUATION FOR THE QUARTER POINT TEMPERATURES TQ(CRACK,TIME) - TWATER(TIME) - SUM(I) - SUM(2) - SUM(3) - SUM(4) C CONTROL FOR THE CONSTANT KAPPA OPTION IF (CONSTK-.EQ. 1) GO TO 150 C TESTFOR THE ACCURACY OF KAPPA FOR THE GIVEN METAL TEMPERATURE, C IF THE DESIRED ACCURACY IS NOT OBTAINED, ITERATE ON KAPPA C FOR THIS CRACK DEPTH AND TIME. KTEST = 1.030 - (5.97E-7)*((T(CRACKTIME))**2) IF ((ABS(KTEST-K)) JLE. 0.0001) GO TO 150 K=KIEST GOTO 110 150 CONTINUE B-4 RETURN END SUBROUTINE TEXP C Calculate WATER TEMPERATURES Using an 'Exponential Decay" Model REAL B. KTEST, K, KO. SUM(4) INTEGER CRACK, TIME, CONSTK, CONSTE INTEGER Q C EXPONENTIAL DECAY MODEL OF THE WATER TEMPERATURE DO 130 TIME = 1, 10 TT = TMAX*TIMFl0. C EQUATION FOR THE TEMPERATURE OF WATER TWATER(TIME) = TO + DT * (I-EXP(-BE*TT)) DO 130 CRACK i 1,5 K=KO 100 WSQ = BE*TH*THK TAU - K*TT/(TH*TH) DO 120N=1,4 ALNQ = AL(N.Q) B = -DT*((2*SIN(ALNQ)/(ALNQ+(SIN(ALNQ))*(COS(ALNQ)))) I *(EXP(-(ALNQ**2*TAU))-EXP(-WSQ*TAU))I((ALNQ**2/WSQ)-1)) X =i ZQ(CRACKYTH SUM(N) = B * COS(ALNQ*(I -X)) 120 CONTINUE C EQUATON FOR THE -QUARTER POINTS" TEMPERATURE VALUES TQ(CRACKTIME) = TWATER(TIME) - SUM(1) - SUM(2) - SUM(3) - SUM(4) C CONTROL FOR THE CONSTANT KAPPA OPTION IF (CONSTK.EQ. 1) GO TO 130 C TEST FOR KAPPA ACCURACY AND CONTROL OF KAPPA OPTION KTEST = 1.030 -