2.8.2.11. Block 51 — INPCHN — Channel-Dependent Options and Integer Input

1

IDBUGV

Channel-dependent debug flag. (See IDBUG0 and IERSTP).

= 0, For no debugging output.

= 1, For steady-state debugging output.

Time step control output for 2 ≤ IDBUG0 ≤ 6.

= 2, Time step, also steady-state fuel behavior output if IDBUG0 = 0.

= 3, Coolant calculation output. (See IERSTP)

= 4, For more coolant debugging output. (See IERSTP )

= 5, For transient heat-transfer debugging output

= 6, Print out all reactivity changes

= 7, Primary loop calculation debug

= 8, Coolant-cladding temperature calculations only

2

IERSTP

Main time step number when debugging output starts. Relevant only if IDBUGV > 0. If IERSTP is negative, debug output starts on step |IERSTP| of a null transient (ISSNUL > 0).

3

IRHOK

Fuel density, heat capacity and thermal conductivity selection parameter.

≤ 0, Temperature dependent tabular forms. See RHOTAB, XKTAB, and CPFTAB for tables.

> 0, Correlations that depend on temperature, and in case of binary/ ternary fuel, also on composition.

Correlations for oxide fuel (IMETAL=0):

> 0, Temperature and porosity dependent correlations for fuel density. Inputs COEFDS, COEFDL are required.

= 1, Correlation for fuel conductivity from GEAP-13582.

= 2, Correlation for fuel conductivity from WARD-FTM-FI-RJS-17. Input COEFK is required.

= 3, Correlation for fuel conductivity based on COMETH3J-FBR.

Note: The parameter values associated with IRHOK=1 or 3 need not be input. For oxide fuel, tabular heat capacity is used even if IRHOK > 0.

Correlations for U-Fs fuel (IMETAL=1):

> 0, Correlation for fuel density uses input REFDEN with thermal expansion function subprogram ALPHF. Correlation for fuel conductivity depends on temperature. Porosity and sodium logging is accounted for by the method of ANL/RAS 85-19.

Note: For U-Fs fuel, tabular heat capacity is used even if IRHOK > 0.

Correlations for ternary/binary fuel (IMETAL=2,3):

Note: Input IFUELM is required for this fuel.

IFUELM:

= 0, IFR handbook data interpolations;

= 1, ANL/RAS 85-19 correlations;

= 2, Mark-V fuel or U-10Zr fuel according to input IMETAL = 2 or 3.

For IFUELM=0:

> 0, Correlation for fuel density depends on temperature. Input RHOZN is required. Correlation for fuel heat capacity depends on alloy composition/type and temperature. Inputs PUZRTP , TFSOL , TFLIQ are required. Correlation for fuel conductivity depends on composition/ type, temperature, porosity and sodium logging. Inputs PRSTY and XLOGNA are required.

For IFUELM=1:

> 0, Correlation for fuel density depends on alloy composition and temperature. Correlation for fuel heat capacity depends on alloy composition and temperature. Correlation for fuel conductivity depends on composition, temperature, porosity and sodium logging. Inputs POROSS , PORMSS , PORCSS , FPORNA are required.

For IFUELM=2:

> 0, Correlation for Mark-V & U-10Zr fuel densities. Input composition from array PUZRTP is used with thermal expansion. Correlation for Mark-V & U-10Zr fuel heat capacities depend on temperature. Correlation for Mark-V & U-10Zr fuel conductivities depend on composition, temperature and burnup. Input BURNFU is required, based on memos of Billone to Briggs dated 1/14/91, 3/8/91 & 10/21/91.

4

NPLN

Number of segments in gas plenum.

1 ≤ NPLN ≤ 6.

5

NREFB

Number of reflector zones below the pin section.

1 ≤ NREFB ≤ 5.

6

NREFT

Number of reflector zones above the pin section.

1 ≤ NREFT ≤ 5.

NREFB + NREFT ≤ 6.

7-13

NZNODE (KZ)

Number of segments in zone KZ. (KZ ≤ 7).

Total number of all segments in all zones ≤ 48. Only one segment per node is necessary, but if the LEVITATE or PLUTO2 region can extend into a node the segments there should be in the range from 0.03 meters to 0.1 meters. Neighboring cells for PLUTO2 and LEVITATE should not differ much in length. A length ratio of 1.5 is still reasonable.

14

NT

Number of radial temperature nodes within the fuel,

3 < NT < 12.

15

IFUELV

Table number of property value to be used for core fuel.

0 < IFUELV ≤ IFUEL1.

16

IFUELB

Table number of property value to be used for blanket fuel. 0 < IFUELB ≤ IFUEL1.

17

ICLADV

Table number of property value to be used for cladding table.

18

NGRDSP

Number of spacer grids in pin. 2 ≤ NGRDSP ≤ 10.

= 0, No spacer grids.

19

KTING

Fission-gas release model option.

= 0, Isotropic release model.

= 1, Weisman model.

20

NAXOP

Model selection for axial expansion and crack volume.

NAXOP	Crack Volume	Axial Plane Strain	Axial Swelling
0	No	No	No
1	No	Yes	No
2	No	No	Yes
3	No	Yes	Yes
4	Yes	No	No
5	Yes	Yes	No
6	Yes	No	Yes
7	Yes	Yes	Yes

If NAXOP = X with X above, mixed plane strain.

= 3X with X above, cladding controlled axial plane strain.

= 2X with X above, then free fuel axial plane strain.

= 1X constrained plane strain.

21

MSTEP

≤ 8, Number of main steady-state power change and constant power intervals. One required for each power change. One required for each constant power interval.

> 8, Table ID for table defining normalized power and flow per pin during pre-transient characterization.

An example table can be found here.

22

ITAU

Irradiation induced cladding swelling incubation parameter options.

=-2, No cladding swelling.

=-1, Lower limit.

= 0, Nominal value.

> 0, Upper limit.

23

IRATE

Irradiation induced cladding swelling rate options.

=-2, No cladding swelling.

=-1, Lower limit.

= 0, Nominal value.

> 0, Upper limit.

24

IHGAP

Fuel-cladding gap conductance selection.

= 0, SAS3D method for calculating HB (See AHBPAR, BHBPAR, CHBPAR, HBMAX, HBMIN, and HBPAR).

= 1, SAS4A method (Ross-Stoute model).

Note: When using the simple bond gap conduction model (ISSFU2 = -1), set IHGAP =0.

25

NPIN

Number of pins per subassembly.

26

NSUBAS

Number of subassemblies in channel.

27

MZUB

Number of segments in upper blanket.

28

MZLB

Number of segments in lower blanket.

29

IHEX

= 0, No hexadecimal printout of the sum of the coolant temperatures.

> 0, Hexadecimal printout of the sum of the coolant temperatures.

30

IRELAX

Stress relaxation options. Use IRELAX = 0.

= 0, No relaxation.

= 1, Creep law estimate, not operational.

= 2, Exact solution, not operational.

31

NGRAIN

Model selection for grain growth theory.

= 0, Limited grain growth theory.

> 0, Unlimited grain growth theory, where NGRAIN is the grain diameter exponent in Eq. (8.3-9).

Suggested value: 4.

32

ISSFUE

= 0, Bypass dynamical model calculations in steady-state, no fuel restructuring and no deformation of fuel and cladding.

= 1, Use dynamic calculation of DEFORM in steady-state.

33

IRAD

Not currently used.

34

ILAG

= 0, Use Eulerian coolant temperature calculation until flow reversal.

= 1, Use Lagrangian coolant temperature calculation from the start.

Suggested value: 0.

35

NOSTRN

Option to avoid radial strain in the cladding even if conditions would produce strain.

= 0, Strain cladding if conditions warrant.

= 1, Allow no cladding strain.

36

JRUPT

Not currently used.

37-44

NPLIN (M)

Number of subdivisions in each of the MSTEP divisions. There should be enough subdivisions to allow for reasonable feedback of the restructuring into the thermal calculation. At least 4 subdivisions for each power change and enough constant power subdivisions so that each does not exceed about 10 days (preferably 2 days during early irradiation).

(M = 1, MSTEP)

45

IROR

Controls assumption used when molten cavity extends to the cracked fuel zone.

= 0, Cavity pressure reduced by R/R before acting on cladding, no crack volume included in cavity.

= 1, Cavity pressure acts on the cladding, all crack volume included in cavity if molten to cracked region.

46

JPRNT1

The lowest axial node for which debug output is produced from DEFORM.

47

JPRNT2

The highest axial node for which debug output is produced from DEFORM.

Note: The debug output is produced for all nodes from JPRNT1 to JPRNT2 inclusive. If these two values are set to 0 but the time step controls and IDBUGF are activated, then the debug output is from the molten cavity routine and the axial expansion and feedback calculation only (the axial node independent part of the calculation).

48

NNBUG1

The time step at which to start the debug output from DEFORM.

49

NNBUG2

The time step which is the last time step for debug output from DEFORM

50

IDBUGF

The control for the type of debug output desired from DEFORM.

= 3, NNBUG1 and NNBUG2 refer to the steady-state time step.

=13, NNBUG1 and NNBUG2 refer to the transient time step.

=23, NNBUG1 and NNBUG2 refer to both the steady-state and the transient time steps.

=-1, Initiates a print option that writes out when DEFORM (the transient DEFORM driver) is entered and left beginning at transient time step NNBUG1.

=-2, This is a special option that produces standard DEFORM output at every time step from NNBUG1 to NNBUG2 inclusive. NNBUG1 and NNBUG2 must be entered as negative numbers. This is active only in the transient calculation.

51-58

NSKIP (M)

Print control for each of the MSTEP divisions in DEFORM. There is one value for each NPLIN value.

< NPLIN, the DEFORM results are printed after each NSKIP subdivisions.

= NPLIN, the DEFORM results are printed at the end of the MSTEP division.

> NPLIN, no DEFORM results are printed.

(M = 1, MSTEP)

59

MPL1

Plots channel pressure history at node MPL1. The nodes MPL1 - MPL7 must be a node in the zone KZPIN.

60

MPL2

Plots channel pressure history at node MPL2.

61

MPL3

Plots channel pressure history at node MPL3.

62

MPL4

Plots channel temperature history at node MPL4.

63

MPL5

Plots channel temperature history at node MPL5.

64

MPL6

Plots channel temperature history at node MPL6.

65

MPL7

Plots channel coolant volumetric flow rate at node MPL7. (MPL7 should correspond to the location of upper flowmeter in experiments, the volumetric flowrate of lower coolant slug is plotted with it.)

66

MPL8

Not currently used.

67

MPL9

Not currently used.

68

KKSBTP

Not currently used.

69

KKSBRI

Not currently used.

70

NRPI

Not currently used.

71

NRPI1

Number of pins per assembly that are assumed to fail (number of unfailed pins is NPIN - NRPI1) when the fuel-pin failure option MFAIL is satisfied.

72

NRPI2

Not currently used.

73

NRPI3

Not currently used.

74

IPSIZE

Fuel particle size option in PLUTO2 and LEVITATE.

= 1, RAFPLA is used from the initial fuel ejection to time TIFP.

= 2, RAFPSM is used all the way.

75

IBUGPL

Debug flag, should currently always be 0.

76

ICFINE

= 0, Automatic time step selection in LEVITATE and PLUTO2 using DTPLIN as the initial and later on the minimum time step.

= 1, Initial time step DTPLIN is used all the time if the main and primary loop time steps remain multiples of DTPLIN. If this is not the case, then the PLUTO2 and LEVITATE time steps are temporarily smaller than DTPLIN.

77

IPRINT

Should currently always be 0.

78

IPLOT

> 0, Plotting file (Logical units 12, 13 and and 14) for PLUTO2 is written every IPLOT milliseconds if IPLOT is greater than zero.

= 0, No plotting data is saved.

79

IBGO

Debug printout starting at LEVITATE and PLUTO2 cycle IBGO.

80

IBSTOP

PLUTO2 or LEVITATE printout debug ends at cycle IBSTOP.

81

IBNEW

Debug levels:

= 0, No debug output.

= 1, Time step debug.

= 2, Some debug output.

= 3, More debug output.

= 4, Lots of debug output.

82

IPGO

Between cycles IPGO and IPSTOP regular full LEVITATE or PLUTO2 output at every IPNEW cycles.

83

IPSTOP

See IPGO.

84

IPNEW

See IPGO.

85

ICLADB

= 0, Cladding is ignored after melting.

= 1, Cladding motion occurs.

= 2, Cladding motion does not occur, but heat transfer to molten cladding does occur.

86

MFAIL

Fuel-pin failure option. (See FSPEC).

= 0, No pin failure.

= 1, FSPEC is failure time.

= 2, FSPEC is fuel failure temperature.

= 3, FSPEC is fuel mass melt fraction at failure.

= 4, FSPEC is cavity pressure at failure.

= 5, FSPEC is cladding yield stress at failure. Not yet operational.

= 6, Eutectic penetration correlation for U-5fs fuel. FSPEC is eutectic temperature.

= 7, Failure criteria consistent with PLUTO2/LEVITATE rip propagation model. Functional ultimate tensile strength used and fully cracked fuel assumed.

= 8, Eutectic penetration and stress-based failure for metal fuel.

= 9, Failure time based on melt fraction at FSPEC with location based on maximum MFAIL = 7 failure criterion.

87

IFAIL

Relevant for MFAIL = 2, the radial node to test.

88

JFAIL

Relevant for MFAIL = 1,2,3,4,5, the axial fuel-pin node to test.

Note: If IFAIL or JFAIL are not specified, the peak value of the pertinent failure quantity is used. If MFAIL = 1, JFAIL must be specified.

89

ISUBAS

Subassembly number, only required for the detailed coolant sub-channel model.

90

JCLN

Axial heat-transfer segment, JCLN, output on the plotting unit for fuel information, 1 ≤ JCLN ≤ MZ.

91

JNEN

Axial heat-transfer segment, JNEN, output on the plotting unit for the cladding temperature, 1 ≤ JNEN ≤ MZ.

92

JNCN

Axial coolant node, JNCN, output on the plotting unit for coolant temperature, 1 ≤ JNCN ≤ MZC.

93

JNSN

Axial coolant node, JNSN, output on the plotting unit for structure temperature, 1 ≤ JNSN ≤ MZC.

94

JRPRO

Axial heat-transfer segment, JRPRO, 1 ≤ JRPRO ≤ MZ, output on the plotting unit for total radial profile (written each time the mass averaged fuel temperature, TBAR, increases by DTFUEL degrees). Not currently operational.

95

IPSIG

Hydrostatic pressure for fuel swelling:

= 1, Use SIGR.

= 2, Use (SIGR + SIGC)/ 2.

= 3, Use (SIGR + SIGC + SIGZ)/ 3.

96

IHTPRS

= 1, Hot pressing of fuel to PRSMIN, not yet operational.

= 0, No hot pressing.

(Recommended IHTPRS = 0)

97

IPRD

Controls the amount of DEFORM output in the transient calculation.

(For steady-state control see NSKIP)

= 0, Only short form output.

= 1, Radial stress, total porosity, crack volume + short form.

= 2, All above + circumferential stresses retained fission-gas distribution, fission-gas porosity.

= 3, All above + axial stresses, grain size distribution, radial mesh locations.

98-107

IDBFLG

IDBFLG(4) > 0 gives reentry temperature debug print. (10) IDBFLG(5) > 2 gives TSCA debug print.

IDBFLG(6) > 0 gives subassembly-subassembly heat transfer prints.

IDBFLG(7) > 0 for film motion of debugs.

IDBFLG(8) > 0 for Wallis flooding correlation debugs.

IDBFLG(9) > 0 for sub-channel analysis coolant debugs

IDBFLG(10) > 0 for sub-channel analysis heat transfer debugs

108-117

IDBSTP

TSCA debug print starts at coolant step IDBSTP(5).

118

IEQMAS

Radial fuel mesh size assumption.

= 0, Equal radial distances between points at which temperatures are calculated.

> 0, Equal cylindrical areas associated with each radial temperature node.

119

IBLPRN

Number of coolant dynamics time steps between boiling printouts.

(See IPO, IPOBOI, and IBLPRT).

Default: 100.

120

IDBGBL

IDBGBL > 3 gives boiling debug print.

121

IDBLST

Boiling debug print starts at coolant step IDBLST.

122

ISSFU2

= 0, Omit DEFORM calculation during the transient. (See IAXEXP).

= 1, Use DEFORM calculation during the transient.

=-1, Simple bond gap conduction model.

123

IHEALC

= 0, No crack healing.

= 1, Crack healing based on fuel swelling rate parameters. (Not operational)

= 2, Cracking healing based on 100% healed if temperature above FTMPCH*TMF(IFUEL).

124

IAXTHF

Determine components active in the axial expansion calculation in DEFORM.

= 0, Only thermal effects.

= 1, Thermal and force effects.

125

IDCLGO

Value of ICOUNT (number of cladding time steps) when cladding debug print begins.

126

IDCLSP

Value of ICOUNT when CLAP debug print ends.

127

IDCLDE

CLAP debug is printed after every IDCLDE call to the CLAP module.

128

IFILM

Number of nodes dried out before switch from wet (a few boiling segments) to dry (larger boiling length) minimum film thickness.

Suggested value: 3 or 4. Do not input 0.

129-131

NZONF

Not currently used.

132-155

IFUELI

Not currently used.

156-179

NODSUM

Not currently used.

180

IFUOPT

Not currently used.

181

IAXEXP

Simple axial expansion reactivity feedback model. Can not be used with DEFORM-4. See ISSFUE and ISSFU2. See MODEEX for model choice.

= 0, No simple axial expansion feedback.

= 1, Calculate feedback.

= 2, Calculate and print simple axial expansion feedback for each channel.

182

IMOMEN

LEVITATE option referring to the convective momentum flux formulation.

= 0, Central formulation (recommended).

= 1, Upstream formulation.

183

JSTRDX

Axial node number in structure corresponding to the above core load pad. Use only if IRADEX = 1,2,3 or -1,-2,-3.

184

IFAE

Fuel adjacency effect in Kramer-DiMelfi cladding failure model.

= 0, No

= 1, Yes

185

ICLADK

Cladding thermal conductivity option.

= 0, Table lookup.

= 1, Functional form cladding conductivity option.

(Not currently operational, use = 0)

186

IFRFAC

= 0, Turbulent and laminar friction factors only.

= 1, Add turbulent and laminar friction factors.

= 2, Turbulent, transition, and laminar friction factors.

AFR, BFR, AFLAM, RELAM, and RETRAN are required when IFRFAC = 2. Boiling, multiple pin, and subchannel models cannot be used with IFRFAC = 2.

187

IRDEXP

Option to use or skip this channel in the radial expansion feedback calculation. Only used when IRADEX \(\neq\) 0.

= 0, Use this channel.

= 1, Skip this channel.

188

IBUGPN

Debug flag, should currently always be 0.

189

IMETAL

Indicates fuel type.

= 0, Oxide fuel.

> 0, Metal fuel.

= 1, Uranium-fissium metal fuel.

= 2, U-Pu-Zr ternary alloy fuel.

= 3, U-Zr binary alloy fuel.

Note that MFUEL (IFUELO= 2) requires IMETAL= 2 or 3.

190

IPNPLT

=1, The PINACLE module produces a printer plot of the axial fuel distribution with every full printout.

191

IFUELO

Option for annular zone formation model in U-Pu-Zr alloy fuel (for the case of IMETAL = 2 only).

= 0, User input zonal compositions and radii (See IFUELC, IZNC, IZNM, MFTZN, RIZNC, and RIZNM. See also PUZRTP, RHOZN, XLOGNA, PRSTY, TFSOL, and TFLIQ).

= 1, (Not functional) Zonal composition and radii computed using the SSCOMP physical model (See TTRANM, TTRANC, POROSS, PORMSS, PORCSS, FPORNA, RHOREF, WUREF, WPUREF, WZRREF, PUBYU, CPCM, CPMO, CZCM, CZMO, CUCM, CUMO, EPSMS, EPSCOM, IDSSC, and BURNFU).

= 2, Fuel composition and geometry computed by MFUEL physical model.

192

IFUELM

Option for the thermal properties of the U-Pu-Zr alloy fuel (for the case of IMETAL > 1 only).

= 0, Use properties interpolated from the IFR metallic fuels handbook data.

= 1, Use properties in the report ANL/RAS 85-19.

= 2, Use Mark-V fuel or U-10Zr fuel properties (depending on IMETAL = 2 or 3) based on the metallic fuels handbook data.

= 3, Use thermal conductivity correlation recommended for use with MFUEL (MFUEL must be active, IFUELO= 2).

193

IFUELC

User input fuel zone specification flag.

= 0, Single radial fuel zone (See IFUELV, IFUELB, MZLB, MZUB).

= 1, Multiple radial fuel zones (See IZNC, IZNM, MFTZN, RIZNC, RIZNM).

= 2, Each radial mesh interval is a unique fuel zone for IMETAL > 1, with input mesh point porosities and compositions determined in interface routine LIFEIF.

194

IPNGO

Between cycles IPNGO and IPNSTP a full PINACLE output is obtained at cycle intervals containing IPNNEW cycles.

195

IPNSTP

See IPNGO.

196

IPNNEW

See IPNGO.

197

IEUTOPT

MFUEL eutectic penetration option.

= 0, Mechanistic eutectic penetration, default axial expansion model

= 1, Empirical eutectic penetration, default axial expansion model

= 2, Mechanistic eutectic penetration, conservative axial expansion model

= 3, Empirical eutectic penetration, conservative axial expansion model

198-202

IDM51

Not currently used.

203

IDKCRV

Power or decay heat curve for this channel.

204

ITP20

Fuel, cladding, and coolant temperatures and voiding data are written on unit 20 every ITP20 time steps. If ITP20 ≤ 0, no unit 20 output.

CHANNEL-TO-CHANNEL HEAT TRANSFER

205

NCHCH

Number of other channels that this channel is in contact with for duct wall-to-duct wall heat transfer. Maximum 8.

If NCHCH < 0, Q(ICH to JCH) = -Q(JCH to ICH)

206-213

ICHCH

Channel number of the K-th channel that this channel is in contact with. See also HACHCH.

If ICHCH < 0, then -ICHCH is a bypass channel number.

If ICHCH(1) < -8, then -ICHCH(K) is the temperature of a constant temperature heat sink in axial zone K.

If ICHCH > 500000, transfer heat from structure of ICH to coolant of ICHCH - 500000.

If ICHCH > 750000, transfer heat from coolant of ICH to coolant of ICHCH - 750000.

MULTIPLE PIN OPTION

214

JJMLTP

Multiple pin option.

=0, No multiple pin treatment for this subassembly. Use single-pin model.

>0, This is the first of JJMLTP channels used to represent the subassembly.

<0, This is one of the additional channels used to represent the subassembly.

Notes:

1) JJMLTP and ICHCH refer to different phenomena: Intra-subassembly heat transfer and inter-subassembly heat transfer.

2) JJMLTP is used for intra-subassembly coolant-to-coolant heat transfer from channel I to I + 1 and from I to I - 1. UACH1 and UACH2 must be supplied to determine the heat transfer coefficients.

3) If M channels are used to represent a subassembly, then they must be consecutive channels, starting with channel ICHN and going to channel ICHN + M - 1. For channel ICHN, JJMLTP = M. For channels ICHN + 1 through ICHN + M - 1, JJMLTP is negative.

4) The maximum value of M is 56 (no more than 56 channels can be used to represent a subassembly).

5) The axial zones outside the pin section go only with the first channel (channel ICHN). In these zones, the coolant flow area per pin and the reflector and structure perimeters per pin must be based on the number of pins in channel ICHN. The reflector zones are ignored for channels ICHN + 1 to ICHN + M - 1.

6) ICHCH refers to subassembly-to-subassembly heat transfer from the current channel to channel ICHCH. This heat transfer is from outer surface node to outer surface node. ICHCH > 1000 is an exception. This exception was included for modeling the thimble flow region of the XX09 subassembly in EBR-II. If ICHCH ≠ 0 then the heat transfer coefficient times area per unit height is specified by HACHCH.

DETAILED COOLANT SUB-CHANNEL MODEL

215-218

JCHMPN (K)

Other channel this coolant sub-channel is in contact with. Used only if ISCH = 1, JJMLTP ≠ 0

Reserved for the detailed coolant sub-channel model.

219

NUMKLT

Number of other channels this coolant sub-channel is in contact with. Used only if ISCH = 1, JJMLTP ≠ 0

Reserved for the detailed coolant sub-channel model.

220

KSWIRL

Value of K in JCHMPN(K) for swirl flow. Normally KSWIRL=0 except in an edge or corner sub-channel. Used only if ISCH = 1, JJMLTP ≠ 0

Reserved for the detailed coolant sub-channel model.

221

NULST1

Number of time steps in a null transient for this subassembly or group of channels. Used only if ISCH = 1, JJMLTP ≠ 0.

Reserved for the detailed coolant sub-channel model.

SSCOMP METAL FUEL BEHAVIOR MODEL

222

IPORC

Options to control open porosity considerations in metal fuel axial swelling during transient calculation.

= 0, Close open porosity by fraction given in FCLOP before axial movement takes place.

= 1, Treat open porosity like closed porosity.

223

IDSSC

MFUEL (IFUELO= 2) detailed print option.

= 0, print only at end of steady-state.

= 1, print at every coarse time step specified by the user.

END SSCOMP METAL FUEL BEHAVIOR MODEL

224

IDRY

= 1, For film motion.

225

ICTYPE

Cladding material indicator.

= 1, 316

= 2, D9

= 3, HT9

= 4, 15-15Ti (only if IDEFOPT > 0)

ICTYPE = 1, 2, and 3 is used for metal fuel models (DEFORM-5 and FPIN). In DEFORM-4, the default cladding is 316. If IDEFOPT > 0 and ICTYPE = 4, 15-15Ti is selected.

SUBASSEMBLY-TO-SUBASSEMBLY HEAT TRANSFER

226-233

IOPCHC (K)

Subassembly-to-subassembly heat transfer option.

= 0, Q(JCH,ICH) = - Q(ICH,JCH)

= 1, Q(JCH,ICH) not set by Q(ICH,JCH)

234-249

KTRANC

Not currently used.

250-273

KTRANM

Not currently used.

274

IPRSKP

Skip the prints for channel ICH if IPRSKP is nonzero. To be used to reduce the volume of printed output from TSPRNT.

275

ITREAT

For use with oxide fuel TREAT experiment analysis.

= 1, Completely crack fuel at start of transient.

276

IOPPL

Not currently used.

277

IBUBND

= 0, Form bubbles at node mid-points.

= 1, Form bubbles at node boundaries.

278

IGASRL

= 0, No gas release.

= 1, Gas release in the boiling model, requires DEFORM-5 life fraction calculation and FRUPT(1).

= 2, Gas release in the boiling model by pin group at the times given by TMFAIL at axial node INDFAL.

279

IGRLTM

Not currently used.

280

IRAPEN

Rapid eutectic formation rate assumption.

= 0, Use exponential eutectic formation rate formula for all temperatures (i.e. no eutectic formation temperature threshold and no rapid eutectic penetration allowed).

= 1, Allow rapid eutectic penetration if molten fuel is available (i.e. cladding temperature above fuel solidus temperature) and temperatures are in rapid penetration region (1353K to 1506K).

> 1, Use exponential eutectic formation rate formula above 923K (650 deg. C) and allow rapid eutectic penetration in the range from 1353K to 1506K (EBR-II Mark-V safety case).

281

LCHTYP

Core channel designator.

= 1, Fuel channel.

= 2, Reflector channel.

= 3, Control rod channel.

282

IGSPRS

Controls fission gas assumptions used with metal fuel pins.

= 0, Equilibrium maintained between plenum and axial segments of fuel by instantaneous gas relocation.

= 1, After initiation of equilibrated pressure during first transient time stop, gas remains in axial segments and may locally pressurize.

283

INDFAL

Failure node number on the fuel pin axial mesh for pin failure and gas release if IGASRL=2.

284

IAXCON

Axial coolant heat conduction option.

= 0, No axial heat conduction in the coolant.

= 1, Axial heat conduction in the coolant. Only operative with the multiple pin option (JJMLTP not equal to 0).

FPIN-2 INPUT

285

IFPIN2

= 0, Do not use FPIN2 metal fuel model.

= 1, Use FPIN2 metal fuel model.

No other data required when IFPIN2=0.

286

IFPI01

= 0, Use FPIN2 in interfaced mode.

= 1, Use FPIN2 in stand-alone mode.

No other data required when IFPI01=1.

287

IHTFLG

= 0, Bypass FPIN2 heat transfer calculation.

= 1, Include FPIN2 heat transfer calculation (for debugging purposes only).

288

LHTOPT

= 0, Perform heat transfer calculation including coolant and structure.

= 1, Perform heat transfer calculation with input values of cladding outer surface temperature.

This input is required only for IHTFLG=1.

289

LCRACK

Fuel cracking option.

= 0, No cracking.

= 1, Radial cracks included.

290

LFPLAS

Option for creep-plastic strains in fuel.

= 0, Allow creep-plastic strains.

= 1, Suppress creep-plastic strains.

291

LCPLAS

Option for creep-plastic strains in cladding.

= 0, Allow creep-plastic strains.

= 1, Suppress creep-plastic strains.

292

LFSWEL

Option for swelling-hotpressing strains in fuel.

= 0, Allow swelling-hotpressing strains.

= 1, Suppress swelling-hotpressing strains.

293

LCSWEL

Option for swelling strains in cladding.

= 0, Allow swelling strains.

= 1, Suppress swelling strains.

294

LLRGST

Option for strain analysis.

= 0, Large strain analysis.

= 1, Small perturbation analysis.

295

LFCSLP

= 0, Fuel-cladding locked when gap is closed.

= 1, Independent fuel-cladding axial displacement.

296

LOUTSW

Print option.

= 0, No detailed printing of results - summary only.

= 1, Normal detailed printout under LFREQA, MFREQA, and LFREQB control.

297

LFREQA

Initial print frequency, number of time steps between normal detailed printout.

298

MFREQA

Total number of time steps under LFREQA control.

299

LFREQB

Final print frequency.

300

LGRAPH

Graphics file option.

= 0, Do not write graphics file.

= 1, Write a graphics datafile.

301

LDBOUT

= 0, No debug output.

= 1, Add debug output to regular LOUTSW=2 output.

302

LDBSTP

= 0, FPIN2 calculation stops when molten cavity freezes.

= 1, Ignore this program stop.

303

LDBFPL

= 0, Use recommended fuel flow stress.

= 1, Simple power law fuel creep (EPSDOT=C0*SIGE**C1).

See XFPLC0 and XFPLC1.

304

LDBFDV

= 0, Use recommended fuel swelling - hot pressing. (Fuel swelling option for metal fuel is the simple grain boundary swelling model, ANL-IFR-27 and ANL/RAS 83-33).

= 1, Use equilibrium swelling model (ANL-IFR-6 & -23).

= 2, Simple power law fuel swelling (EPVDOT=C0*SIGM**C1).

See XFDVC0 and XFDVC1 .

305

LDBCPL

= 0, Use recommended cladding flow stress.

= 1, Ideal plastic flow for cladding (SIGY=C0+C1*EPBAR).

See XCIPL0 and XCIPL1 .

= 2, Use high-temperature power-law creep.

= 3, Use simple power law cladding creep

(EPSDOT=C0*SIGE**C1).

See XCIPL0 and XCIPL1 .

306

LGPRES

Not currently used.

307

LGAPCL

= 0, Use fuel-cladding opening/closure model.

= 1, Fuel-cladding gap always closed.

308

LCPROP

= 0, Use material property correlations.

= 1, Use temperature independent material properties.

(SAS thermal properties are used for IFPI01=0).

309

LSKIPM

= 0, Perform mechanical calculations.

= 1, Bypass mechanical calculations.

310

LGCLOS

= 0, Use gap closure routine at 100% fuel melting.

= 1, Do not close gap (if open) at 100% fuel melting.

311-334

LDBOTA (J)

Axial debug print vector.

= 0, No print.

= 1, Print.

335-345

LDBOTF (IF)

Fuel radial debug print vector.

= 0, No print.

= 1, Print.

346-348

LDBOTC (IC)

Clad radial debug print vector.

= 0, No print.

= 1, Print.

349-360

NONE

Reserved.

PRIMAR-4 MULTIPLE INLET/OUTLET PLENA

361

NSEGMP

PRIMAR-4 segment number to which this channel is assigned in the multiple inlet/outlet plena model (See IFMIOP).

CONTROL ROD DRIVE FEEDBACK

362

ICHUIS

= 0, Coolant from channel ICH is included in the upper internal structure temperature calculation for control rod drive expansion reactivity.

= 1, This channel is not used.

PINACLE

363

LQSLTP

= 0, Wall friction is not considered in calculating the velocity of the sodium slug above the fuel.

= 1, Wall friction is considered.

POWER AND REACTIVITY MESH

364

IPOWRZ

Axial power shape input option.

= 0, Enter channel axial power shape in array PSHAPE on the MZ axial mesh; JMAX=24.

= 1, Enter channel axial power shapes for core and axial blanket fuel in arrays PSHAPC (core fuel power shape) and PSHAPB (blanket fuel power shape) on the MZC axial mesh; JMAX=48.

365

IREACZ

Axial reactivity worth input option.

= 0, Enter channel axial reactivity worth for Doppler, coolant void, cladding motion, fuel motion, and structure motion in arrays WDOPA, VOIDRA, CLADRA, FUELRA, and STRCRA on the MZ axial mesh; JMAX=24.

= 1, Enter channel axial reactivity worth for Doppler, coolant void, cladding motion, fuel motion, and structure motion in arrays WDOPA, VOIDRA, CLADRA, FUELRA, and STRCRA on the MZC axial mesh; JMAX=48.

FUEL ZONE TYPE ASSIGNMENT

366-389

IZNC (J)

Outermost radial mesh interval of central zone at axial segment J. Maximum value = NT. See also RIZNC. Used only for IFUELC= 1.

390-413

IZNM (J)

Outermost radial mesh interval of middle (intermediate) zone at axial segment J. Maximum value = NT. See also RIZNM. Used only for IFUELC= 1.

414-485

MFTZN (L,J)

Fuel type (IFUEL) assignment to radial zones at axial segment J. Maximum value = 8. Used only for IFUELC= 1.

L=1, fuel type assigned to central (inner) zone. L=2, fuel type assigned to middle intermediate) zone. L=3, fuel type assigned to outer zone.

Note: A maximum of three radial zones may be specified at each axial level. Fewer than three zones may also be specified, with a minimum of a single zone assigned to all the fuel in the pin at a given axial location. Zones are assigned assuming azimuthal symmetry; the central (inner) zone begins at the fuel centerline and extends outward radially through radial mesh interval IZNC(J). The middle (intermediate) zone begins at radial mesh interval IZNC(J)+1 and extends outward radially through radial mesh interval IZNM(J). The outer zone begins in radial mesh interval IZNM(J)+1 and extends outward radially to the fuel surface (radial temperature node NT). The central zone may be eliminated by setting IZNC(J) = 0. The middle zone may be eliminated by setting IZNM(J) = IZNC(J). The outer zone may be eliminated by setting IZNM(J) = NT. Both inner and middle zones may be eliminated by setting IZNC(J) = IZNM(J) = 0. The central zone is present only if IZNC(J) > 0. The middle zone is present only if IZNM(J) > IZNC(J). The outer zone is present only if IZNM(J) < NT. The MFTZN array assigns a fuel type to a zone; fuel types and zones have a one- to-one correspondence at each axial level. The fuel types assigned here specify the fuel thermo-physical properties to be used in the solution of the fuel pin heat transfer equations. See IFUELC.

486

IPINFG

Metal fuel fission gas model flag.

= 0, Use DEFORM-5 fission gas formulation.

> 0, Use PINACLE fission gas formulation.

487

IPINRE

PINACLE ejected fuel re-entry flag.

= 0, Re-entry of fuel ejected into above-core region not permitted.

> 0, Allow re-entry into active core region of fuel ejected into above-core region.

488

IPORFG

Density correction to porosity in DEFORM-5 fission gas model.

= 0, Neglect transient temperature/density impact on open porosity volume.

> 0, Adjust open porosity volume accounting for transient temperature/density changes.

489

IPRSS1

Steady state fuels characterization initial values print in subroutine SSINC1.

= 0, Print initial values for pointwise fuel mass and porosity; Pu content for IMETAL = 2; Zr content for IMETAL = 2 or 3; Fe and Ni content for IMETAL = 2 and IFUELO = 1; Na content for IMETAL > 0.

= 1, No prints.

490

IMKVPL

EBR-II Mark-V safety case plotting data.

= 0, Do not save data for post-processing.

> 0, Compute and save maximum fuel temperature in low Zr zone, maximum inner cladding temperature, and minimum coolant saturation temperature in the channel on each main time step (Entry TSPLT1 in subroutine TSPLOT). Compute and save reactor power, flow, and power-to-flow ratio on each main time step (Entry TSPLT3 in subroutine TSPLOT).

491

MZCHCH

See FED.

492

KZEMFM

See FED.

493

MTREAT

TREAT fuel channel modeling flag.

= 0, Use standard fuel channel models.

> 0, Use special TREAT fuel channel models. See also CFLAT and FFLAT. RBR, RER, and ROUTFP contain TREAT fuel assembly half-thickness dimensions. In fuel and cladding heat transfer calculations, the correct periphery will be used in place of the circumference. The structure field will be eliminated. The fission gas plenum will be eliminated. Air properties will be used for the coolant. The coolant pressure drop calculation will be eliminated. One pin per channel will be used. The DEFORM, PRIMAR, and multiple pin models are not allowed. Only equally-spaced radial heat transfer mesh (IEQMAS=0) is allowed.

494

IFLOOD

= 0, Use FVAPM to determine the friction factor multiplier for the vapor.

= 1, Use Wallis flooding correlation to determine whether to use the Wallis friction factor multiplier for the vapor. This option should only be used with the film motion model, IDRY = 1.

DETAILED COOLANT SUB-CHANNEL MODEL

ISCH = 1 and JJMLTP ≠ 0

495

IFT24

Output detailed coolant sub-channel model variables on fort.24 every IFT24 time steps. If IFT24 = 0, no output on fort.24.

Reserved for the detailed coolant sub-channel model

496

ILATF

Always use ILATF = 0 to include lateral flow terms in the momentum equation for the detailed coolant sub-channel model.

Reserved for the detailed coolant sub-channel model.

497

IDEFOPT

= 0, use default DEFORM-4 models

> 0, use extended DEFORM-4 models.

Note that the extended DEFORM-4 models will be fully available at later versions. Currently, only 15-15Ti cladding models (ICTYPE = 4) are supported.

498

IDEFSTFAL

= 0, do not use stochastic clad damage evaluation model

> 0, use stochastic clad damage evaluation model as a sidebar analysis. Only applicable when IDEFOPT > 0 and ICTYPE = 4.

499

IDMICH

Not currently used.

500

LDETL

Visualization Data Detail Flag:

= 0, Do not write visualization data for this channel

= 1, Write coolant temperatures only

= 2, Coolant and duct wall (structure) temperatures

= 3, Coolant, structure, and bulk fuel temperatures

= 4, Coolant, structure, average fuel, and detailed radial fuel temperatures.

501

NFT24

Number of items to be output on fort.24

Reserved for the detailed coolant sub-channel model.

502-521

JCFT24 (K)

Axial node for the Kth output on fort.24.

Reserved for the detailed coolant sub-channel model.

522-541

ITYP24 (K)

Variable type for the Kth output on fort.24

= 1, PCBAR2 (coolant pressure, middle of the node)

= 2, W2RT2 (coolant flow rate)

= 3, TCOOL2 (coolant temperature at the bottom of the node)

= 4, TCBAR2 (coolant temperature at the middle of the node)

= 5, PCOOL2 (coolant pressure, bottom of the node)

= 6, WLAT2(1) (coolant lateral flow rate to first adjacent subchannel

= 7, WLAT2(2) (coolant lateral flow rate to second adjacent subchannel

= 8, WLAT2(3) (coolant lateral flow rate to third adjacent subchannel

= 9, WLAT2(4) (coolant lateral flow rate to fourth adjacent subchannel

=10, NITER (number of iterations for the last solution in SOLVIT

Reserved for the detailed coolant sub-channel model.

542

NULPT1

Print results every NULPT1 steps for the null transient (see NULST1)

Reserved for the detailed coolant sub-channel model.

545

ZoneLossCoefTableID

Table ID referencing channel-specific anisotropic Re-dependent loss coefficient data.

An example table can be found here.

546-600

IDMICH

Not currently used.