ML18086B208
| ML18086B208 | |
| Person / Time | |
|---|---|
| Site: | University of Lowell |
| Issue date: | 03/31/2017 |
| From: | Edward Helvenston Research and Test Reactors Licensing Projects Branch |
| To: | |
| Helvenston E, NRR/DLP, 415-4067 | |
| References | |
| Download: ML18086B208 (1) | |
Text
NOTES:This file focuses on the flow distribution calculation - all axial power info in separate file Here we do calcs for both the UMLRR and WPI fuel elementsFor the bypasses, all the calcs have 5 RBs, 4 control blades , and 1 RegBladeWe also do a calc based on only flow area for comparison purposesBasic geometry parameters and assumption are tabulated at the top and then the PLTEMP input is prepared belowThis file is based on the original work done by JRWhite back in 1991 for the HEU to LEU conversion.Added the power-to-flow plot on 7/30/15Reviewed in March 2017 -- Fixed a couple of typos in the comments and changed the # channels from 19 to 18 in the min flow area calc for the fuel channels... This was modified to account for the fact that the two end channels are only approximately 1/2 the size of the interior channels (19 full channels was simply not correct). This change only affected the flow area distribution for the no-friction case (Tables 8 and 9 in the safety analysis report), but this was only used for comparison purposes anyway. Thus, this change does not affect any of the results and conclusions from the PLTEMP full-core analysis.Fuel Assembly Info (units ==> centimeters)UMLRR WPIfueltop channel area (assume 2.6" x 2.6")43.61343.613top channel length (about 4.5" for top end box)11.43011.430# fuel assyfuel5 rad baskets4 large blades1 regbladetotalplate height 63.50062.54820704.43588.595106.81910.068909.918channel width6.6046.70621739.65788.595106.81910.068945.140channel thickness0.2960.27122774.87988.595106.81910.068980.362bottom channel area (assume 2.6" x 2.6")43.61343.61323810.10188.595106.81910.0681015.584bottom channel length (about 6" for grid plate)15.24015.24024845.32388.595106.81910.0681050.805min flow area35.22232.69825880.54488.595106.81910.0681086.02726915.76688.595106.81910.0681121.249rad basketflow area (1.87" diameter hole)17.719length (35")88.900min flow area17.719# fuel assyfuel5 rad baskets4 large blades1 regbladetotal200.7740.0970.1170.0111.000control bladeshroud area (10.844 x 0.75)52.471210.7830.0940.1130.0111.000blade area (10.65 x 0.375)25.766220.7900.0900.1090.0101.000top control in flow area (shroud - blade)26.705230.7980.0870.1050.0101.000top length (30")76.200240.8040.0840.1020.0101.000middle control out flow area52.471250.8110.0820.0980.0091.000middle length (21")53.340260.8170.0790.0950.0091.000bottom (10.844 x.47) -- very conservative since ignore holes32.882bottom length (6")15.240min flow area26.705# fuel assyfuel5 rad baskets4 large blades1 regbladetotal20653.95788.595106.81910.068859.440regbladeshroud area (2.465 x 2.465)39.20121686.65588.595106.81910.068892.138blade area (2.125 x 2.125)29.13322719.35388.595106.81910.068924.835top control in (shroud - blade)10.06823752.05088.595106.81910.068957.533top length (26")66.04024784.74888.595106.81910.068990.231middle control out 39.20125817.44688.595106.81910.0681022.929middle length (25")63.50026850.14488.595106.81910.0681055.627bottom (grid box)43.613bottom length (6")15.240min flow area10.068# fuel assyfuel5 rad baskets4 large blades1 regbladetotal200.7610.1030.1240.0121.000convert cm to m100210.7700.0990.1200.0111.000convert cm^2 to m^210000220.7780.0960.1160.0111.000230.7850.0930.1120.0111.000PLTEMP Data for LEU Core Flow Distribution Calculation240.7920.0890.1080.0101.000Card 0100250.7990.0870.1040.0101.000ANAME case descriptionumlrr_flow21wpi_flow21260.8050.0840.1010.0101.000Card 0200IHheat transfer correlation option0IBboiling correlation option0Summary ResultsICHFCHF correlation option0NFTYPnumber of fuel types1# fuel assyflow area calcPLTEMP calcflow area calcPLTEMP calcNCTYPnumber of bypass types3200.7830.6880.7710.677NEDITnumber of Del P intervals between edits20210.7920.7050.7790.689NAXDISaxial temp calc flag (0/1 no/yes)0220.7990.7140.7870.699NELPRTelement edit flag (ULowell mod - usually zero)0230.8060.7220.7940.708Card 0300 (fuel type 1)240.8130.7320.8010.716NELFnumber of fuel elements21250.8190.7400.8080.724NFnumber of axial regions3260.8250.7460.8140.730WFGESguess for flow rate (kg/s)1.00FBHCF for bulk water temp rise1.00FQHCF for heat flux1.00FHHCF for heat transfer coeff1.00Card 0301 (fuel type 1)FZ(J)axial peaking factors (NELF entries) -- no effect here1.50Card 0302 (fuel type 1) (axial region 1) (top)AFflow area (m^2)4.361E-034.361E-03DFhydraulic diameter (m)6.604E-026.604E-02LFlength of region (m)1.143E-011.143E-01ZFentrance coeff (reentrant pipe) -- all fric losses based on Mott 5th Ed1.000E+001.000E+00Card 0302 (fuel type 1) (axial region 2 - must be fuel)AFflow area (m^2) (calc from card 305 if zero)0.000E+000.000E+00DFhydraulic diameter (m) (calc internally if 0)0.000E+000.000E+00LFlength of region (m)6.350E-016.255E-01ZFentrance coeff (reentrant pipe)1.000E+001.000E+00Card 0302 (fuel type 1) (axial region 3)AFflow area (m^2) (calc from card 305 if zero)4.361E-034.361E-03DFhydraulic diameter (m) (calc internally if 0)6.604E-026.604E-02LFlength of region (m)1.524E-011.524E-01ZFentrance+exit coeff (2 sudden expansions)2.000E+002.000E+00Note: need a conversion factor to go from assembly mass flow rate (kg/s) in PLTEMPNote: In the documentation for this work I decided to treat this conversion factor Card 0303to overall pump volumetric flow rate in gal/min (gpm)a little differently. I wrote it as wassy = 0.0630 (f/N) Qpump (gpm)FCOEFfriction factor coeff0.316where N = Nassy and f = 0.50 + 0.0088NFEXPFexponent in friction factor correlation0.250# fuel assyConv. FactorLet's use these to compute some of the same data as before (just as a check)Card 030420469.8Ptot (kW)1000Qpump (gpm)1700NCHNFnumber of coolant channels1921487.0# fuel assykW/assyfgpm/assykW/gpmwchan (kg/s)UMLRRWPIIDFfuel material option (03/01 UMLRR U3Si2-AL /WPI UAlx-Al)3122503.72050.00.67657.50.8700.2011.028E+031.107E+03IDCclad type (3 - Al 6061)323520.02147.60.68555.40.8590.1949.916E+021.068E+03UDENuraniun density (g/cc) (w/o if IDF = 0)3.420E+001.700E+0024535.92245.50.69453.60.8480.1889.587E+021.033E+03UNFUELwidth of unfueled end of fuel plate (m)2.595E-033.103E-0325551.42343.50.70251.90.8370.1829.286E+021.000E+03Llength of fuel region (m)5.969E-015.969E-0126566.52441.70.71150.40.8270.1769.011E+029.706E+02CLADclad thickness (m)3.800E-043.810E-042540.00.72049.00.8170.1718.757E+029.433E+02TAEMmeat thickness (m)5.100E-047.620E-042638.50.72947.70.8070.1678.523E+029.181E+02Card 0305 (central channels)AFFflow area (m^2)1.957E-041.817E-04Ptot (kW)1000Qpump (gpm)1650DFFhydraulic diameter (m) (calc internally if zero)0.000E+000.000E+00# fuel assykW/assyfgpm/assykW/gpmwchan (kg/s)UMLRRWPIPERFwetted perimeter (m)1.380E-011.395E-012050.00.67655.80.8970.1959.975E+021.075E+03XIFheated perimeter (m) (calc internally if zero)0.000E+000.000E+002147.60.68553.80.8850.1889.624E+021.037E+03Card 0305 (end channels)2245.50.69452.00.8740.1829.305E+021.002E+03AFFflow area (m^2)9.784E-059.083E-052343.50.70250.40.8630.1769.013E+029.709E+02DFFhydraulic diameter (m) (calc internally if zero)0.000E+000.000E+002441.70.71148.90.8520.1718.746E+029.421E+02PERFwetted perimeter (m)6.900E-026.977E-022540.00.72047.50.8420.1668.500E+029.156E+02XIFheated perimeter (m) (calc internally if zero)0.000E+000.000E+002638.50.72946.30.8320.1628.273E+028.911E+02Card 0306CIRCFwidth of fuel plate (repeat 18 times)6.604E-026.706E-02Ptot (kW)1000Qpump (gpm)1600Card 0307(don't forget dummy endplates - RP = 0 for UMLRR only)# fuel assykW/assyfgpm/assykW/gpmwchan (kg/s)UMLRRWPIFACTFRP factor by plate (NELF*18 entries) 1.001.002050.00.67654.10.9250.1899.673E+021.042E+03loop over 3 types of bypasses (if present)Note: the power to flow ratio vs Nassy is also of interest. As the # assy increases, the power2147.60.68552.20.9130.1839.332E+021.005E+03Card 0400 (bypass 1 - radiation baskets)per assy decreases (for constant Ptot), but the flow rate per assy aslo decreases. Thus, 2245.50.69450.40.9010.1779.023E+029.719E+02NCRSnumber of bypasses of this type5addressing how the power to flow (kW/gpm) changes is important.2343.50.70248.90.8900.1718.740E+029.415E+02NCnumber of axial regions12441.70.71147.40.8790.1668.481E+029.135E+02WCGESguess for flow rate (kg/s)1.00Ptot (kW)1000Qpump (gpm)17002540.00.72046.10.8680.1618.242E+028.878E+02Card 04012638.50.72944.80.8580.1578.022E+028.641E+02ACflow area (m^2)1.772E-03# fuel assykW/assygpm/assykW/gpmDChydraulic diameter (m)4.750E-022050.057.50.869Ptot (kW)1000Qpump (gpm)1400LClength of region (m)8.890E-012147.655.80.854# fuel assykW/assyfgpm/assykW/gpmwchan (kg/s)UMLRRWPIZCentrance+exit loss coeff (contract+expand)1.500E+002245.554.00.8422050.00.67647.31.0570.1668.464E+029.117E+02Card 04022343.552.30.8312147.60.68545.71.0430.1608.166E+028.796E+02FCOECfriction factor coeff3.160E-012441.750.70.8222245.50.69444.11.0300.1547.895E+028.504E+02FEXPCexponent in friction factor correlation2.500E-012540.049.20.8122343.50.70242.81.0170.1507.647E+028.238E+02Card 0400 (bypass 2 - control blades)2638.547.70.8062441.70.71141.51.0040.1457.421E+027.993E+02NCRSnumber of bypasses of this type42540.00.72040.30.9920.1417.212E+027.769E+02NCnumber of axial regions32638.50.72939.20.9800.1377.019E+027.561E+02WCGESguess for flow rate (kg/s)1.00Card 0401 (top) -- control in -- Ptot (kW)1000Qpump (gpm)1300ACflow area (m^2)2.670E-03# fuel assykW/assyfgpm/assykW/gpmwchan (kg/s)UMLRRWPIDChydraulic diameter (m)9.296E-032050.00.67643.91.1380.1547.859E+028.466E+02LClength of region (m)7.620E-012147.60.68542.41.1230.1487.583E+028.168E+02ZCentrance loss coeff (reentrant pipe)1.000E+002245.50.69441.01.1090.1437.331E+027.897E+02Card 0401 (middle) -- control out -- 2343.50.70239.71.0950.1397.101E+027.649E+02ACflow area (m^2)5.247E-032441.70.71138.51.0820.1356.891E+027.422E+02DChydraulic diameter (m)3.564E-022540.00.72037.41.0680.1316.697E+027.214E+02LClength of region (m)5.334E-012638.50.72936.41.0550.1286.518E+027.021E+02ZCentrance loss coeff (sudden expansion)2.411E-01Card 0401 (bottom) -- grid region -- ACflow area (m^2)3.288E-03DChydraulic diameter (m)2.288E-02LClength of region (m)1.524E-01ZCloss coeffs (use two sudden expansions)2.000E+00Card 0402FCOECfriction factor coeff3.160E-01FEXPCexponent in friction factor correlation2.500E-01Card 0400 (bypass 3 - regulating rod)NCRSnumber of bypasses of this type1NCnumber of axial regions3WCGESguess for flow rate (kg/s)1.00Card 0401 (top) -- control in -- ACflow area (m^2)1.007E-03DChydraulic diameter (m)8.636E-03LClength of region (m)6.604E-01ZCentrance loss coeff (reentrant pipe)1.000E+00Card 0401 (middle) -- control out -- ACflow area (m^2)3.920E-03DChydraulic diameter (m)6.261E-02LClength of region (m)6.350E-01ZCentrance loss coeff (sudden expansion)5.523E-01Card 0401 (bottom) -- grid box -- ACflow area (m^2)4.361E-03DChydraulic diameter (m)6.554E-02LClength of region (m)1.524E-01ZCentrance loss coeff (sudden expansion)1.000E+00Card 0402FCOECfriction factor coeff3.160E-01FEXPCexponent in friction factor correlation2.500E-01Card 0500DPOinitial Del P (MPa)1.000E-03DDPDel P increment (MPa)2.500E-04DPMAXmax Del P (MPa)4.000E-03POWERtotal power (MW) (for portion modelled)1.000E+00TINinlet water temp (C)4.300E+01Pinlet pressure (MPa) (about 25 psi converted to Mpa)1.724E-01Card 0600ITERmax number of iterations100CONVconvergence ctiterion0.0001ETAflow instability parameter (default=25)25.00Card 0700NNnumber of axial interfaces (nodes+1)rest not neededCard 0701for flowZRdistance of interface from inlet (normalized)distributionQVZrelative power density at interface(fz = 1.5 for hot assy)Card 0702NIJKnumber of channels for axial editsCard 0703 (repeat NIJK times)NIIelement typeNJJelement numberNKKchannel numberPrepared in June 2015 with latest update by JRWhite in July 2015 (and again in March 2017)UMLRR Flow Areas (cm^2)UMLRR Normalized Flow Areas WPI Normalized Flow Areas WPI Flow Areas (cm^2)mass flux (kg/s/m^2)mass flux (kg/s/m^2)UMLRR AssemblyWPI Assemblymass flux (kg/s/m^2)mass flux (kg/s/m^2)mass flux (kg/s/m^2)y = 8.786E-03x + 5.041E-010.6400.6600.6800.7000.7200.7400.76020212223242526Series2Series4Linear (Series4)y = 16.111x + 148.62450.0475.0500.0525.0550.0575.020212223242526Series1Linear (Series1)0.8000.8100.8200.8300.8400.8500.8600.8700.8802021222324252627Series1