Thermal Properties for Zry Cladding

Calculates the thermal conductivity and the specific heat, under constant pressure, for zirconium alloy cladding based on either the MATPRO or IAEA models.

Description

The ThermalZry material model computes the thermal conductivity and specific heat capacity for zirconium alloy cladding. Two different sets of thermal properties models are available: one based on the MATPRO relations (Allison et al., 1993) and a second set based on IAEA relations (IAEA, 2006). The MATPRO models are the default behavior in Bison.

Thermal Conductivity

Both the MATPRO model and the IAEA model for thermal conductivity used experimental data from Zircaloy-2 and Zircaloy-4.

MATPRO Model

The correlation for thermal conductivity of zirconium alloy (Allison et al., 1993) is given by: (1) where the temperature in Kelvin. This correlation is valid up to 2098 K.

IAEA Model

The IAEA thermal conductivity of zirconium alloy is (IAEA, 2006): (2) where is the temperature in K. The IAEA model is valid for temperatures between 300K and 1800K.

Specific Heat Capacity

The MATPRO and the IAEA models for specific heat capacity for zirconium alloy are both based on experimental data for Zircaloy-2 and applied to Zircaloy-4; no information on the specific heat capacity of Zircaloy-4 is available in the open literature.

MATPRO Specific Heat Capacity Model

The specific heat capacity for Zircaloy alloys is based upon the tabulated data show in Table 1. Linear interpolation is used between the tabulated values.

Table 1: Temperature dependent specific heat capacity of Zircaloy (Allison et al., 1993)

Temperature (K)Specific Heat Capacity (J/kg-K)
300.0281
400.0302
640.0331
1090.0375
1093.0502
1113.0590
1133.0615
1153.0719
1173.0816
1193.0770
1213.0619
1233.0469
1248.0356

IAEA Specific Heat Capacity Model

The IAEA model (IAEA, 2006) for specific heat capacity of zirconium alloy is a piece wise function that explicitly accounts for mixed and phases as the zirconium alloy undergoes a phase transition with temperature. The piecewise function is constructed to weight the specific heat capacity calculation based on whether the or phase is more previlant. (3) where is the temperature in K, and the calculated specific heat capacity has units of (J/kg-K). The IAEA heat capacity model is valid for temperatures between 273K and 2000K.

Example Input Syntax


[./thermal]
  type = ThermalZry
  block = 1
  temp = temp
[../]
(test/tests/tensor_mechanics/zry_oxidation_cladding/oxidation_cladding_zry_tm.i)

Input Parameters

  • thermal_conductivity_scale_factor1The scaling factor for zircaloy thermal conductivity

    Default:1

    C++ Type:double

    Description:The scaling factor for zircaloy thermal conductivity

  • computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the Material via MaterialPropertyInterface::getMaterial(). Non-computed Materials are not sorted for dependencies.

    Default:True

    C++ Type:bool

    Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the Material via MaterialPropertyInterface::getMaterial(). Non-computed Materials are not sorted for dependencies.

  • tempCoupled Temperature

    C++ Type:std::vector

    Description:Coupled Temperature

  • zry_thermal_properties_modelMATPROReference source for the set of thermal property calculations to use: MATPRO IAEA

    Default:MATPRO

    C++ Type:MooseEnum

    Description:Reference source for the set of thermal property calculations to use: MATPRO IAEA

  • boundaryThe list of boundary IDs from the mesh where this boundary condition applies

    C++ Type:std::vector

    Description:The list of boundary IDs from the mesh where this boundary condition applies

  • blockThe list of block ids (SubdomainID) that this object will be applied

    C++ Type:std::vector

    Description:The list of block ids (SubdomainID) that this object will be applied

Optional Parameters

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Description:Set the enabled status of the MooseObject.

  • use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

    Default:False

    C++ Type:bool

    Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • seed0The seed for the master random number generator

    Default:0

    C++ Type:unsigned int

    Description:The seed for the master random number generator

  • implicitTrueDetermines whether this object is calculated using an implicit or explicit form

    Default:True

    C++ Type:bool

    Description:Determines whether this object is calculated using an implicit or explicit form

  • constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

    Default:NONE

    C++ Type:MooseEnum

    Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

Advanced Parameters

  • output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)

    C++ Type:std::vector

    Description:List of material properties, from this material, to output (outputs must also be defined to an output type)

  • outputsnone Vector of output names were you would like to restrict the output of variables(s) associated with this object

    Default:none

    C++ Type:std::vector

    Description:Vector of output names were you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

Input Files

References

  1. C. M. Allison, G. A. Berna, R. Chambers, E. W. Coryell, K. L. Davis, D. L. Hagrman, D. T. Hagrman, N. L. Hampton, J. K. Hohorst, R. E. Mason, M. L. McComas, K. A. McNeil, R. L. Miller, C. S. Olsen, G. A. Reymann, and L. J. Siefken. SCDAP/RELAP5/MOD3.1 code manual, volume IV: MATPRO–A library of materials properties for light-water-reactor accident analysis. Technical Report NUREG/CR-6150, EGG-2720, Idaho National Engineering Laboratory, 1993.[BibTeX]
  2. IAEA. Thermophysical properties database of materials for light water reactors and heavy water reactors: final report of a coordinated research project 1999-2005. Technical Report IAEA-TECDOC-1496, IAEA, 2006.[BibTeX]