MetallicFuelCoolantWastage

Compute wastage thickness on the cladding-coolant interface.

Description of Original Correlations

Although the majority of the cladding wastage in sodium-cooled fast reactors (SFRs) occurs at the cladding-fuel interface, known as fuel-cladding chemical interaction (FCCI), which is covered by MetallicFuelWastage, a minor but non-negligible part of the cladding wastage forms on the outer surface of the cladding due to chemical interaction between cladding and primary sodium coolant. The cladding-coolant interaction on HT9 and SS316 cladding is described by this MetallicFuelCoolantWastage Materials module. The model was developed at Argonne National Laboratory as part of the Integral Fast Reactor (IFR) program (Briggs et al., 1995)) As an alternative option, a sodium corrosion thinning correlation developed based on experiments performed by Hanford Engineering Development Laboratory (HEDL) is also given (Anantatmula and Brehm, 1985; Nam et al., 1998).

In the original cladding-coolant interaction wastage correlation, the thickness of the wastage ( in ) is dependent on temperature ( in ) and irradiation time ( in ). The wastage thickness correlation is comprised of three terms: where the first term is the damage depth due to mass loss (corrosion thinning), the second term is ferrite layer growth, and the third term is intergranular penetration. Here, is the damage factor in , is the mass loss rate in , is the ferrite layer growth coefficient in , is intergranular penetration depth in .In the current model, is only applicable to austenitic cladding (i.e., SS316).

Here, is an activation energy controlled process, where and are respectively pre-exponential coefficient () and activation energy (), is gas constant ().

is a linear function of temperature with coefficients and and is only meaningful when the value is non-negative,

has similar form as the , with and as coefficients,

Hence, the correlation can be written as follows,

Direct Differentiation of the Correlation

In order to use the original correlation in BISON time iteration, the expression must be written in a differential form. A direct approach is by adopting the time derivative of the correlation,

This approach may lead to large discrepancies from the original correlation prediction when the time step is large and/or the temperature changes significantly during operation.

Effective Time Approach

To make the differential form of the correlation more consistent with the original correlation, an effective time is calculated every time step as if the temperature has been maintained at the current value since the beginning. where is the wastage thickness of the last time step, is the current temperature, can be calculated by solving the quadratic equation. Therefore, the wastage thickness of this time step is,

This approach ensures numerical consistency with the original form of the wastage thickness correlation.

Example Input Syntax

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [wastage]
    type = MetallicFuelCoolantWastage<<<{"description": "Compute wastage thickness on the cladding-coolant interface.", "href": "MetallicFuelCoolantWastage.html"}>>>
    clad_material<<<{"description": "Fuel cladding material choice. Options: HT9 SS316 HT9_HEDL custom"}>>> = HT9
    use_effective_method<<<{"description": "Flag to use effective method that gives better consistency with original correlation"}>>> = true
    temperature<<<{"description": "The coupled temperature (K)"}>>> = temp
    scale_factor<<<{"description": "Scale factor used for fine tuning"}>>> = 1.0
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
[]
(test/tests/metallic_fuel_coolant_wastage/cc_wastage_ht9_eff.i)

Input Parameters

  • temperatureThe coupled temperature (K)

    C++ Type:std::vector<VariableName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The coupled temperature (K)

Required Parameters

  • A0.55771Pre-preexponential coefficient (kg m^-2 s^-1)

    Default:0.55771

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Pre-preexponential coefficient (kg m^-2 s^-1)

  • B13.026e-11Ferrite layer growth term factor 1 (m/s^0.5/K)

    Default:3.026e-11

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Ferrite layer growth term factor 1 (m/s^0.5/K)

  • B22.5529e-08Ferrite layer growth term factor 2 (m/s^0.5)

    Default:2.5529e-08

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Ferrite layer growth term factor 2 (m/s^0.5)

  • F0.00021Damage factor (m kg^-1 m^2)

    Default:0.00021

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Damage factor (m kg^-1 m^2)

  • QR18117.8Activtion energy divided by gas constant (K)

    Default:18117.8

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activtion energy divided by gas constant (K)

  • X13.916e-11Intergranular penetrtion term factor 1 (m/s^0.5/K)

    Default:3.916e-11

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Intergranular penetrtion term factor 1 (m/s^0.5/K)

  • X23.419e-08Intergranular penetrtion term factor 2 (m/s^0.5)

    Default:3.419e-08

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Intergranular penetrtion term factor 2 (m/s^0.5)

  • blockThe list of blocks (ids or names) that this object will be applied

    C++ Type:std::vector<SubdomainName>

    Controllable:No

    Description:The list of blocks (ids or names) that this object will be applied

  • boundaryThe list of boundaries (ids or names) from the mesh where this object applies

    C++ Type:std::vector<BoundaryName>

    Controllable:No

    Description:The list of boundaries (ids or names) from the mesh where this object applies

  • clad_materialHT9Fuel cladding material choice. Options: HT9 SS316 HT9_HEDL custom

    Default:HT9

    C++ Type:MooseEnum

    Options:HT9, SS316, HT9_HEDL, custom

    Controllable:No

    Description:Fuel cladding material choice. Options: HT9 SS316 HT9_HEDL custom

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

    Default:True

    C++ Type:bool

    Controllable:No

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

  • 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 computeQpProperties() 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

    Options:NONE, ELEMENT, SUBDOMAIN

    Controllable:No

    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 computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

  • declare_suffixAn optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

  • scale_factor1Scale factor used for fine tuning

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor used for fine tuning

  • use_avg_tempTrueFlag to use avarage temperature of current and last timestepsfor better accuracy

    Default:True

    C++ Type:bool

    Controllable:No

    Description:Flag to use avarage temperature of current and last timestepsfor better accuracy

  • use_effective_methodFalseFlag to use effective method that gives better consistency with original correlation

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Flag to use effective method that gives better consistency with original correlation

Optional Parameters

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

    C++ Type:std::vector<std::string>

    Controllable:No

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

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Controllable:Yes

    Description:Set the enabled status of the MooseObject.

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

    Default:True

    C++ Type:bool

    Controllable:No

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

  • seed0The seed for the master random number generator

    Default:0

    C++ Type:unsigned int

    Controllable:No

    Description:The seed for the master random number generator

  • 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

    Controllable:No

    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.

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<std::string>

    Controllable:No

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

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

    Default:none

    C++ Type:std::vector<OutputName>

    Controllable:No

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

Outputs Parameters

  • prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

  • use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

    Default:False

    C++ Type:bool

    Controllable:No

    Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

References

  1. RP Anantatmula and WF Brehm. Sodium compatibility of ht-9 and fe-9cr-1mo steels. Technical Report, Hanford Engineering Development Lab., 1985.[BibTeX]
  2. L. L. Briggs, L. K. Chang, and D. J. Hill. Safety Analysis and Technical Basis for Establishing an Interim Burnup Limit for Mark-V and Mark-VA Fueled Subassemblies in EBR-II. Technical Report ANL-NSE-1, Argonne National Laboratory, 1995.[BibTeX]
  3. Cheol Nam, Woan Hwang, and Dong-Seong Sohn. Statistical failure analysis of metallic u-10zr/ht9 fast reactor fuel pin by considering the weibull distribution and cumulative damage fraction. Annals of Nuclear Energy, 25(17):1441–1453, 1998.[BibTeX]