ML092650370
ML092650370 | |
Person / Time | |
---|---|
Site: | Seabrook ![]() |
Issue date: | 08/12/2009 |
From: | Westinghouse |
To: | Office of Nuclear Reactor Regulation |
References | |
SBK-L-09168 LTR-SGMP-09-100 NP | |
Download: ML092650370 (158) | |
Text
Attachment 7 Westinghouse Electric Company LLC LTR-SGMP-09-100 NP-Attachment, "Response to NRC Request for Additional Information on H*; Model F and Model D5 Steam Generators"
WESTINGHOUSE NON-PROPRIETARY CLASS 3 LTR-SGMP-09-100 NP-Attachment Westinghouse Electric Company Response to NRC Request for Additional Information on H*;
Model F and Model D5 Steam Generators August 12, 2009 Westinghouse Electric Company LLC P.O. Box 158 Madison, PA 15663
© 2009 Westinghouse Electric Company LLC All Rights Reserved
LTR-SGMP-09-100 NP-Attachment Response to NRC Request for Additional Information on H*;
Model F and Model D5 Steam Generators
References:
- 1. NL-09-0547, Vogtle Electric Generating Plant License Amendment Request to Revise Technical Specification(TS) Sections 5.5.9, "Steam Generator (SG) Program" and TS 5.6.10, "Steam Generator Tube Inspection Report for Permanent Alternate Repair Criteria," Southern Company, May 19, 2009.
- 2. RS-09-071, " License Amendment Request to Revise Technical Specifications (TS) for Permanent Alternate Repair Criteria," Exelon Nuclear, June 24, 2009.
Docket Nos. 50-445 and 50-446, License Amendment Request 09-007, Model D5 Steam Generator Alternate Repair Criteria," Luminant, June 8, 2009.
- 4. SBK-L-09118, "Seabrook Station: License Amendment Request 09-03; Revision to Technical Specification 6.7.6.k, "Steam Generator (SG) Program," for Permanent Alternate Repair Criteria (H*)," May 28, 2009.
- 5. Vogtle Electric Generating Plant, Units 1 and 2, Request for Additional Information Regarding Steam Generator Program (TAC Nos. ME1339 and ME1340)," United States Nuclear Regulatory Commission, July 10, 2009.
- 6. Braidwood Station, Units 1 and 2, and Byron Station, Unit Nos. 1 and 2 - Request for Additional Information Related to Steam Generator Permanent Alternate Repair Criteria (TAC Nos. ME1613, ME1614, ME1615, and ME1616)," United States Nuclear Regulatory Commission, July 20, 2009.
- 7. Comanche Peak Steam Electric Station, Units 1 and 2 - Request for Additional Information Regarding the Permanent Alternate Repair Criteria License Amendment Request (TAC Nos. ME1446 and ME1447)," United States Nuclear Regulatory Commission, July 23, 2009.
- 8. WCAP-17071-P, "H*: Alternate Repair Criteria for the Tubesheet Expansion Region in Steam Generators with Hydraulically Expanded Tubes (Model F)," Westinghouse Electric LLC, April 2009.
- 9. WCAP-17072-P, "H*: Alternate Repair Criteria for the Tubesheet Expansion Region in Steam Generators with Hydraulically Expanded Tubes (Model D5)," Westinghouse Electric LLC, May 2009.
- 10. "Vogtle Electric Generating Plant, Units 1 and 2, Request for Additional Information Regarding Steam Generator Program (TAC Nos. ME1339 and ME1340)," United States Nuclear Regulatory Commission, August 5, 2009 2
LTR-SGMP-09-100 NP-Attachment Introduction In response to formal requests for technical specification amendments, References 1, 2, 3 and 4, the USNRC formally requested additional information in References 5, 6 and 7. A preliminary request was received in response to Reference 4. This document provides responses to NRC RAI on the Vogtle, Byron/Braidwood and Comanche Peak requests for a permanent license amendment to implement H*. These plants represent the Model F and Model D5 steam generators for which the H* technical justification is provided in Reference 8 and 9. It is anticipated that similar RAIs may be issued as other LARs are submitted that include other models of SG, specifically Models 44F and 51 F. The intent of these responses is to provide a generic response for all applicable models of SGs to the extent possible, recognizing that there may be specific numerical differences for the models of SG not yet addressed (Model 44F and Model 51 F). If necessary, a second issue of these responses will be provided to specifically address the Model 44F and Model 51 RAI when they are received.
The NRC questions are repeated verbatim for each of the plants who received formal or draft RAI in the tables preceding the response to each question. The current NRC RAIs are specifically in regard to WCAP-17071-P (Model F H*) and WCAP-17072-P (Model D5 H*);
responses are provided primarily for these reports, but additional information is provided as appropriate for the H* reports for the other models of SGs , WCAP-1 7091-P for the Model 44F and WCAP-17092-P for the Model 51 F. Because the reports all utilize the same methodology, model-specific information provided will generally be different in only the numerical information.
Subsequent to the initial issue of the RAI (References 5, 6 and 7), the NRC issued follow-up questions (Reference 10) to question numbers 4, 20 and 24 and an additional question regarding a TS commitment for applying the leakage factors. The responses to the follow-up questions to original question numbers 20 and 24 are included directly in the response to these questions below. The response to RAI#4 and the follow-up question to RAI#4 will be provided under separate cover.
Where references are made in a response to other responses included in this document, the basis for the reference is the RAI received by Vogtle. For example, the Vogtle RAI#20 response applies to the Byron/Braidwood RAI#21 as noted in the tabularization of the questions preceding each response.
3
LTR-SGMP-09-100 NP-Attachment Vogtle 1. Reference 1, page 6-21, Table 6-6: This table contains a RAI number of undefined parametersand some apparent inconsistencies with Table 5-2 on page 5-6. Please define the input parametersin Table 6-6.
WCGS 1. Reference 1, page 6-21, Table 6-6: This table contains a number of undefined parametersand some apparent inconsistencies with Table 5-2 on page 5-6. Please define the input parametersin Table 6-6.
B/B 1. Reference 1, Page 6-21, Table 6-6: This table contains a number of undefined parametersand some apparent inconsistencieswith Table 5-2 on page 5- 6. Please define the input parametersin Table 6-6.
CPSES 1. Reference 1, page 6-21, Table 6-6: This table contains a number of undefined parametersand some apparent inconsistencies with Table 5-2 on page 5-6. Please define the input parametersin Table 6-6.
Seabrook 1. Reference 1, Page 6-21, Table 6-6: This table contains a number of undefined parametersand some apparent inconsistencies with Table 5-2 on page 5-6. Please define the input parametersin Table 6-6.
Response
Table 6-6 in WCAP-17071-P and WCAP-17072-P is provided principally as a reference to provide a bridge to the source of basic design data maintained by Westinghouse and as a historical reference from prior H* reports. Although many of the entries in Table 6-6 are not used in the H* analysis, the table was provided to show traceability to the principal sources of the design data, the Westinghouse Power Capability Working Group (PCWG) sheets and the Systems Standards 1.3F and 1.3, which provide transient response data for component design.
The references to Millstone Unit 3 in WCAP-17071-P and to Byron Unit 2 in WCAP-17072-P reflect that these plants are the limiting plants for the Model F and Model D5 SGs that are candidates for application of H*.
Updated Tables 6-6 for the Model F and Model D5 are provided as Tables RAIl-2 and RAI1-3.
The references in the tables have been updated from those contained in Revision 0 of WCAP-17071-P and WCAP-17072-P.
4
LTR-SGMP-09-100 NP-Attachment Table RAIl-1 Updated Table 6-6 of WCAP-17071-P: Summary of H* Millstone Unit 3 Analysis Mean Input Properties Plant Name Millstone Unit 3 Plant Alpha NEU Plant Analysis Type Hot Leg SG Type F Input Value I Unit Reference Accident and Normal Temperature Inputs a,c,e NOP Thot OF PCWG-06-9 NOP T1ow OF PCWG-06-9 SLB TS AT OF __1.3F SLB CH AT OF 1.3F Shell AT OF 1.3F FLB Primary AT Hi OF 1.3F FLB Primary AT Low OF 1.3F SLB Primary AT OF I_1.3F SLB Secondary AT OF 1.3F Secondary Shell AT Hi OF PCWG-06-9 Secondary Shell AT Low OF PCWG-06-9 Cold Leg AT OF PCWG-06-9 Hot Standby Temperature F oF PCWG-06-9 Operating Pressure Input Faulted SLB Primary Pressure 2560.0 psig 1.3F Faulted FLB Primary Pressure 2642.01 psig 1.3F Normal Primary Pressure 2235.0 psig PCWG-06-9 Cold Leg AP ac,e NSD-RMW 1 '1ps_g 070 NOP Secondary Pressure - Low psig PCWG-06-9 NOP Secondary Pressure - Hi psig PCWG-06-9 Faulted FLB Secondary Pressure 0.0 psig 1.3F Faulted SLB Secondary Pressure 0.0 psig 1.3F Notes. 1. The value for Faulted FLB Primary Pressure used in the H*analysis is 2650 psi which conservatively bounds the limiting value of 2642 psi as identified in SSDC 1.3F for the Model F SGs. The value of 2642 psig for peak primary-secondary pressure differential differs from the value provided in Table 5-3 (2657 psig) reported in WCAP-1 7071-P.
5
LTR-SGMP-09-100 NP-Attachment Table RAIl-2 Updated Table 6-6 of WCAP-17072-P: Summary of H* Byron Unit 2 Analysis Mean Input Properties Plant Name Byron 2 Plant Alpha CBE Plant Analysis Type Hot Leg SG Type D5 Input Value I Unit Reference Accident and Normal Temperature Inputs a,c,e NOP Thot OF PCWG-2741 NOP Tjo, OF PCWG-2741 SLB TS AT OF 1.3, Rev. 2 SLB CH AT OF 1.3, Rev. 2 Shell AT OF 1.3, Rev. 2 FLB Primary AT OF 1.3, Rev. 2 SLB Primary AT OF 1.3, Rev. 2 SLB Secondary AT OF 1.3, Rev. 2 Secondary Shell AT Hi OF PCWG-2741 Secondary Shell AT Low OF PCWG-2741 Cold Leg AT OF PCWG-2741 Hot Standby Temperature* F oF PCWG-2741 Operating Pressure Input Faulted SLB Primary Pressure 2560.0 psig 1.3, Rev. 2 Faulted FLB Primary Pressure 2560.0 psig 1.3, Rev. 2 Normal Primary Pressure 2235.0 psig PCWG-2741 Cold Leg AP a,c,e. NSD-RMW Ipsig 070 NOP Secondary Pressure - Low psig PCWG-2741 NOP Secondary Pressure - Hi psig PCWG-2741 Faulted FLB Secondary Pressure 0.0 psig 1.3, Rev. 2 Faulted SLB Secondary Pressure 0.0 psig 1.3, Rev. 2 Much of the data provided in Table 6-6 is not utilized in the final H* analysis. Table RAI1-3 provides a summary of whether the data is utilized in the reference analysis of H* and in which analysis model it is used (See Figure 1-1 in the respective reports). It is emphasized that changqes made in Tables RAIl-1 and RAI1-2 do not affect the H* results provided in References 7 and 8 of this document.
6
LTR-SGMP-09-100 NP-Attachment Table RAIl-3 Utilization of Data from Table 6-6 Input [ Where Used Accident and Normal Temperature Inputs NOP Thot H* Integrator Spreadsheet NOP T1ow H* Integrator Spreadsheet SLB TS AT Not Used SLB CH AT Not Used Shell AT Not Used FLB Primary AT Hi Not Used FLB Primary AT Low Not Used SLB Primary AT Not Used SLB Secondary AT Not Used Secondary Shell AT Hi H* Integrator Spreadsheet; same as Secondary Fluid Temperature at NOP High Tv, Conditions Secondary Shell AT Low H* Integrator Spreadsheet; same as Secondary Fluid Temperature at NOP Low T.v, Conditions Cold Leg AT Not Used Hot Standby Temperature H* Integrator Spreadsheet Operating Pressure Input Faulted SLB Primary Pressure H* Integrator Spreadsheet Faulted FLB Primary Pressure H* Integrator Spreadsheet Normal Primary Pressure H* Integrator Spreadsheet Cold Leg AP Not Used NOP Secondary Pressure - Low H* Integrator Spreadsheet NOP Secondary Pressure - Hi H* Integrator Spreadsheet Faulted FLB Secondary Pressure Not Used Faulted SLB Secondary Pressure Not Used 7
LTR-SGMP-09-100 NP-Attachment The definitions of the entries in the Table 6-6 of WCAP-17071-P and WCAP17072-P are presented below. Also, discussion is provided regarding the consistency of the values in Table 6-6 of the respective reports with Tables 5-1 through 5-6 of the reports.
NOP Thot The steam generator hot leg temperature at high Tavg normal operating conditions at 100%
power (considered to be the same as the reactor vessel outlet temperature).
Model F: [ ]a,c,e OF at the inlet of the tubes at high Tavg normal operating conditions at 100%
power for Millstone Unit 3 is consistent with the value provided in Table 5-1 (WCAP-17071-P).
Model D5: [ ]a,c,e oF at the inlet of the tubes at high Tavg normal operating conditions at 100% power Byron Unit 2 is consistent with the value provided in Table 5-1 (WCAP-1 7072-P).
NOP TIow The steam generator hot leg temperature at the inlet of the tubes at low Tavg normal operating conditions at 100% power (considered to be the same as the reactor vessel outlet temperature).
Model F: [ ]a,c.e oF is consistent with the value provided in Table 5-1.
Model D5: [ ]acOF is consistent with the value provided in Table 5-1.
SLB TS AT Model F: [ ]a'ce °F, ([ ]ace °F - 70°F ) [ ]a'c'e °F: The steam generator hot and cold leg temperature difference that occurs during a postulated, steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state tubesheet metal temperature and the ambient temperature surrounding the steam generator (assumed to be 700 F). The value of [ ]a,c,e oF is not used in the analysis.
Model D5: [ ]ace OF , ([ a,c,e OF - 70'F)= [ ]aoe OF: The steam generator hot and cold leg temperature difference that occurs during a postulated steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state tubesheet metal temperature and the ambient temperature surrounding the steam generator (assumed to be 700 F). The value of [ ]a,c,, oF is not used in the H* analysis.
SLB CH AT-Model F: 348°F, []a,c,e OF: The steam generator hot and cold leg temperature difference that occurs during a postulated steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state channelhead metal temperature and the ambient temperature surrounding the steam generator (assumed to be 700 F). The value of [ ]a,c,e OF is not used in the H* analysis.
8
LTR-SGMP-09-100 NP-Attachment Model D5: 227OF, ([ ]a,c,e OF: The steam generator hot and cold leg temperature difference that occurs during a postulated steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state channelhead metal temperature and the ambient temperature surrounding the steam generator (assumed to be 70'F). The value of [ ],,e oF is not used in the H* analysis.
Shell AT Model F: ]ac,.e OF -701F) = [ ]a,c,e OF : The steam generator secondary side temperature difference that occurs during a postulated steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state secondary side shell metal temperature and the ambient temperature surrounding the steam generator (assumed to be 70'F). The [ ]a'c°e OF value is not used in the H* analysis.
Model D5: [ ]ac.,e OF. The steam generator secondary side temperature difference that occurs during a postulated steam line break event during the maximum pressure difference across the tubesheet of 2560 psi between the steady-state secondary side shell metal temperature and the ambient temperature surrounding the steam generator (assumed to be 70'F). The [ ]ae-F value is not used in the H* analysis.
The secondary side temperature during a postulated SLB is used in the H* analysis for both the Model F ([ ]ace OF) and Model D5 ([ ]a.c'e OF) SGs.
FLB Primary ATHI The reduction in NOP Thot temperature that occurs during a postulated feedwater line break during the maximum pressure difference across the tubesheet of 2642 psi (Model F), 2560 psi (Model D5) corresponding to the high Tavg plant condition.
Model F: [ ]a,c,e OF is consistent with the value provided in Table 5-3 ([ ]ace OF- 54OF =
[ ]a,,e OF). The - 54 0 F value is not used in the H* analysis.
Model D5: [ ]a'ce °F is consistent with the value provided in Table 5-3 ([]a'ceOF)-
[ ] a'c'eOF [ ace OF). The [ ]ac'e OF value is not used in the H* analysis.
The primary side temperature that occurs during a postulated FLB initiating from the high Tavg plant condition, [ ]ace OF is used in the H* analysis for the Model F SG. The no load temperature of [ ]a,c,e oF is used for the Model D5 SGs.
FLB Primary ATLow The reduction in NOP TIow temperature that occurs during a postulated feedwater line break during the maximum pressure difference across the tubesheet of 2642 psi.
9
LTR-SGMP-09-100 NP-Attachment Model F: [ ]ace IF is consistent with the value provided in Table 5-3 ([ ]a,c,e OF [ ]a,c,e OF =
[]a,c, OF). The -54°F value is not used in the H* analysis.
Model D5: [ 1a'ce°F is not included in WCAP-17072-P. The [ ],,Ce°F value is not used in the H* analysis.
The primary side temperature that occurs during a postulated FLB initiating from the low Tvg plant condition, [ ]ace OF, is used in the H* analysis for the Model F. The no load temperature of [ ]a'ce oF is used for the Model D5 SGs.
SLB Primary AT Model F: The reduction in no load temperature of [ ],,e IF [ ],ce OF) to [ ]a,,,e OF that occurs in the reactor coolant system during a postulated SLB during the maximum pressure difference across the tubesheet of 2560 psi. The value in Table 6-6 should be [ ]a,ce oF to be consistent with SSDC 1.3F and Table 5-2. The [ I ,c,e oF value is not used in the H* analysis.
Model D5: The reduction in no load temperature of [ ]a,,ce OF ([ ]a,c,e OF) to [ ]ac-e OF that occurs in the reactor coolant system during a postulated SLB during the maximum pressure difference across the tubesheet of 2560 psi. The value in Table 6-6 is consistent with SSDC 1.3, Rev 2 and Table!,5-2. The [ ]ac' OF value is not used in the H* analysis.
The primary side temperature that occurs during a postulated SLB, [ ]-,c-e OF, is used in the H*
analysis for the Model F. The primary side temperature, [ ]O'c'e °F, is used for the Model D5 SGs.
SLB Secondary AT Model F: The reduction in no load temperature of [ ]ac*e OF ([ ]a,c~e OF) to [ ]a*,e F that occurs on the secondary side of the steam generator during a postulated SLB during the maximum pressure difference across the tubesheet of 2560 psi. The value in Table 6-6 should be [ ]a,c,eoF to be consistent with Table 5-2.
Model D5: The reduction in no load temperature of [ ]ace OF ([ ]a,,ce OF) to [ ))ac OF that occurs on the secondary side of the steam generator during a postulated SLB during the maximum pressure difference across the tubesheet of 2560 psi. The value in Table 6-6 should be [ ]a'ceoF to be consistent with SSDC 1.3, Rev. 2 and Table 5-2.
As noted above, the secondary side temperature during a postulated SLB is used in the H*
analysis for both the Model F ([ ]a,.,e OF) and Model D5 ([ ]a,ce OF) SGs.
10
LTR-SGMP-09-100 NP-Attachment Secondary Shell ATHi For the Model F SG, [ ]jac,e OF (should be [ lace oF) is the average temperature between the secondary side steam temperature and the feedwater temperature during NOP Hi Tavg operation ([ ]a,c,e OF + [ ]a,c,e 'F). This value is the same as the secondary fluid temperature during high Tavg normal operating conditions. The same value calculated for the Model D5 SGs is [ ]a,c,e OF.
Secondary Shell ATow For the Model F SGs, [ ]a,c,e OF is the average temperature between the secondary side steam temperature and the feedwater temperature during NOP Low Tavg operation ([ Ia,c, OF
+ [ ]a,.,e OF)/2= [ ]a~c,e OF). This value is the same as the secondary fluid temperature during low Tavg normal operating conditions. The same value calculated for the Model D5 SGs is [ ]a,c,e OF.
Cold Leg AT Model F: The temperature difference between the hot and cold leg of the Millstone 3 SGs during NOP Low Tavg is 66.6 0 F. This value is not used in the H* analysis.
Model D5: The temperature difference between the hot and cold leg of the Byron 2/Braidwood 2 SGs during NOP Low Tavg is 63 0 F. This value is not used in the H* analysis.
Hot Standby Temperature The zero load temperature, ],,c,e OF.
This value is used in the H* analysis for both the Model F and Model D5 SGs.
Faulted SLB Primary Pressure The maximum pressure difference that occurs across the tubesheet during a postulated SLB.
Model F: 2560 psig is consistent with the value reported in Table 5-2.
Model D5: 2560 psig is consistent with the value reported in Table 5-2.
Faulted FLB Primary Pressure The maximum pressure difference that occurs across the tubesheet during a postulated FLB.
Model F: The value (2650 psig) used for the Model F SG in the H* analysis bounds the actual FLB pressure differential, 2642 psi identified in SSDC 1.3F. As noted above, Table 5-3 of 11
LTR-SGMP-09-100 NP-Attachment WCAP-1 7071-P should be corrected to 2642 psig for the entries for peak primary-to-secondary pressure.
Model D5: The maximum FLB pressure differential for the Model D5 SGs is 2560 psi.
Normal Primary Pressure The primary side pressure during normal operation.
Model F: 2235 psig is consistent with the absolute primary pressure reported in Table 5-1 of 2250 psia.
Model D5: 2235 psig is consistent with the absolute primary pressure reported in Table 5-1 of 2250 psia.
Cold Leg AP The overall pressure drop that occurs in a steam generator tube as fluid flows through the tube from hot leg to cold leg.
Model F: [ ]a,c,e psig (Millstone 3). This value is not used in the H* analysis.
Model D5: [ ]a~c~e psig (Byron 2/Braidwood 2). This value is not used in the H* analysis.
NOP Secondary Pressure Low The steam pressure on the secondary side of the steam generators for NOP Low Tavg.
Model F: [ ]aoe psig is consistent with the value reported in Table 5-1 as [ ]ace psia.
Model D5: [ ],,ce psig is consistent with the value reported in Table 5-1 as [ ]a,c,e psia.
NOP Secondary Pressure Hi The steam pressure on the secondary side of the steam generators for NOP Hi Tavg.
Model F: [ ]ac, psig is consistent with the value reported in Table 5-1 as [ 1',c,' psia.
Model D5: [ ]a,c,e psig is consistent with the value reported in Table 5-1 as [ ]ace psia.
Faulted FLB Secondary Pressure 0 psig, for the Model F and Model D5 SGs, the steam pressure on the secondary side of a steam generator during a postulated FLB. This value is not used in the H* analysis.
12
LTR-SGMP-09-100 NP-Attachment Vogtle 2. Reference 1, page 6-23, Section 6.2.2.2: Why was the finite RAI element analysis not run directly with the modified temperature distributionrather than running with the linear distribution and scaling the results?
WCGS 2. Reference 1, Section 6.2.2.2: Please explain why the finite element analysis not run directly with the modified temperature distribution ratherthan running with the linear distributionand scaling the results?
B/B 2. Reference 1, Section 6.2.2.2: Why was the finite element analysis not run directly with the modified temperature distributionratherthan running with the lineardistributionand scaling the results?
CPSES 2. Reference 1, Section 6.2.2.2: Please explain why the finite element analysis was not run directly with the modified temperature distribution ratherthan running with the linear distribution and scaling the results?
Seabrook 2. Reference 1, Section 6.2.2.2: Why was the FEA analysisnot run directly with the modified temperature distributionrather than running with the linear distributionand scaling the results?
Response
The finite element analysis was run with the. modified temperature distribution as described in section 6.2.2.2.5 of WCAP 17071-P (Model F) and WCAP-17072-P (Model D5). However, since the modified temperature distribution required a different meshing scheme, the displacement results could not initially be used as inputs to the H* contact pressure analysis. The difficulty in applying the results for the modified temperature distribution is what led to the development of Figures 6-21 through Figure 6-23.
Additional tools were developed to accommodate the displacement results from the modified temperature distribution during the course of refining the analysis procedures to accommodate other steam generator designs. When the actual tubesheet displacements from the modified thermal distribution are used, instead of the results from the linear temperature distribution being scaled, the actual change in H* distance is much less than the [ ]ac.e inches reported in Section 6.2.2.2.5. The value of the adjustment for the thermal distribution effect in H* for the different models of SG in the H* fleet is given in Table RAI2-1. All of the values in Table RAI2-1 assume zero residual contact pressure from tube installation effects. The results in column (3) come from using the scaled linear tubesheet displacements in the H* contact pressure analysis.
The results in column (4) come from using the tubesheet displacements from the modified thermal distribution.
13
LTR-SGMP-09-1 00 NP-Attachment Table RAI2-1 Updated NOP Thermal Offset Factors SG Model Re Thermal Offset Thermal Offset SM eeport (Scaled Result) (Applied)
(1) (2) (3) (4)
Model F WCAP 17071-P a,c,e in. ]a,c,e in.
Model D5 WCAP 17072-P [ a,c,e in. [ ]a in.
Model 44F WCAP 17091-P 0.00 in. 0.00 in.
Model 51F WCAP 17092-P 0.00 in. 0.00 in.
Vogtle 3. Reference 1, page 6-38, Section 6.2.3: Why is radial RAI displacement the "figure of merit" for determining the bounding segment? Does circumferentialdisplacement not enter into this? Why is the change in the tube hole diameter not the "figure of merit?"
WCGS 3. Reference 1, Section 6.2.3: Please explain why radial displacementis the "figure of merit"for determining the bounding segment ? Does circumferentialdisplacementnot enter into this? Why is the change in tube hole diameter not the "figure of merit?"
B/B 3. Reference 1, Section 6.2.3: Why is radialdisplacement the "figure of merit" for determining the bounding segment? D oes circumferentialdisplacement not enter into this? Why is the change in tube hole diameter not the "figure of merit?"
CPSES 3. Reference 1, Section 6.2.3: Please explain why radial displacementis the "figure of merit" for determining the bounding segment. Does circumferentialdisplacement not enter into this? Why is the change in tube hole diameternot the "figure of merit?"
Seabrook 3. Reference 1, Section 6.2.3: Why is radialdisplacement the "figureof merit" for determining the bounding segment?
Does circumferential displacement not enter into this? Why is the change in tube hole diameter not the "figure of merit?"
Response
Radial displacement is calculated in two different ways in the H* analysis: the global scale and the local scale.
On the scale of the steam generator itself, otherwise referred to as the global scale, the radial displacement of the entire tubesheet is calculated. At this level, the tubes are not included in the structural model and there is no direct way to calculate the change in the tube hole diameter. It is not possible to calculate the change in the tube hole diameter at the global scale because the tube holes physically do not exist but are represented by 14
LTR-SGMP-09-100 NP-Attachment the effective anisotropic material properties of the tubesheet. Therefore, from the global perspective, it is not possible to use the change in hole diameter as a "figure of merit."
On the local scale, the displacements of the tube and tubesheet collar are calculated in the radial and circumferential directions. As described in Section 6.3 of WCAP-17071-P (Model F) and WCAP -17072-P (Model D5), the expansion of a hole of diameter D in the tubesheet at a radius R is given by:
Radial: AD = D {dUR(R)/dR}
Circumferential: AD = D {UR(R)/R}
UR is available directly from the finite element results as the global radial displacement for a given point in the tubesheet. The value for dUR(R)/dR is obtained by numerical differentiation of the combined displacement field. The maximum expansion of a hole in the tubesheet is in either the radial or circumferential direction. Typically, these two values are within [ ]a,c.e% of each other. However, it is clear from the relationship described in Section 6.3 that maximizing the radial displacement at the global scale (i.e.,
increasing UR) results in maximizing the circumferential and radial displacement of the tubesheet material at the local scale.
The connection between the local and global scales is the global radial displacement of the tubesheet. This is because the applied boundary conditions and the structures attached to the tubesheet have the greatest effect on the displacement in the radial direction. The tubesheet displacement in the circumferential direction due to the applied pressure loading is typically constant at a small negative value on the order of [ ]a,ce inch or less. Therefore, the radial displacement is the best indicator, or "figure of merit,"
of the effect of different operating conditions on tubesheet displacement due to pressure loading. Radial displacement is also a good "figure of merit" for the change in tube hole diameter because maximizing the global radial displacement leads to the maximum calculated circumferential and radial tubesheet displacements at the local level.
Therefore, the global radial displacement of the tubesheet as described in Section 6.2.3 is the appropriate choice for determining the bounding segment of the tubesheet with respect to the contact pressure analysis.
15
LTR-SGMP-09-100 NP-Attachment Vogtle 4. Reference 1, page 6-69: In Section 6.2.5.3, it is concluded RAI that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact bet ween the tube and tubesheet, over a portion of the circumference, for a distance above H *. Is the conclusion in Section 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
WCGS 4. Reference 1, page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintaincontact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H *. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
B/B 4. Reference 1, Page 6-7: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wallcracks at the H* distance, there may be little or no net pressure acting on the tube fo r some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps prior to the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H *. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location.
16
LTR-SGMP-09-100 NP-Attachment CPSES 4. Reference 1, page 6-70: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H *. Is the conclusion in Section 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
Seabrook 4. Reference 1, Page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H *. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
to Reference 10 provided follow-up questions to RAI#4. In Reference 10, the follow questions Reference 10 provided follow-up questions to RAI#4. In Reference 10, the follow questions to RAI#4 were titled RAI#1. The follow-up questions from Reference 10 are reproduced below:
RAI#1 (Reference 10)
- 1. Address following questions as partof response to RAI 4 (Vogtle):
- a. Clarify the nature of the finite element model ("slice" model versus axisymmetric SG assembly model) used to generate the specific information in Tables 6-1, 2, and 3 (and accompanying graph entitled "EllipticalHole Factors') of Reference 6-15. What loads were applied? How was the eccentricity produced in the model? (By modeling the eccentricity as part of the geometry? By applying an axisymmetric pressure the inside of the bore?) Explain why this model is not scalable to lower temperatures.
- b. Provide table showing maximum delta diameters (total diameter distortion) and maximum eccentricities (maximum diameter minus minimum diameter) from the 3 dimensional (3-D) finite element analysis for normal operating and steam line break (SLB), for model F and D5.
- c. In Figure 2 of the White Paper, add plot for original relationship between reductions in contact pressure and eccentricity as given in Reference 6-15 in the graph 17
LTR-SGMP-09-100 NP-Attachment accompanying Table 6-3. Explain why this originalrelationship remains conservative in light of the new relationship. Explain the reasons for the differences between the curves.
- d. When establishing whether contact pressure increases when going from normal operating to steam line break conditions, how can a valid and conservative comparison be made if the normal operating case is based on the original delta contact pressure versus eccentricity curve and the SLB case is based on the new curve?
Response
The responses to RAI#4 of References 5, 6 and 7 and to the follow-up question, RAI#1 of Reference 10, will be provided under separate cover.
18
LTR-SGMP-09-1 00 NP-Attachment
- 5. Reference 1, Page 6-87: -Are the previously calculated scale RAI Vogtle factors and delta D factors in Section 6.3 conservative for steam line break and fee dwater line break? Are they conservative for an intact dividerplate assumption? Are they conservative for all values of primarypressure minus crevice pressure that may exist along the H* distance for intact tubes and tubes with through-wall cracks at the H* distance? How is tube temperature (TT) on page 6-87 determined? For normal operating conditions, how is the TT assumed to vary as function of elevation?
WCGS 5. Reference 1, Section 6.3: Please verify if the previously calculatedscale factors and delta D factors in Section 6.3 conservative for (1) steam line break (SLB) and a feedwater line break (FLB); (2) an intact dividerplate assumption; and (3) all values of primary pressureminus crevice pressure that may exist along the H* distance for intact tubes and tubes with through-wall cracks at the H* distance.
B/B 5. Reference 1, Page 6-86, Section 6: Are the previously calculatedscale factors and delta D factors in Section 6.3 conservative for steam line break and feedwater line break (FLB) ? Are they conservative for an intact divider plate assumption? Are they conservative for all values of primary pressure minus crevice pressure that may exist along the H*
distance for intact tubes and tubes with through-wall cracks at the H* distance?
CPSES 5. Reference 1, Section 6.3, Page 6-86: Please verify if the previously calculatedscale factors and delta D factors in Section 6.3 are conservative for (1) a steam line break (SLB) and a feedwater line break (FLB); (2) an intact divider plate assumption; and (3) all values of primary pressure minus crevice pressure that may exist along the H* distance for intact tubes and tubes with through-wall cracks at the H*
distance.
Seabrook 5. Reference 1, Section 6.3: Are the previously calculatedscale factors and delta D factors in Section 6.3 conservative for steam line break (SLB) and feed line break (FLB) ? Are they conservative for an intact divider plate assumption? Are they conservative for all values of primary pressure minus crevice pressure that may exist along the H
- distance for intacttubes and tubes with throuqh-wall cracks at the H* distance?
19
LTR-SGMP-09-100 NP-Attachment
Response
Note: The page reference, 6-87 (Model F) appearto be incorrect in the question. Section 6.3 begins on page 6-83 (Model F). The page reference for the Model D5 is correct as stated in the Byron/BraidwoodRAI.
- 1) The previously calculated scale factors and delta D factors in Section 6.3 are conservative for all of the analyzed Model F and Model D5 conditions, including normal operating, steam line break and feedwater line break, as appropriate. Use of the contact pressure data described in Reference RAI5-1 would increase the tube-to-tubesheet contact pressure in the Model F H* analysis.
- 2) The previously calculated scale factors and delta D factors in Section 6.3 are conservative for an intact divider plate assumption. The results on page 6-87 assume that a greater level of weld and divider plate degradation exists in the SG (DPF = [ ]ace) than in the rest of the H* structural analysis (DPF = [ ]ac,e (DPF = Divider Plate Factor).
- 3) The previously calculated scale factors and delta D factors in Section 6.3 are conservative for all values of primary pressure minus crevice pressure regardless of their location within the tubesheet. This is because the calculated scale factors and delta D factors applied unit pressure loads to either side of the tube and weld structure in the model such that either the primary side of the tube and tubesheet were pressurized or the secondary side of the tube and tubesheet (including the crevice) were pressurized. In the reference elliptical hole study, the gap elements that were selected for use in the two dimensional study also penalized the tube-tubesheet contact pressure by preventing line on line contact between the tube outside diameter (OD) and the tubesheet/sleeve inside diameter (ID) which results in a lower estimate of the tube-to-tubesheet contact pressure.
- 4) The tube temperature (TT) is assumed to be equal to the primary fluid temperature for the operating condition of interest. The tube temperature is assumed to not vary as a function of elevation within the tubesheet.
RAI5
References:
RAI5-1. LTR-NRC-09-26, "LTR-SGMP-09-66 P-Attachment, "White Paper: Low Temperature Steam Line Break Contact Pressure and Local Tube Bore Deformation Analysis for H*"(Proprietary)," May 13, 2009 20
LTR-SGMP-09-100 NP-Attachment WCGS 6. Reference 1, page 6-87: Please provide information on how RAI the tube temperature (TT) on page 6-87 was determined?
For normal operating conditions, how is the TT assumed to vary as function of elevation?
B/B 6. Reference 1, Page 6-9: How is tube temperature(TT) on page 6-96 determined? For normal operatingconditions, how is the TT assumed to vary as a function of elevation?
CPSES 6. Reference 1, page 6-96: Pleaseprovide information on how the tube temperature (TT) on page 6-96 was determined. For normal operatingconditions, please explain how the TT is assumed to vary as function of elevation.
Seabrook 6. Reference 1, Page 6-8: How is tube temperature (TT) on page 6-87 determined? For normal operating conditions, how is the TT assumed to vary as function of elevation?
Response
The tube temperature (TT) is assumed to be equal to the primary fluid temperature for the operating condition of interest. The tube temperature is assumed to not vary as a function of elevation within the tubesheet.
21
LTR-SGMP-09-100 NP-Attachment Vogtle 6. Reference 1, page 6-97, Figure 6-75:-Contactpressures for RAI nuclear plants with Model F SGs are plotted in Figure 6- 75, but it is not clear what operating conditions are represented in the plotted data. olease clarify.
WCGS 7. Reference 1, page 6-97, Figure 6-75: Contact pressures for nuclear plants with Model F SGs are plotted in Figure 6- 75, but it is not clear what operating conditions are represented in the plotted data. Please clarify.
B/B 7. Reference 1, Page 6-104, Figure 6-77: Contactpressures for nuclear plants with Model D5 steam generatorsare plotted in Figure 6-77, but it is not clear what operating conditions are represented for the plants shown in the plotted data; please clarify.
CPSES 7. Reference 1, page 6-104, Figure 6-77: Contactpressures for nuclearplants with Model D5 SGs are plotted in Figure 6-77, but it is not clear what operatingconditionsare represented for the plants shown in the plotted data. PI ease clarify.
Seabrook 7. Reference 1, Page 6-97, Figure 6-75: Contact pressuresfor nuclear plants with Model F steam generatorsare plotted in Figure 6-75, but it is not clear what operating conditions are represented in the plotted data, please clarify.
Response
Figure 6-75 (WCAP-17071-P) shows the contact pressure results for the fleet of Model F steam generators for the main feedwater line break (FLB), main steam line break (SLB) and normal operating low average temperature (NOP Low Tavg) conditions. Figure RAI6-1 provides an update of Figure 6-75 in WCAP-17071-P with an expanded legend that describes each curve in the figure.
Figure 6-77 (WCAP-17072-P) shows the contact pressure results for the fleet of Model D5 steam generators for the main feed line break (FLB), main steam line break (SLB) and normal operating low average temperature (NOP Low Tavg) conditions. Figure RAI6-2 provides an update of Figure 6-77 in WCAP-17072-P with an expanded legend that describes each curve in the figure.
22
LTR-SGMP-09-100 NP-Attachment Figure RAI-6-1 Revised Figure 6-75 (WCAP-17071-P) a,c,e Figure RAI6-2 Revised Figure 6-77 (WCAP-17072-P) a,c,e 23
LTR-SGMP-09-1 00 NP-Attachment Vogtle 7. Reference 1, page 6-113, Reference 6-5: This reference RAI seems to be incomplete; please provide a complete reference.
WCGS 8. Reference 1, page 6-112, Reference 6-5: This reference seems to be incomplete. Please provide a complete reference.
B/B 8. Reference 1, Page 6-120, Reference 6-5: This reference seems to be incomplete; please provide a complete reference.
CPSES 8. Reference 1, page 6-120, Reference 6-5: This reference appearsto be incomplete. Pleaseprovide a complete reference.
Seabrook 8. Reference 1, Page 6-112, Reference 6-5: This reference seems to be incomplete; please provide a complete reference.
Response
The complete reference is:
Slot, Thomas, "Stress Analysis of Thick Perforated Plates," TECHNOMIC Publishing Company, Inc., Westport, Connecticut, 1972.
Vogtle 8. Reference 1, page 6-113, Reference 6-15: Table 6-3 in RAI Reference 6-15 (SM-94-58, Rev 1) appearsinconsistent with Table 6-2 in the same reference. Explain how the analysis procgresses from Table 6-2 to Table 6-3.
WCGS 9. Reference 1, page 6-113, Reference 6-15: Table 6-3 in Reference 6-15 (SM-94-58, Revision 1) appears inconsistent with Table 6-2 in the same reference. Please explain how the analysis progressesfrom Table 6-2 to Table 6-3.
B/B 9. Reference 1, Page 6-121, Reference 6-15: Table 6-3 in Reference 6-15 (SM-94-58, Rev. 1) appears inconsistent with Table 6-2 in the same reference. Explain how the analysis progressesfrom Table 6-2 to Table 6-3.
CPSES 9. Reference 1, Page 6-121, Reference 6-15: Table 6-3 in Reference 6-15 (SM-94-58, Revision 1) appears to be inconsistent with Table 6-2 in the same reference. Please explain how the analysis progressesfrom Table 6-2 to Table 6-3.
Seabrook 9. Reference 1, Page 6-113, Reference 6-15: Table 6-3 in Reference 6-15 (SM-94-58, Rev 1) appearsinconsistent with Table 6-2 in the same reference. Explain how the analysis orocyresses from Table 6-2 to Table 6-3 24
LTR-SGMP-09-100 NP-Attachment
Response
The values for initial and final eccentricity for the contact pressure ratio of 0.91 listed in Table 6-3 of Reference 6-15 (SM-94-58, Rev. 1) are calculated as follows using the values from Table 6-2:
Initial Eccentricity = (Dmax-Dmin)/ ]a,,,e inch Tube Hole ID = [
a,c,e Final Eccentricity = ((Hole Delta D (900) - Hole Delta D (00))/ ]ace inch Tube Hole ID) =[ ]a,c,e The values for eccentricity in Table 6-2 of the reference should have been divided by the nominal diameter of the tubesheet hole [ ]a,c,e inch) to be consistent with Table 6-3.
Vogtle 9. Reference 1, page 8-9, Figure 8: -There is an apparent RAI discontinuityin the plotted data of the adjustment to H* for distributedcrevice pressure,please provide any insight you may have as to why this apparentdiscontinuity exists.
WCGS 10. Reference 1, page 8-9, Figure 8-1: There is an apparent discontinuity in the plotted data of the adjustment to H* for distributedcrevice pressure. Please provide any insight yoL may have as to why this apparentdiscontinuity exists.
B/B 10. Reference 1, page 8-9, Figure 8-1: There is an apparent discontinuity in the plotted data of the adjustment to H* for distributedcrevice pressure. Please provide any insight you may have as to why this apparentdiscontinuity exists.
CPSES 10. Reference 1, page 8-9, Figure 8-1: There is an apparent discontinuity in the plotted data of the adjustment to H* for distributed crevice pressure. Please provide any insight you may have as to why this apparentdiscontinuity exists.
Seabrook 10. Reference 1, Page 8-9, Figure 8-1: There is an apparent discontinuity in the plotted data of the adjustment to H* for distributed crevice pressure, please provide any insight you may have as to why this apparentdiscontinuity exists.
Response
Figure 8-1 (WCAP-17071-P, WCAP-17072-P) summarizes the variability cases run to determine the H* value response to variation of the input parameters (ctT, cLTS, ET, ETS) individually or in combination. The values of the variables were chosen to provide sufficient data to define the potential surface of interactions between the variables. No attempt was made to bias the variables in a manner that would yield H* values in the range between 3.8 inches and 4.2 inches; therefore this gap is coincidental.
25
LTR-SGMP-09-100 NP-Attachment Figure RAI9-1 shows a composite of the Pcrev corrections for all of the models of SGs considered, Models F, D5, 44F and 51F SGs under H* (Ref: WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP-17092-P). Figure RAI9-1 shows the same characteristic shape of the Pcrev correction but also shows that the H* responses are different for the different structures. The "apparent discontinuity" in the curve for the Model F is much less pronounced for the Model D5 and other models of SG and, in the case of the Model 44F, is populated by calculated data points. Because the same analysis methods are employed for all of the Model-specific structures, it isconcluded that the apparent discontinuity in Figure 8-1 of WCAP-17071-P and WCAP-17072-P is related principally to the structural response of the specific SG model being addressed, and does not imply a potential calculation error.
Figure RAI9-1 also shows that in each of the structures considered, there are steps in the Pcrev correction curves (e.g., between 3.8 and 4.2 inches in the Model F, at about 6.6 inches in the Model D5, at about 3.5 inches and 4.5 inches for the Model 44F and 51 F).
To investigate the step in the curve between initial predictions- of H* and the Pcrev correction, several cases were considered for the Model F SGs for H* values between 3.8 inches and 4.2 inches as a typical case to evaluate the issue generically. These cases were synthesized by adjusting the values of the four influencing parameters (aT, cLTS, ET and ETS), based on interpolation among existing variabilities, in an attempt to yield H* values in this range. Each of the four parameters was adjusted in at least one case to meet this objective.
The following are the additional cases that were examined:
Input Parameters H*(raw) Pcrev Comment aTS ETS aT ET f- _ _c_ _
1 -1 -2 -2 Original Case 5 4 0 0
-1 0 -3.25 0
-1 0 -3 -5 4.5 0 0 0 5 4 0 -1 4.5 0 0 -1 5 0 0 0 Original Case Figure RAI9-2 shows the results of this study. The Pcv correction values are essentially constant within the narrow range of initial H* predictions that define the step in the overall curve, Figure RAI9-1, except for a single point at approximately [ ]a,c,e inches.
As discussed below, the interpolation between the limited number of points representing the crevice pressure distribution and the fixed number of points representing the thickness of the tubesheet leads to isolated conditions at which the integration scheme cannot converge to a single value. A minor departure (less than about 0.005 inch) in either direction results in convergence of the integration. The point at [ ]a,,,e inches is at such a condition. It does not suggest that the crevice pressure correction is undefined at that location.
26
LTR-SGMP-09-100 NP-Attachment As described in each of the H* WCAP reports, the correction for Pcrev is an iterative process. Following the initial prediction of H*, which assumes that a tube separation is located at the primary face of the tubesheet and, therefore, assumes the crevice pressure is distributed over the entire thickness of the tubesheet, the calculation process depicted in Figure 1-1 of the report is repeated but with the crevice pressure distributed over the length of the initial prediction of H*. The resulting prediction of H* will exceed the initial prediction. This process is iterated until the input values and output values of H* converge to the same number. The convergence criteria are set to 2 decimals because the H* distance cannot practically be measured to the second decimal. In some instances, depending on the specific combination of input parameters that lead to the initial prediction of H*, the variation of H* is less than the convergence criteria. In that case, the default is at the larger value of the Pcrev.
The H* integrator model utilizes discrete, dimensionally fixed points through the thickness of the tubesheet to represent the tube to tubesheet contact pressure. The representation of the distributed crevice pressure as discussed in Section 6 of the report utilizes a discrete number of points whose axial dimensions vary according to the assumed position of the flaw. Thus, the same number of points describes the crevice pressure profile regardless if the flaw is assumed at the bottom of the tubesheet or at some other location within the tubesheet. Only the slope of the distribution between the points changes. Because of a mismatch between the crevice pressure axial definition and the tubesheet contact pressure axial definition, the integration model cannot converge to a single value at certain discrete points, depending on the model of SG under consideration. For the Model F SG, this point occurs at approximately 4 inches from the top of the tubesheet. The axial range within which this occurs is extremely narrow, less than [ . ]a,ce inch (see Figure RAI9-2), and the non-convergence results in a very limited range of the axial crevice pressure correction factor, less than [ ]."
inch. For the Model F SG, a variation of initial H* prediction of approximately 0.005 inch from the critical axial length results in the model converging again at the lower value of Pcrev correction as also shown on Figure RAI9-2. This result applies generically to the Model D5, 44F and 51 F SGs as well.
For practical application in determining the final value of H*, it is noted that when the adjustments for BET and NOP thermal distribution are included, the predicted values of H* are far removed from the points in the Pcrev correction curves where the model does not converge for all models of SGs. The recommended values of H*, prior to the correction for Pcrev, for the different models of SG are:
Model F: 9.81 inches (Ref: WCAP-17071-P)
Model D5 12.11 inches (Ref: WCAP-17072-P)
Model 44F 11.06 inches (Ref: WCAP-17091-P)
Model 51F 11.14 inches (Ref: WCAP-17092-P)
In all cases, the point of non-convergence of the model does not affect the final recommended value of H*.
27
LTR-SGMP-09-100 NP-Attachment Figure RAI9-1 Puve Correction Profiles for Models F, D5, 44F and 51F SGs
-- ac,e 28
LTR-SGMP-09-100 NP-Attachment Figure RAI9-2 a,c,e 29
LTR-SGMP-09-100 NP-Attachment Vogtle 10. Reference 1, Page 8-6, Section 8.1.4: Clarify whether the RAI "biased"H* distributionsfor each of the four input variables are sampled from both sides of the mean H* value during the Monte Carlo process, or only on the side of the mean H*
value yielding an increased value of H*.
WCGS 11. Reference 1, Page 8-6, Section 8.1.4: Please clarify whether the "biased"H* distributionsfor each of the four input variables are sampled from both sides of the mean H* value during the Monte Carloprocess, or only on the side of the mean H* value yielding an increasedvalue of H*.
B/B 11. Reference 1, Page 8-6, Section 8.1.4: Clarify whether the "biased"H* distributionsfor each of the four input variables are sampled from both sides of the mean H* value during the Monte Carlo process, or only on the side of the mean H*
value yielding an increased value of H*.
CPSES 11. Reference 1,Page 8-6, Section 8.1.4: Please clarify whether the "biased"H* distributionsfor each of the four input variables are sampled from both sides of the mean H* value during the Monte Carloprocess, or only on the side of the mean H* value yielding an increasedvalue of H*.
Seabrook 11. Reference 1, Page 8-6, Section 8.1.4: Clarify whether the "biased"H* distributionsfor each of the four input variables are sampled from both sides of the mean H* value during the Monte Carlo process, or only on the side of the m ean H*
value yielding an increased value of H*
Response
As shown in Figure 8-11 of the report (WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP-17092-P), the variation of the parameters that resulted in the greatest increase in the value of H* were chosen as the "biased" influence factors from which to sample in the Monte Carlo (MC) process. These distributions were normal distributions determined from the mean H* and greatest H* variation resulting from equal valued positive and negative variations of the respective parameters. Note that for the case of coefficient of thermal expansion of the tube, a decrease in the coefficient results in an increase in the H* value and also reflects the broadest distribution. For the coefficient of thermal expansion of the tubesheet, an increase in the coefficient results in increasing H* and also results in the broadest distribution.
Both sides of the biased influence factors were sampled during the Monte Carlo analysis. Sampling from the broadest distributions results in the broadest H* distribution and the largest values of H* corresponding to the desired probabilistic goal, in this case, 95/50.
Figure RAI10-1 shows the results of the Monte Carlo sampling from the interaction surface (see RAI#20) for the resulting values of H* between the upper 93% and 98% of the simulations. (The 98% upper limit was chosen for convenience). The highest values 30
LTR-SGMP-09-100 NP-Attachment of H* are concentrated in a well defined region bounded approximately by the tube coefficient of thermal expansion (COT)between [ ]a,c,e and tubesheet coefficient of thermal expansion (cLTS) between [ ]a,c,e The conclusion that the maximum values of H* are produced from samples in approximately the center of the interaction surface defined by Figure 8-5 in the report applies to both the Model F and Model D5 SGs. Consequently, the use of the broadest distributions that increase the value of H* will tend to focus on the region in question because the broadest H*
distributions are defined by negative variations of orT and by positive variations of cCTS.
Selections from the negative sides of the broadest distributions will not result in maximum values of H*. If picks are made from both distributions on the negative side of the biased influence distributions, the result will be an over-prediction of the lower tail of the H* distribution. This is noted in WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP17092-P and is of no consequence because only the maximum value of H* is of concern.
31
LTR-SGMP-09-100 NP-Attachment Figure RAIIO-1 a,c,e 32
LTR-SGMP-09-100 NP-Attachment Figure RAIIO-2 a,c,e 33
LTR-SGMP-09-100 NP-Attachment Vogtle 11. Reference 1, page 8-14, Figure 8-6: The legend for one of the RAI interactionsshown between a-Ts and ETs appearsto have a typo in it, please review and verify that all values shown in the legend are correct.
WCGS 12. Reference 1, page 8-14, Figure 8-6: The legend for one of the interactionsshown between the coefficient of thermal expansion of the tube (aTs) and tubesheet (ETs) appears to contain typographicalerror. Please review and verify that all values shown in the legend are correct.
B/B 12. Reference 1, Page 8-14, Figure 8-6: The legend for one of the interactionsshown between a-Ts and ETs appears to have a typo in it. Please review and verify that all values shown in the legend are correct.
CPSES 12. Reference 1, page 8-14, Figure 8-6: The legend for one of the interactionsshown between the coefficient of thermal expansion of the tubesheet (aTS) and Young's modulus of the tubesheet (ETS) appearsto contain a typographicalerror. Please review and verify that all values shown in the legend are correct.
Seabrook 12. Reference 1, Page 8-14, Figure 8-6: The legend for one of the interactionsshown between aTs and ETs appears to have a typo in it, please review and verify that all values shown in the legend are correct.
Response
The uppermost curve, defined by the star point, which is currently labeled as aTS=- 3 should be labeled as ats=+3. All other values in the legend are correct.
34
LTR-SGMP-09-100 NP-Attachment Vogtle 12. Reference 1, page 8-20, Case S-4: Why does the RAI assumption of a 2-sigma value for the coefficient of thermal expansion of the tube (aT) and the tubesheet (atTS) to determine a "very conservative biased mean value of H*"
conservatively bound the interactioneffects between aT and aTS? Describe the specifics of how the "very conservative biased mean value of H*," as shown in Table 8-4, was determined.
WCGS 13. Reference 1, page 8-20, Case S-4: Why does the assumption of a 2-sigma value for the coefficient of thermal expansion of the tube (aT) and the tubesheet (aTS) to determine a "very conservative biased mean value of H*"
conservatively bound the interactioneffects between aT and aTS? Pleasedescribe the specifics of how the "very conservative biased mean value of H*" as shown in Table 8-4, was determined.
B/B 13. Reference 1, Page 8-20, Case S-4: Why does the assumption of a 2-sigma value for the coefficient of thermal expansion of the tube (aT) and the tubesheet (aTs) to determine a "very conservative biased mean value of H*"
conservatively bound the interactioneffects between aT and aTs? Describe the specifics of how the "very conservative biased mean value of H*" as shown in Table 8-4, was determined.
CPSES 13. Reference 1, page 8-20, Case S-4: Why does the assumption of a 2-sigma value for the coefficient of thermal expansion of the tube (aT) and tubesheet (aTS) to determine a "very conservative biased mean value of H*" conservatively bound the interaction effects between aT and aTS? Please describe how the "very conservative biased mean value of H*__"as shown in Table 8-4, was determined.
Seabrook 13. Reference 1, Page 8-20, Case S-4: Why does the assumption of a 2-sigma value for the coefficient of thermal expansion of the tube (aT) and the tubesheet (aTS) to determine a "very conservative biased mean value of H*"
conservatively bound the interactioneffects between aT and aTS? Describe the specifics of how the "very conservative biased mean value of H*, as shown in Table 8-4, was determined.
Response
The very conservative mean value of H*, [ ]a.ce inches (Model F), [ ]a.ce inches, (Model D5), is determined by arbitrarily assuming that the 2-sigma values of all variables defines the mean value of H*. To determine these values, it was assumed that the input variables to the structural evaluation (i.e, the entire H* calculation process as shown in 35
LTR-SGMP-09-100 NP-Attachment Figure 1-1 of the report) were set at their 2-sigma values, and the resulting H* was termed the "conservative mean." Table RAI12-1 illustrates the input values that define the mean value of H* and the "very conservative mean" value of H*. The SRSS approach was then applied using the influence factors from Table 8-2 in the report for the 95/50 whole-bundle value appropriate to the model SG being considered. The result is essentially equivalent to the 5-sigma variation case, Case S4 on Table 8-3 of the report. Note that because the 2-sigma input parameter value of H* was determined by the entire calculation process shown in Figure 1-1 of WCAP, the interaction effects of the variables at the 2-sigma level are included in this calculation.
Table RAI12-1 Definition of "Conservative Mean" H*
Definition Analysis Input Parameters and their Values
_aT atTS ET ETS Mean H* mean mean mean mean 1
Conservative Mean-2a( 1 ) Mean+2a( ) Mean-2cy(k) Mean-2cy1 )
Mean H* I (1) Values chosen in direction of increasing H*
36
LTR-SGMP-09-1 00 NP-Attachment Vogtle 13. Reference 1, page 8-22, Case M-5: The description for this RAI case seems to'correspondto a single tube H* estimate rather than a whole bundle H* estimate. How is the analysis performed for a whole bundle H* estimate?
WCGS 14. Reference 1, page 8-22, Case M-5: The description for this case seems to correspondto a single tube H* estimate rather than a whole bundle H* estimate. Please explain how is the analysis performed for a whole bundle H* estimate?
B/B 14. Reference 1, Page 8-22, Case M-5: The descriptionfor this case seems to correspond to a single tube H* estimate rather than a whole bundle H* estimate. How is the analysis performed for a whole bundle H* estimate?
CPSES 14. Reference 1, page 8-22, Case M-5: The description for this case seems to correspondto a single tube H* estimate rather than a whole bundle H* estimate. Please explain how the analysis is performed for a whole bundle H* estimate?
Seabrook 14. Reference 1, Page 8-22, Case M-5: The description for this case seems to correspondto a single tube H* estimate rather than a whole bundle H* estimate. How is the analysis performed a whole bundle H
- estimate?
Response
Case M-5 is the Monte Carlo (MC) sampling analogy to Case S-2. A single tube analysis would sample from the la influence distributions to determine the overall distribution of H*, and from the resulting H* distribution, choose the 95% probability value of the upper tail. Case M-5 pre-biases the influence factor distributions by choosing the influence factor distributions at the 4.285a (Model F) (4.237a Model D5) values divided by 4.285 (Model F) (4.237, Model D5). Thus, the input distributions are pseudo-ludistributions that are already biased by the number of standard deviations required to represent a whole bundle analysis as was done in Case S-2. The use of the greater value influence functions results in a broader final H* distribution from which the 95/50 value represents the whole bundle. The basis for the 4.285.y (Model F.)(4.237a, Model D5) value to represent the whole bundle case is discussed in the report.
It was recognized that the assumption of normality of the influence factor distribution could influence the results from the MC approach included in the report. Nevertheless, the MC cases were included in the report to provide a basis for evaluating multiple variability cases that could not be considered using the SRSS approach. The response to RAI#20 provides a comprehensive analysis based on the interaction surface of Figure 8-5 and utilization of the Monte Carlo technique.
37
LTR-SGMP-09-100 NP-Attachment Vogtle 14. Reference 1, page 8-2: Case M-5 states, "Interactioneffects RAI are included because the 4.285 sigma variations were used that already include the effective interactionsamong the variables." Case M-5 also states that the 4.285 sigma variationscome from Table 8-2. However, Table 8-2 does not appear to include interactions among the variables. Explain how the 4.285 sigma variations include the effect of interactionsamong the variables.
WCGS 15. Reference 1, page 8-22: Case M-5 states, "Interactioneffects are included because the 4.285 sigma variations were used that alreadyinclude the effective interactionsamong the variables."Case M-5 also states that the 4.285 sigma variationscome from Table 8-2; however, Table 8-2 does not appearto include interactions among the variables. Please explain how the 4.285 sigma variationsinclude the effect of interactionsamong the variables.
B/B 15. Reference 1, Page 8-22: Case M-5 states, "Interactioneffects are included because the 4.237 sigma variationswere used that alreadyinclude the effective interactionsamong the variables." Case M-5 also states that the 4.237 sigma variations come from Table 8-2. However, Table 8-2 does not appearto include interactions among the variables. Explain how the 4.237 sigma variations include the effect of interactionsamong the variables.
CPSES 15. Reference 1, page 8-22: Case M-5 states, "Interactioneffects are included because the 4.237 sigma variations were used that already include the effective interactionsamong the variables." Case M-5 also states that the 4.237 sigma variations come from Table 8-2. however, Table 8-2 does not appearto include interactionsamong the variables. Please explain how the 4.237 sigma variationsinclude the effect of interactionsamong the variables.
Seabrook 15. Reference 1, Page 8-22: Case M-5 state,: "Interactioneffects are included because the 4.285 sigma variations were used that alreadyinclude the effective interactionsamong the variables." Case M-5 also states that the 4.285 sigma variations come from Table 8-2. However, Table 8-2 does not appearto include interactionsamong the variables.
Explain how the 4.285 sigma variationsinclude the effect of interactionsamong the variables
Response
Because the 4.285a (Model F), 4.237a (Model D5) variations were calculated using the complete calculation process depicted in Figure 1-1 of the report (WCAP-17071-P, WCAP-1 7072-P), the variations include the structural interaction effects for each variable assuming that all other variables are at their mean value. If multiple variables were 38
LTR-SGMP-09-100 NP-Attachment perturbed simultaneously, a greater effect on H* would be expected. The Monte Carlo sampling scheme used did not support the use of compound parameter variations.
The response to RAI#20 provides an in-depth analysis of the interaction effects among the significant variables using the Monte Carlo method and sampling from the interaction surface of Figure 8-5.
Vogtle 15. Reference 1, page 8-22, Case M-6, first bullet: Should the RAI words "divided by 4.285" appear at the end of the sentence?
WCGS 16. Reference 1, page 8-22, Case M-6, first bullet: Please verify if the words "dividedby 4.285" appearat the end of the sentence.
B/B 16. Reference 1, Page 8-22, Case M-6, first bullet: Should the words "divided by 4.237" appear at the end of the sentence?
CPSES 16. Reference 1, page 8-22, Case M -6, first bullet: Please verify if the words "dividedby 4.237"should appearat the end of the sentence?
Seabrook 16. Reference 1, Page 8-22, Case M-6, first bullet: Should the words "divided by 4.285" appearat the end of the sentence?
Response
The first bullet under Case M-6 on page 8-22 is clarified by adding the phrase "divided by 4.285" (Model F) ("4.237"-Model D5) between "4.285a" (Model F) ("4.237a"-Model D5) and "from."
39
LTR-SGMP-09-100 NP-Attachment Vogtle 16. Reference 1, page 8-23, Case M-7: Was the "2 sigma RAI variation of all variables"divided by a factor of 2?
WCGS 17. Reference 1, page 8-23, Case M-7: Please verify if the "2 sigma variationof all variables"divided by a factor of 2.
B/B 17. Reference 1, Page 8-23, Case M-7: Was the "2 sigma variation of all variables"divided by a factor of 2?
CPSES 17. Reference 1, page 8-23, Case M-7: Please verify if the "2 sigma variationof all variables"was divided by a factor of 2.
Seabrook 17. Reference 1, Page 8-23, Case M-7: Was the "2 sigma variation of all variables"divided by a factor of 2?
Response
For case M-7, the 2-sigma variation was treated as if it were 1-sigma variation. This assumption is somewhat arbitrary and intended only as a hypothetical case to show the effect on H* if it were assumed that the calculated standard deviation are much larger.
Therefore, the 2-sigma variation was NOT divided by 2.
This case is an arbitrary sensitivity study that addresses the H* result if the 1 influence factors were more than doubled. Starting from the basic mean structural prediction of H*, [ ]a,,,e for the Model F ([ ]a.c.e, for the Model D5) inches, it was assumed that the 2a influence distributions applied instead of the 1c influence distributions, and the MC sampling was from the 2cy distributions. The principal objective of this case was to show that very conservative assumptions do not lead to a major impact on the value of H*.
40
LTR-SGMP-09-100 NP-Attachment Vogtle 17. Reference 1, page 8-23, Case M-7: Explain how this case RAI includes the interactioneffects between the two principle variables, acT and aTS.
WCGS 18. Reference 1, page 8-23, Case M-7: Please explain how thi.
case includes the interaction effects between the two principle variables, aT and aTS.
B/B 18. Reference 1, Page 8-23, Case M-7: Explain how this case includes the interactioneffects between the two principle variables, aT and aCTS.
CPSES 18. Reference 1, page 8-23, Case M-7: Please explain how thi.s case includes the interaction effects between the two principal variables, aT and aTS.
Seabrook 18. Reference 1, Page 8-23, Case M-7: Explain how this case includes the interaction effects between the two principle variables, aT and aTS.
Response
Case M-7 assumes that the 1oc variability of H* in the parameters is based on the 2a influence factors calculated for each parameters. Because the influence factors are calculated using the entire calculation flow depicted in Figure 1-1 of the report, the interactive effect of the key parameters at the 2a is reflected. The calculations were performed by perturbing one parameter at a time; therefore, the combined interaction of perturbing multiple parameters is not reflected. However, the assumption that the 2a variation in the direction of increasing H* represent one standard deviation of the H*
influence factors and the extreme value calculation process provide a very conservative estimate of H*.
The response to RAI#20 provides an in-depth analysis of the interaction effects among the significant variables.
41
LTR-SGMP-09-100 NP-Attachment Vogtle 18. Reference 1, page 8-25, Table 8-4: Explain why the mean H*
RAI calculated in the fifth case does not require the same adjustments, as noted by the footnotes, that all other cases in the table require.
WCGS 19. Reference 1, page 8-25, Table 8-4: Please explain why the mean H* calculated in the fifth case does not require the same adjustments, as noted by the footnotes, th at all other cases in the table require.
BPB 19. Reference 1, Page 8-25, Table 8-4: Explain why the mean H*
calculatedin the fifth case does not require the same adjustments, as noted by the footnotes, that all other cases in the table require.
CPSES 19. Reference 1, page 8-25. Table 8-4: Please explain why the mean H* calculated in the fifth case does not require the same adjustments, as noted by the footnotes, that all other cases in the table require.
Seabrook 19. Reference 1, Page 8-25, Table 8-4: Explain why the mean H* calculatedin the fifth case does not require the same adjustments, as noted by the footnotes, that all other cases in the table require.
Response
The superscripts referring to the notes were inadvertently omitted from the mean H*
value for Case S-4 in Table 8-4. The mean value of H* shown, [ ]ac,e inches (Model F) ([ ]ace inches, Model D5), includes the adjustment for BET and for the NOP thermal distribution.
This omission exists also in WCAP-17072-P for the Model D5 SGs, but has been corrected on subsequent H* reports for the Model 44F and 51F SGs (WCAP-17091-P and WCAP17092-P).
42
LTR-SGMP-09-100 NP-Attachment Vogtle 19. Reference 1, page 8-25, Table 8-4: Verify the mean H*
RAI shown in the last case in the table.
WCGS 20. Reference 1, page 8-25, Table 8-4: Please verify the mean H* shown in the last case in the table.
B/B 20. Reference 1, Page 8-25, Table 8-4: Verify the mean H*
shown in the last case in the table.
CPSES 20. Reference 1, page 8-25, Table 8-4: Please verify the mean H* shown in the last case in the table.
Seabrook 20. Reference 1, Page 8-25, Table 8-4: Verify the mean H*
shown in the last case in the table
Response
(Please also see the response to Question 16.)
The mean value of H*for Case M7 is correct as shown on Table 8-4.
The purpose of this case was to determine the whole bundle, extreme value of H* and to show the effect on H* if the uncertainties were doubled at the same time as discussed in the response to Question 16. The process to achieve this was to calculate the mean H*
as for all other cases, except case S4, on Table 8-4, but then to arbitrarily assume that the 2 G variations as the values for the 1a influence distributions of H*. The intent of this case was to show that extreme assumptions of variability do not invalidate the H*
concept.
43
LTR-SGMP-09-100 NP-Attachment Vogtle 20. Section 8 of Reference 1: The variabilityof H* with all relevant RAI parametersis shown in Figure 8-3. The interactionbetween caT and aTs are shown in Figure 8-5. Please explain why the direct relationshipsshown in these two figures were not sampled directly in the Monte Carlo analysis, instead of the sampling method that was chosen. Also, please explain why the sampling method chosen led to a m ore conservative analysis than directly sampling the relationshipsin Figures 8-3 and 8-5.
WCGS 21. Section 8 of Reference 1: The variabilityof H* with all relevant parametersis shown in Figure 8-3. The interactionbetween aIT and a'TS are shown in Figure 8-5. Please explain why the direct relationshipsshown in these two figures were not sampled directly in the Monte Carlo analysis, instead of the sampling method that was chosen. Also, please explain why the sampling method chosen led to a m ore conservative analysis than directly sampling the relationshipsin Figures 8-3 and 8-5.
B/B 21. Section 8 of Reference 1: The variabilityof H* with all relevant parametersis shown in Figure 8-3. The interactionbetween aT and aZTS are shown in Figure 8-5. Please explain why the direct relationshipsshown in these two figures were not sampled directly in the Monte Carlo analysis, instead of the sampling method that was chosen. Also, please explain why the sampling method chosen led to a m ore conservative analysis than directly sampling the relationshipsin Figures 8-3 and 8-5.
CPSES 21. Section 8 of Reference 1: The variabilityof H* with all relevant parametersis shown in Figure 8-3. The interactionbetween aT and aTS are shown in Figure 8-5. Please explain why the direct relationshipsshown in these two figures were not sampled directly in the Monte Carlo analysis, instead of the sampling method that was chosen. Also, please explain why the sampling method chosen led to a m ore conservative analysis than directly samplingthe relationshipsin Figures 8-3 and 8-5.
Seabrook 21. Section 8 of Reference 1: The variability of H* with all relevant parametersis shown in Figure 8-3. The interactionbetween aT and aTS are shown in Figure 8-5. Please explain why the direct relationshipsshown in these two figures were not sampled directly in the Monte Carlo analysis, instead of the sampling method that was chosen. Also, please explain why the sampling method chosen led to a m ore conservative analysis than directly sampling the relationshipsin Figures 8-3 and 8-5.
44
LTR-SGMP-09-100 NP-Attachment
Response
General The recommended value of H* is based on the square root of the sum of the squares (SRSS) approach to combining the uncertainties for H*. The Monte Carlo cases included in the report were included as a vehicle to study different sensitivities to H*
parameters variations and were provided as support for the SRSS recommendation.
The peer review (Expert Panel's) conclusions were that the SRSS approach was a suitably conservative approach given the many conservatisms built into the H* analysis.
The significant conservatisms included in the H* analysis are summarized in Section 1 of the report(s) and again identified in Section 10 of the report(s).
Figures 8-3 and 8-5 were developed during, and immediately after, the peer review of the H* project, which was followed in close order by publishing the report. The staffs
,observation that Figures 8-3 and 8-5 reasonably define an interaction surface, which could be utilized directly for a Monte Carlo sampling assessment, is correct. Therefore, a Monte Carlo analysis based on the interaction surface defined by Figure 8-4 in the respective WCAP reports for the different models of SGs was completed. This analysis provided the opportunity to quantify some of the conservatisms that are included in the technical justification of H*. The approach to this issue was to consider the most significant conservatisms in the overall H* analysis and quantify their effects on the recommended value of H* to show that the recommended value of H* is conservative.
The sequence of the analysis was as follows:
- 1. Application of the Monte Carlo methodology discussed in the H* reports, except for case M-6, assumes that each simulation of H* includes a different value of the properties of the tubesheet. Thus, if 100,000 simulations are performed, each simulation includes a different random pick of tubesheet properties. Among the population of H* candidate plants, there are 60 steam generators; therefore, the actual population of tubesheets is limited to 60. To better address the limited population of tubesheets, the reference MC sampling is a staged process corresponding to the simulation of one steam generator tubesheet/tube bundle combination. A set of tubesheet properties is selected, and for that set, the corresponding tube properties are sampled 5626 times for the Model F SG tube population (4570 times for the model D5 SG tube population), and as appropriate for the other models of SGs. The above process is repeated 10,000 times. This provides a more accurate simulation reflecting the limited number of tubesheets in the population.
- 2. The probabilistic analysis in Section 8 of the report(s) assumes that the entire tube bundle consists of tubes located at the worst case location in the tube bundle (e.g.,
the most limiting radius in the most limiting sector of the tube bundle as shown in Section 6.2.3). As shown in Figure 6-1 of the report, the worst tube is defined by a very narrow segment of tubes, while all other tubes are shown to have a lower value 45
LTR-SGMP-09-100 NP-Attachment of H*. Therefore, the bundle was divided into a number of sectors as discussed below, and the 0.95 probability at 50% confidence value of H* was defined on the combined probability of the sector probability for all tubes. This analysis is still quite conservative because all tubes are still assumed to be in the limiting azimuthal sector of the tube bundle (the sector perpendicular to the divider plate including about 50 from the centerline of the tubesheet. See Section 6.2.3 of the report). Tubes more than about 5 pitches removed from the centerline perpendicular to the divider plate have been shown to have lower values of H*.
- 3. The current analysis for the Model F and Model D5 SGs includes a correction factor for the NOP thermal distribution through the tubesheet. The factor was developed very conservatively using a ratio technique (see report Section 6.2.2.2.5). This correction factor was re-evaluated (see the response to RAI#2) and a more realistic value of it is applied in this analysis. The correction factor does not apply equally to all models of SG; therefore, the effect of the correction factor is SG model-specific and has already been included in the reports for the Model 44F and Model 51F SGs (WCAP-17091-P and WCAP-17092-P). The analysis results below identify the SG models where this improved analysis applies.
- 4. The H* analysis assumes no contribution from residual contact pressure (RCP). All test data to date, including data from tests performed prior to 2008, has shown that a positive value of RCP exists after hydraulic expansion of the tubes. Tests were performed during the current H* program that confirmed a significant level of RCP, and also showed that within a short distance of motion, the forces required to continue to move the tube by far exceeded the maximum pull out forces that could be generated under very conservative assumptions. The analysis quantifies the effect of RCP on the calculated value of H* and benchmarks the RCP to the tests that were performed during the H* development.
A. Sector Analysis Based on the profile of the predicted mean H*, the tube bundle is divided into 9 annular sectors as shown in Figures RAI20-1, -2, -3 and -4 for the models of SG included in the H* population (Reference RAI20-1). (In Figure RAI20-4, for the Model 51F SGs, the appropriate sector division results in only 7 sectors; however, additional sectors with 1 tube, each, were added at both ends for convenience of the calculations.) The normalized H* is determined from the raw H* calculation results, prior to adjustment of the H* value by the addition of correction factors for the BET and NOP tubesheet thermal profile. This is done to obtain a true normalization, unaffected by any constants.
However, the final value of H*, after the MC sampling for AH, is based on the adjusted maximum mean value of H* as shown in the appropriate sector in Tables 1, 2, 3 and 4.
The adjustment for crevice pressure referenced to the predicted H* is made after the all other factors have been accounted for. Thus, for each sector:
46
LTR-SGMP-09-100 NP-Attachment E ] a,c,e
- Where, F is the sector normalization factor from Figures RAI20-1, -2, -3 and -4, H*(BET+ Tnop) is the raw H* value adjusted for BET and NOP thermal distribution AH*uncert is the adjustment for interaction effects from the MC analysis AHpcrev is the adjustment for crevice pressure The normalized value of H* in each sector is based on the maximum value of H* in that sector; thus, the sector evaluation is inherently conservative.
The number of tubes in each sector is determined from the row and column numbers and the model-specific pitch of the tubes. Tables RAI20-1, -2, -3, and -4 summarize the sector populations for each of the models of SGs.
B. Interaction Surface and Monte Carlo Sampling A simulation model was developed to evaluate limiting values of H* for specific classes of steam generators. The Monte-Carlo based model evaluates extreme values of H* on a single steam generator basis, repeating the process to construct a distribution of maximum H* values. The final output of the model is the 95/50 estimate of extreme H*
within any one steam generator.
The components of variance included in the model are the coefficients of thermal expansion (CTE's) for the tubesheet and the individual tubes. These have been shown in the H* reports to be, by far, the most significant contributors to variations in H* for the tubesheet/tube bundle combinations. The essential function describing H* variation for specific value pairs of the thermal expansion coefficients has been developed and is shown in Figure 8.5 of the H* reports. It should be noted that full interaction effects are included.
The basic structure of the simulation is shown in Figure RAI20-5 and represents one Monte Carlo trial. The process shown produces one realization of the extreme H* for a given steam generator. Repetition, involving 10,000 trials produces a distribution from which a 95/50 estimate of H*can be obtained by robust nonparametric means. As shown in Figure RAI20-5, the core process involves a random selection of one value of tubesheet CTE and N values of tube CTE, where N is the number of tubes in the steam generator or specific region of interest. The resulting N pairs are propagated through the 47
LTR-SGMP-09-1 00 NP-Attachment fitted surface to produce N values of H* which are then sorted to identify the maximum (extreme) value of H* which is stored for further use.
The above process can be easily applied on a regional (SG sector) basis by running the simulation for each region separately based on region-specific values on tube population size and average H*. The composite H* for the entire steam generator can be obtained by the following equation:
H*c = [ 1 ac,e for M Region model It is most important that the H* values for the individual regions are not sorted prior to application of the above post-processing because of the need to maintain tubesheet identity between regions.
C. Sector Application of Interaction Effects The interaction data shown in Figure 8-5 of the H* WCAPs were developed for the limiting tube radius (i.e., the tubesheet radius in Figures 1, 2, 3 and 4) where the normalized value of H* is 1. Because of the complex nature of the H* analysis, it was necessary to determine if the interaction effects at the limiting H* radius adequately represented the interactions at other tubesheet radii. Two radii were selected to represent the most probable locations where significant effects, if they exist, might materialize: 1) A tubesheet sector near the limiting radius, and 2) A tubesheet sector far removed from the limiting radius.
It was shown in the reports that the influence of Young's modulus on the final values of H* is negligible and that there was no significant interaction between the Young's moduli of the materials and the coefficient of thermal expansion of the materials. The existing interactions are limited to the coefficients of thermal expansion of the tube and tubesheet materials. Therefore, the same matrix of sensitivity cases that defined Figure 8-5 in the reports was run for each of the two tubesheet sector chosen as noted above.
In all cases it was determined that the interaction effects defined in Figure 8-5 of the report(s) for the location of the maximum mean H* value bounded the interaction effects of the other sectors considered. Therefore, for conservatism and simplicity, the range of interaction effects (i.e., AH* = f(aT, CLTS)) for the maximum mean value of H* shown on Figure 8-5 was applied for all sectors of the tubesheet.
Figures RAI20-6 and RAI20-7 show the results of this evaluation for the tubesheet sectors selected for the Model F and Model D5 SGs. The interaction profile for the mean, 3a and 5a variation of tubesheet coefficient of thermal expansion are shown to cover the significant range of variability. In all cases, the variability of the location of the maximum value of H* is greater than, or equal to, the variability at other radial locations on the tubesheet. Therefore, the application of the variability for the radial location of the maximum value of H* for all other radial locations is justified and conservative.
48
LTR-SGMP-09-100 NP-Attachment D. Results from Sector Based Samplingq from the Interaction Surfaces Table RAI20-5 (a) summarizes the recommended values of H* from the H* reports for all of the affected Model SGs together with the results of the Monte Carlo (MC) sampling from the interaction surface defined in Figure 8-5 of each report. The MC sampling was based on the sector approach described above and also the approach shown in Figure RAI20-5 to limit the number of tubesheet simulation. The result from this sampling must be adjusted for the crevice pressure distribution referenced to the location of the initially predicted value of H*. The correction for crevice pressure is taken from Figure 8-1 of the respective reports. After the adjustments are made for the crevice pressure reference location, the values of H* are slightly greater than the recommended values of H* from the respective reports.
Table RAI20-5(b) extends the evaluation of the conservatism of the recommended SRSS-based values of H* by adjusting the Monte Carlo sampling results for the updated values of the adder for the NOP thermal distribution in the tubesheet for the Model F and Model D5 SGs. The updated NOP thermal distribution factor for the Model 44F and Model 51F SGs are already included in the respective reports (WCAP-17091-P and WCAP-1 7092-P); consequently there is no adjustment made for these models of SG.
The original NOP thermal distribution adjustment factor was developed on a very conservative basis, using the scaling method described in Section 6.2.2.2.5 of WCAP-17071-P and WCAP-17072-P. As the analysis for H* evolved, a direct method of applying the tubesheet NOP thermal distribution in the structural analysis was developed; this method is describe in Section 6.2.2.2.5 of WCAP-17072-P (Model D5 report). For the Model D5 SG, the necessary correction based on the updated method was [ ]a~c,e inch compared to [ ,2,c'e inch based on the scaling technique. A similar analysis was subsequently performed for the Model F SG and it was determined that the appropriate correction for the NOP thermal distribution is [ ]a*c'e inch instead of the
[]jace inches included in the recommended value of H* in WCAP-17071-P.
When the updated correction for the NOP thermal distributions are applied, and the necessary correction for crevice pressure reference location is applied, the final value of H* for the Model F SG is [ ]a,ce inches and, for the Model D5, is [ ]a,ce inches (see Table RAI20-5(b)). Both of these values are less than the recommended values of H*,
respectively, for the Model F and Model D5 SGs. Thus, it is concluded that the recommended values of H*, based on the SRSS approach as shown in the respective reports for the Model F and Model D5 SGs, are conservative.
It should be noted that the adjustment of the NOP thermal distribution correction factor does not impact which operating condition, NOP or SLB, is the limiting condition. The limiting value of H* is determined by three times normal operating pressure before and after the adjustment for the NOP thermal distribution. Section 6.4.5 of the Model F report, WCAP-17071-P, and the Model D5 report, WCAP-17072-P, discusses the determination of the H* values. When the NOP thermal distribution is directly included in 49
LTR-SGMP-09-100 NP-Attachment the structural analysis to determine tubesheet displacements, the NOP condition remains the limiting condition for H*.
E. Determination of Residual Contact Loads from Pull Out Tests In prior analyses for H*, pull out test data has been used to calculate a residual contact pressure, which is distributed over the length of the tubesheet and included in the
,integration of pull out force over length to determine the length at which the pull out and resisting forces are equal. However, the pull out resistance can also be used to offset the pull out forces. Both methods were studied and it was determined that the same result was achieved, regardless of which method was applied. Because offsetting the applied loads requires fewer assumption (i.e., coefficient of friction) and results in more conservative values of H*, this approach was selected to determine the effect of the hydraulic expansion only on the calculated value of H*.
Reference RAI20-2, provided as Appendix A to this document, summarizes the pull out test program performed in support of the H* development. The data from the pull out tests and Monte Carlo simulation were used to determine a conservative value of end cap load reduction. As in prior pull out tests, there was considerable scatter in the pull test data. The highest pull force recorded at 0.25 inch cross head displacement was S]a,c,e Ibf, and the lowest pull force recorded at 0.25 inch cross head displacement was [ ]a,c,e Ibf. Monte Carlo simulation was then used to determine a 5/50 value (i.e.,
the lower 95% bound) of the pull test data.
The Monte Carlo simulations used the pull test data to establish means and standard deviations for the pull forces that were observed. Two sets of data for each of three tube diameters (0.688 inch, 0.750 inch, and 0.875 inch) were provided: One considered the 13 in. expansion lengths only and the other considered all expansion lengths (13, 15 and 17 inches) combined. Seven distributions were used: 1) A truncated (at 0) normal distribution, 2) a lognormal distribution, 3) an Erlang distribution, 4) a Gamma distribution, 5) an inverse Gaussian distribution, 6) a Pearson Type V distribution, and 7) a Weibull distribution. All except the truncated normal were chosen because their domains range from 0 to + infinity, their domains are continuous, and their fitting parameters for the means and standard deviations used were within their allowable values. One hundred thousand iterations were run for each simulation, and the 5/50 values of pull force recorded for each distribution. The most conservative result, Ibf]ace, came from the simulation that used the Weibull distribution, and this number is yery consistent with the lowest observed pull test datum. Note that the Weibull distribution is widely recommended to model distributions in lieu of a truncated normal distribution. The figure below illustrates the results of the Monte Carlo sampling based on the Weibull distribution of the test data. Complete details of the above analysis can be found in Reference RAI20-2.
50
LTR-SGMP-09-100 NP-Attachment The recommended end car load reduction is I r.ace Ibf.
a,c,e (Figure corresponding to the Monte Carlo simulation using a Weibull distribution for the Model F SG data, using the 13 inch expansion length only.
The 5/50 value of pull force is [ ]a.ce lbf.)
F. Application of Residual Contact Load The H* results in Figure 8-5 of WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP-1 7092-P show that H* is sensitive to the variations in the coefficient of thermal expansion (CTE) of the tube (caT) and the tubesheet (caTs). The reports also show that H*
is not significantly sensitive to variations in the Young's modulus (E) of the tube or the tubesheet. The results in Figure 8-5 in WCAP-17071-P, WCAP-17072-P, WCAP-1 7091-P and WCAP-1 7092-P also demonstrate that the worst case trend in the variation of the thermal expansion coefficients is when the aT is decreasing and CLTS is increasing. In other words, H* increases the most when the coefficients of thermal expansion are varied to reduce the contact pressure between the tube and the tubesheet due to thermal growth.
It is possible to reduce the order of the problem (i.e., reduce the number of dimensions involved in the sensitivity study) given the knowledge of which values and directions of variation in CTE are most important to the problem. Figures RAI10-1 and RAI10-2 show the combinations of aT and ctTS that are most likely to produce a worst case H*value. The values of CTE standard deviations for both the tube and tubesheet are combined into an effective variable using the following relationship:
51
LTR-SGMP-09-100 NP-Attachment asrss = A__0T7 _+(UTS )Y The possible variation in sign of either CTE standard deviations is not included in this equation because the only values of interest occur when the tube CTE variation is negative relative to the mean and the tubesheet CTE variation is positive relative to the mean. This reduced form of variation in CTE is then used to compare the change in H*
due to the application of residual pre-load between the tube and the tubesheet due to the installation and hydraulic expansion of the tube.
There are multiple ways to achieve the same value of Oxsrss. For example, a TS CTE variation of +5a about the mean and a tube CTE variation of -5a about the mean are each equal to a combined a of 5 (assuming only one is non-zero). Likewise, a combination of tube and TS CTE variations of -3/+4 and -4/+3 will also yield an a of 5.
However, the net change in H* with respect to the material properties are very similar for a single value of asss regardless of values of its component parts. In cases where there are multiple possibilities for a unique value of ctsrss, the combination of TS and tube CTE that produced the smallest reduction in H* was used. Figure RAI20-8 shows the multiple curves that were used to create the surface seen in RAI20-9.
Hydraulic expansion of the tube into contact with the tubesheet tube bore introduces a pre-load that must be overcome before the tube can translate within the tubesheet tube bore. This means that in addition to the pull out resistance that a tube develops due to internal pressure, thermal growth, etc., the pull out resistance of the tube due to the hydraulic expansion must also be overcome in order for the tube to freely translate within the tube bore. However, the hydraulic expansion process has only a small effect on the development of contact pressure between the tube and the tubesheet compared to that developed due to operating pressure and temperature. Therefore, the installation effect, termed residual contact load (RCL), is included as a reduction of the applied end cap load. Recall that the end cap loads are based on the mean +2a tubesheet bore diameter and are thus very conservative.
The reduction in end cap load, for the th value of pull out resistance due to installation effects is:
Pi = End Cap Load = nAp gr,' - DL - RCLi
- Where, n is the applicable safety factor for the SG operating condition based on the SIPC, zip is the primary to secondary pressure differential, r, is the outside tube radius, 52
LTR-SGMP-09-100 NP-Attachment DL is the dead load of the tube straight leg above the top of the tubesheet and RCL is the value of installation pre-load determined from test results.
The minimum pull out force from section F above, [ ],o,e lbf, was used. The dead weight of the straight leg portion of the tube above the tubesheet was also included because it also provides a resistance to tube pull out. The dead weight of the straight legs of the tubes varies between [ ]ace and [ ]a,,ce lbs depending on the length of the tube straight leg; an average value of [ ]ace lbs was used.
As an example, for the NOP Low Tavg condition at Millstone Unit 3, the end cap load due to the pressure acting on the tube is [ ],,c lbf. Assuming the minimum value of pull out force from the test data and an average dead weight of the tube straight leg, the applied end cap load that must be balanced by the distribution of contact pressure between the tube and the tubesheet is equal to [ ]c..e lbf - [ a,c,e lbf - [ ]a,ce lbf, or [ ]a,c,e lbf.
Using the RCL to reduce the end cap load on the tube has been shown to be conservative in a direct comparison with the alternative method, that is, converting the pull out force to a residual contact pressure and including it in the integration for H*.
Further, reduction of the applied load does not affect the contact pressure distribution between the tube and the tubesheet. For instance, if there was a combination of material properties and operating conditions that resulted in a very small or zero value of contact pressure for some portion of the tube below the top of the tubesheet, the application of RCL as a reduction of applied load does not change the predicted contact pressure. The first point of positive contact between the tube and the tubesheet is still determined based on the structural analysis of the tubesheet. An additional benefit from applying the RCL as a reduction to the applied end cap load is that there is no need to develop a distribution of the residual effect of the tube installation as a function of elevation in the tubesheet. The test results can be directly used to determine the pre-load on the tube.
A value of H* is determined for any value of RCL for the limiting SG operating condition at the limiting TS radius and sector in the bundle. The process for determining the H*
value is shown in the following flow chart.
53
LTR-SGMP-09-100 NP-Attachment Step I Step 2 Step 3 Dekrmn~e, Det,' '*06e:
Allijlpýd Ejid CAip Load = hficeaste PR siiIilaion.
- Lovad l fU) He ,]
RCP Effect(f)Y BETOff + Ses.epeatet eandLoadof va Tubee e f)Hrenice DoL chageissu Of set ttevuo Determine New H'ui faPinc The result of this process is a surface of the response in H* to changes in RaL andc srss the square root of the sum of squares of the specific variations in CTE from one MC simulation). If the values for RCL are normalized to an assumed value, say [ ]ace lbf, and the values of H* are taken as the change in H* relative to the value of H* with an RCL of zero, the result is a non-dimensional surface that can be used in conjunction with a Monte Carlo analysis to determine the reduction in H* due to the inclusion of ROL.
Figure RA120-9 is a surface plot of the change in the Model F H* values as a function of RCL and U srss. Figure RA120-10 is a surface plot of the change in the Model D5 H*
values as a function of RCL and a,,,s.
Figure RAI20-9 and RAI20-10 illustrate that the effect of including the RCL as a reduction in the applied tube end cap load is dependent on both the H* value and the material parameters. This is a logical result because if H* is small (correlated to a small value of asrss) then the effect of RCL should also be small because there is enough contact pressure to maintain equilibrium with the load on the tube regardless of the value of RCL. However, if H* is large, because of some combination of material parameters or operating conditions that produce less contact pressure between the tube and the tubesheet, then the presence of any value of RCL has a much larger effect on H*. For example, in Figure RAI20-8, assuming an RCL ratio of 1 (RCL -1 ] a,c,e lbf) with an asrss of 0 results in a very small correction to the final H* distance on the order of
[
],,,e inch. However, ifthe RCL ratio is equal to 1 and srss is equal to 5, the change in H* is 2 or more inches, or a factor of 4 greater.
The effects of residual contact pressure (RCP) are implemented in the extreme-value simulation model using a functional representation of the developed steam generator-specific data described above. The function describes the correction term ( AH* ) in terms of two variables:
54
LTR-SGMP-09-100 NP-Attachment AH* = G( RCL, Alpha)
Where:
RCL = Residual contact load Alpha = Effective thermal expansion coefficient A typical description of this surface is shown in Figure RAI20-8. As can be seen from the figure, the behavior of the function is somewhat complex. The value of the function generally increases with both independent variables which makes some simplification possible based on a conservatively low estimate of one of the variables.
A lower limit constant value of RCL was chosen, in part to assure a more robust computational behavior in the implementation of the RCL effects modeling. The value cited in the response to part F of this RAI corresponds to a RCL ratio of approximately 1.0. Figure RAI20-11 shows the resulting AH* as a function only of Alpha. This and corresponding functions for each steam generator class, were implemented in the full simulation model.
The actual implementation into the simulation model was straightforward. Since the RCL correction is subtractive, the computation of Alpha and AH* is performed directly after the computation of H* within the simulation. The computation is performed for all tube/tubesheet combinations in the entire simulation process. The reduction in the computed extreme values of H* is typically on the order of 1-2 inches, and is steam generator-specific.
It is important to note that the change in H* due to the crevice pressure adjustment, thermal offset and BET is already included in the analysis. The distribution of the crevice pressure adjustment shown in Figure 8-1 of the H* reports is not required in this instance.
That is because the reduction of the end cap load changes how the H* value will react to a change in contact pressure distribution. So it is necessary to incorporate the change in H* due to the RCL reduction of the end cap load with the crevice pressure adjustment to produce a net change in H* using consistent methods. Therefore, the result of using the RCL surface to determine the change in H* is the net effect of all adjustments to H* and no further corrections are required.
Table RAI2-05(c) summarizes the effects of the application of residual pull out load (RCL) on the value of H*. When the 5/50 pull out force from the test data is applied using the Monte Carlo approach that samples from Figure 8-5 in the reports and also from the RCL correction surface discussed above, the values of H* are reduced approximately 1 to 2 inches for all affected models of SGs. The resulting values of H* for the Model F and Model D5 SGs are further reduced by application of the updated NOP temperature distribution correction factor. As can be seen from Table RAI20-5, the recommended values of H* for the respective SGs in the applicable reports exceeds the values determined when the conservative factors inherent in the recommended values are considered in the analysis.
55
LTR-SGMP-09-100 NP-Attachment G. Summary and Conclusions The recommended values of H* for the different models of SGs as provided in the respective reports (WCAP-17071-P [Model F], WCAP-17072-P [Model D5], WCAP-17091-P [Model 44F] and WCAP-17092-P [Model 51F]) were based on very conservative assumptions. Additional analysis, using Monte Carlo techniques and the variables interaction surfaces defined in Figure 8.5 of the reports, was performed to quantify the conservatism of these assumption with regard to the recommended values of H* for the different models of SGs. Four principal conservatisms were evaluated:
- 1. The number of tubesheets was limited to a number less than the number of tubes in the bundles of the respective SG models. The total population of SGs among the H*
candidate plants is 60 including 4 different models of SGs. The number of tubesheets simulated for each SG was limited to 10,000.
- 2. Instead of assuming that all tubes in the bundle are located at the single worst case position that defines the recommended value of H*, the bundles were divided into sectors. This approach retains significant conservatism because the maximum value of H* in each sector was used for the analysis and the limiting interaction variances were applied to all sectors. It is noted that all sectors considered are located in the limiting azimuthal sector of the tubesheet as discussed in Section 6.2.3 of the reports.
- 3. The conservative adder for NOP tubesheet thermal distribution was re-evaluated by including the thermal distribution directly in the structural analysis. The resulting adders to H* are realistic values that reflect the actual response of the tubesheet structure to the applied thermal distribution. This applies only for the Model F and Model D5 SGs because the updated thermal correction factor is already included in the recommended H* Values for the Model 44F and 51F SG. Modification of the thermal distribution factors does not change that the NOP conditions are the limiting conditions that determine the value of H*.
- 4. Based on pull out tests performed during the H* development, the effect of the minimum measured pull out forces at 0.25 inch of tube travel on the values of H*
were evaluated. The pull out force data was applied directly as a reduction of the applied loading instead of utilizing an intermediate conversion of pull out force to contact pressure. This approach is more direct, and its specific application is conservative because the 5/50 value of pull out force was used. In reality, much greater values of pull out force were demonstrated in the tests at 0.25 inch travel.
Still greater pull out forces were observed during the tests for greater values of tube travel, even exceeding the limiting applied design loads for H*. Therefore, the application of the 5/50 value of pull out force from the tests is conservative.
After addressing the above factors, the final values of H* are significantly less than the values recommended for all affected models of SGs. Therefore, the recommended values of H* for each of the models of SG are shown to be conservative.
56
LTR-SGMP-09-100 NP-Attachment RAI#20
References:
RAI20-1 LTR-SGMP-09-92;"Tubesheet Sector Definition for H* Revised Probabilistic Analysis," July 10, 2009.
RAI20-2 LTR-SGMP-09-98, "H* Pull Test Program Summary," July 27, 2009.
57
LTR-SGMP-09-100 NP-Attachment Figure RAI20-1 Model F
-- I a,c,e Figure RAI20-2 Model D5 a,c,e 58
LTR-SGMP-09-100 NP-Attachment Figure RAI20-3 Model 44F a,c,e Figure RAI20-4 a,c,e Model 51F 59
LTR-SGMP-09-100 NP-Attachment Figure RAI20-5 Monte Carlo Simulation Process SELECT RýANDOM TUBE
,GTE()
Repeat for number of tubes in sector or bundle as appropriate IN SG OR 60
LTR-SGMP-09-100 NP-Attachment Figure RAI20-6 Model F: Interaction Profiles for Sector-Base Sampling ac,e 61
LTR-SGMP-09-100 NP-Attachment Figure RAI20-7 Model D: Interaction Profiles for Sector-Base Sampling a,c,e 62
LTR-SGMP-09-100 NP-Attachment Figure RAI20-8 AH* for Various Values of Oxsrss and RCL Ratio (a -cxsrss, RCLref = 8001bf) a,c,e 63
LTR-SGMP-09-100 NP-Attachment Figure RAI20-9 Model F Response Surface for the Change in H* as a Function of RCL and (srss a,c,e 64
LTR-SGMP-09-100 NP-Attachment Figure RAI20-10 Model D5 Response Surface for the Change in H* as a Function of RCL and Usrss a,c,e 65
LTR-SGMP-09-100 NP-Attachment Figure RAI20-11 Change in H* as a Function of Usrss a,c,e 66
LTR-SGMP-09-100 NP-Attachment Table RAI20-1 Model F SG Sector Populations Model F TS Radius 0-11 11-17 17-23 23-29 29-35 35-41 41-47 47-53 53-60 Max Mean H* ]a,c,e Max Mean a,c,e H*Factor Number of Tubes Table RAI20-2 Model D5 SG Sector Populations Model D5 TS Radius 0-6 6-12 12-18 18-24 24-30 30- 36-42 42-48 >48
_________ 36
],c,e Max Mean H*
a,c,e Max Mean H*Factor Number of Tubes Table RAI20-3 Model 44F SG Sector Populations Model 44F TS Radius <9 9-15 15-21 21-27 27-33 33-39 39-45 45-51 >51
]a,c,e Max Mean H*
a,c,e Max Mean H*Factor Number of Tubes Table RAI20-4 Model 51F SG Sector Populations TS Radius <9 9-17 17-24 24-32 32-41.85 41.85-52.52 >52.52 Max Mean H* a,c,e a,c,e Max Mean H*Factor Number of Tubes 67
LTR-SGMP-09-100 NP-Attachment Table RAI20-5 Results of Monte Carlo Sam ling and Valuation of Conservatism Surface Sampling from Figure 8-5 of the Report(s) with SG Model Report Case S-2 Limited Number of Tubesheets and Sector Based Approach 95/50 (inch) Pcrev (inch) ",':Final H*ýi(inlch).- 95/50 (inch) Pcrev Final H*
a,c,e D5 13.8+- ____
44F *'3 >~ _____
51F , 3.14 :_._,__,__:
-4a) Sampling from Intera ,on Surface Figure 8-5 H* After Correction for NOP Surface Sampling, Limited Tubesheets Corrected for SG Model Surface Sami Thermal Distribution NOP Thermal Offset Sampling 95/50 H*(inch) Original (inch) Revised 95/50 H*(inch) Pcrev (inch) Final H*(inch)
(inch c'e a,c,e F
D5___________
44F _ _ _=
51F =
Surface Sampling Final H*
oH* for Minimum Final H*After Correction for NOP (inch)
SG from Figure 8-5 of the Report(s) with Limited Out Force (inch) Including Minimum Thermal Distribution Number of Pull Out Force Tubesheets 95/50 (inch) NA (95/50) (inch) (inch) (95/50) (inch ,c,e F -
D5 : ._-- i_- "_____
44F (c) Adjustment for Residual Contact Pressure Notes:
- 1. The value of H*before correction for Pcrev is used because the interaction surface is based on the H*value without the Pcre, adjustment.
68
LTR-SGMP-09-100 NP-Attachment There are a number of utility specific RAIs with numbers in the range of RAI 21 through 23, depending on the specific utility. The responses to utility-specific RAIs are provided under separate cover by the utilities.
Vogtle 24. Reference 1, Page 9-6, Section 9.2.3.1: The feedwater line RAI break heat-up transientis part of the plant design and licensing basis. Thus, it is the NRC staffs position that H* and the "leakage factors," as discussed in Section 9.4, should include considerationof this transient. Explain why the proposed H* and leakage factor values are conservative, even with considerationof the feedwater line break heat-up transient.
WCGS 24. Reference 1, page 9-6, Section 9.2.3.1: The FLB heat-up transientis part of the plant design an d licensing basis.
Thus, it is the staffs position that H* and the "leakage factors," as discussed in Section 9.4, should include consideration of this transient. Explain why the proposed H*
and leakage factor values are conservative, even with consideration of the FLB heat-up transient.
B/B 23. Reference 1, Page 9-6, Section 9.2.3.1: The FLB heat-up transientis part of the plant design and licensing basis. Thus, it is the NRC staffs position that H
- and the "leakage factors,"
as discussed in Section 9.4, should include considerationof this transient. Explain why the proposed H* and leakage factor values are conservative, even with considerationof the FLB heat-up transient.
CPSES 24. Reference 1, page 9-6, Section 9.2.3.1: The FLB heat-up transientis part of the plant design and licens ing basis.
Thus, it is the NRC staff's position that H* and the "leakage factors," as discussed in Section 9.4, should include considerationof this transient.Pleaseexplain why the proposed H* and leakage factor values are conservative, even with consideration of the FLB heat-up transient.
Seabrook 24.Reference 1, Page 9-6, Section 9.2.3.1: The feedwater line break heat-up transientis part of the plant design and licensing basis. Thus, it is the staff's position that H* and the "leakagefactors," as discussed in Section 9.4, should include considerationof this transient. Explain why the proposedH*
and leakage factor values are conservative, even with considerationof the feedwater line break heat-up transient.
Response
Radiological consequences are a function of the source term and activity transport.
- Source term refers to the activity available for release. This is controlled by the Tech Specs on primary and secondary activity and the iodine spike 69
LTR-SGMP-09-1 00 NP-Attachment considerations required by the NRC. Fuel damage is not expected for either the SLB or the FLB. As a result, the source term would be the same. This is the case for the H* plants under consideration.
- Activity transport is dependent upon initial locations of the activity and the mechanism for transport that are applicable. For both the SLB and the FLB events, the dose calculation would use the Tech Spec leakage rate for tube leakage. For both the SLB and FLB events, the secondary break would be assumed to occur outside containment such that the faulted SG releases would occur directly to the atmosphere. As a result, the activity transport would be the same.
Therefore, calculation of the dose consequences for the SLB event would be identical to the calculations that would be made for a FLB event. In this subject, the item that is addressed by the H* program is to define a criterion such that the Tech Spec tube leakage is adequately bounding for both the FLB and SLB. The approach that is used is to define single values for a conservative temperature and conservative pressure differential for the determination of the leakage rate. For the purposes of the dose calculation, these single values are effectively assumed both to be simultaneously occurring and to be continuous for the duration of the calculation. In the dose calculations, these leakage conditions are assumed to last anywhere from multiple hours to multiple days. In some cases, the timeframe is based on allowing for plant cooldown to 212°F which would also require depressurization of the RCS and, therefore reducing the severity of conditions contributing to the tube leakage. However, as noted above, the dose calculations do not consider a more realistic set of conditions for tube leakage.
With respect to temperature transients, a review of the steam generator design transients and the FSAR safety analyses determined that they are not appropriate for defining a temperature basis for the leakage calculations. Those calculations are focused on different criteria and include assumptions which would be overly conservative and operationally limiting for the H* program. For secondary side breaks occurring outside containment, reasonable assumptions would result in a much greater cooling capability of the steam generator secondary side inventory in comparison to FSAR safety analyses. Moreover, based on engineering judgment, a more realistic time-dependent leakage for these events, would result in dose consequences that are less than those reported in the UFSAR under the current licensing basis for both a postulated SLB and FLB event (including a FLB heatup event) due to the reduction in pressure across the tubes and tubesheet that occurs over the long term duration of the accident that is not currently accounted for in the dose analysis. This effect has not been quantified nor does it need to be. It is simpler to define a reasonable peak temperature to use as a basis for the duration of the dose calculation. As identified in WCAP-17071-P and WCAP-17072-P, Westinghouse believes that, with assumptions consistent with an outside containment break, and considering operator actions that are consistent with current Emergency Operating Procedures (EOPs), the no-load condition for a plant is a reasonable condition to use as the basis of the primary to secondary 70
LTR-SGMP-09-100 NP-Attachment leakage. The EOPs, for an event which results in a Safety Injection, provide for reduction/termination of safety injection flow and for initiation of cooldown and depressurization of the reactor coolant system to the point that the RHR system can be placed in operation for continued cooling of the RCS. These actions will significantly reduce the pressure and temperature of the RCS to conditions less than the limiting conditions proposed for use in the H* calculations. The cooldown is initiated by releasing steam from the steam generators by using the systems that are available at the time. The steaming of the steam generators provides the additional benefit of increasing the available AFW flow injection due to a reduced pressure that the AFW pumps have to overcome. Either by considering a reduction in dose consequences due to a more realistic time dependent leakage for these events or by considering that the FLB event is best represented as a cooldown event, it is concluded that no change is required to the leakage factors for the Model F and D5 SGs as reported in WCAP-17071-P and WCAP-17072-P.
However, the NRC staff has pointed out that the "figure of merit" in the technical specification performance criterion is "the leakage rate assumed in the accident analysis" and that a FLB heatup event is part of the current licensing basis for certain plants in the H* fleet. Therefore, to ensure that there is sufficient margin between the accident leakage and operational leakage during a postulated FLB as required by the plant Technical Specifications and to ensure that the implementation of the H* criterion remains within the current licensing basis, an adjustment to the leakage factors provided in Table 9-7 of WCAP-17071-P and WCAP-17072-P has been made that accommodates the design specification FLB heatup event. As noted above, the use of temperatures from this transient is judged to be non-realistic and overly conservative.
As described in WCAP-17071-P and WCAP-17072-P, for the Model F SGs, the FLB design transient represents a double-ended rupture of the main feedwater line concurrent with both a Station Blackout (loss of main feedwater and reactor coolant pump coastdown) and Turbine Trip. For the Model D5 SGs, the maximum RCS temperature of 670°F exceeds the saturation temperature which is not predicted to occur by the worst case Chapter 15 Safety Analysis Transient response.
Because a FLB heat-up event would result in an increase in primary-to-secondary leakage due to the reduction in viscosity of the reactor coolant, the extent of temperature increase must be quantified to address the impact on radiological consequences for the H* plants with Model F and Model D5 steam generators. Referring to references 9-12 and 9-13 of WCAP-17071-P and WCAP-17072-P, the maximum temperature rise for a Model F SG is less than 6.5°F above the normal operating hot leg temperature (approximately 630'F). For the Model D5 steam generator, the maximum increase in temperature is 50°F above the normal operating hot leg temperature (approximately 670°F). This would require a negligible increase in leakage factor for the Model F SGs reported in Table 9-7 (to a maximum of 2.05 from 2.03) and a slight increase in leakage factors reported for the Model D5 steam generators (to a maximum of 2.31 from 2.03).
71
LTR-SGMP-09-100 NP-Attachment The maximum temperature rise for a Model F SG is 66°F above the normal operating cold leg temperature (to approximately 620'F). For the Model D5 steam generator, the maximum increase in cold leg temperature is 120°F above the normal operating cold leg temperature (to approximately 670'F). This would require a maximum increase in leakage factor of 1.23 times the factors provided for the individual Model F SGs in Table 9-7 (to a maximum value of 2.50) and a maximum increase in leakage factor of 1.63 times the factors provided for the individual Model D5 steam generators to a maximum value of 3.16. The leakage factor for the cold leg is limiting for a FLB heat-up event and should be incorporated into the reporting requirements for the plant technical specifications.
Revised versions of Tables 9-1 and 9-7 from WCAP-17071-P and WCAP-17072 -P are provided in this RAI response to reflect the potential increase in temperature that may occur during a postulated FLB event.
Finally, the feedwater line break heat-up transient definition is not a concern for the H*
structural analysis. The SIPC requirement for calculating the end cap load during the faulted condition (1.4DP) during the feedwater line break condition does not exceed the end cap load applied to the tubes during NOP (3.ODP). In fact, the applied end cap load during feedwater line break, regardless of whether it is a heat-up or cool down feedwater line break, is several hundred pounds less than the end cap load applied during normal operating conditions. Therefore, normal operating conditions are bounding for the structural determination of H* in all cases. Please refer to Section 5 and Section 6 in the WCAP 17071-P and WCAP-17072-P for a discussion of the calculated end cap loads and the contact pressure results for the feedwater line break condition.
72
LTR-SGMP-09-100 NP-Attachment Table RA124-1 Revised Table 9-1 Reactor Coolant System Temperature Increase Above Normal Operating Temperature Associated With Design Basis Accidents (References 9-12 and 9-13) a,c,e (1) The postulated SLB does not result in a temperature increase above normal operating conditions as the SLB is a cooldown transient, only the postulated FLB can result in a heatup event dependent upon accident analysis assumptions. The postulated FLB is not part of the licensing basis for plants with Model 44F and Model 51 F SGs.
73
Table RA124-2 Revised Table 9-7 Final H* Leakage Analysis Leak Rate Factors (Revised)
Transient SLB/FLB Locked Rotor Control Rod Ejection FLB- VR 3 ,4 LR/NOP Leak VR3 Leak Adjusted SLB/NOP SLB/FLB AP Ratio VR Leak VR3 CRE LRFA Plant Name AP Ratio @ Leak Rate Adjusted CRE/NOR @Rate
@ Rate a,c,e 2642 @ Factor LR LRF(1" AP Ratio 3030 Factor (High Tavq) p Factor(LRF) 2711 psia (LRF) psia (LRF)
Byron Unit 2 and 1.93 1.61 3.11 Braidwood Unit 2 Salem Unit 1 1.79 1.21 2.16 Robinson Unit 2 1.82 1 1.82 Vogtle Unit 1 and 2 2.02 1.23 2.48 Millstone Unit 3 2.02 1.23 2.49 Catawba Unit 2 1.75 1.52 2.65 Comanche Peak 1.94 1.63 3.16 Unit 2 Vandellos Unit 2 1.97 1.22 2.41 Seabrook Unit 1 2.02 1.23 2.49 Turkey Point Units 1.82 1 1.82 3 and 4 Wolf Creek 2.03 1.23 2.50 Surry Units 1 and 2 1.80 1 1.80 Indian Point Unit 2 1.75 1 1.75 Point Beach Unit 1 1.73 1 1.73
- 1. Includes time integration leak rate adjustment discussed in Section 9.5.
- 2. The larger of the AP's for SLB or FLB is used.
- 3. VR - Viscosity Ratio
- 4. VR - Viscosity Ratio in SG cold leg during a postulated FLB heatup event
LTR-SGMP-09-100 NP-Attachment APPENDIX A
SUMMARY
OF 2008-2009 PULL OUT TEST PROGRAM IN SUPPORT OF H*
75
LTR-SGMP-09-100 NP-Attachment Abstract Steam generator tubes made of Alloy 600 (A600) were hydraulically expanded in AISI 1018 cold rolled, carbon steel, cylindrical collars, which simulate the steam generator tubesheet, and then pulled by an MTS machine out of the collars in order that tube-to-tubesheet joint (hereafter referred to as "joint" or "the joint") strength might be measured. Nine tubes from each of Model F, Model D5, and Model 44F tubes were tested for pull out resistance, three at each expansion length (13 inches, 15 inches, and 17 inches). The pull out test parameters were established so that the results can be considered to be prototypic of the as-built condition of the steam generators within the H* fleet (i.e., the test specimens were designed and manufactured to be within the manufacturing tolerances for dimensional variations, material properties, and process control parameters for the H* fleet steam generator tube joints).
The pull force capacity associated with 0.25 inch tube displacement relative to the tubesheet ranged from approximately [ ],,b,, lbf to approximately [ ]a,bc Ibf. The values for the maximum pull force ranged from approximately [ ]a,bc Ibf to approximately [ ]a,b,c lbf within a maximum relative displacement of 2.02 inches, regardless of the tube outside diameter or hydraulic expansion length[10].
Monte Carlo simulations were performed in order to better define a 5/50 value of pull force, which is based on the presence of residual contact pressure, for use in the H* analysis. The minimum 5/50 value of the pull force has been observed to be [ ]a,bc Ibf, and this corresponds very well to the lowest recorded pull force from the testing.
Introduction H* (pronounced "H star") is the length of hydraulically expanded steam generator tube that must remain intact within the tubesheet in order for the joint to resist pull out and leakage due to normal operating or accident conditions. The basis of the H* program is such that residual contact pressure between the tube and the tubesheet is not considered in the structural or leakage calculations. Hence, any indication of joint strength from test program data is a measure of conservatism contained in the H* analysis.
Westinghouse commenced a test program in which steam generator tubes were hydraulically expanded in cylindrical collars representing the tubesheet and pulled to measure joint strength.
There were no tack expansions, hard rolled expansions, or welds to consider. Initially, the H*
program applied to Model F steam generators, but it has been expanded to include Model D5, Model 44F, and Model 51 F steam generators. The following sections of this document summarize the results of this test program.
76
LTR-SGMP-09-100 NP-Attachment Experimental Materials Alloy 600 tubes representing those from Model F, Model D5, and Model 44F steam generators were cut to seventeen, nineteen, and twenty-one inch lengths. The Model F steam generator tube was taken from Heat NX7368 and is believed to be mill annealed. The Model D5 and Model 44F steam generator tubes were taken from Heats 2645 and 752570, respectively, and both are in the thermally treated condition. The chemical analyses for these materials are contained in Table 1 and the mechanical properties are contained in Table 2. Note that the mechanical properties listed in Table 2 are from the providers' certifications and from testing done at Industrial Testing Laboratory Services (ITLS). The latter tests were done according to ASTM E8-08 [1].
The cylindrical collars representing the tubesheet were cut to fifteen, seventeen, and nineteen inch lengths from AISI 1018 cold-rolled, carbon steel. The chemical analysis and mechanical properties of Heat 777553 are contained in Tables 3 and 4, respectively. It should be noted that the outer diameters of the collars were chosen to be [ ]a,ce times the outer diameter of the tubes so that the stiffness of the actual tubesheet plate is correctly represented. This ratio is based on the work of Middlebrooks et al. [2].
The list of tube and collar pairings is presented in Table 5. The indices are to be read as follows: the first two indices refer to the overall length of the tube or collar, the second three indices refer to the nominal OD of the tube or the nominal ID of the collar, and the last two indices refer to the sample number. The "A" suffix refers to a second manufacture of the same sample. It should be noted that two of the tests done were originally planned to be diagnostic in nature. The collars were rebored so that they would contain an inner diameter surface finish of 250 micro-inch rms max. vice an engineered finish of 250 micro-inch rms. These collars were from Heat 730492, and its properties are also contained in Tables 3 and 4.
Pre-ExpansionMeasurements The inner diameters of the collars were measured by the vendor (Tooling Specialists, Inc.,
Latrobe, PA) at distances corresponding to 25%, 50%, and 75% of the length of the collar, relative to the serialized end. Two measurements, ninety degrees apart, were made with an intramic at each location, and the two values at each location were then averaged. Surface roughness measurements were also made by the vendor at the 25% and 75% distances using a profilometer. Lack of an extension device for the profilometer did not permit roughness measurements at the 50% distance.
After being cut, the inner and outer diameters of the A600 tubes were measured by an intramic and the surface roughness of the outer diameters were measured with a profilometer at Westinghouse RRAS (R. Fetter). The diameter measurements were made, relative to the non serialized end, at distances that overlap those made in the collars. Thus, the 25%, 50%, and 75% distances correspond to those percents of the collar's length, not the tube's length. The 77
LTR-SGMP-09-100 NP-Attachment inner and outer diameters were also measured at two points for each distance, ninety degrees apart, and the two values were averaged.
Hydraulic Expansion The tubes were inserted in the collars such that the non serialized end of the tube was flush with the serialized end of the collar. Thus, the serialized end of the tubes protruded from the collars by two inches. The tube/collar assemblies were then inserted on an O-ring mandrel, which was connected to a screw drive pressurizing system. The tube/collar assemblies were pressurized to a nominal pressure of [ ]ac~e psi per Process Specification 81013RM, Revs. 4 through 10 applicable [ ]. The nominal expansion pressure was typically exceeded, but the excess was less than [ ]ace psi, which is within the tolerance of the equipment ([ ]a'c,' psi). This work was performed at Westinghouse's Waltz Mill facility by M. Gallik and A. Stett. The details of the tube expansion test plan are contained in [4].
Post-ExpansionMeasurements After the hydraulic expansions were completed, measurements of the tubes' inner diameters were again made with an intramic by Westinghouse RRAS (R. Fetter), and eddy current measurements of individual tube/collar assemblies were performed by Westinghouse with the 3-coil +point and standard bobbin coil probes (R. Pocratsky). Once the measurements were complete, the end caps were welded to the tubes at their serialized ends. The tube/collar assemblies with the end caps welded on are shown in Figures 1 through 9.
Heat Treatment Real tubesheet Z-channels are given a post-weld heat treatment (PWHT) with an electric "belt" wrapped around the channel. In order to simulate that PWHT, the tube/collar assemblies were heat treated in air at nominally [ ]a*c'e OF for nominally 3 hours3.472222e-5 days <br />8.333333e-4 hours <br />4.960317e-6 weeks <br />1.1415e-6 months <br /> in a Blue M furnace, Model B-2730-Q. This was accomplished at Westinghouse's Churchill site (A. Neville). The actual PWHT temperature applied to the Z-channel is 11501F. However, it was determined [5] that
]a,ce OF is higher than the vast majority of the tubes will experience by the PWHT.
Instrumentation Prior to testing, the exposed ends of the tubes were fitted with two 350 ohm, quarter bridge strain gauges. Additionally, two linear variable differential transformers (LVDTs) were used in order to model the displacement of each end of the tube relative to the collar. All electronic readouts (load cell, cross-head displacement, displacements of the LVDTs, and strains) were recorded on a Strainbook data acquisition system.
Pull Tests The pull tests were performed according to the test program described in [6] in air and at room temperature. The mechanical operation of the MTS system was performed by M. Gallik and 78
LTR-SGMP-09-100 NP-Attachment A. Stett, while the electronic recording of the data was done by A. Roslund, all of whom work at Westinghouse's Waltz Mill facility. The sequence of the testing was activation of the Strainbook and confirmation of its recording, initiation of the pull test, continuation of the pull test until approximately two inches of cross-head displacement were achieved, and finally, the cessation of pull testing and electronic recording of data.
Post-Test Evaluations After the pull testing was complete, another set of eddy current measurements were made at Westinghouse's Waltz Mill facility by R. Pocratsky. Again, both the 3-coil +point and standard bobbin coil probes were used.
Monte Carlo Analysis In support of the test program, Monte Carlo simulations were run, based on means and standard deviations from the test data, in order to determine a 5/50 bound on pull force. The simulations were performed in two ways and on a tube OD basis: one simulation considered the thirteen inch expansions only, and the other simulation considered all nine tests together, ignoring expansion length difference. Seven distributions were chosen and the fitting parameters set so that the resulting distribution has the mean and standard deviation of the test data. The first was a truncated normal distribution. The other six were chosen so that their domains span zero to positive infinity, their domains are continuous, and fitting parameters are within their allowable ranges. They were lognormal, Erlang, Gamma, inverse Gaussian, Pearson Type V, and Weibull. In each simulation, 100,000 iterations were run.
Discussion of Key Parameters Tube pull out force capacity (based on residual contact pressure) can be derived from the measured pull out forces from the test that simulate the as-manufactured condition of the steam generators. All of the tests performed to date have demonstrated that a positive value of residual contact pressure exists after the hydraulic expansion process. However, the results from these tests depend on a number of factors including dimensional variations of the tubes and tube collars, surface finish variations, potential manufacturing artifacts on the tubesheet (collar) bore, tube joint process variables, and material properties of the test specimens. The key items identified are addressed below.
The NRC staff has raised the concern that sufficient information must be provided to adequately characterize the potential range in values of residual contact pressure between the tube and the tubesheet (due to the hydraulic expansion process) which may be encountered within the whole plant [7]. At that time, only limited pull out data existed upon which the residual contact pressure was estimated. The staff pointed out [7] that the residual contact pressure, and thus the residual load capacity, is highly sensitive to several parameters including hydraulic expansion pressure, tube yield strength, tube material strain hardening properties, and initial (pre-expansion) gap between the tube and the tubesheet. The NRC staff further pointed out in
[7] that additional information was necessary to establish whether the pull out test specimens 79
LTR-SGMP-09-100 NP-Attachment adequately envelop the range of values of those parameters that may be encountered in the as-built steam generators.
Consequently, two actions have been taken to address the NRC staff concerns. First, an analysis was performed to identify the key parameters that affect the residual contact pressure and to quantify the effects of uncertainties. Secondly, a new pull out test program was initiated to provide test results that can be directly compared to the key parameters as identified by analysis in support of the development of the H* criterion.
The analysis model used to evaluate the residual contact pressure was a two-dimensional, plane strain, finite element model using the ANSYS computer code as described in Section 7 of
[8]. Based on a review of Table 7-3 of [8], the key parameters impacting pull out force capacity are:
" Initial tube gap
- Tube yield strength
" Tube joint expansion pressure
" Strain hardening.
Other parameters important to pull out force capacity not considered in the analytical model are surface roughness and variations in the diameter of the tubesheet bore (waviness).
Table 6 provides a comparison of the as-built to as-tested parameters in the new test program.
Based on a review of Table 6, several points can be made regarding the key parameters of the pull out testing.
- It is expected that standard gun drilling practices used in the manufacture of steam generators would typically result in nominal gaps between the tube and tubesheet. No special controls were placed on the initial gap size as the test program was meant to be as prototypic as possible.
- The yield strength of the tubes used for the test specimens simulating the as-built configuration of the Model F and Model D5 steam generators was conservatively high compared to the as-built mean values (f ]a.b.c ksi vs. [ ]a,b,c ksi), because higher yield strengths result in less tube deformation for a given expansion pressure. The yield strength for the tubes used for the Model 44F steam generators was slightly less than the as-built mean yield strength ([ ]ab.C ksi vs. [ ],b-c ksi).
- The expansion pressure used in the manufacture of the test specimen was consistent with what is specified in [3] and is directly applicable to the as-built conditions of the steam generators in the H* fleet.
- The surface roughness of the tubes outer diameters and the collars inner diameters was well within the tolerances of the as-built conditions of the steam generators in the H*
fleet.
- The mechanical properties of the materials used for the test specimens are within ASME Code specifications for the respective materials. Thus, use of the ASME Code values for the key parameters of the H* study is valid.
80
LTR-SGMP-09-100 NP-Attachment Other differences between the materials used for the test specimens for the pull out tests are addressed below.
The use of mill annealed tube vice thermally treated tubing for the Model F specimens has been evaluated and found to be acceptable. For room temperature testing, the key material property affecting the residual contact pressure is yield strength. The difference between yield strength of mill annealed material and thermally treated material is presented in Table 7 and further discussed below. Based on the similarity of mechanical properties between the two materials, it is concluded that there is no adverse effect on the test results. The yield strength value used for the Model F test specimens was [ ]a,b,c ksi, which would result in a reduction of residual contact pressure [9].
The test specimen collar is manufactured from AISI 1018 cold rolled, carbon steel. The material used in the H* fleet is actually A508 Class 2a carbon steel. The use of the different material does not adversely affect the pull test results since the primary property of the material in this case is elastic flexural rigidity of the tubesheet (i.e., elastic modulus), and since the tube expansion operation does not produce significant yielding of the tubesheet (the yield strength of the AISI 1018 cold rolled, carbon steel at room temperature is -83 ksi), the use of higher strength material for the collar is acceptable (see pp. 8-9 of [8]). Thus, it is concluded that the pull out testing is representative of the as-built condition of the steam generators in the H* fleet.
Results Table 8 shows the results of the pull tests, while Table 9 through Table 15 shows the results of the Monte Carlo simulations. The latter results are calculated for the pull out force at 0.25 inch displacement.
Discussion Discussion between Westinghouse and the NRC staff has led to the decision that the pull out force of record should be the pull out force at 0.25 inch cross head displacement. The following discussion and analysis will, therefore, be based on that quantity.
Figure 10 plots the pull out force as a function of the collar ID surface roughness. The graph also provides information on the tube expansion lengths and the tube diameters that were tested. Intuitively, it would be expected that tube pull out force would increase with increasing tube diameter (which provides greater surface area in contact), increasing tube expansion length (which does the same thing), and increasing surface roughness. However, the results in Figure 10 do not necessarily support these assumptions. The highest pull out force for 0.25 inch cross head displacement (approximately [ ]a,bc kips) occurred for both a test specimen with the largest tube OD, the largest collar ID surface roughness, and the smallest expansion length, as well as for a test specimen having the largest tube OD, one of the lowest collar ID surface roughness values, and the smallest expansion length. The next highest pull out force ([ ]a,b-c kips) occurred for tubes with varying degrees of collar ID surface roughness, for all tube ODs, and for all expansion lengths. The lowest pull out force ([ ]a,b,c lbf) occurred for a test 81
LTR-SGMP-09-100 NP-Attachment specimen with a 0.75 inch tube OD, a collar ID surface roughness of - 50 micro-inch (rms), and an expansion length of 15 inches. The lowest pull out force for a Model F test specimen was less than [ ]a,b-c kips. This specimen had a collar ID surface roughness less than 40 micro-inch (rms) and an expansion length of 13 inches. The lowest pull out force for a Model 44F specimen was less than [ ]ab,c kips. This specimen had a collar ID surface roughness of less than 40 micro-inch (rms) and a tube expansion length of 15 inches.
Similarly, the pull test results are shown as a function of tube expansion length in Figure 11.
These results also show the lack of correlation between pull out force and tube OD and expansion length.
The pull force necessary to move a tube in the collar is a consequence of three main factors:
the residual contact pressure due to the hydraulic expansion, the surface roughness of the tube and the collar, and any geometric irregularities due to machining of the tube and collar, which are then subject to hydraulic expansion. As shown by analysis, the initial gap between the outer diameter of the unexpanded tube and the tubesheet bore hole can adversely affect the resulting residual contact pressure. Small variations along the length of the collar ID (waviness) due to the gun drilling process are significant contributors to the pull out resistance. Geometric irregularities are present as initial gaps between the tube and the collar and as bulges in the tubes. One possible explanation for the significant variation in the test results may be that the waviness was not well profiled due to the difficulty of quantifying this variable. Nonetheless, the pull out test results do appear to be consistent with the expected as-built condition.
Recall that nine pull out tests were performed for each tube OD. Analysis of variance (ANOVA in statistics) is a collection of statistical models and their associated procedures in which the observed variance is partitioned into components due to different explanatory variables. In its simplest form, ANOVA provides a statistical test of whether the means of several groups of data are all equal. One such method is called the F-test. Therefore, the F-test was conducted on the pull out test capabilities comparing the variance of each set of 9 samples for each tube diameter using Microsoft EXCEL. The F-test was used to determine whether or not there was any statistical difference between tube OD and pull test results. The answer was that it cannot be concluded that there is any difference in the variance between each sample set and that the means for tube pull out force for each of the outer diameters may be equal. Therefore, it is judged that all of the data can be considered to be one data set.
However, the NRC staff stated in [71 that there is a need to adjust the pull out data so as to produce an estimate of the residual contact pressure that is conservative for the range of H*
values that are being proposed. In order to address this concern for the new pull out test data (i.e., the expansion length of some of the pull out test data exceed the calculated H* values), the sample sets for the different tube ODs were not combined. They were separated by expansion length, even though the F-test results suggest that the mean values of the tube pull out capacity are the same for different tube ODs and considering variations in expansion length and surface roughness.
To investigate this further, the Monte Carlo simulations were performed. Each tube OD was broken up into two sets (13 inch expansion length only and all expansion lengths) and 82
LTR-SGMP-09-100 NP-Attachment distributions were chosen based on the criteria previously defined. Using the calculated means and standard deviations from each data set, the fitting parameters for the seven distributions chosen were calculated. Note that the fitting parameters for the normal and lognormal distributions are simply the mean and standard deviation. In each case the 5/50 value was recorded, and the lowest of these corresponded to a pull force of [ ]a,,,e Ibf. This was calculated for the Model F tube, 13 inch expansion length only, and using the Weibull distribution (see Table 15). This value is very consistent with the lowest actual pull force from the test data ([ ]a,b,c lbf).
83
LTR-SGMP-09-100 NP-Attachment Conclusion Based on the results of the pull tests and Monte Carlo analyses, it is concluded that the end cap load used in the H* analysis can be conservatively reduced by [ ] lbf. H* can then be IV recalculated accordingly.
References (1) ASTM E8/E8M-08, "Standard Test Method for Tension Testing of Metallic Materials,"
West Conshohocken: ASTM International, 2008.
(2) W. B. Middlebrooks, D. L. Harrod, and R. E. Gold, Nuclear Engineeringand Design 143, 1993, pp. 159-169.
(3) Process Specification 81013RM, "Hydraulic Tube Expansion," Rev.4 through Rev. 10, February 1, 1979 through July 24, 1981.
(4) J. T. Kandra, TP-CDME-08-3, "Test Procedure for Tube Expansion for H*," August 25, 2008.
(5) D. L. Harrod, WNEP-9725, "The Westinghouse Tube-to-Tubesheet Joint Hydraulic Expansion Process," July 1997.
(6) J. T. Kandra, TP-CDME-08-1, "Pull-Out Test Program for H*," August 25, 2008.
(7) NRC Letter, "Wolf Creek Generating Station - Withdrawal of License Amendment on Steam Generator Tube Inspections (TAC No. MD0197)," United States Nuclear Regulatory Commission, Washington, D.C., February 28, 2008.
(8) WCAP-17071-P, "H*: Alternate Repair Criteria for the Tubesheet Expansion Region in Steam Generators with Hydraulically Expanded Tubes (Model F);" April 2009.
(9) DP-SGDA-05-2, "Data Package for H-Star Pull Test of 7/8 Inch Tubing form Simulated Tubesheet, PA-MSC-0199 WOG Program for Steam Generator Models 44F and 51F,"
November 2005.
(10) LTR-SGMP-09-98, "H* Pull Test Program Summary," Westinghouse Electric Company LLC, July 27, 2009.
84
LTR-SGMP-09-100 NP-Attachment Table 1 Chemical Analyses of the A600 Materials Used in This Test Program Steam Generator Model MdlF D5 44F Chemical Analysis Plymouth Tube Co. Huntington Alloys, Sourcemia Anliuyt Tb C Inc. AB Sandvik Steel Source Salisbury, MD HnigoW Huntington, WVV Heat NX7368 2645 752570 Element (W/o)
C 0.04 0.033 0.025 Mn 0.41 0.34 0.79 P N/A 0.007 0.009 S 0.001 0.001 0.002 Si 0.30 0.09 0.33 Cr 14.87 15.44 16.60 Ni 76.21 75.45 72.45 Cu 0.15 0.23 0.010 Co 0.04 0.04 0.011 Fe 7.98 8.42 9.29 B N/A 0.003 N/A Table 2 Mechanical Properties of the A600 Materials Used in This Test Program Steam Mechanical Generator Heat Property j v (psi) CUT (psi) Elongation (%)
Model Source NX7368 Vendor 59,700 106,600 39 F
ITLS 58,000 108,000 32 D5 2645 Vendor 43,000 97,000 41.5 ITLS 54,000 110,000 35 44F 752570 Vendor 47,500 101,700 45.5
_ _ ITLS 46,000 101,000 40 85
LTR-SGMP-09-100 NP-Attachment Table 3 Chemical Analyses of the 1018 Cold-Rolled, Carbon Steel Used in This Test Program I ~~ i - _______Element (w/o____
IChemical Steel Analysis Heat C Mn Si S P
__ _ Source Steel Bar AISI Corp. 777553 0.17 0.84 0.27 0.030 0.005 1018 Greensboro, NC Steel Bar AISI Corp. 730492 0.18 0.79 0.22 0,030 0.010 1018 Greensboro, NC Table 4 Mechanical Properties of the 1018 Cold-Rolled, Carbon Steel Used in This Test Program Mechanical Steel Heat Property (7y (ksi) GUT (ksi) Elongation (%)
Source DuBose National AISI 1018 777553 Energy 83.0 90.0 18 Services, Inc.
Clinton, NC DuBose National AISI 1018 730492 Energy 67.5 79.3 25 Services, Inc.
Clinton, NC 86
LTR-SGMP-09-100 NP-Attachment Table 5 Steam Generator Tube and Collar Pairings Used in This Test Program Tube Heat Collar Heat 17-688-01A NX7368 15-699-01A 777553 17-688-02A NX7368 15-699-02A 777553 17-688-03 NX7368 15-699-03A 777553 19-688-01 NX7368 17-699-01A 777553 19-688-02 NX7368 17-699-02A 777553 19-688-03 NX7368 17-699-03A 777553 21-688-01 NX7368 19-699-OIA 777553 21-688-02 NX7368 19-699-02A 777553 21-688-03 NX7368 19-699-03A 777553 17-750-01A 2645 15-762-01A 777553 17-750-02A 2645 15-762-02A 777553 17-750-03 2645 15-699-03 730492 19-750-01 2645 17-762-01A 777553 19-750-02 2645 17-762-02A 777553 19-750-03 2645 17-762-03A 777553 21-750-01 2645 19-762-OIA 777553 21-750-02 2645 19-762-02A 777553 21-750-03 2645 19-762-03A 777553 17-875-01A 752570 15-888-01A 777553 17-875-02A 752570 15-888-02A 777553 17-875-03 752570 15-762-03 730492 19-875-01 752570 17-888-01A 777553 19-875-02 752570 17-888-02A 777553 19-875-03 752570 17-888-03A 777553 21-875-01 752570 19-888-01A 777553 21-875-02 752570 19-888-02A 777553 21-875-03 752570 19-888-03A 777553 87
LTR-SGMP-09-100 NP-Attachment Table 6 Residual Contact Pressure Critical Parameter Comparison Key Parameters Model F Model D5 Models 44F -51F As-Built As-Tested As-Built As-Tested As-Built As-Tested Average Initial -a,c,e Gap (inches)
Tube Yield Strength (ksi)
Expansion Pressure (ksi)
Tube Outer Diameter Surface Roughness p in. rms Collar Inner Diameter Surface Roughness p in. rms Tube OD (in)
Collar ID (in)
Table 7 Comparison of Yield Strength Between Mill Annealed and Thermally Treated Alloy 600 AlloAs:TestMi Anne' ed `led a Trae Alloy 6OQTherfal W....
i..
a,c,e Minimum Mean Maximum Standard Deviation Number of Tests 361 307 Tube Size (OD) 7/8 inch 7/8 inch Data Reference [1] Reference [1]
Yield Strength values are in units of ksi.
88
LTR-SGMP-09-100 NP-Attachment Table 8 Results of the Pull Testing Load at 0.25" Max. Displacement Displacement Tube ID Heat Collar ID Heat (kip) Load (kip) at Max. Load (in) a,b,c 17-688-01A NX7368 15-699-01A 777553 F 17-688-02A NX7368 15-699-02A 777553 17-688-03 NX7368 15-699-03A 777553 19-688-01 NX7368 17-699-01A 777553 19-688-02 NX7368 17-699-02A 777553 19-688-03 NX7368 17-699-03A 777553 21-688-01 NX7368 19-699-01A 777553 21-688-02 NX7368 19-699-02A 777553 21-688-03 NX7368 19-699-03A 777553 17-750-01A 2645 15-762-01A 777553 17-750-02A 2645 15-762-02A 777553 17-750-03 2645 15-699-03 730492 19-750-01 2645 17-762-01A 777553 19-750-02 2645 17-762-02A 777553 19-750-03 2645 17-762-03A 777553 21-750-01 2645 19-762-01A 777553 21-750-02 2645 19-762-02A 777553 21-750-03 2645 19-762-03A 777553 17-875-01A 752570 15-888-01A 777553 17-875-02A 752570 15-888-02A 777553 17-875-03 752570 15-762-03 730492 19-875-01 752570 17-888-01A 777553 19-875-02 752570 17-888-02A 777553 19-875-03 752570 17-888-03A 777553 21-875-01 752570 19-888-01A 777553 21-875-02 752570 19-888-02A 777553 21-875-03 752570 19-888-03A 777553 89
LTR-SGMP-09-100 NP-Attachment Table 9 Monte Carlo Results for the Truncated Normal Distribution Distribution Normal Distribution (truncated at 0) 5/50 Pull Out Case 1 Parameters to Define the Distribution Force (kip)
Model F Name Symbol ________ Mean Stand. Dev. [ ] a~c~e 13" Expansion Valu G c 5/50 Pull Out Case 2 Parameters to Define the Distribution Force( kip)
ModelFSybli Name Mean Stand. Dev.
_______ace All Expansions VSymbol Value[ ] a,[,,[ ] a,,, ] c,e Case 3 Parameters to Define the Distribution 5/50 Pull Out Force (kip)
Model D5 Sy Name *-mb'ol Mean Stand. Dev.
________ ]
[_______ a'c'e ac 13" Expansion Valu Value [ ai
] a,c,e I a,c,e 5/50 Pull Out Case 4 Parameters to Define the Distribution Force (ki )
Model D5 Name Mean Stand. Dev.
All Expansions Symbol ___ ]______ G] a ace Value [ ,ce a,,
5/50 Pull Out Case 5 Parameters to Define the Distribution Force (kip)
Model 44F Name Smo _______
Mean Stand. Dev. ~~
13" Expansion Symbol Value[ i. j
]ace[ G]aoe I, [~iL 5/50 Pull Out Case 6 Parameters to Define the Distribution Forc P)
Force (ki)
Name Mean Stand. Dev.
Model 44F S[ ace All Expansions VSymbol Value[ -3"[]jaaC,e ]ace[ ]ace 90
LTR-SGMP-09-100 NP-Attachment Table 10 Monte Carlo Results for the LogNormal Distribution Distribution Lognormal Distribution 5/50 Pull Out Parameters to Define the Distribution Forc P)
Case 1 Force (kip Name Mean Stand. Dev.
Model F Symbol __ a [ ] a,c,e 13" Expaso ... Value ]a,c,e [ aJc,e Distriution5/50 Pull Out Parameters to Define the Distribution P Case 2 Force (kip)
Model F Name Mean Stand. Dev.
Symbol aace All Expansions. Value ace [I] a,c,e 5/50 Pull Out Case 3 Parameters to Define the Distribution Force (ki )
Model D5 Name Mean Stand. Dev.
Symbol ______[ a c,e 13" Expansion Value Case 4 Parameters to Define the Distribution 5/50 Pull Out Force (kip)
Model D5 Name Symbol Mean
________[_]___ Stand. Dev. a,*
All Expansions Value [ ][
5/50 Pull Out Case 5 Parameters to Define the Distribution Forc P)
Force (kip)
Model 44F Name Mean Stand. Dev.
13" Expansion Symbol a,c,e a,c,e 5/50 Pull Out Case 6 Parameters to Define the Distribution Force ki Name Mean Stand. Dev.
Symbol a,c,e Model 44F __
All Expansions Value _ ____
91
LTR-SGMP-09-100 NP-Attachment Table 11 Monte Carlo Results for the Erlang Distribution Distribution Erlang Distribution 5/50 Pull Out Case 1 Parameters to Define the Distribution Forc P)
Force (kip)
Model F Name Cont. Shape Par. Cont. Scale Par.
Symbol m [ lace 13" Expansion Value ] a,c,e ] a,c,e 5/50 Pull Out Case 2 Parameters to Define the Distribution Force (kip)
Model F All F All Expansions Name Symbol Value Cont. Shape Par.
[ m] ae]_a,c,e_[
1Cont. Scale P[
Par.
]a cIe 5/50 Pull Out Case 3 Parameters to Define the Distribution Forc ki Force (kip)
Model D5 Name Cont. Shape Par. Cont. Scale Par.
13" Exao Symbol m P[ Ia~c, 13" Expansion Value ] a,j,e ],j,e I 5150 Pull Out Case 4 Parameters to Define the Distribution Forc P)
Force (kip)
Model D5 Name Cont. Shape Par. Cont. Scale Par.
All Expansions Symbol m ac,e
] a,c,e ] a,c,e Value Case 5 Parameters to Define the Distribution 5/50 Pull Out Force (kip)
Model 44F Name Cont. Shape Par. Cont. Scale Par.
13" Symbol m P I ac Expansion Value [ ] a,c,e I a,c,e 5/50 Pull Out Case 6 Parameters to Define the Distribution Forc P Force (kip)
Model 44F Name Cont. Shape Par. Cont. Scale Par.
All Expansions Symbol m I [ ]ea cce Value I ,c,[ _
92
LTR-SGMP-09-100 NP-Attachment Table 12 Monte Carlo Results for the Gamma Distribution Distribution Gamma Distribution Parameters to Define the Distribution 5/50 Pull Out Case 1 Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
Model F Symbol U 13 [ I a,c,e 13" Expansion Value j ai2e [ j a,c,e 5/50 Pull Out Case 2 Parameters to Define the Distribution Force (kip)
Model F Name Symbol Cont. Shape Par.
_________ Cont. Scale 3 Par. [ ] ac~e All Expansions Value [ c] a,,e [ ]a,o,e 5/50 Pull Out Case 3 Parameters to Define the Distribution Forc P)
Force (ki )
Model D5 Name Symbol Cont. Shape Par.
_________1 Cont. Scale Par. [ ]a~c~e 13" Expansion Value [ ] a,c,e a,c e Case 4 Parameters to Define the Distribution 5/50 Pull Out Force (kip)
Model D5 Name Cont. Shape Par. Cont. Scale Par.
Moexaso Symbol aX ]PacIe All Expansions Value I[~ ] a,c,e 5/50 Pull Out Case 5 Parameters to Define the Distribution Force (kip)
Model 44F Name Cont. Shape Par. Cont. Scale Par.
Expansion Symbol ace a,c,e [ ]ac~e 13" Expansion Value [ ]jjj 5/50 Pull Out Case 6 Parameters to Define the Distribution Forc P)
Force (kip)
Model 44F Sybodel_____44F Name Cont. Shape Par. Cont. Scale Par. ac All Expansions Symbol Value I Cac] a,c,e_[]a 93
LTR-SGMP-09-100 NP-Attachment Table 13 Monte Carlo Results for the Inverse Gaussian Distribution Distribution ] Inverse Gaussian Distribution 5/50 Pull Parameters to Define the Distribution Forc P Out )
Case 1 Force (kip)
Name Model FSybl________xjce Cont. Par. Cont. Par.
13" Expansion VSymbol Iuj [ ]ecc, Defin t5/50 Pull Out Case 2 Parameters to Define the Distributionki)
Model FSybl_______[Name Cont. Par. Cont. Par.
Value [ X All Expansions 5/50 Pull Out Case 3 Parameters to Define the Distribution Forc P)
Force (kip)
Name Cont. Par. Cont. Par.
Model D5 Symbol _ __ _[ ]a,c,e 13" Expansion Value Case 4 Parameters to Define the Distribution 5/50 Pull Out Cont. Par. Cont. Par.
Name Model D5 Symbol _ _ [_I ] a,c,e All Expansions Value a,,e 5/50 Pull Out Parameters to Define the Distribution Forc P)
Case 5 Force (kip)
Name Model 44FSybl________[ Cont. Par. Cont. Par. ce 13" Expansion Symbol Value Ia]ace
] [ ace [,c,e 5/50 Pull Out Case 6 Parameters to Define the Distribution Forc P Force (kip)
Model 44F All Expansions A__lExpansions_
Name Symbol Value
]-,[ Cont. Par.
a,c,e Cont. Par.
[_]_a___e
]a,c,e 94
LTR-SGMP-09-100 NP-Attachment Table 14 Monte Carlo Results for the Pearson Type V Distribution Distribution [Pearson Type V Distribution 5/50 Pull Out Parameters to Define the Distribution Forc P)
Case 1 Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
Model F Symbol Cac 13 [ ]acc e 13" Expansion Value [ace Distriution5/50 Pull Out Parameters to Define the Distribution P Case 2 Model F Name Symbol Cont. Shape oa Par. 1 Par.
Cont. Scale [ ] a~ce All Expansions Value a,c,e [,c,e 5/50 Pull Out Case 3 Parameters to Define the Distribution Force (kip)
Model D5 Name Symbol 1Cont. Shape Par. 13 Par.
Cont. Scale [ ]ace 13" Expansion Value a,c,e ]a,c,e Case 4 Parameters to Define the Distribution 5/50 Pull Out Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
Model D5 Symbol a* 1 [ ] a~c~e All Expansions Value ac,,e ],c,e 5/50 Pull Out Case 5 Parameters to Define the DistributionFoc(kp Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
13" Expansion Valuembl a,c,e [a,c,e 5/50 Pull Out Case 6 Parameters to Define the Distribution Forc Pu Force (kip)
Model 44F Name Cont. Shape Par. Cont. Scale Par.
All Expansions Symbol a,c, 13 ]ac.e Value [_[
95
LTR-SGMP-09-100 NP-Attachment Table 15 Monte Carlo Results for the Weibull Distribution Distribution Weibull Distribution 5/50 Pull Out Parameters to Define the Distribution Forc P)
Case 1 Force (kip)
Model F Name Cont. Shape Par. Cont. Scale Par.
[
Symbol Ct _ Pac~e 13" Expansion Value] a,c,e ]a,c,e 5/50 PullO 5/50 Pull Out Case 2 Parameters to Define the Distribution Force kip)
Name Cont. Shape Par. Cont. Scale Par.
Model F Symbol o_
___ _[ ]PacIe All Expansions Value ] a,c,e I a,c,e 5/50 Pull Out Case 3 Parameters to Define the Distribution Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
Model D5 Symbol a 13 [ ]Pace 13" Expansion Value a,c, ace 5/50 Pull Out Case 4 Parameters to Define the Distribution Forc P)
Force (kip)
Model D5 Symbol Name Cont. Shape ace Par. 3ao~ Par.
Cont. Scale [ ] acoe All Expansions Value ace 5/50 Pull Out Parameters to Define the Distribution Forc P Case 5 Force (kip)
Model 44F Name Symbol Cont. Shape a Par. 13[ Par.
Cont. Scale ] a~ce 13" Expansion Value [ ]_a,c,e [ a,c,e 5/50 Pull Out Parameters to Define the Distribution Forc P)
Case 6 Force (kip)
Name Cont. Shape Par. Cont. Scale Par.
Model 44F Symbol cjac ]Pace All Expansions Value [ ____I___Ic 96
LTR-SGMP-09-100 NP-Attachment Figure 1 The Model F 13 Inch Expansion Tube/Collar Assembly a,c,e I
97
LTR-SGMP-09-100 NP-Attachment Figure 2 The Model F 15 Inch Expansion Tube/Collar Assembly a,c,e 98
LTR-SGMP-09-100 NP-Attachment Figure 3 The Model F 17 Inch Expansion Tube/Collar Assembly a,c,e 99
LTR-SGMP-09-100 NP-Attachment Figure 4 The Model D5 13 Inch Expansion Tube/Collar Assembly ac,e 100
LTR-SGMP-09-100 NP-Attachment Figure 5 The Model D5 15 Inch Expansion Tube/Collar Assembly ace 101
LTR-SGMP-09-100 NP-Attachment Figure 6 The Model D5 17 Inch Expansion Tube/Collar Assembly ace 102
LTR-SGMP-09-100 NP-Attachment Figure 7 The Model 44F 13 Inch Expansion Tube/Collar Assembly _
- a,c,e 103
LTR-SGMP-09-1 00 NP-Attachment Figure 8 The Model 44F 15 Inch Expansion Tube/Collar Assembly 1 a,c,e 104
LTR-SGMP-09-100 NP-Attachment Figure 9 The Model 44F 17 Inch Expansion Tube/Collar Assembly a,c,e 105
LTR-SGMP-09-100 NP-Attachment Figure 10 A Plot of Pull Out Force vs. Surface Roughness a,c,e 106
LTR-SGMP-09-100 NP-Attachment Figure 11 The Pull Out Force vs. Expansion Length for a Given Tube OD a,c,e 107
Attachment 8 Westinghouse Electric Company LLC LTR-SGMP-09-109 NP-Attachment, "Response to NRC Request for Additional Information on H*; RAI #4; Model F and Model D5 steam Generators"
WESTINGHOUSE NON-PROPRIETARY CLASS 3 LTR-SGMP-09-109 NP-Attachment Westinghouse Electric Company Response to NRC Request for Additional Information on H*; RAI #4; Model F and Model D5 Steam Generators August 25, 2009 Westinghouse Electric Company LLC P.O. Box 158 Madison, PA 15663
© 2009 Westinghouse Electric Company LLC All Rights Reserved
- Electronicallyapproved records are authenticated in the Electronic Document Management System 1
LTR-SGMP-09-109 NP-Attachment Response to NRC Request for Additional Information on H*; RAI #4; Model F and Model D5 Steam Generators
References:
- 1. NL-09-0547, Vogtle Electric Generating Plant License Amendment Request to Revise Technical Specification(TS) Sections 5.5.9, "Steam Generator (SG) Program" and TS 5.6.10, "Steam Generator Tube Inspection Report for Permanent Alternate Repair Criteria,"
Southern Company, May 19, 2009.
- 2. RS-09-071, "License Amendment Request to Revise Technical Specifications (TS) for Permanent Alternate Repair Criteria," Exelon Nuclear, June 24, 2009.
Docket Nos. 50-445 and 50-446, License Amendment Request 09-007, Model D5 Steam Generator Alternate Repair Criteria," Luminant, June 8, 2009.
- 4. SBK-L-09118, "Seabrook Station: License Amendment Request 09-03; Revision to Technical Specification 6.7.6.k, "Steam Generator (SG) Program," for Permanent Alternate Repair Criteria (H*)," May 28, 2009.
- 5. Vogtle Electric Generating Plant, Units 1 and 2, Request for Additional Information Regarding Steam Generator Program (TAC Nos. ME1339 and ME1340)," United States Nuclear Regulatory Commission, July 10, 2009.
- 6. Braidwood Station, Units 1 and 2, and Byron Station, Unit Nos. 1 and 2 - Request for Additional Information Related to Steam Generator Permanent Alternate Repair Criteria (TAC Nos. ME1613, ME1614, ME1615, and ME1616)," United States Nuclear Regulatory Commission, July 20, 2009.
- 7. Comanche Peak Steam Electric Station, Units 1 and 2 - Request for Additional Information Regarding the Permanent Alternate Repair Criteria License Amendment Request (TAC Nos.
ME1446 and ME1447)," United States Nuclear Regulatory Commission, July 23, 2009.
- 8. WCAP-17071-P, "H*: Alternate Repair Criteria for the Tubesheet Expansion Region in Steam Generators with Hydraulically Expanded Tubes (Model F)," Westinghouse Electric LLC, April 2009.
- 9. WCAP-17072-P, "H*: Alternate Repair Criteria for the Tubesheet Expansion Region in Steam Generators with Hydraulically Expanded Tubes (Model D5)," Westinghouse Electric LLC, May 2009.
- 10. "Vogtle Electric Generating Plant, Units I and 2, Request for Additional Information Regarding Steam Generator Program (TAC Nos. ME1339 and ME1340)," United States Nuclear Regulatory Commission, August 5, 2009 2
LTR-SGMP-09-109 NP-Attachment
- 11. LTR-SGMP-09-100, "Response to NRC Request for Additional Information on H*; Model F and Model D5 Steam Generators," August 2009
- 12. LTR-NRC-09-26, "LTR-SGMP-09-66 P-Attachment, "White Paper: Low Temperature Steam Line Break Contact Pressure and Local Tube Bore Deformation Analysis for H*
(Proprietary)," Westinghouse Electric Company LLC, May 13, 2009.
- 13. Seabrook Station, Unit No. 1- Request for Additional Information Regarding Steam Generator Program (TAC Nos. ME1 386)," United States Nuclear Regulatory Commission, August 13, 2009 Introduction In response to formal requests for technical specification amendments, References 1, 2, 3 and 4, the USNRC formally requested additional information in References 5, 6, 7 and 13. The Vogtle, Seabrook, Byron/Braidwood and Comanche Peak requests for a permanent license amendment to implement H* represent the Model F and Model D5 steam generators for which the H* technical justification is provided in References 8 and 9.
Subsequent to the initial issue of the RAI (References 5, 6, 7 and 13), the NRC issued follow-up questions (Reference 10) to questions numbers 4, 20 and 24 and an additional request regarding a technical specification (TS) commitment for applying the leakage factors. Except for RAI#4, responses to all of the RAIs, including the follow-up questions in Reference 10, were provided in Reference 11. The affected licensees provided separate responses in regard to the commitment for applying leakage factors.
The response to RAI#4 required additional explanation as discussed with the NRC staff on August 11, 2009 and was, therefore, not included in Reference 11. The additional questions related to RAI#4 that were identified during the August 11, 2009 telephone conference were summarized by Westinghouse and were the basis of the discussion at a meeting among the NRC, several licensees and Westinghouse on August 17 and 18, 2009. These additional questions are reproduced in the response to RAI#4, below. Specific discussion is included in the response to address the additional questions.
To summarize, this document provides the response to the initial RAI#4 as included in References 5, 6, 7 and 13, response to the follow-up question relating to RAI#4 in Reference 10 and response to the additional questions raised during the conference call on August 11, 2009.
Utilities, other than referenced in this document, have requested amendments to their licensees in parallel with the response to these RAI's. The technical RAIs are generic in nature because the analysis methods are the same for all affected plants. Therefore, this response to RAI#4 is generic for all Models of SGs that are candidates for application of H*. However, this letter 3
LTR-SGMP-09-109 NP-Attachment specifically augments Reference 11 to complete the responses to NRC RAIs for WCAP-17071-P (Model F H*) and WCAP-17072-P (Model D5 H*).
Vogtle 4. Reference 1, page 6-69: In Section 6.2.5.3, it is concluded RAI that the tube outside diameter and the tubesheet tube bore PartA inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps prior to the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in Section 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
WCGS 4. Reference 1, page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
B/B 4. Reference 1, Page 6-7: In Section 6.2.5.3, it is concluded that the tube outside diameterand the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location.
4
LTR-SGMP-09-109 NP-Attachment CPSES 4. Reference 1, page 6-70: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in Section 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
Seabrook 4. Reference 1, page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameteralways maintain contact in the predicted range of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps priorto the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location?
PartB: The additionalquestions relatingto RAI#4 as provided in Reference 10 are:
Address following questions as part of response to RAI#4 (Vogtle):
- a. Clarify the nature of the finite element model ("slice" model versus axisymmetric SG assembly model) used to generate the specific information in Tables 6-1, 2, and 3 (and accompanying graph entitled "EllipticalHole Factors")of Reference 6-15. What loads were applied? How was the eccentricity producedin the model? (By modeling the eccentricity as part of the geometry? By applying an axisymmetric pressure the inside of the bore?) Explain why this model is not scalable to lower temperatures.
- b. Provide table showing maximum delta diameters (total diameter distortion) and maximum eccentricities (maximum diameter minus minimum diameter) from the 3 dimensional (3-D) finite element analysis for normal operating and steam line break (SLB), for model F and D5.
- c. In Figure 2 of the White Paper,add plot for original relationshipbetween reductions in contact pressure and eccentricity as given in Reference 6-15 in the graph accompanying Table 6-3. Explain why this original relationshipremains conservative 5
LTR-SGMP-09-109 NP-Attachment in light of the new relationship. Explain the reasons for the differences between the curves.
- d. When establishing whether contact pressure increases when going from normal operating to steam line break conditions, how can a valid and conservative comparison be made if the normal operating case is based on the original delta contact pressure versus eccentricity curve and the SLB case is based on the new curve?
PartC: Additional Questions Providedin the August 11, 2009 telephone conference:
- a. Overall Hiqh Level Question
- 1. Discuss if the eccentricity effect on contact pressureis occurring as described. It is the opinion of the NRC staff that the eccentricity effect may not be as significant as being reportedby Westinghouse.
- b. Other Key Questions
- 1. The eccentricitiesincluded in Table RAI 4-4 appearlargerthan anticipated. Need to confirm that positive contactpressure exists around the entire circumference of the tube and state this clearly in the response.
- 2. The difference between initialand final eccentricity included in Table RAI4-2 needs to be explained. In particular,the exclusive use of the relationship between initial eccentricity and scale factor in calculating contactpressure needs to be justified.
- 3. The basis for applying the correlation for scale factor outside an "eccentricity"range of between IE-3 to IE-4 inch in the calculationof contact pressure needs to be further explained. Values for displacements included in Table 6-18 (of WCAP-1 7071) suggest that contact pressure may be lost at displacementranging between IE-3 in to IE-4 inch.
- 4. Provide the calculation basis for the upperand lower curves provided in Figure RAI 4-2
- 5. Resolve the apparentinconsistency between Item 4 on page 25 and the statement below Figure RAI4-1 regardinghow the model in Figure RAI4-1 is loaded.
- c. Key Remaininq Issues
- 1. Provide the basis for why the ADhole adjustment for contact pressure made using the old model remains conservative.
- 2. Provide an appropriatebasis for demonstrating thatjoints tighten during a postulated SLB event. Why is it acceptable to compare the contact pressures calculatedusing the originalmodel for NOP to the contact pressures calculated using the new model for SLB for the Model D5 SGs?
- 3. If both old and new models are conservative, is there an appropriatebasis to show the relative conservatism of the methods?
6
LTR-SGMP-09-109 NP-Attachment To facilitate a continuous response to the total RAI#4 questions, the questions received originally (Part A), those received as follow-up questions (Part B) and those identified during the 8/11/09 telephone conference (Part C) are re-arranged as noted below. The location of responses to specific questions is shown in bold type after the question. Also, in the responses, the specific questions addressed by the responses are repeated in bold type in the box at the start of the response.
PartC: Sub a.
Discuss if the eccentricity effect on contact pressure is occurring as described. It is the opinion of the NRC staff that the eccentricity effect may not be as significant as being reportedby Westinghouse. (See Section 1.0)
PartB Address following questions as part of response to RAI#4 (Vogtle):
- a. Clarify the nature of the finite element model ("slice" model versus axisymmetric SG assembly model) used to generate the specific information in Tables 6-1, 2, and 3 (and accompanying graph entitled "EllipticalHole Factors') of Reference 6-15. What loads were applied? How was the eccentricity produced in the model? (By modeling the eccentricity as part of the geometry? By applying an axisymmetric pressure the inside of the bore?) Explain why this model is not scalable to lower temperatures. (See Section 1.2)
- b. Provide table showing maximum delta diameters (total diameter distortion) and maximum eccentricities (maximum diameter minus minimum diameter) from the 3 dimensional (3-D) finite element analysis for normal operating and steam line break (SLB), for model F and D5. (See Section 1.1)
- c. In Figure 2 of the White Paper,add plot for originalrelationship between reductionsin contact pressure and eccentricity as given in Reference 6-15 in the graph accompanying Table 6-3. Explain why this originalrelationshipremains conservative in light of the new relationship. Explain the reasons for the differences between the curves. (See Section 4.1)
- d. When establishing whether contact pressure increases when going from normal operating to steam line break conditions, how can a valid and conservative comparison be made if the normal operating case is based on the original delta contact pressure versus eccentricity curve and the SLB case is based on the new curve? (See Section 4.2) 7
LTR-SGMP-09-109 NP-Attachment PartC: Sub b. Other Key Questions
- 1. The eccentricitiesincluded in Table RAI 4-4 appearlargerthan anticipated. Need to confirm that positive contact pressure exists around the entire circumference of the tube and state this clearly in the response. (See Section 3)
- 2. The difference between initial and final eccentricity included in Table RAI4-2 needs to be explained. In particular,the exclusive use of the relationshipbetween initialeccentricity and scale factor in calculating contactpressure needs to be justified. (See Section 1.2)
- 3. The basis for applying the correlationfor scale factor outside an "eccentricity"range of between IE-3 to IE-4 inch in the calculation of contact pressure needs to be further explained. Values for displacements included in Table 6-18 (of WCAP-17071) suggest that contact pressure may be lost at displacement ranging between 1E-3 in to IE-4 inch.
(See Section 2.0)
- 4. Provide the calculation basis for the upper and lower curves provided in Figure RAI 4-2.
(See Section 2.1)
- 5. Resolve the apparent inconsistency between Item 4 on page 25 and the statement below Figure RAI4-1 regardinghow the model in Figure RAI4-1 is loaded. (See Section 1.2)
Part C: Sub c. Key Remaining Issues
- 1. Provide the basis for why the ADhole adjustment for contactpressure made using the old model remains conservative. (See Section 2.2)
- 2. Provide an appropriatebasis for demonstratingthat joints tighten during a postulated SLB event. Why is it acceptable to compare the contact pressurescalculated using the original model for NOP to the contactpressures calculatedusing the new model for SLB for the Model D5 SGs? (See Section 2.3)
- 3. If both old and new models are conservative, is there an appropriatebasis to show the relative conservatism of the methods? (See Section 2.4)
PartA: (OriginalRAI#4 from Reference 5)
Reference 1, Page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameter always maintain contact in the predictedrange of tubesheet displacements. However, for tubes with through-wall cracks at the H* distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps prior to the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in 6.2.5.3 valid for the entire H* distance, given the possibility that the tubes may contain through-wall cracks at that location? (See Section 5.0) 8
LTR-SGMP-09-109 NP-Attachment 1.0 General Background on Approach and Models Discuss if the eccentricity effect on contactpressure is occurring as described. It is the opinion of the NRC staff that the eccentricity effect may not be as significantas being reportedby Westinghouse.
Response
The reference structural model for the H* calculation as described in References 8 and 9 is a 3D FEA model that utilizes the equivalent properties approach for perforated plates in accordance with Reference 6-15 of the H* WCAP reports. This model provides the tubesheet displacements that are utilized in the calculation of H*. Included in the displacement output from the 3D FEA model are the radius and depth dependent x- and y- axis displacements for the tubesheet. These displacements are the input to the H* integrator model that uses the inputs to calculate contact pressures based on thick-shell equations. The tubesheet displacements from the FEA model indicate that the tubesheet bores become eccentric after application of all thermal and pressure loads. The displacement results from the 3D FEA model are the difference between the completely unloaded case and the fully loaded case for the conditions of interest (i.e., NOP, SLB).
The information from the 3D FEA model, that the tubesheet bores become eccentric, led to a question regarding continued tube-to-tubesheet contact in the eccentric tubesheet bore. The impact of tubesheet bore hole out-of-roundness (eccentricity) on the calculation of tube to-tubesheet contact pressures was originally addressed using a scale factor approach as described below and in Reference 6-15 of the H* WCAP reports. The fit developed in Reference 6-15 , a third order polynomial, was appropriate for the conditions for which it was developed but it provided physically impossible results when extrapolated significantly outside its data basis such as was the case for the SLB conditions for the Model D5 SGs.
To resolve this issue, a separate model, was developed as described in Section 6.2.5 and shown in Figure 6-48 of Reference 8 and 9, to assess tube-to-tubesheet contact under the fully loaded condition (e.g., AP and thermal loading) for the small eccentricities that were calculated during the much "colder" temperature postulated SLB conditions for the Model D5 SGs than for the Model F SGs. To properly represent the tube in tubesheet condition, this model considered a tubesheet equivalent cell (the local TS material around a tubesheet bore) and a tube. To address the question if continued contact would exist between the tube and tubesheet after the tubesheet bore becomes eccentric, the tube expansion was analytically simulated to provide a condition of tube to tubesheet contact in a non-eccentric tubesheet bore. This condition was the reference condition for the subsequent loading of the model by pressure loads (thermal loads were not included) and by applying displacement boundary conditions (e-bar) to simulate the expected range of tubesheet bore eccentricity. The unloaded, post-tube expansion simulation conditions of the model was the reference condition for the displacements provided in Tables 6-18 and 6-19 of the H* reports, References 8 and 9.
9
LTR-SGMP-09-109 NP-Attachment While eccentricity was the specific focus of this study because of the question raised about continued tube to tubesheet contact in an eccentric condition, the analytical model naturally also provided information on tubesheet bore dilation, the diametral growth of the tubesheet bore represented by the average of the maximum and minimum diameters of the eccentric tubesheet bore. Examination of the results from this model, as is discussed further below, resulted in two significant conclusions:
- 1. For the tubesheet bore eccentricities and dilation due to the applied loading in the limiting plants in the models of SG considered, the tube remains in contact with the tubesheet bore.
.2. While tubesheet bore eccentricity contributes to the reduction in contact pressure between the tube and the tubesheet, tubesheet bore dilation appears to be the principal cause of reduction of contact pressure between the tube and the tubesheet.
1.1 Discussion of 3D FEA Model for H* Analysis Provide table showing maximum delta diameters (total diameter distortion) and maximum eccentricities (maximum diameter minus minimum diameter) from the 3 dimensional (3-D) finite element analysis for normal operating and steam line break (SLB), for model F and D5.
Response
The 3D FEA Model and its application for determining the tubesheet displacements are extensively described in Section 6 of the H* WCAP reports (References 8 and 9). It is important to note that the 3D FEA model includes the entire tubesheet complex (i.e., tubesheet, stub barrel, channelhead and divider plate) but excludes the tubes. The model utilizes an equivalent material approach from Reference 6-5 in the WCAP reports to represent the deformation of the tubesheet under the applied loading conditions (NOP, SLB/FLB). Displacements in Cartesian coordinates are calculated for these conditions at any location on the tubesheet. The displacements calculated are the changes from an unstressed, room temperature condition, after all thermal and pressure loads appropriate to the operating conditions are applied. Application of a uniform temperature increase causes uniform dilation at each tubesheet bore. Application of pressure loads causes distortions in the structure due to bending. The 3D FEA model provides integrated total displacements of each tubesheet bore location.
Table RAI4-1 is a summary of the maximum eccentricities and ADs for the Model F and Model D5 limiting plants as calculated based on the UR (tubesheet radial displacement) results from the 3-D lower SG complex model.
10
LTR-SGMP-09-109 NP-Attachment Table RAI4-1: Summary of Model D5 and Model F NOP and SLB Eccentricity Results SG Model Elev. Avg. Eccentricity Data Max. Eccentricity Data Avg. A D Max. A D Above BTS"11 NOP SLB NOP SLB NOP SLB NOP SLB
,c,e
- in in/in in/in in/in in/in in in in in F
F F
D5 D5 D5 F
D5 Eccentricity, e A D,02 A D, 902 Notes:
- 1. BTS is Bottom of the Tubesheet Plant Condition Value inch/inch inch inch a,c,e Byron SLB MAX Byron SLB MIN Byron SLB AVG The original Table RAI4-4 is Millstone SLB MAX provided here for convenience Millstone SLB MIN Millstone SLB AVG Byron NOP MAX Byron NOP MIN Byron NOP AVG Millstone NOP MAX Millstone NOP MIN Millstone NOP AVG 11
LTR-SGMP-09-109 NP-Attachment 1.2 Discussion of the "Slice" Model Clarify the nature of the finite element model ("slice" model versus axisymmetric SG assembly model) used to generate the specific information in Tables 6-1, 2, and 3 (and accompanying graph entitled "EllipticalHole Factors")of Reference 6-15. What loads were applied? How was the eccentricityproduced in the model? (By modeling the eccentricityas part of the geometry? By applying an axisymmetric pressure the inside of the bore?) Explain why this model is not scalable to lower temperatures.
IThe difference between initial and final eccentricityincluded in Table RAI4-2 needs to be explained. In particular,the exclusive use of the relationshipbetween initialeccentricity and scale factor in calculatingcontactpressureneeds to be justified.
Resolve the apparentinconsistency between Item 4 on page 25 and the statement below Figure RA14-1 regardinghow the model in Figure RA14-1 is loaded.
Response
The "slice model" is shown in Figure 6-9 of Reference 6-15 in WCAP-17071-P, WCAP-17072-P, WCAP-17091-P, and WCAP-17092-P.
The data in Tables 6-1, 6-2, and 6-3 of Reference 6-15 of the H* WCAP reports, are derived from this plane stress model ("slice model") developed in WECAN/PLUS and the contact pressure equation identified on page 6-87 of WCAP-17071-P, page 6-95 of WCAP-17072-P, page 6-91 of WCAP-17091-P and page 6-84 of WCAP-17092-P as described below.
For convenience Tables 6-1, 6-2, and 6-3 of Reference 6-15 are replicated below and re-named as follows:. Table 6-1 is renamed as Table RAI4-2, Table 6-2 is renamed as Table RAI 4-3, and Table 6-3 is renamed as Table RAI4-4.
The "initial" eccentricities (defined as DMAX - DMIN) applied in the "slice" model in Table RAI4-3 and Table RAI4-4 are directly incorporated into the model geometry. That is, the initial eccentricity is built into the model geometry. The eccentricity values in the model were assumed values for tubesheet tube bore deformation based on engineering judgment and prior experience.
In the "slice" model analysis, the tubesheet is assumed to have a thermal expansion coefficient of zero (0) in/in/"F and the tube material is assumed to have the appropriate ASME Code thermal expansion coefficient values. (The TS coefficient of thermal expansion is set to zero to provide a loading mechanism for the model. When a temperature is applied, the tube "grows" into the tubesheet collar. The temperature difference applied to the tube in the "slice" model was 500'F, for a total tube temperature of 570°F. [Applied 500'F + 70'F assumed room temperature]). The sole purpose of the development of the "slice" model was to provide a sensitivity study to relate the effects of assumed eccentricity (DMAx - DMIN) conditions to contact pressures from which the contact pressure ratios were developed. No attempt was made to 12
LTR-SGMP-09-109 NP-Attachment reproduce the contact pressures that would be calculated by the 2-D axisymmetric model that was previously used to develop the tubesheet displacements.
The "final" eccentricity (DMAx - DMIN) values in Table RAI4-3 and Table RAI4-4 were also determined using the "slice model": The final eccentricity values are the (DMAX - DMIN) results of applying the loading conditions on the slice model: The loads applied to the "slice" model were thermal loads only as follows:
0 psig - Primary Side Pressure 0 psig - Secondary Side Pressure 500 OF- Tubesheet AT 500 OF- Channel Head AT 500 OF- Shell AT As discussed in Reference 6-15, Table RAI4-3 was constructed using the displacement results from the plane stress model analysis for the elliptical holes along with the contact pressure equations. The effective change in hole diameter was calculated as follows using a series of assumed scale factors:
a,c,e The ADMAx and ADMIN were taken from the radial and circumfential change in tube bore diameter in the "slice" model.
The corresponding contact pressure for each scale factor was then determined as follows:
a,c,e (RA14-2)
Equation RAI4-2 is a generic representation of how tube to tubesheet contact pressure is calculated in the H* integrator spreadsheet analysis. The equation is equivalent to the equation for P2 shown on page 6-87 in WCAP-17071-P, page 6-95 in WCAP-17072-P, page 6-91 in WCAP-1 0791-P and page 6-84 in WCAP-1 7092-P.
The scale factors for a given input eccentricity in Table RAI 4-3 result in contact pressure ratios using the thick shell equations that are equal to the contact pressure ratios calculated using the "slice" model for initial eccentricities (defined as DMAX - DMIN) equivalent to 0.0002, 0.0004, 0.0006 and 0.0008 inches, respectively, compared to the contact pressures for a circular hole (DMAx - DMIN =0). These scale factors are identified in bold print in Table RAI4-3. The data for the scale factors as a function of "initial" eccentricity was fit by a third order polynomial equation provided on page 6-85 of WCAP-1 7071-P and page 6-86 of WCAP-1 7072-P.
13
LTR-SGMP-09-109 NP-Attachment Based on a review of Table RAI4-3 and Table RAI4-4, the scale factor [ ]ace is the appropriate scale factor for calculating a reduction factor for contact pressure of [ ] a,c,e associated with an initial eccentricity of [ I a,c,e ((DMAx - DMIN)/ [ ]ac.e inch) from the "slice" model. The scale factor of [ ]ac,e relates to a contact pressure reduction factor of Ia..re and corresponds to an initial eccentricity of [ ]ace inch, and so forth.
The "final eccentricity" values corresponding to the same scale factors highlighted in bold in Table RAI 4-3 (and Table RAI4-4) are not used in determining the reduction in contact pressure because the resulting third order polynomial relationship between scale factor and eccentricity is bounded by the relationship for "initial eccentricity", i.e., the resultant scale factors, and hence the reduction in contact pressure due to eccentricity, would be less using the third order fit resulting from the "final" eccentricity values from Table RAI 4-3. For example, for an eccentricity of 1 E-3 in/in, the scale factor is [ ] 5'c,' as compared to [ ]a,c,e for the trend line associated with the "initial" eccentricity results. Figure RAI 4-1 illustrates this. This figure shows a comparison of the trend line analysis for "initial" eccentricity and "final" eccentricity.
Referring to Equation RAI 4-1, larger scale factors result in a greater reduction in contact pressure due to eccentricity.
Table RAI4-2 Reproduced Table 6-1 of Reference 6-15 Sleeve O.D. Tube O.D.
1 DeltaF1MT Eccentricity 2r (ineS AwarLn a0) i_(3) J lfA(1) (2) Awnrn (1) - pi_,(3) ace 0.0000 0.0002 0.0004 0.0006 0.0008 Notes: This table-is developed from the model shown in Figure RAI4-1, below.
- 1. The units of these columns are stress in psi.
- 2. The "delta" in this table refers to the maximum deviation from a constant value of the mean linearized radial stress around the tube bore.
- 3. The ratio is calculated by dividing the contact pressure between the tube and the tubesheet at a given eccentricity by the contact pressure between the tube and the tubesheet in a round tube bore (e=0.0). For example, the ratio of [ Ic e calculated in Table 6-1 is a ratio of the average contact pressure at an eccentricity of 0.0002 in of[ ]aG.epsi divided by the average contact pressure at an eccentricity of f 1.... psi.
14
LTR-SGMP-09-109 NP-Attachment Table RAI4-3 Reproduction of Table 6-2 of Reference 6-15 Primary Pressure 0 psig Secondary Pressure 0 p'Sig Tubesheet Delta T 500 oF Shell Delta T 500 F Channel Head Delta T 500 oF Sleeve OD Delta D [ ]..e in Tube ID Delta D [ iin Tube OD Delta D (Thermal) [ itn1~
Sleeve/Tube Interaction Coefficients [ ],¢ Tube/Tubesheet Interaction Coefficients [ ].
Eccentricity (1) (2) (3) (4) (5) (6)
I Final Hole Delta Hole Delta S/T T/TS Initial (inch) Max/Min D (0 Deg) D (90 Contract Contact Ratio (inch),c,e o n Deg) Pressure Pressure 0.0000 Minimum a,c,e Average Maximum 0.0002 Minimum Averae a,e Maximum 0.0004 Minimum Avera e a, LMaximum 15
LTR-SGMP-09-109 NP-Attachment Table RAI4-3 (Cont'd.)
Eccentricity a,c,e (1) (2) (3) (4) (5) (6) a,c,e 0.0006 F Minimum Averageace
-Maximum__-_ _
0.0008 Minimum Average ac Maximum Note: The values in Bol-d identify the source data for Table RAI4-3 Table RAI4-4 Reproduction of Table 6-3 of Reference 6-15 Nominal Hole Diameter I a,c,e 11 Fcce.ntricitv )
Initial Initial Pressure Delta (in/in) Ratio Dia Final Max/Min (in) (in/in) Factor a,c,e I
0.0000 0.0002 0.0004 0.0006 0.0008 tI) I nese values are the values Tor initiai and Tinal eccentricity from Table RAI4-2 are divided by the nominal tubesheet hole diameter [ ]a,c,e I
]a,c,e 16
LTR-SGMP-09-109 NP-Attachment a,c,e Figure RAI4-2: Scale Factor Comparison (Initial versus Final Eccentricity)
The method for calculating the contact pressure for using the "old" method for the Model F SGs (all plant conditions) and the Model D5 SGs (NOP and FLB conditions) and the "new" method for calculating the contact pressure the Model D5 SGs only (SLB conditions) are described below:
Old Method (Reference 6-15):
- 1. The UR used in the calculation of the circumferential and radial AD is based on the linearly scaled 2D axisymmetric FEA model (3-D model for the current H* analysis) of the lower SG complex
- 2. The circumferential and radial AD's are used in the scale factor (SF) equation to determine the ADhole (see equation RAI4-1) that is used to determine the reduction in contact pressure as a function of eccentricity (e), equation RAI4-2.
- 3. The relationship between 11D and e is based on the 2-D plane model shown in Figure 6-9 of SM-94-58, Rev.1.
17
LTR-SGMP-09-109 NP-Attachment .
- 4. The model in Figure 6-9 of SM-94-58, Rev.1 includes the initial applied eccentricities (DMAx -
DMIN) geometry definition of the model.
- 5. The "slice" model provides the input for using the SF relationship (Eqn. RAI4-1). The SF is determined by comparing the "slice" model results to the axisymmetric model results for a TS collar and tube model at a given radius in the TS over the full thickness of the TS.
- 6. The result is then used to calculate the reduction in contact pressure as a function of TS elevation and radius due to TS displacement and tube bore eccentricity. This is appropriate because the conditions for the Model F SG and Model D5 SG (NOP and FLB conditions) are within the range of data for which the scale factor relationship is applicable.
New Method (WCAP-1 7071 -P, WCAP-1 7072-P):
- 1. The UR used in the calculation of the circumferential and radial AD comes from a 3-D FEA model of the lower SG complex with condition-specific inputs applied.
- 2. The circumferential and radial AD's are compared to determine the maximum AD that will give the maximum reduction in contact pressure as a function of eccentricity (e).
- 3. The relationship between AD and e is based on the 2-D [ face model shown in WCAP-17071-P and WCAP-17072-P, section 6.2.5. The model is shown in Figure 6-49 of the WCAP reports. The range of eccentricity used in this study conservatively exceeds the values of tube bore eccentricity calculated from the perforated TS model in Section 6.2.4.
- 4. The model in Figure 6-49 of the H* WCAP reports applies boundary conditions to the outer edge of the tube pitch material and does not directly affect the material that is deforming in the tube and tubesheet cell.
- 5. The TS deformations and tube to tubesheet contact pressure results that produce the maximum reduction in contact pressure at the minimum value of TS tube bore eccentricity are then fit with a linear relationship.
- 6. The result of the linear relationship is used to determine the reduction in contact pressure between the tube and the tubesheet directly. There are no intermediate equations or results.
A correct prediction of contact pressure loss requires the knowledge of both the proper values of DMAX and DMIN associated with the different pressure and temperature conditions at a given tubesheet radius and elevation as well as the value of eccentricity. The values of DMAX and DMIN are a function of the radial deflection of the tubesheet, UR, as determined by the finite element analysis model (which previously was a 2-D axisymmetric model of the SG lower 18
LTR-SGMP-09-109 NP-Attachment assembly and at present, is a 3-D model of the SG lower assembly). The results from the
'slice" model cannot be linearly scaled to lower temperatures because the method of super-position has been shown during the development of the current H* analysis to not apply to the non-linear combination of materials and loading in the lower SG complex. This conclusion led to the development of the 3D FEA model that is the reference model for the H* analysis. A discussion of this is provided in Section 6.1.2 of WCAP-17071-P and WCAP-17072-P.
1.3 Discussion of the Unit Cell Model to Calculate Contact Pressures The "Unit Cell" model is extensively discussed in Section 6.2.5 of the H* WCAPs (References 8 and 9). The specific goal of this model was to determine if tube to tube contact would remain when the tubesheet is deformed due to operating loads. An equivalent tubesheet cell is modeled, that is, a tubesheet bore with surrounding tubesheet material, and a tube in the tubesheet bore (see Figure 6-48 of the H* WCAPs). For the primary purpose of this model - to study if tube-to-tubesheet contact is present during the limiting tubesheet deformations - the model was initialized by simulating the tube expansion process. The expansion process was conservatively simulated by applying a low value of expansion pressure [ Ia,c,e inside the tube, resulting in initial tube to tubesheet contact, and then removing the tube expansion internal pressure. The calculated dilation of the tubesheet bore due to the simulation of the tube expansion is [ I a.c.e inch for all models of SG considered.
As discussed in Section 6.2.5 of the H* WCAP reports, the operating pressure loads, were applied to the initialized model in a sequential manner, and the resulting contact pressures were calculated when a range of displacements (termed "E-bar") were applied as boundary conditions to the model. Figure RAI4-2 shows the updated sequential loading (includes application of thermal loads) of the model and relates it to the steps discussed in Section 6.2.5 and Tables 6-18 and 6-19 of the H* WCAPs. The "E-bar" values shown as the displacement inputs on Tables 6-18 and 6-19 in the H* WCAP reports are uni-directional displacements (in inches) that are NOT the same as eccentricity and also not the same as AD. (Eccentricity is defined as the difference between the maximum and minimum diameters of a bore divided by the nominal diameter of the bore. The units of eccentricity are inch/inch.) The displacement inputs applied to the unit cell model are assumed values that based on prior analyses that envelope the expected tubesheet displacement for all of the applicable operating conditions. It is important to note that the unit cell model as described in Section 6.2.5 of the H* WCAP reports utilizes boundary conditions chosen to minimize the tube-to-tubesheet contact pressures for the applied relative displacements.
To interpret the results from the unit cell model properly, the following must be observed:
- To address if tube to tubesheet contact continues for all the assumed tubesheet displacements, the appropriate reference condition is the initialized condition (after Step 4) of the model that simulates a tube expanded in the tubesheet bore.
19
LTR-SGMP-09-109 NP-Attachment
- To compare the results of the unit cell model with the 3D FEA model, the appropriate reference condition of the unit cell model is the initial model (Step 0) without the tube expansion simulated and thermal loads must be included.
Figures RAI4-3 and RAI4-4 show the average tubesheet bore dilation (AD) as a function of tubesheet relative displacement (E-bar) for the Model F and Model D5. The average tube bore dilation at zero E-bar input is the result of the temperature and pressure loading of the unit cell model. Initially, application of the displacement input "E-bar" results in more significant hole dilation, but rapidly takes on a shallower slope as the applied displacement increases. The curves are characteristically the same for the Model F and Model D5 steam generators and also for the different operating conditions, NOP and SLB, for the different models of SGs.
Similarly, Figures RAI4-5 and RAI4-6 show the tubesheet bore eccentricity "e" as a function of tubesheet relative displacement (E-bar) for the Model F and Model D5. Eccentricity initially increases with application of the displacement boundary condition (E-bar) simulating the load due to pressure differential across the tubesheet, but the rate of increase decays with increasing E-bar. A significant difference is noted between NOP and SLB conditions at large values of E-bar. This difference reflects the fact that the uniform growth of the tube bore hole due to increased temperature overwhelms the effect of application of the displacement boundary condition (E-bar) on tubesheet bore eccentricity. During the SLB event, the temperature is decreased and the differences in DMAX and DMIN remain more significant as the displacement boundary condition is increased, although the rate of increase in the difference between DMAX and DMIN is reduced at some point.. Eventually, at NOP conditions, the difference between DMAX and DMIN tends to become decrease even though a greater displacement (E-bar) is applied, leading to a reduction of eccentricity "e."
Figures RAI4-7 and RAI4-8 show the contact pressure as a function of tubesheet relative displacement (E-bar) for the Model F and Model D5 for both NOP and SLB conditions based on the unit cell model. As expected, both NOP and SLB contact pressure decrease with increasing displacement inputs, ultimately going to zero at a very large value of applied displacements. It is to be noted that the maximum displacement assumed is significantly greater than would be predicted by the 3D FEA model. Over the entire range of assumed displacement conditions, the SLB contact pressure exceeds that for NOP conditions.
Table RAI4-5 summarizes the eccentricity, AD and predicted contact pressure using the unit cell model for various values of applied displacement (E-bar) for both the model F and Model D5 SGs. The true eccentricity ([Dmax-Dmin]/Dnom) is shown for the applied displacement, E-bar.
Table RAI4-5 also provides a comparison of the AD predicted by the unit cell model for the two reference conditions noted above, that is, for the total AD from the model without the simulated tube expansion (reference step 0 in Table 6-18) and for the initialized case with the tube expansion simulated (reference step 4 in Table 6-18).
20
LTR-SGMP-09-109 NP-Attachment Further, Table RAI4-5 provides a summary of contact pressures between the tube and the tubesheet for various applied values of E-bar for the Model F and Model D5 SGs. The "Modified Contact Pressure" is the "Raw Contact Pressure" from the unit cell model adjusted for the actual tube expansion process ([ ]a,,,e psi compared to the simulation at [ ] a,c,e psi) real Model F and Model D5 geometry and more realistic operating conditions of pressures and temperatures. For all cases of applied displacement, positive contact pressure remains between the tube and tubesheet. It should be noted that the largest value of applied displacement (E-bar) is well in excess of the displacement predicted by the 3D FEA model.
Table RAI4-6 provides similar data to that in Table RAI4-5, except that the data is based on the 3D FEA model.
Comparison of Tables RAI4-5 and RAI4-6 leads to the following observations:
1.The ADs from the 3D FEA model are significantly less than the corresponding ADs from the unit cell model from the unloaded to the fully loaded condition (i.e., from step 0 to step 9) for both NOP and SLB conditions. This leads to the conclusion that the unit cell model displacement results and contact pressure predictions conservatively represent the reference 3D FEA model results.
- 2. The eccentricities from the unit cell model are generally comparable to those from the 3D FEA model. A more exact comparison is difficult based on the available data; however, it is clear that the actual range of eccentricities from the 3D FEA model was adequately addressed by the unit cell model.
- 3. The method of Reference 6-15 of the H* WCAP report for adjusting contact pressure provides acceptable results for all conditions except the SLB condition for the Model D5 SGs The method of Reference 6-15 significantly under-predicts contact pressure for the Model D5 SLB conditions. Referring to Figure RAI4-6, the method for calculating the reduction in contact pressure defined by the White Paper, when adjusted for temperature effects, shows that SLB contact pressure is increased relative to normal operating conditions.
21
LTR-SGMP-09-109 NP-Attachment Table RAI4-5 Eccentricity, Contact Pressure and AD Results from Unit Cell Model Square Cell Results Square Cell Results i Square Cell Results Square Cell - Average Delta D Raw Contact Modified Contact 2 Fnn-, ntrirnifv P re_.ssure.(1) pressure(1) Sten 0( ) - Sten 9(3 ) Steo 4(4) - Steo 9(3)
SG "E bar" NOP SLB NOP SLB NOP SLB NOP SLB NOP SLB a,c,e Model in in/in in/in psi psi psi psi in in in in F
F F
F D5 D5 D5 D5 Notes:
- 1. Accounts for expansion pressure and geometry.
- 2. See Section 6.2.5 H*WCAP. Step 0 is the condition of the unit-cell model prior to any modifications for tube expansion, loading, etc.
- 3. See Section 6.2.5 H*WCAP. Step 9 is the condition of the unit cell model after all loading conditions have been applied.
- 4. See Section 6.2.5 H*WCAP. Step 4 is the initialized condition of the Unit Cell model after tube expansion has been simulated.
22
LTR-SGMP-09-109 NP-Attachment Table RAI4-6 Eccentricity, Bore Dilation and Contact Pressure from 3D FEA Model Hstar Analysis I Hstar Analysis Hstar Analysis - Avg. A D Eccentricity NOP SLB NOP SLB NOP SLB SG Model and Contact Pressure a,c,e Reduction Model in/in in/in psi psi in in F - Ref. 6-15 Limiting Radius - F - Ref. 6-15 D5 - Ref. 6-15 D5 - White Paper Limiting Radius - D5 - Ref. 6-15 Limiting Radius - D5 - White Paper F - Updated Model (
D5 - Updated Model (1) i): upoate ivMoael Results based on estimates Trom approximate values in Tinite element analysis ana do not reflect ine result of a regression analysis.
23
LTR-SGMP-09-109 NP-Attachment Figure RAI4-2 Unit Cell Model and Loading Sequence (1) (2) (3)
Loading Steps:
- 0. Initial Model
- 1. Initial Tube to TS gap
- 2. Pressurize tube to 16ksi
- 3. Pressurize tube to 28ksi
- 4. Release Pressure on Tube
- 5. Apply AT(1 )
- 6. Apply "E-bar"
- 7. Apply AP=[ ]a~c,e psi
- 8. Apply AP=[ ]a,c,e psi
- 9. Apply AP=[ ]a.ce psi Step 0 Step 5 Step 6-9 Notes: (1) The application of the unit AP= 0 tAP= 0 AP> 0 cell model in support of Tables 6-18 AT>>0 and 6-19 does not include application AT=0 AT=0 of AT.
Unexpanded Tube Expanded Tube Expanded Tube e=0 e=0 e>0 e bar = 0 ebar=0 e bar> 0 24
LTR-SGMP-09-109 NP-Attachment a,c,e Figure RAI4-3 Relationship between "E-bar" and AD; Model F a,c,e Figure RAI4-4 Relationship between "E-bar" and AD; Model D5 25
LTR-SGMP-09-109 NP-Attachment a,ce Figure RAI4-5 Relationship between "E-bar" and Eccentricity "e"; Model F a,c,e Figure RAI4-6 Relationship between "E-bar" and Eccentricity "e"; Model D5 26
LTR-SGMP-09-109 NP-Attachment a,c,e F
Figure RAI4-7 Relationship between "E-bar" and Contact Pressure; Model F a,c,e Figure RAI4-8 Relationship between "E-bar" and Contact Pressure; Model D5 27
LTR-SGMP-09-109 NP-Attachment 2.0 Comparison of Slice Model and Unit Cell Model Results The basis for applying the correlationfor scale factor outside an "eccentricity"range of between 1E-3 to 1E-4 inch in the calculationof contact pressureneeds to be further explained. Values for displacements included in Table 6-18 (of WCAP-17071) suggest that contact pressure may be lost at displacementrangingbetween 1E-3 in to 1E-4 in.
Response
Interpretation of the displacements noted in Table 6-18 of the WCAP reports was clarified in the prior response, Section 1.3. The values noted in the column titled "Displacement Total" refer to the condition of the unit cell model after Step 4 of the loading sequence (See Figure RAI4-2).
When the true reference condition (Step 0) for total displacement is considered, the values of total displacement are significantly larger as noted previously.
Westinghouse agrees that the derivation of the fit in Reference 6-15 is non-intuitive and limited in its application. However, the results of applying the fit described in reference 6-15 are acceptable relative to a best case finite element model (unit cell with thermal and AP loading) for the reasons described below.
Westinghouse also agrees that the fit that describes the reduction in contact pressure for the steam line break condition in the Model D5 White Paper does not account for the reduction in contact pressure due to tube bore dilation in the same manner as the fit described in Reference 6-15. The results of using the fit described in Reference 6-15 also match the expected trend from a best case finite element model. See the response to b.4 below for more details.
A series of tubesheet tube bore eccentricities were applied to the tubesheet cell model and combined with different pressure and temperature loads. The average, maximum and minimum values of the tube-to-tubesheet (T/TS) contact pressures around the circumference of the tube were reported. The values of tubesheet relative displacement, pressure and temperature that were used in the analysis are summarized in the table below.
Input Conditions for Unit Cell Model (no correlation implied)
Internal Temperature
- Pressure Difference in AP, psi AT, -F ace 0.00 2.OE-04 4.OE-04 28
LTR-SGMP-09-109 NP-Attachment Normal operating (NOP) conditions in the Model D5 and Model F steam generators are represented by a AP of [ ]a,,,e psi and a AT of [ ] ace OF. Main steam line break (SLB) conditions in the Model D5 are represented by a AP of [ ] a,c,e psi and a AT of [ ] a,c,e OF.
The value of AP in the tubesheet cell can change as a function of elevation in the tubesheet due to the distribution of crevice pressure. The results of the study include the data for a depth ratio of 0.9 which is an elevation roughly 2 inches below the top of the tubesheet. The values of AP represented in this study account for the region of interest near the top of the tubesheet where the maximum eccentricity in the tubesheet is expected and where the crevice fluid is transitioning from the crevice conditions to the secondary side fluid conditions. The region roughly 2 inches below the top of the tubesheet is also where a significant portion of the T/TS contact pressure develops so it is a good indicator of trends in the effect that different operating conditions have on the contact pressure.
The original results in section 6.2.5 of WCAP-1 7071-P were used to verify that the reduction in T/TS contact pressure as a function of tubesheet tube bore eccentricity was appropriate for the Model F SG. The original relationship that is used to define the reduction in T/TS contact pressure as a function of eccentricity is described in section 6.3 of WCAP 17071-P and WCAP 17072-P. However, the result of applying the fit described in section 6.3 to the Model D5 SG during SLB was shown to be inconsistent with the expected trend from the more detailed analysis described in section 6.2.5. The results of section 6.2.5 were then used to define a new relationship between the reduction in T/TS contact pressure and tube bore eccentricity. This new relationship is described in the Model D5 White Paper (Reference 12). Figure RAI4-8 shows the result of applying the new relationship to the Model D5 SLB conditions (i.e., White Paper results, Reference 12) in comparison with the results from the old 3 rd order polynomial relationship. Because the tubesheet temperature induced hole dilation, potentially the most significant factor in contact pressure reduction, was not considered in the Model D5 condition results, a third curve was added to the figure titled "Model D5 FEA trend." This curve represents the most accurate calculation of the contact pressure ratio.
Figure RAI4-9 shows the contact pressure ratio (PCSLB/PCNoP) as a function of tubesheet relative displacement, E-bar. It is clear from Figure RAI4-9 that the results of using the old fit for the Model D5 SLB are inconsistent with the more detailed analysis. At SLB conditions, the tubesheet bore dilation is relatively larger than at NOP conditions due to the increased bending of the TS and decreased thermal expansion. Therefore, it is expected that the TITS contact pressure ratio should increase by a factor of at least [ ]a,c,e (see Figure RAI4-9) when going from NOP to SLB. It is also expected that the tube to tubesheet contact pressure should decrease with increasing tube bore eccentricity. The H* results using the old fit for the Model D5 clearly do not follow either expectation from the detailed analysis. However, when the new fit results are applied to the H* calculation process the relationship between T/TS contact pressure in the Model D5 is much more reasonable and follows the expected trend from the more detailed analysis.
The Model F H* contact pressure results, using the old fit, are well within the range predicted by the more detailed analysis in section 6.2.5 and the additional work described in this RAI response. See Figure RAI 4-10 below. This means that the old fit is appropriate to use for the 29
LTR-SGMP-09-109 NP-Attachment Model F NOP and SLB conditions and the NOP condition in the Model D5 SG. The results of using the fit described in Reference 6-15 match the expected trend from a best case finite element model for the NOP and SLB conditions for the Model F SGs and NOP conditions for the Model D5 SG.
To further address the concern that contact pressure may be lost at displacements ranging between 1E-3 in and 1E-4 in, the "Unit Cell" model is extensively discussed in Section 1.3 of this response above.
30
LTR-SGMP-09-109 NP-Attachment a,c,e Figure RAI 4-9 a,c,e Figure RAI 4-10 31
LTR-SGMP-09-109 NP-Attachment 2.1 Calculation Basis for Contact Pressure Reduction Factors Provide the calculationbasis for the upper and lower curves provided in Figure RAI 4-2 I
Response
The original figure RAI4-2, referred to in the question, is reproduced here as RAI4-10 to provide the foundation for the question and the response. Note that the scale of the y-axis has been corrected as discussed in the meeting on August 17, 2009.
a,c,e Figure RAI4-10 (original Figure RAI4-2)
The upper curve in the figure above is based on the data from the following table:
32
LTR-SGMP-09-109 NP-Attachment Eccentricity (ADmax- Reduction in Contact Normalization Contact Pressure ADmin) (in) Pressure (psi)(1) Basis Reduction (psi) Factor(psi/psi) 0 0 0 2E-4 ]ace psi 1200 ]ac,e 4E-4[ a,c,e psi[ ]ace 5E-4_[ ]a,,e psi[ ]ace 6E-4 [ ] a,c,e [ ]a,c,e Notes: (1) Contact stress reductions are based on the values on Table RAI4-3 Referring to Table RAI 4-3, the contact pressure for a round tube bore hole is calculated to be I I ",ce psi (Ratio = 1.0). The contact pressure for a tube bore hole that results in a contact pressure ratio reduction of [ ]a,c,e (Ratio = [ ],.c,e), which corresponds to an eccentricity of 2E-4 inch, is [ ]a,c,e psi. The absolute reduction in contact pressure is
[ ] a,.,e psi.
The total reduction in contact pressure using the new model is approximately [ ] a,c,e psi (see Figure 6-69 of WCAP-17072-P). To plot the absolute reduction in contact pressure of
] ace psi for an eccentricity of 2E-4 on Figure RAl4-10, the value is normalized by the total reduction in contact pressure of [ ] a,c, psi from the new method. This value represents a reduction in contact pressure of [ ]a~c~e Again, referring to Table RAI 4-3, the contact pressure for a round tube hole is calculated to be I ] a,,,e psi. The contact pressure for a tube bore hole that results in a contact pressure ratio reduction of [ ]a,c,e (Ratio = [ ] ace), which corresponds to be eccentricity of 4E-4 inch, is [ a,,e psi. The absolute reduction in contact pressure is [ a,c,e psi.
Again, the total reduction in contact pressure using the new model is approximately [ ] a,c,e psi (see Figure 6-69 of WCAP-17072-P). To plot the absolute reduction in contact pressure of I I a,c,e psi for an eccentricity of 4E-4 on Figure RAI4-10, the value is normalized by the total reduction in contact pressure of [ ],c, psi from the new method. This value represents a reduction in contact pressure of [ ] a*,e.
The same calculation was completed for an eccentricity of 6E-4 in. The value for 5E-4 in is an interpolated value between 4 E-4 in and 6E-4 in.
The bottom curve in the figure above is generated using the 3 rd order polynomial fit. The results are summarized in the following table:
33
LTR-SGMP-09-109 NP-Attachment E, eccentricity (in) T/TS Contact Normalized Contact Pressure Reduction Pressure Reduction (psi) a,c,e a,c,e 6.36E-07 a,c,e ac,e 5.53E-05 a.c,e a,c,e 3.16E-04 a,c,e a,c,e 5.69E-04 a,c,e a,c,e 9.07E-04 2.2 Conservatism of 3 rd Order Polynomial Fit from WCAP Reference 6-15 Provide the basis for why the ADhoWe adjustment for contact pressure made using the old model remains conservative. I
Response
The key conclusions from the comparison of the Reference 6-15 analysis, the WCAP results and the results of the square cell tubesheet model are:
1.) The fit described in Reference 6-15 of the H* WCAP reports is conservative when applied to the NOP condition in both the Model D5 and Model F SG. The fit tends to under-estimate the contact pressure during NOP by as much as [ ]a,c,e psi to [ ]ace psi) for the Model F SG and as much as [ ]acoe % for the Model D5 SG ([ ]ace psi to
]a,,c.e psi) (see Table RAI4-6).
2.) The fit described in Reference 6-15 of the H* WCAP reports is comparable when applied to the SLB condition in the Model F SG. The fit described in the Model D5 White Paper tends to over-estimate the contact pressure, by as much as [ ]ac,,e %, during SLB
([ ]aec'e psi to [ ]ac,, psi) because the White Paper does not fully account for the change in tube bore diameter during the transient.
3.) The fit described in Reference 6-15 of the H* WCAP reports significantly under-estimates the contact pressure, by as much as [ ]a'c'e %, during the D5 SLB condition (from
]a,c,e psi to [ ]ac.e psi).
4.) The square cell tubesheet finite element model shows an increase in contact pressure when going from NOP to SLB conditions in both the D5 and F SGs.
34
LTR-SGMP-09-109 NP-Attachment 5.) Using the results from the square cell model to estimate the magnitude of the contact pressure reduction from the change in tube bore diameter calculated using the 3D finite element results from the lower SG tubesheet complex model show that the contact pressure still increases when going from NOP to SLB conditions in both the Model F and Model D5 SG.
The results of this analysis show that NOP contact pressures that define H* in the Model F and Model D5 SG are conservative and that a more realistic model of contact pressure reduction as a function of tube bore deformation (including both dilation and eccentricity) would predict an increase in tube to tubesheet contact pressure at SLB conditions compared to NOP conditions.
(See also Section 2.3) 2.3 SLB vs. NOP Contact Pressures Provide an appropriatebasis for demonstratingthatjoints tighten duringa postulated SLB event. Why is it acceptableto compare the contact pressures calculatedusing the originalmodel for NOP to the contactpressures calculatedusing the new model for SLB for the Model D5 SGs?
Response
Table RAI4-5 provides a summary of contact pressures between the tube and the tubesheet for various applied values of E-bar for the Model F and Model D5 SGs. Comparison of the eccentricity values calculated using the unit cell model (see Table RAI4-5) with the eccentricity values calculated from the 3D FEA model (see Table RAI4-1) shows that the eccentricities from both models are comparable. It is not reasonable to expect exact matches of numbers between the two models, however, the order of magnitude of the calculated eccentricities is the same.
Given that the two structural models provide similar eccentricities, the unit cell model shows that for these eccentricities, positive contact pressure exists between the tubes and the tubesheet for the-entire range of displacements considered. Further, the results show that the contact pressures at SLB conditions exceed those at NOP conditions (See Table RAI4-6). See also the discussion in Section 2.4 below.
35
LTR-SGMP-09-109 NP-Attachment 2.4 Relative Conservatism of "Old" and "New" Fit If both old and new models are conservative, is there an appropriatebasis to show the relative conservatism of the methods?
Response
As noted above in Section 1.3 of this response, tube bore dilation is a more significant factor in determining tube-to-tubesheet contact pressure at higher temperatures and the effect of eccentricity on contact pressure is reduced at higher temperatures. The methodology for addressing the effect of eccentricity on contact pressure discussed in Reference 6-15 and utilized in WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP-17092-P reflects this fact and it, therefore, provides acceptably accurate contact pressure results at higher temperatures (i.e., for all conditions except the "colder" SLB condition). This includes NOP, SLB (higher temperature, > 400 0 F, and FLB, where appropriate).
Also, as noted in Section 1.3 of this report, the effect of eccentricity on contact pressure loss is a more significant factor at the lower SLB temperatures for the Model D5 SG, but tube bore dilation due to temperature and pressure needs to be considered (which was not addressed in the "new" method, a.k.a the White Paper method discussed in WCAP-17072-P or 17091-P).
Moreover, the original 3rd order polynomial fit significantly over-predicts contact pressure loss during the "colder" Model D5 SLB transient (and Model 44F two loop plant SLB).
Therefore, a more detailed model for contact pressure during a postulated SLB was developed.
Referring to Table RAI4-6, it shows that contact pressure increases during a SLB event
([ ]a,c,e psi) relative to NOP ([ ] a.c,e psi) with primary and secondary side temperatures as low as 212°F when comparing contact pressures for NOP conditions for the unit cell to contact pressures for SLB for the unit cell.
Again, referring to Table RAI4-6, it has been shown when comparing contact pressures for NOP conditions for the unit cell to contact pressures for SLB for the unit cell for the Model F SG (higher temperature SLB conditions), that contact pressure increases during a postulated SLB (from [ ]a,c,e psi at NOP to [ ]a,c,e psi at SLB).
36
LTR-SGMP-09-109 NP-Attachment 3.0 Comparison of 3D FEA and Unit Cell Model Results The eccentricitiesincluded in Table RAI 4-4 appearlarger than anticipated. Need to confirm that positive contactpressure exists around the entire circumference of the tube and state this clearly in the response.
Response
Comparison of the eccentricity values calculated using the unit cell model (see Table RAI4-5) with the eccentricity values calculated from the 3D FEA model (see Table RAI4-1) shows that the eccentricities from both models are comparable. It is not reasonable to expect exact matches of numbers between the two models, however, the order of magnitude of the eccentricities calculated is the same. Given that the two structural models provide similar eccentricities, the unit cell model shows that for these eccentricities, positive contact pressure exists between the tubes and the tubesheet for the entire range of displacements considered.
Further, the results show that the contact pressures at SLB conditions exceed those at NOP conditions.
4.0 Additional Background Information For Key Questions and Issues RAI#4 evolved in several stages, each stage building on the prior stage. Reference 10 provided additional questions to augment those that were provided by Reference 5. Responses were prepared and were discussed in a telephone conference on August 11, 2009. During this telephone conference, additional questions were raised as identified in the introduction of this document. The following are responses that were provided in response to Reference 10 that were discussed in the August 11, 2009 telephone conference. They are historical in nature and are provided to complete the record of information provided in response to the NRC request for additional information.
4.1 Comparison of "Old and New" Relationship for Reduction in Contact Pressure and Eccentricity In Figure 2 of the White Paper,add a plot for originalrelationshipbetween reductions in contactpressure and eccentricity as given in Reference 6-15 in the graph accompanying Table 6-3. Explain why this original relationshipremains conservative in light of the new relationship. Explain the reasons for the differences between the curves.
In order to superimpose the results of the "old" and "new" analyses for reduction in contact pressure related to eccentricity, the data for the "old" method must be normalized in the same fashion that Figure 2 has been normalized. The plot of contact pressure reduction included in 37
LTR-SGMP-09-109 NP-Attachment Figure 2 of the White Paper represents the total reduction in contact pressure associated with a given eccentricity. The information from Table 6-3 represents the ratio of the contact pressure calculated at a given eccentricity divided by the contact pressure calculated for a tubesheet bore with no eccentricity. For the new analysis, the total reduction in contact pressure for the eccentricities (DMAx - DMIN) for a range of up to [ ]ac.e inch is determined to be [ ]ac,e psi. For the old analysis, the total reduction in contact pressure for eccentricities in the same range is calculated to be [ ]a,c~e psi. The normalization basis is the same for both curves on the figure.
Figure RAI4-11, showing the normalized results as discussed during the August 17, 2009 meeting, is provided below. (Figure RAI4-11 is the same as Figure RAI4-10 in Section 2.1 of this document, except that the values of the "Old Polynomial Results" have been corrected on Figure RAI4-10 by a factor of 2 as discussed in the August 17, 2009 meeting.) The curve labeled "Old" Model Results is based on the data from Table RAI4-3 (Table 6-2 of Reference 15 of the WCAP report). The curve labeled "New" Model reproduces Figure 2 in the White Paper (Reference 12). The curve labeled "D5 SLB Polynomial Fit" are the results when the eccentricity data and ADhole for the Model D5 SLB condition are applied directly to the polynomial fit, equation 6-8 in WCAP-17072-P and similar equation on page 6-85 in WCAP-17071-P. The latter curve is based on the maximum displacement conditions at the top of the tubesheet for the Model D5.
The curve labeled "Old Model Results" (top curve on Figure RAI4-11) is misleading relative to making an assessment of the conservatism of the new analysis method compared to the old analysis method. Unlike the new analysis method, which is only applied to the SLB case for the Model D5 SGs, the old analysis method has not been applied as a linear function as represented in the figure as the uppermost curve (solid squares). In reality, the old data fit (top curve on Figure RAI4-11), which is a 3rd order polynomial fit, when extrapolated significantly outside its supported data range (i.e., at temperatures either significantly above or below 500 0 F), provides physically unrealistic results as shown on Figure RAI4-11 (bottom curve, A-symbols). The Model D5 SLB condition puts the tubesheet at a nearly uniform temperature of less than 3001F, which is far outside of the range for which the eccentricity relationship was developed in Reference 6-15 in the WCAP reports.
The original relationship remains conservative because it predicts greater reduction of tube to tubesheet contact pressure than the new method for all operating conditions. However, the original relationship is only valid when ADmin and ADmax are within [ ]a,c,e % and eccentricity is within [ ]a.c.e inch to [ ]*,,,e inch range, (i.e., the basis of the original fit).
The maximum tube bore distortions occur at the top of the tubesheet. The results from applying the old fit for the relationship versus the new fit for the relationship for the Model D5 SLB tubesheet displacements and contact pressures are shown in Table RAI4-7. The tube-to-tubesheet (T/TS) contact pressure result due to thermal expansion of the tube and the pressure expansion of the tube including the effect of the crevice pressure distribution, is the same in the both the "old" and "new" cases in the Table RAI4-7.
38
LTR-SGMP-09-109 NP-Attachment Table RAI4-7 Summary of Model D5 SLB Contact Pressure Results for Different Eccentricity Fit Relationships Model D5 T/TS PCON Reduction T/TS PCON Condition Value Eccentricity Old New Old New a,ce SLB Avg SLB Max SLB Min The results in Section 6.2.4 of WCAP-17071-P and WCAP-17072-P show that the average expected tubesheet-tube-bore eccentricity is on the order of [ ]ac.e inch. The results in Table RAI4-7 show that the old method of calculating the reduction in contact pressure due to tubesheet-tube-bore eccentricity and change in diameter is conservative for larger values of eccentricity and AD (predicts greater decrease in contact pressure) than the new fit. However, it is inappropriate to use the old method at smaller values of eccentricity and AD because it provides physically impossible results (see Table RAI4-7). For example, the "old" method predicts a larger decrease in contact pressure for a smaller eccentricity on the order of 10-7 inch than for a larger eccentricity on the order of 103 inch. The "new" method, by comparison, predicts a slightly positive increase in contact pressure for an eccentricity of 10-7 inch and a large reduction in contact pressure for an eccentricity of 10-4 inch or greater, a physically realistic result. The reason that the "old" method predicts such a different reduction in contact pressure for small values of eccentricity is that these small eccentricity values are well outside the range of the data upon which the "old" relationship was developed. However, when used within its intended range of eccentricities and tubesheet bore displacement, the "old" method provides valid and conservative results. The "new" method of calculating the reduction in T/TS contact pressure is linear and directly accommodates small calculated values of eccentricity. It is also clear from the results in Table RAI4-7 that the results from the old method when used in its supported eccentricity range are highly conservative compared to the "new" method.
39
LTR-SGMP-09-109 NP-Attachment a,c,e Figure RAI4-11 Original Figure RAI4-2 Discussed at the August 17, 2009 Meeting 40
LTR-SGMP-09-109 NP-Attachment 4.2 Use of Both "Old" and "New" Fit When establishing whether contact pressure increases when going from normal operating to steam line break conditions, how can a valid and conservative comparison be made if the normal operating case is based on the original delta contact pressure versus eccentricity curve and the SLB case is based on the new curve?
Response
It is important to note than the new analysis method is only used for the SLB condition for the Model D5 steam generators. Comparison of contact pressures between the normal operating condition and the SLB condition is made for the Model F steam generators in the H* fleet in WCAP-1 7071-P on a consistent basis.
It is Westinghouse's engineering judgment that the old methodology provides an accurate determination of contact pressures during normal operating conditions and postulated accident conditions (FLB and SLB) when peak temperatures range between [ ]a,c,e OF and eccentricities are between [ 1a,c,e inch and [ ]c 0 inch and Dmanx and Dmin are within I ]a~c,e % of each other.
Application of the new method to calculate eccentricities and values of Drax and Dmin that fall outside the above noted range provides conservative results because the plane strain model upon which it based over-estimates the stiffness of the tube and tubesheet structure leading to lower contact pressure results as a function of eccentricity. The new method also excluded the effect of temperature and therefore, conservatively bounds the lower temperatures of the Model D5 SLB transient. The T/TS contact pressure results during SLB are still expected to bound the TITS contact pressure results during NOP because, even though the tube bore eccentricity during SLB is generally greater than that during NOP, the overall growth of the tube bore during NOP is greater than that during SLB. Larger magnitudes of tube bore growth are directly related to decreasing tube-tubesheet contact pressure regardless of the value of calculated tube bore eccentricity.
It is appropriate to compare the Model D5 SLB contact pressure results from the "new" method to the Model D5 NOP results from the "old" method because each condition uses the appropriate fitto conservatively determine the reduction in T/TS contact pressure due to tube bore eccentricity and tube bore growth.
The sole purpose of the new methodology was to develop a more accurate way of calculating contact pressures during a postulated SLB for the Model D5 steam generators. The comparison provided in Figure 6-83 of WCAP-17072-P remains a valid comparison.
41
LTR-SGMP-09-109 NP-Attachment 5.0 Part A (Original RAI#4)
Reference 1, Page 6-69: In Section 6.2.5.3, it is concluded that the tube outside diameter and the tubesheet tube bore inside diameteralways maintain contactin the predicted range of tubesheet displacements. However, for tubes with through-wallcracks at the H*
distance, there may be little or no net pressure acting on the tube for some distance above H*. In Tables 6-18 and 6-19, the fourth increment in the step that occurs two steps prior to the last step suggests that there may be no contact between the tube and tubesheet, over a portion of the circumference, for a distance above H*. Is the conclusion in 6.2.5.3 valid for the entire H* distance,given the possibility that the tubes may contain through-wall cracks at that location?
The following response to RAI#4 was included in Reference 11. The same response is included here to complete the record of information provided in regard to RAI#4 of References 5, 6 and 7.
Response
The conclusions reached in Section 6.2.5.3 of WCAP-17071-P are valid for the entire H*
distance because of the following considerations:
- 1. The primary source of contact pressure between the tube and the tubesheet is differential thermal expansion between the tubes and the tubesheet. The analysis in Section 6.2.5.3 specifically excludes the effect of thermal expansion of the tube from the analysis. The tubesheet is assumed to deform due to the combination of pressure and thermal loads which produces the tube bore ovalization and leads to the displacements applied in this model. Only the residual effects from installation are considered for the tube in steps 1 through 5. The tube internal pressure applied in these steps only simulates the hydraulic expansion pressure to establish the initial conditions for the following step. The conditions assumed for this study are not possible during any operating condition in the steam generator but are conservative relative to actual SG conditions. (Note: Residual contact pressure is not used in the calculation of H* values in Section 6. The residual effects of installation are included in the results of Section 6.2.5.3 so that the sensitivity of a strain hardened tube to tubesheet tube bore deformation can be studied.)
- 2. Step 5 on Tables 6-18 and 6-19 is not representative of any condition in the steam generator because it assumes that the tubesheet is at operating temperature with an applied primary-to-secondary pressure differential while the tubes remain at room temperature and are not pressurized. That is why Steps 1 through 5 are described as "initializing" steps in the process. It is physically impossible for these conditions to occur simultaneously in the same steam generator.
42
LTR-SGMP-09-109 NP-Attachment
- 3. Because no pressure loading is applied to the tube in Step 5 of the analysis discussed in section 6.2.5, the results presented in Tables 6-18 and 6-19 are applicable regardless of whether, or not, a through-wall crack exists at the H* location. The more representative case is Step 6 shown on Tables 6-18 and 6-19 in which tube internal pressure is included. For that case, the potential point of zero contact pressure is at an applied displacement a factor of 5 greater than for Step 5, and far in excess of what is reasonably predicted for the actual tubesheet deformation. The factor of 5 difference in required displacement to cause the contact pressure to reduce to zero more than adequately covers the postulated potential local reduction in crevice pressure due to a circumferential separation at the location of H*. Recall also, that no thermal expansion of the tube is considered in this analysis.
It is also noted that tables 6-18 and 6-19 are the results of a sensitivity study that is not intended to represent the integrated calculation for H*. The integrated H* analysis is a complex process that combines the effects of several types of loading and deformation into an integrated estimate of the tube-to-tubesheet contact pressure. Therefore, it is not appropriate to consider a sensitivity study out of the context of the greater analysis. The integrated analysis presented in the complete Section 6 shows that for the combined case of the thermal effects, pressure effects, and tubesheet displacement there is tube-to-tubesheet contact pressure throughout the tubesheet.
It is acknowledged that the cut end of a tube is radially less stiff than a tube that is radially loaded at a point away from the tube end, and that the presumption of a tube sever at the H*
distance may represent the case of a tube end. The decreased tube-end stiffness is referred to as "compliance." In other words, a tube that is loaded at the cut end provides less resistance to the load than a tube with equal load applied a distance removed from the tube-end. Thus, conceptually, a local "end effect" could be expected to occur due to the increased compliance of the tube-end.
The calculation process for H* shown in Figure 1-1 of the H* WCAP reports and discussed in several places in the report notes that an adjustment is made to the initial prediction of H* to account for the distributed crevice pressure referenced to the predicted H* position. Thus, the greatest crevice pressure is always located at the final value of H*. Increased tube compliance cannot result in a higher local crevice pressure than is already included in the analysis because, at the point of sever, the primary side pressure is the crevice pressure.
It may be postulated that the increased tube compliance results in reduced contact pressure because the net differential pressure across the tube wall is zero. At the tube-end, the current analysis already includes a zero differential pressure due to the adjustment process for distributed crevice pressure. Therefore, the net reduction in contact pressure would be limited to the axial length of the local effect and would further depend on the slope of the decrease in crevice pressure.
For the Model F and Model D5 SGs, the bounding value of isolation distance above the tube end is 0.6 inch based on classical solutions for the design of pressure vessels (Timoshenko).
43
LTR-SGMP-09-109 NP-Attachment The isolation distance is the generically applicable minimum separation distance from an applied load to a point of interest in order to safely assume that the load is in the far field relative to the point of interest. Specific structures and load cases may have different isolation distances but the classical result by Timoshenko for a pressure vessel will conservatively bound any specific cases. For this length, the slope of the contact pressure curve would have to decrease by a factor of at least [ ]a,ce before the value of H* is affected by more than [ ]a,cC inch. If the tube is conservatively modeled as a center-loaded beam on an elastic foundation compared to an end-loaded beam on an elastic foundation, the resulting worst case change in structural compliance and the resulting contact pressure slope could be a factor of up to 2. Alternatively, similar analyses for the cross sections of curved beams suggest that the change in compliance of the structure could be as high as a factor of 6. Neither case approaches the factor of [ ],,c,e required based on classical pressure vessel analysis to impact the value of H*; therefore, no additional adjustments to H* are necessary to address the potential end effects.
6.0 Summary of the Response to RAI #4 A summary of the response to the original RAI# 4 and additional questions related to RAI 4 are provided below:
- 1. No additional adjustment to the value for H* is necessary to address the potential for end effects. This is because the greatest crevice pressure is always located at the final value of H*. At the H* distance, the current analysis already includes a zero pressure differential due to the adjustment process for the distributed crevice pressure. Therefore, the net reduction in contact pressure would be limited to the axial length of the local effect and would further depend on the slope of the decrease in crevice pressure. It is judged that the slope of the contact pressure curve would not decrease at a rate such that the value of H*
would be affected.
- 2. Tube bore dilation is a more significant factor in determining tube-to-tubesheet contact pressure at higher temperatures and the effect of eccentricity on contact pressure is reduced at higher temperatures. The methodology for addressing the effect of eccentricity on contact pressure discussed in Reference 6-15 and utilized in WCAP-17071-P, WCAP-17072-P, WCAP-17091-P and WCAP-17092-P reflects this fact and, therefore, it provides acceptably accurate contact pressure results at higher temperatures (i.e., for all conditions except the "colder" SLB condition). This includes NOP, SLB (higher temperature, > 400 0 F, and FLB, where appropriate).
- 3. The results of using the fit described in Reference 6-15 match the expected trend from a best case finite element model for the NOP and SLB conditions for the Model F SGs and NOP conditions for the Model D5 SG.
- 4. The ADs from the 3D FEA model are significantly less than the corresponding ADs from the unit cell model from the unloaded to the fully loaded condition (i.e.,
44
LTR-SGMP-09-109 NP-Attachment from step 0 to step 9) for both NOP and SLB conditions. This leads to the conclusion that the unit cell model displacement results and contact pressure predictions conservatively represent the reference 3D FEA model results.
- 5. The eccentricities from the unit cell model are generally comparable to those from the 3D FEA model. A more exact comparison is difficult based on the available data; however, it is clear that the actual range of eccentricities from the 3D FEA model was adequately addressed by the unit cell model.
- 6. Based on items 4) and 5) which demonstrate the acceptability of the use of the unit cell model for benchmarking the 3-D FEA model, the method for calculating the reduction in contact pressure defined by the unit cell model, when adjusted for temperature effects, shows that SLB contact pressure is increased relative to normal operating conditions for the Model D5 steam generators.
- 7. It has also been shown when comparing contact pressures for NOP conditions for the unit cell to contact pressures for SLB for the unit cell for the Model F SG (higher temperature SLB conditions), that contact pressure increases during a postulated SLB.
- 8. Given that the two structural models provide similar eccentricities, the unit cell model shows that for these eccentricities, positive contact pressure exists between the tubes and the tubesheet for the entire range of displacements considered.
Based on the above, it is concluded that the NOP contact pressures that define H* in the Model F and Model D5 SG are conservative and that a more realistic model of contact pressure reduction as a function of tube bore deformation (including both dilation and eccentricity) would predict positive contact pressure around the entire circumference of the tube and an-increase in tube to tubesheet contact pressure at SLB conditions compared to NOP conditions.
The conclusions reached in the response to RAI#4 apply equally for the Model 44F and Model 51F SGs.
45
Attachment 9 List of Regulatory Commitments
List of Regulatory Commitments The following table identifies those actions committed to by NextEra Energy Seabrook, LLC Seabrook in this document. Any other statements in this submittal are provided for information purposes and are not considered to be regulatory commitments. Please direct questions regarding these commitments to Mr. Michael O'Keefe, Licensing Manager.
Regulatory Commitment Due Date / Event
- 1. NextEra Energy Seabrook, LLC commits to Prior to entering Mode 4 perform a one-time verification of tube during startup following expansion locations to determine if any refueling outage 13 in the significant deviations exist from the top of fall of 2009 tubesheet to the bottom of expansion transition (BET). If any significant deviations are found, the condition will be entered into the plant's corrective action program and dispositioned.
Additionally, if any significant deviations are found, those deviations will be reported to the NRC via the Steam Generator Tube Inspection Report or separate timely correspondence.
- 2. NextEra Energy Seabrook, LLC commits to Each inspection of the monitor for tube slippage as part of the steam Seabrook Station Steam generator tube inspection program. Slippage Generators monitoring will occur for each inspection of the Seabrook Station Steam Generators.
- 3. For the Condition Monitoring assessment, the During each inspection component of operational leakage from the prior required by Technical cycle from below the H* distance will be Specification 6.7.6.k, "Steam multiplied by a factor of 2.50 and added to the Generator (SG) Program" total accident leakage from any other source and compared to the allowable accident induced leakage limit. For the Operational Assessment, the difference between the allowable accident induced leakage and the accident induced leakage from sources other than the tubesheet expansion region will be divided by 2.50 and compared to the observed operational leakage.
An administrative operational leakage limit will be established to not exceed the calculated value.
Attachment 9 List of Regulatory Commitments
List of Regulatory Commitments The following table identifies those actions committed to by NextEra Energy Seabrook, LLC Seabrook in this document. Any other statements in this submittal are provided for information purposes and are not considered to be regulatory commitments; Please direct questions regarding these commitments to Mr. Michael O'Keefe, Licensing Manager.
Regulatory Commitment Due Date / Event
- 1. NextEra Energy Seabrook, LLC commits to Prior to entering Mode 4 perform a one-time verification of tube during startup following expansion locations to determine if any refueling outage 13 in the significant deviations exist from the top of fall of 2009 tubesheet to the bottom of expansion transition (BET). If any significant deviations are found, the condition will be entered into the plant's corrective action program and dispositioned.
Additionally, if any significant deviations are found, those deviations will be reported to the NRC via the Steam Generator Tube Inspection Report or separate timely correspondence.
- 2. NextEra Energy Seabrook, LLC commits to Each inspection of the monitor for tube slippage as part of the steam Seabrook Station Steam generator tube inspection program. Slippage Generators monitoring will occur for each inspection of the Seabrook Station Steam Generators.
- 3. For the Condition Monitoring assessment, the During each inspection component of operational leakage from the prior required by Technical cycle from below the H* distance will be Specification 6.7.6.k, "Steam multiplied by a factor of 2.50 and added to the Generator (SG) Program" total accident leakage from any other source and compared to the allowable accident induced leakage limit. For the Operational Assessment, the difference between the allowable accident induced leakage and the accident induced leakage from sources other than the tubesheet expansion region will be divided by 2.50 and compared to the observed operational leakage.
An administrative operational leakage limit will be established to not exceed the calculated value.