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)

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


block = 1
fast_neutron_fluence = fast_neutron_fluence
zircaloy_material_model_type = stress_relief_annealed
eigenstrain_name = irr_eigenstrain
[../]

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
[../]

## 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

• 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

## 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]