HT9FailureClad

Failure model for HT-9 cladding. Contains multiple models for steady state (burnup calculations) and transient operations.

Description

HT9FailureClad is the model for HT9 cladding failure during both steady state and transient irradiation conditions. Here the steady state and transient are the terms used by reactor physicists: steady state is mainly referred to those irradiation conditions during normal reactor operation, while transient is mainly referred to those irradiation conditions during accident scenarios. Unlike the name implies, the steady state conditions may also involve some time-dependent irradiation conditions.

This cladding failure model contains two types of cladding failure models: One is based on the cumulative damage fraction (CDF) and the other is based on the constrained cavity growth (CCG) with Diffusion and Creep with Sliding (D&CS) mechanism. The CDF model is mainly used for steady state (normal reactor operation) simulations while the CCG model is used for power transients (accidents scenarios). Some CDF models may also be used for transient scenarios.

Both types of the creep-fracture failure models are available in BISON similar to the approach used in Karahan and Buongiorno (2010).

For the CDF-style approach, a variety of options are provided based on different time-to-rupture correlations. First, a pair of long and short time-scale CDF options, which are reported in Karahan and Buongiorno (2010), are available as the cdf_long and cdf_short option. Steady-state simulations are usually handled by the Cumulative Damage Fraction (CDF) model using the long time-scale Dorn parameter correlation (Nam et al., 1998). Transient scenarios are handled with the short time-scale Dorn parameter correlation (Nam et al., 1998). A cdf_both option is also available to adopt the longer predicted time-to-rupture of the two correlations. Second, the time-to-rupture correlation set reported in the Metallic Fuel Handbook (MFH)(Hofman et al., 2019) is provided as the cdf_mfh option. Last but not least, the transient specific time-to-rupture correlation in the Westinghouse Hanford Corporate report is provided as the cdf_whc option (Huang, 1994).

The advanced Constrained Cavity Growth (CCG) with Diffusion and Creep with Sliding (D&CS) model (Tvergaard, 1985) may also be used for short time-scale transients if all parameters are known.

Cumulative Damage Fraction

The CDF model compares the time-to-rupture value with experimentally obtained results as a function of stress and absolute temperature. When this value (CDF) equals one, the material has failed.

Thus, the key for CDF prediction is an appropriate time-to-rupture correlation, which is usually a function of temperature and hoop stress. In literature, there exist a series of time-to-rupture correlations for HT9. Each of these correlations has its own features and applicable scenarios, which will be introduced in detail as follows.

DiMelfi/Nam Long/Short Time-Scale Correlation Set

The time of rupture function in hours is found in where is 154 kcal/mole for long time-scale data and 70.17 kcal/mole for short time-scale data, is the Boltzmann constant, and is the absolute temperature in K. is the Dorn parameter coefficient ranging from 3.915 10 at 650 C to 1 at 600 C for long time-scale data and 2.778 10 for all short time-scale data. The Dorn parameter comes from the curve fits for long and short time-scale correlations, respectively, where is the hoop stress in Pa.

For long time-scale correlations with temperatures outside of the range of the Dorn parameter coefficient , the value is set to the appropriate end value to avoid extrapolation. Interpolation of the value is conducted within the temperature range. If both correlations are chosen, then the longer rupture time result is taken as the actual time of rupture (DiMelfi et al., 1993). When choosing both correlations, erroneous values are possible. Both correlations are recommended when simulation values fall within the range where the long and short time-scale correlations overlap.

MFH Correlation Set

The Metallic Fuel Handbook (MFH) provides another set of correlations for the creep rupture time for HT9. Similar to the DiMelfi's approach, the MFH creep rupture correlation set also adopts a two-stage approach to account for high-stress and low-stress scenarios. For both stages, the creep rupture time has the same form as follows, where is time-to-rupture in hour, is temperature in Kelvin (K), and is hoop stress in MPa. On the other hand, , , and are the three fitting coefficients that are different for high-stress and low-stress scenarios. The values of these coefficients can be found in the following table.

ABC
High Stress-32.49057781-11800
Low Stress-35.17345858-5563.1

The choice of which set of coefficients to use should be made to ensure that the correlation predicts the shortest time-to-rupture.

WHC Correlation Set

In 1994, Westinghouse Hanford Corporate (WHC) reported an updated creep rupture time correlation for HT9 under transient conditions. In the WHC correlation, the time-to-rupture () has the following form, where is the Dorn parameter, is the activation energy fitted to be 70.17 kcal/mole, is the gas constant, and is temperature in Kelvin (K).

The Dorn parameter is a function of applied hoop stress (), where = 730 MPa is the nominal hardness of the material. Fitting coefficient = 12.47 is assumed to be independent from temperature, while is a function of temperature as follows,

Summary of Time-to-Rupture Correlation Selection

All the four time-to-rupture correlations are compared in Figure 1, showing that use of different CDF options would lead to very dissimilar CDF prediction. Generally speaking, cdf_long and cdf_mfh should be used in those steady state simulation, while cdf_short and cdf_whc should be used in those short-term transient scenarios.

Figure 1: Comparison between different time-to-rupture correlations available to select for BISON HT9 CDF calculation.

It is worth mentioning that Briggs et al. (1995) provides a cladding failure probability assessment based on CDF values, which is based on the cdf_mfh correlation. Therefore, cdf_mfh is the ideal selection if the cladding failure probability needs to be evaluated accordingly. Briggs et al. (1995) also suggests use of a time-to-rupture correlation for transient/accident scenarios reported back in 1985 by WHC, which was later updated to be the cdf_whc correlation in 1994 (Huang, 1994).

Constrained Cavity Growth

The CCG with D&CS model calculates the crack radius of periodic cavities along grain boundaries. The cavity centers are spaced equally at a distance of . Failure occurs when . The correlation is very sensitive to the user supplied value of . For most situations, the value of is either known from experiment data or fitted for a known failure of an experiment which is being simulated. The short time-scale CDF model may also be used for short time-scale failure simulations where the value of is unknown or in question.

The crack radius growth rate is related to the cavity volume growth rate by with being defined as where being the ratio of the grain boundary free energy to twice the grain surface free energy.

The volume growth rate is the sum of the rigid grain growth rate and power law creeping material growth rate for .

The average normal stress is . The sintering stress can be calculated as . The grain boundary diffusion parameter is where is the boundary diffusivity, is the atomic volume, is the Boltzmann constant, and is the absolute temperature in K. The area fraction of the grain boundary is determined with where . The Von Mises stress is . The hydrostatic (mean) stress is . The effective creep strain- rate is . With the assumption that the material follows a power law creep, is the value of the power, , and .

The crack length is assumed to begin at a minimum value. The crack length is never allowed to fall below this value. The crack may shrink after it has grown. However, after failure has occurred, the crack is assumed to be permanent and can no longer shorten. The crack length is found by taking the calculated growth rate and multiplying by the current time increment.

Example Input Syntax

HT9FailureClad has tests for both CDF and CCG models. The tests also verify errors are produced when required input is missing or undefined. As both models depend on other models, investigating the entire test input file is recommended.

Cumulative Damage Fraction

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [FailureClad]
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 1
    type = HT9FailureClad<<<{"description": "Failure model for HT-9 cladding. Contains multiple models for steady state (burnup calculations) and transient operations.", "href": "HT9FailureClad.html"}>>>
    hoop_stress<<<{"description": "Hoop stress in cladding"}>>> = hoop_stress
    temperature<<<{"description": "Temperature of cladding"}>>> = temp
    method<<<{"description": "Failure method choice. Options: ccg_dcs cdf_long cdf_short cdf_both cdf_mfh cdf_whc"}>>> = cdf_long
  []
[]
(test/tests/solid_mechanics/failurecladHT9/cdf_based/axisymm_thin_cyl_q4_cdf.i)

Constrained Cavity Growth

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [FailureClad]
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 0
    type = HT9FailureClad<<<{"description": "Failure model for HT-9 cladding. Contains multiple models for steady state (burnup calculations) and transient operations.", "href": "HT9FailureClad.html"}>>>
    hoop_stress<<<{"description": "Hoop stress in cladding"}>>> = hoop_stress
    temperature<<<{"description": "Temperature of cladding"}>>> = temp
    method<<<{"description": "Failure method choice. Options: ccg_dcs cdf_long cdf_short cdf_both cdf_mfh cdf_whc"}>>> = ccg_dcs
    # For CCG with D&CS
    eff_strain_rate_creep<<<{"description": "Effective creep strain rate in cladding (1/s)"}>>> = creep_rate_aux
    hydrostatic_stress<<<{"description": "Hydrostatic (mean) stress in cladding"}>>> = hydrostatic_stress
    von_mises_stress<<<{"description": "Von Mises stress in the cladding"}>>> = vonmises
  []
[]
(test/tests/solid_mechanics/failurecladHT9/ccg_short/ccgtest.i)

Input Parameters

  • hoop_stressHoop stress in cladding

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Hoop stress in cladding

  • temperatureTemperature of cladding

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Temperature of cladding

Required Parameters

  • a_initial7e-08Initial cavity radius (m)

    Default:7e-08

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Initial cavity radius (m)

  • atomic_volume1.18e-29Atomic Volume (m^3)

    Default:1.18e-29

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Atomic Volume (m^3)

  • avogadro_num6.02214e+23Avogadro's Number (1/mol)

    Default:6.02214e+23

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Avogadro's Number (1/mol)

  • b3.5e-06Distance from cavity center to midplane between cavities (m)

    Default:3.5e-06

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Distance from cavity center to midplane between cavities (m)

  • 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

  • boltzmann1.38065e-23Boltzmann constant (J/K)

    Default:1.38065e-23

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Boltzmann constant (J/K)

  • 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

  • boundary_diffusivity1.1e-12Grain boundary diffusivity (m^3/s)

    Default:1.1e-12

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Grain boundary diffusivity (m^3/s)

  • boundary_free_energy0.85Grain boundary free energy (J/m^2)

    Default:0.85

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Grain boundary free energy (J/m^2)

  • cdf_failure_namecdf_failureName of created CDF failure material

    Default:cdf_failure

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:Name of created CDF failure material

  • comparedless_equalOptions for variable _compared_ to criteria: greater_than greater_equal less_equal less_than

    Default:less_equal

    C++ Type:MooseEnum

    Options:greater_than, greater_equal, less_equal, less_than

    Controllable:No

    Description:Options for variable _compared_ to criteria: greater_than greater_equal less_equal less_than

  • 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_criteria0Numerical value providing criteria value.

    Default:0

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Numerical value providing criteria value.

  • 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

  • creep_n_power5Power law creeping material constant

    Default:5

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Power law creeping material constant

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

  • eff_strain_rate_creepEffective creep strain rate in cladding (1/s)

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Effective creep strain rate in cladding (1/s)

  • function_criteriaFunction name providing criteria value.

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Function name providing criteria value.

  • hydrostatic_stressHydrostatic (mean) stress in cladding

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Hydrostatic (mean) stress in cladding

  • methodcdf_longFailure method choice. Options: ccg_dcs cdf_long cdf_short cdf_both cdf_mfh cdf_whc

    Default:cdf_long

    C++ Type:MooseEnum

    Options:ccg_dcs, cdf_long, cdf_short, cdf_both, cdf_mfh, cdf_whc

    Controllable:No

    Description:Failure method choice. Options: ccg_dcs cdf_long cdf_short cdf_both cdf_mfh cdf_whc

  • surface_free_energy2.1Grain surface free energy (J/m^2)

    Default:2.1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Grain surface free energy (J/m^2)

  • variable_checkVariable name which is compared to criteria. Example: Var < 0, true=failed

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Variable name which is compared to criteria. Example: Var < 0, true=failed

  • von_mises_stressVon Mises stress in the cladding

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Von Mises stress in the cladding

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. 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]
  2. R. J. DiMelfi, E. E. Gruber, J. M. Kramer, and T. H. Hughes. Strength and rupture-life transitions caused by secondary carbide precipitation in HT-9 during high-temperature low-rate mechanical testing. Technical Report ANL/RE/CP–75828, Argonne National Laboratory, December 1993.[BibTeX]
  3. G. L. Hofman, M. C. Billone, J. F. Koenig, J. M. Kramer, J. D. B. Lambert, L. Leibowitz, Y. Orechwa, D. R. Pedersen, D. L. Porter, H. Tsai, and A. E. Wright. Metallic fuels handbook. Technical Report ANL-NSE-3, Argonne National Laboratory, 2019.[BibTeX]
  4. FH Huang. Transient failure behavior of HT9. Technical Report WHC-SA-2513, Westinghouse Hanford Co., July 1994. URL: https://www.osti.gov/biblio/10173288, doi:10.2172/10173288.[BibTeX]
  5. A. Karahan and J. Buongiorno. A new code for predicting the thermo-mechanical and irradiation behavior of metallic fuels in sodium fast reactors. Journal of Nuclear Materials, 396(2):283–293, 2010. doi:10.1016/j.jnucmat.2009.11.022.[BibTeX]
  6. 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. doi:10.1016/S0306-4549(98)00043-7.[BibTeX]
  7. V. Tvergaard. Effect of grain boundary sliding on creep constrained diffusive cavitation. Journal of the Mechanics and Physics of Solids, 33(5):447–469, 1985. doi:10.1016/0022-5096(85)90009-2.[BibTeX]