Zircaloy Irradiation Growth Eigenstrain

Calculates eigenstrain from irradiation growth in Zircaloy cladding using either the Franklin or ESCORE models

Description

This class, ZryIrradiationGrowthEigenstrain, calculates the zircaloy strains that are not caused by external force loading on the cladding material: here the strain due to irradiation growth. This model for irradiation growth of Zr4 cladding calculates the increment of cladding axial strain due to irradiation growth, and was taken from Franklin (1982). (1) where is a coefficient determined by the zircaloy material type, units of (n/(cm-sec)), is the total fast neutron fluence, and is the fluence exponent (dimensionless). An additional set of material properties, using the same irradiation growth strain relationship, is from Rashid et al. (2004).

Table 1: Material Parameters for Irradiation Growth Strain for Various Zircaloy Microstructures (Franklin, 1982; Rashid et al., 2004)

Zircaloy Material TypeIrradiation Growth Strain Coefficient (A)Irradiation Growth Exponent (n)
Stress relief annealed (Franklin model)
Recrystallization annealed (Franklin model)
Partial recrystallization annealed (Franklin model)
Zirlo (Franklin model)
ESCORE model (Rashid model)

This method builds the irradiation growth strain tensor based on the geometry options set by the user. Since irradiation growth should occur in the axial direction only while being volume conserving, it is necessary to specify a strain increment for the other two directions. This strain increment is given by (2)

Example Input Syntax


[./irradiation_growth_eigenstrain]
  type = ZryIrradiationGrowthEigenstrain
  block = 1
  fast_neutron_fluence = fast_neutron_fluence
  zircaloy_material_model_type = stress_relief_annealed
  eigenstrain_name = irr_eigenstrain
[../]
(test/tests/tensor_mechanics/zry_mechanics/zry_irradiation_growth/test_irradiation_growth_rev2.i)

The eigenstrain_name parameter value must also be set for the strain calculator, and an example parameter setting is shown below:


[./strain]
  type = ComputeAxisymmetricRZFiniteStrain
  block = 1
  eigenstrain_names = irr_eigenstrain
[../]
(test/tests/tensor_mechanics/zry_mechanics/zry_irradiation_growth/test_irradiation_growth_rev2.i)

Input Parameters

  • fast_neutron_fluenceThe fast neutron fluence

    C++ Type:std::vector

    Description:The fast neutron fluence

  • eigenstrain_nameMaterial property name for the eigenstrain tensor computed by this model. IMPORTANT: The name of this property must also be provided to the strain calculator.

    C++ Type:std::string

    Description:Material property name for the eigenstrain tensor computed by this model. IMPORTANT: The name of this property must also be provided to the strain calculator.

Required Parameters

  • zircaloy_material_model_typestress_relief_annealedType of irradiation growth volumetric swelling formulation. Choices are: stress_relief_annealed recrystallization_annealed partial_recrystallization_annealed zirlo ESCORE_IrradiationGrowthZr4

    Default:stress_relief_annealed

    C++ Type:MooseEnum

    Description:Type of irradiation growth volumetric swelling formulation. Choices are: stress_relief_annealed recrystallization_annealed partial_recrystallization_annealed zirlo ESCORE_IrradiationGrowthZr4

  • 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.

  • initial_fast_fluence0The initial fast neutron fluence

    Default:0

    C++ Type:double

    Description:The initial fast neutron fluence

  • base_nameOptional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases

    C++ Type:std::string

    Description:Optional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases

  • model_rthetaFalseSet to true for a plane strain model

    Default:False

    C++ Type:bool

    Description:Set to true for a plane strain model

  • 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

  • growth_direction1The direction you want to apply irradiation growth

    Default:1

    C++ Type:unsigned int

    Description:The direction you want to apply irradiation growth

  • 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. DG Franklin. Zircaloy-4 cladding deformation during power reactor irradiation. In Zirconium in the Nuclear Industry. ASTM International, 1982.[BibTeX]
  2. Y Rashid, R Dunham, and R Montgomery. Fuel Analysis and Licensing Code: FALCON MOD01. Technical Report, Electric Power Research Institute, December 2004.[BibTeX]