ML15075A349

From kanterella
Jump to navigation Jump to search
NUH32PHB-0411, Revision 0, Grid Convergence Study for the Ansys Model for the Nuhoms 32PHB DSC, Non-Proprietary
ML15075A349
Person / Time
Site: Calvert Cliffs  Constellation icon.png
Issue date: 03/19/2015
From:
AREVA
To:
Office of Nuclear Material Safety and Safeguards
Shared Package
ML15075A350 List:
References
NUH32PHB-0411, Rev.0
Download: ML15075A349 (26)
v • d • e


Text

ENCLOSURE 16 Non-Proprietary NUH32PHB-0411, Revision 0, Grid Convergence Study for the ANSYS Model for the NUHOMS 32PHB DSC Calvert Cliffs Nuclear Power Plant March 10, 2015

CONTROLLED COPY E-281 I Form 3.2-1 Calculation No.: NUH32PHB.0411 Calculation Cover Sheet Revision No.: 0 AR EVA Revision 8 Page: I of 21 0

DCR NO (if applicable): NUH32PHB-018 PROJECT NAME: NUHOMS 32PHB System PROJECT NO: 10955 CLIENT: CENG - Calvert Cliff Nuclear Power Plant (CCNPP)

CALCULATION TITLE:

Grid Convergence Study for the ANSYS Model for the NUHOMS 32PHB DSC

SUMMARY

DESCRIPTION:

1) Calculation Summary This calculation performs a grid convergence study of the ANSYS finite element model used for the thermal evaluation of the NUHOMS 32PHB DSC using the procedure specified in ASME V&V 20-2009

[4] and determines the discretization error of the solution.

2) Storage Media Description Secure network server initially, then redundant tape backup If original issue, is licensing review per TIP 3.5 required?

Yes El No 0 (explain below) Licensing Review No.: N/A This calculation is prepared to support a Site Specific License Application by CCNPP that will be reviewed and approved by the NRC. Therefore, a 10CFR72.48 licensing review per TIP 3.5 is not applicable.

Software Utilized: Version:

ANSYS 10.OA1 Calculation is complete:

Originator Name and Signature: Sharad Pachpute Date:

Calculation has been checked for consistency, completeness and correctness:

Checker Name and Signature: Amit Mathur & Date: 0 3)04 11-"

Calculation is approved for use:

PATELGirish

    • o=AREVA GROUP, 2.5.4.45=Tl 1D2D8D413995674D417FCF, cn=PATEL Girish 2015.03.04 09:18:24 -05'00' Project Engineer Name and Signature: Date:

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 2of21 REVISION

SUMMARY

REV. DATE DESCRIPTION AFFECTED AFFECTED PAGES DISKS 0 3/4/15 Initial Issue All All

A Calculation Calculation No.:

Revision No.: 0 NUH32PHB-0411 AR EVA Page: 3of21 TABLE OF CONTENTS Paqe 1.0 P u rp o s e ............................................................................................................................. 4 2.0 Assumptions and Conservatism ................................................................................. 4 3.0 Design Input/D ata ...................................................................................................... . . 4 3.1 Theoretical Order of Accuracy of ANSYS Solver ................................................. 4 4 .0 Me th o d o lo g y ...................................................................................................................... 4 5 .0 R e fe re nce s ........................................................................................................................ 8 6.0 Nom enclature ...................................................................................................... .. ... 9 7 .0 C o m p utatio ns .................................................................................................................. 10 8.0 Results / Conclusion and Recommendations ............................................................ 11 9.0 List of C om puter Files ............................................................................................... . . 13 LIST OF TABLES Page Table 1 Grid Sizes and Refinement Factors ............................................................. 14 Table 2 Maximum Fuel Cladding Temperatures ........................................................ 14 Table 3 Discretization Error of the 32PHB DSC Finite Element Model ...................... 15 Table 4 ANSYS Model Summary .............................................................................. 16 Table 5 ANSYS DSC Shell Profile Map Summary .................................................... 17 Table 6 ANSYS Run Summary ................................................................................. 18 Table 7 Associated Files and Macros ........................................................................ 18 LIST OF FIGURES Page Figure 1 Maximum Fuel Cladding Temperature Vs Number of Elements ................... 19 Figure 2 Axial and Radial Mesh Distribution for Grids # 1 and 2 ................................ 20 Figure 3 Axial and Radial Mesh Distribution for Grids # 3 and 4 ................................ 21

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 4 of 21 1.0 PURPOSE This calculation performs a grid convergence study of the ANSYS finite element model used for the thermal evaluation of the NUHOMS 32PHB DSC in [1] and determines the discretization error of the solution. The discretization error is determined using the five steps specified in Section 2-4.1 of ASME V&V 20-2009 [4]. In addition to the discretization error, this calculation computes the observed order of accuracy (p) and compares it to the theoretical order of accuracy of the ANSYS solution. A total of four grids are considered for this evaluation.

2.0 ASSUMPTIONS AND CONSERVATISM The assumptions and conservatism considered for 32PHB DSC model are the same as those described in Section 3.0 of [1] and remain unchanged.

3.0 DESIGN INPUTIDATA 3.1 Theoretical Order of Accuracy of ANSYS Solver PROPRIETARY 4.0 METHODOLOGY PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 5of21 PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 6 of 21 PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 Page: 7of21 PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 8of21

5.0 REFERENCES

1 Calculation, "Thermal Evaluation of NUHOMS 32PHB DSC for Storage and Transfer Conditions", CCNPP, Calculation No. NUH32PHB-0403, Rev.1.

2 Calculation, "Thermal Evaluation of NUHOMS 32PHB Transfer Cask for Normal, Off-Normal, and Accident Conditions", Transnuclear, Inc., Calculation No. NUH32PHB-0402, Rev.l.

3 ANSYS Computer Code and On-Line User's Manuals, Version 10.OA1.

4 American Society of Mechanical Engineers, "Standard for Verification and Validation in Computational Fluid Dynamics and Heat Transfer", ASME V&V 20-2009, November 30th 2009.

5 Fundamental FEA Concepts and Applications, ANSYS Inc.

6 Examining Spatial (Grid) Convergence, NPARC Alliance CFD Verification and Validation Website, http://www.qrc.nasa..qov/WWW/wind/valid/tutorial/spatconv.html 7 Ismail B. Celik et.al, "Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications," Journal of Fluids Engineering, ASME, July 2008, Vol. 130.

8 USNRC, Office of Nuclear Regulatory Research, "Thermal Analysis of Horizontal Storage Casks for Extended Storage Applications," NUREG/CR-7191, Dec. 2014.

9 USNRC, Office of Nuclear Material Safety and Safeguards, "Standard Review Plan for Dry Cask Storage Systems - Final Report," NUREG-1536, Rev.1, July 2010.

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 9of21 6.0 NOMENCLATURE DSC Dry shielded canister HLZC Heat load zone configuration ASME American society of mechanical engineers GCI Grid convergence index N Number of elements h Representative cell size r Refinement factor N

AVi Volume of the half-symmetric DSC model

r. Outer radius of the 32PHB DSC model L Length of the 32PHB DSC model

) Simulation variable (maximum fuel cladding temperature) p Observed order of accuracy Difference in simulation variable between consecutive mesh refinements 21 (l)2e*Extrapolated simulation variable ea Relative error Fs Factor of safety GCIne Fine grid GCI Unum Discretization error (dimensional)

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 10 of 21 7.0 COMPUTATIONS PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 11 of 21 8.0 RESULTS I CONCLUSION AND RECOMMENDATIONS PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 12of21 PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 13of21 9.0 LIST OF COMPUTER FILES All the runs are performed using ANSYS version 10.OA1 [3] on ATN022 with an Intel(R)

Xeon(R) CPU E5-1620 @ 3.60 GHz 64-bit processor.

The ANSYS models, associated files and macros are listed in Table 4 through Table 7.

A Calculation Calculation No.:

Revision No.:

NUH32PHB-0411 0

AREVA Page: 14 of 21 Table I Grid Sizes and Refinement Factors PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 15of21 Table 3 Discretization Error of the 32PHB DSC Finite Element Model PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 16 of 21 Table 4 ANSYS Model Summary Run Name Description Date and Time Acc48.db Database file used as the input to create the refined mesh Acc48 _____db__ model [1]

__Sd_ dfd A macro is used to generate the nodes and elements in air 12/22/2014 gaps. 05:10 PM Grid # 2 Model Grid # 2 Step 1:

32PHB Mesh_G2.[ext] This step is used to extract the geometry of a cross section

[ext] = inp, out, err, db from the database model of Grid # 1. The mesh in this section 01/09/2015 is refined and then the cross section is extruded to get a sliced 19:57:57 32PHBMESHMAPG2.node model. The Material Types and Element Types are assigned in the subsequent steps.

32PHBMESH_MAPG2.[ext] Grid # 2 Step 2:

[Ext] = inp, out, err, db, mntr, rth This step reads the database model of Grid # 1 simulation and 01/09/2015 32PHBMESHMAP G2 4.bfin assigns arbitrary temperatures to the components which are 20:04:00 32PHBMESHMAPG2_7.bfin required to be segregated. The temperatures are mapped onto 32PHB MESH MAP G2 8.bfin the node file created in Step 1 and three .bfin files are created.

Grid # 2 Step 3:

With the help of the .bfin files created in Step 2, temperatures 32PHBMESH_MATG2.[ext] are mapped on to the database model created in Step 1. 01/09/2015

[ext] = inp, out, err, db Elements are segregated with the help of temperatures 20:04:11 assigned to them and the Material Types and Element Types are assigned to all the components of the sliced model Grid # 2 Step 4: 12/24/2014 32PHBSlice_G2.inp This file is used to fill-in the air gaps in the sliced model created 08:13 PM in Step 3.

32PH Model_G2.[ext] #2 Ste 5: 01/09/2015

[ext] = np,3 out, err, db Sliced model created after filling the air gaps are reflected to 20:33:30 achieve a full length Finite Element Model for Grid # 2.

Grid # 3 Model Grid # 3 Step 1:

32PHBMesh_G3.[ext] This step is used to extract the geometry of a cross section

[ext] = inp, out, err, db from the database model of Grid # 1. The mesh in this section 01/09/2015 is refined and then the cross section is extruded to get a sliced 21:07:09 32PHBMESH_MAPG3.node model. The Material Types and Element Types are assigned in the subsequent steps 32PHBMESH_MAPG3.[ext] Grid # 3 Step 2:

[ext] = inp, out, err, db, mntr, rth This step reads the database model of Grid

  1. 1 simulation and 01/09/2015 32PHBMESHMAP G3 4.bfin assigns arbitrary temperatures to the components which are 21:18:20 32PHBMESHMAP G3_7.bfin required to be segregated. The temperatures are mapped onto 32PHBMESHMAPG3 8.bfin the .node file created in Step 1 and three .bfin files are created.

Grid # 3 Step 3:

With the help of the .bfin files created in Step 2, temperatures 32PHBMESH_MATG3.[ext] are mapped on to the database model created in Step 1. 01/09/2015

[ext] = inp, out, err, db Elements are segregated with the help of temperatures 21:18:46 assigned to them and the Material Types and Element Types are assigned to all the components of the sliced model.

Grid # 3 Step 4: 01/08/2015 32PHB_Slice_G3.inp This file is used to fill-in the air gaps in the sliced model created 11:00 AM in Step 3.

Grid # 3 Step 5:

32PHB

[ext] =BModel_G3.[ext]

inp, out, err, db Grd#3Se :01/09/2015 Sliced model created after filling the air gaps are reflected to 23:14:25 e achieve a full length Finite Element Model for Grid # 3.

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 17of21 Grid # 4 Model Grid # 4 Step 1:

32PHBMesh_G4.[ext] This step is used to extract the geometry of a cross section

[ext] = inp, out, err, db from the database model of Grid # 1. The mesh in this section 01/10/2015 32PHB_MESHMAPG4.node is refined and then the cross section is extruded to get a sliced 00:24:07 model. The Material Types and Element Types are assigned in the subsequent steps 32PHBMESH_MAPG4.[ext] Grid # 4 Step 2:

[ext] = inp, out, err, dib, mntr, rth This step reads the database model of Grid # 1 simulation and 01/10/2015 32PHBMESHMAP G4 4.bfin assigns arbitrary temperatures to the components which are 00:44:47 32PHBMESHMAP G4 7.bfin required to be segregated. The temperatures are mapped onto 32PHB MESH MAP G4 8.bfin the node file created in Step 1 and three .bfin files are created.

Grid # 4 Step 3:

32PHBMESH_MATG4.[ext] With the help of the .bfin files created in Step 2, temperatures

[ext] = inp, out, err, db are mapped on to the database model created in Step 1. 01/10/2015 Elements are segregated with the help of temperatures 00:45:38 assigned to them and the Material Types and Element Types are assigned to all the components of the sliced model.

Grid # 4 Step 4: 01/08/2015 32PHB_Slice_G4.inp This file is used to fill-in the air gaps in the sliced model created 11:05 AM in Step 3.

32PHBModel_G4.[ext] #4 Ste 5: 01/10/2015

[ext]= inp, out, err, db Sliced model created after filling the air gaps are reflected to 06:52:43 achieve a full length Finite Element Model for Grid # 4.

Table 5 ANSYS DSC Shell Profile Map Summary Date Timeand Run Name Description 32PHB TC VERTTRANS_20hr MapG2.[ext] Vertical transfer DSC shell temperature profile 01/09/2015

[ext] = inp, out, err, cbdo @ 20 hrs from Transfer Cask model for Grid # 2 20:35:58 32PHB TC VERTTRANS_20hrMapG3.[ext] Vertical transfer DSC shell temperature profile 01/09/2015

[ext] = inp, out, err, cbdo @ 20 hrs from Transfer Cask model for Grid # 3 23:17:38 32PHBTCOVERTTRANS_20hrMapG4.[ext] Vertical transfer DSC shell temperature profile 01/10/2015

[ext] = inp, out, err, cbdo @ 20 hrs from Transfer Cask model for Grid # 4 06:58:28

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 A REVA Page: 18of21 Table 6 ANSYS Run Summary Run Name Description Date and Time 32PHBTC2M.[ext] Run for Grid # 1 to determine the maximum fuel cladding

[ext] temperature at Vertical Transfer Conditions Table 8-1 in [1]

[ext] = inp, out, err, db, mntr, rth @ 20 hrs, 29.6 kW 32PHBTC2M_G2.[ext] Run for Grid # 2 to determine the maximum fuel cladding 01/09/2015

[ext] = inp, out, err, db, mntr, rth temperature at Vertical Transfer Conditions 21:03:30

@ 20 hrs, 29.6 kW 21:03:30 32PHBTC2MG3.[ext] Run for Gridat#Vertical temperature 3 to determine the maximum fuel cladding Transfer Conditions 01/10/2015

[ext] = inp, out, err, db, mntr, rth @ 20 hrs, 29.6 kW 00:19:11 32PHBTC2M_G4.[ext] Run for Grid # 4 to determine the maximum fuel cladding 01/10/2015

[et]inp, out, err, db, mntr, rth temperature at Vertical Transfer Conditions

@ 20 hrs, 29.6 kW 09:24:47

[ext]

Table 7 Associated Files and Macros FilelMacro Description Date and Time 32PHB TC VERT TRANS.db, Database and Results files to retrieve the DSC shell [2]

32PHBTC VERT-TRANS.rth temperature profile for Vertical Transfer @ 20hrs, 29.6 kW 32PHBHLZC2.MAC Macro for Heat Generation Boundary Condition for all models, Table 8-4 in Grid # 1,2,3 and 4 [1l Macro.mac Macro to get Maximum/Minimum temperatures Table 8-4 in Results.mac Macro to list maximum and average 32PHB DSC component Table 8-4 in temperatures [1]

32PHBMat1.inp Material database for 32PTH2 DSC model Table 8-4 in 01/27/15 05:04/PM NUH32PHB_041 1_RO.xlsx Spreadsheet to compute the discretization errors 05:04 PM 1

Calculation No.: NUH32PHB-0411 Revision No.: 0 Page: 19 of 21 PROPRIETARY

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AR EVA Page: 20of21 PROPRIETARY Figure 2 Axial and Radial Mesh Distribution for Grids # I and 2

Calculation No.: NUH32PHB-0411 Calculation Revision No.: 0 AREVA Page: 21of21 PROPRIETARY Figure 3 Axial and Radial Mesh Distribution for Grids # 3 and 4

ATTACHMENT (2)

MARKED UP TECHNICAL SPECIFICATION PAGE Calvert Cliffs Nuclear Power Plant March 10, 2015

Markup of Technical Specifications CN eu -~

TCdJ%.,,cQ Scoh_4C "V )

3/4.3 TRANSFER CASK 3/4.3.2 TIME LIMIT FOR COMPLETION OF NUHOMS 32 PHB DSC TRANSFER OPERATION LIMITING CONDITION FOR OPERATION 3.32.1 The time limit for completion of transfer of a loaded and welded NUHOMS 32 PHB DSC from the cask handling area to the HSM is as follows:

a. No time limit for a DSC with a total heat load of < 21,12 kW 7l'"°ours for a DSC with a total heat load > 21.12 kW and s 23.04 kW 7b
c. 46 ,hours for a DSC with a total heat load > 23.04 kW and 5 25.6 kW
d. 1-8 20Mhours for a DSC with a total heat load > 25.6 kW and ;29.6 kW APPLICABILITY: This specification is applicable to NUHOMS 32 PHB DSCs only. The time limit is defined as the time elapsed after the initiation of draining the transfer cask/DSC annulus water until completion of insertion of the DSC into the HS3 t < RAI-6-2a ACTION: Initiate one of the following actions within e-1gk4 hours if the specified time limit is exceeded. The chosen action may be temporarily suspended under administrative controls to change from one action to another.
1. Complete the transfer of the DSC to the HSM or,
2. If the transfer cask is in the cask handling area in a vertical orientation fill the transfer casK!DSC annulus with clean water or,
3. Ifthe transfer cask is in a horizontal orientation, initiate air circulation by starting one of the blowers provided on the transfer skid or,
4. eturn the transfer cask to the cask handling area and fill the transfer caslDSC Insert-I for RAI annulus with clean water,-e;-i-ate-appFegoa ernaI-ee-,g-of-tlhi*e-trarfr-6-2b *.

Insert-1

  • If air circulation is initiated, it must be maintained for a minimum duration of 8 hours9.259259e-5 days <br />0.00222 hours <br />1.322751e-5 weeks <br />3.044e-6 months <br /> before it is turned off. Once the air circulation is turned off, a maximum duration of 6 hours6.944444e-5 days <br />0.00167 hours <br />9.920635e-6 weeks <br />2.283e-6 months <br /> is available to complete the transfer to the HSM-HB or re-establish the air circulation.

, If air circulation is initiated and maintained for a minimum duration of 20 hours2.314815e-4 days <br />0.00556 hours <br />3.306878e-5 weeks <br />7.61e-6 months <br /> before if it is turned off, a maximum duration of 8 hours9.259259e-5 days <br />0.00222 hours <br />1.322751e-5 weeks <br />3.044e-6 months <br /> is available to complete the transfer to the HSM-HB or re-establish the air circulation.

ATTACHMENT (3)

AREVA TN AFFIDAVIT Calvert Cliffs Nuclear Power Plant March 10, 2015

E-41330 R-1 Enclosure I AFFIDAVIT PURSUANT TO 10 CFR 2.390 AREVA Inc.

State of Maryland ) SS.

County of Howard )

I, Paul Triska, depose and say that I am a Vice President of AREVA Inc., duly authorized to execute this affidavit, and have reviewed or caused to have reviewed the information which is identified as proprietary and referenced in the paragraph immediately below. I am submitting this affidavit in conformance with the provisions of 10 CFR 2.390 of the Commission's regulations for withholding this information.

The information for which proprietary treatment is sought listed below:

  • AREVA TN Design Calculation NUH32PHB-0203 R-2
  • AREVA TN Design Calculation NUI32PHB-0203 R-1
  • AREVA TN Design Calculation, NUH32PHB-0401 R-2
  • AREVA TN Design Calculation, NUH32PHB-0411 R-0 0 AREVA TN Design Calculation, 10955-TLAAO1 R-0 This documentation has been appropriately designated as proprietary.

I have personal knowledge of the criteria and procedures utilized by AREVA Inc. in designating information as a trade secret, privileged or as confidential commercial or financial information.

Pursuant to the provisions of paragraph (b) (4) of Section 2.390 of the Commission's regulations, the following is furnished for consideration by the Commission in determining whether the information sought to be withheld from public disclosure, included in the document described above, should be withheld.

1) The information, which is owned and has been held in confidence by AREVA Inc., sought to be withheld from public disclosure involves details regarding AREVA Inc.'s approach to Calvert Cliffs Nuclear Power Plant intended use of the 32PHB dry spent fuel storage system.
2) The information is of a type customarily held in confidence by AREVA Inc. and not customarily disclosed to the public. AREVA Inc. has a rational basis for determining the types of information customarily held in confidence by it.
3) Public disclosure of the information is likely to cause substantial harm to the competitive position of AREVA Inc. because the information consists of details regarding AREVA Inc.'s approach to Calvert Cliffs Nuclear Power Plant intended use of the 32PHB dry spent fuel storage system, the application of which provides a competitive economic advantage. The availability of such information to competitors would enable them to modify their product to better compete with AREVA Inc., take marketing or other actions to improve their product's position or impair the position of AREVA Inc.'s product, and avoid developing similar data and analyses in support of their processes, methods or apparatus.

Further the deponent sayeth not. P r.

Paul Triska %I111,,

Vice President, AREVA Inc. . .... *,,.

Sub ribed and sw to e before this 9 th day of March, 2015. "

N PUBLlC Zy Commission Expires lo/ /7 / 2ao, *',"k  !-xa'.

,. COU%4.,-

0 Page 1 of I

Retrieved from "https://kanterella.com/w/index.php?title=ML15075A349&oldid=2449511"