3.13. Subroutine Descriptions and Flowcharts

The subroutines used in the core-channel thermal hydraulics calculations are described below, grouped by the phase of the calculation where they are used.

Table 3.13.1 Steady-state Thermal Hydraulics Routines

Routine	Description
SSTHRM	Driver for the steady-state calculations:
SSCOOL	Calculates steady-state coolant, structure, reflector, and gas plenum temperature for each axial node. Also calculates coolant pressure and saturation temperature at each axial node.
SSCLM1	Calculates steady-state coolant pressures and saturation temperatures for the multiple pin option.
SSHTR	Calculates steady-state fuel and cladding temperatures in the core and blankets.
SSNULL	Driver for the steady-state null transient.

Table 3.13.2 Pre-Voiding Transient Temperature Calculations

Routine	Description
TSCL0	Driver of the pre-voiding thermal hydraulics.
TSHTRN	Driver for the pre-voiding transient temperature calculations. Calls TSHTN1, TSHTN2, TSHTN3, TSHTN4, TSHTM2, and TSHTM3 as appropriate.
TSHTN1	Calculates reflector, coolant, and structure temperatures in a reflector zone.
TSHTN2	Calculates cladding, coolant, structure, and plenum gas temperatures in the gas plenum region.
TSHTN3	Calculates fuel, cladding, coolant, and structure temperatures in the core and axial blankets.
TSHTN4	Calculates the axial node size and the coolant density and specific heat for each axial node.
TSHTN5	Called from TSHTN3 to adjust fuel and cladding temperatures to account for the heat of fusion if melting is occurring.
TSHTM2	Multiple pin version of TSHTN2.
TSHTM3	Multiple pin version of TSHTN3.

Table 3.13.3 Fuel and Cladding Temperatures During Boiling

Routine	Description
TSBOIL	Driver for the boiling module.
TSHTRV	Calculates fuel and cladding temperatures in the core and axial blankets.

Table 3.13.4 Pre-voiding Coolant Flow Rates and Pressures

Routine	Description
TSCNV1	Extrapolates coolant temperatures, computes coolant flow rate and pressures. Also tests for start of boiling. Sets variables for coupling with PRIMAR-4. Called by TSCLO every coolant step.
TSCNV2	Called by TSCNV1 to compute flow-rate coefficients xI1(JC), xI2(JC), xI3L(JC), xI3T(JC), xI4(JC) and xI5(JC).
TSCNV3	Called by TSCLO, only at the end of a heat-transfer step, to get heat-transfer coefficients Herc(JC) and Hsic(JC). Calls TSCNV8 to get liquid heat-transfer coefficients hc(JC).
TSCNV7	Computes \(\rho_{\text{c}}\), \(\mu_{\text{c}}\) for liquid sodium.
TSCNV8	Computes liquid heat-transfer coefficient \(h_{\text{c}} \left( \text{JC} \right)\).
TSCLM1	Multiple pin version of TSCNV1.

Table 3.13.5 Auxiliary Routines

Routine	Description
CFUEL	Calculates the specific heat of the fuel as a function of temperature. For coding efficiency, one call to CFUEL gives the specific heats for all radial nodes at an axial node, rather than using as separate call to CFUEL for each radial node of each axial node.
CCLAD	Calculates the specific heat of the cladding as a function of temperature. One call to CCLAD returns the values for all axial nodes in the core and axial blankets.
KFUEL	Calculates the thermal conductivity of the fuel as a function of temperature. One call to KFUEL provides the values for all radial nodes at an axial node.
KCLAD	Calculates cladding thermal conductivity as a function of temperature. One call to KCLAD returns the values for all axial nodes in the core and axial blankets.
INTRP	General interpolation routine that uses linear interpolation from a table to obtain y as a function of x. A single call to INTRP with an array of x’s will return an array of resulting y’s.
SHAPE	Multiplies the array of axial power shapes by the current time-dependent power to obtain the current power, in watts, at each axial node in a channel.
HBSMPL	Calculates the bond gap conductance at each axial node if the simple bond gap conductance model is used. Otherwise, DEFORM-4 calculates the bond gap conductance; and HBSMPL is not called.
INVRT3	Solves a tri-diagonal matrix equation of arbitrary size.

Figure 3.1.2 is a flowchart for subroutine TSCLO, Figure 3.13.1 is a flowchart for subroutine TSHTRN, and Figure 3.13.2 is a flowchart for subroutine TSHTN3. The logic in the other pre-voiding single pin transient heat-transfer routines is very simple and straightforward. Subroutine TSHTRV is similar to TSHTN3, except that TSHTRV does not calculate coolant and structure temperatures. Also, in TSHTRV the axial nodes are completely de-coupled, so the order in which they are treated is immaterial. For simplicity, the axial node loop in TSHTRV always starts at the bottom and works up.

Figure 3.13.1 Flowchart for Subroutine TSHTRN

Figure 3.13.2 Flowchart for Subroutine TSHTN3