MetallicFuelWastage

Computes wastage thickness on the fuel-cladding interface.

Description

MetallicFuelWastage finds the thickness of the wastage formed on the interface of metallic fuel and steel cladding in fast neutron reactors using a series of methods. Some of these correlations, including flux_ht9, flux_ss316, flux_d9, burnup_ht9_legacy, and flux_burnup_ht9, are used by LIFE-METAL. Others have been optimized to provided improved accuracy based on recently developed EBR-II operating conditions in the FIPD database. The correlations were developed and implemented by Argonne National Laboratory.

The wastage thickness, (), includes three components, fission product wastage (), carbon depletion zone (), and cladding thinning that is not included in the former two components (). Three modes of correlations have been implemented into BISON: the flux dependent mode, the burnup dependent mode, and the flux-burnup dependent mode.

In the flux dependent mode, the growth rate of wastage () is calculated as where is the neutron fast flux (), is the gas constant (), is temperature (), is time (), and (), (), and () are empirical parameters. The flux dependent mode can be used for three cladding materials (HT9, stainless steel 316, and D9, which correspond to methods flux_ht9, flux_ss316, and flux_d9).

In the burnup dependent mode, the growth rate of wastage is calculated as where is the atom percent burnup, is the gas constant (), is temperature (), is time (), and () and () are empirical parameters. The burnup dependent mode is only available for HT9 at this moment. Users may use either method burnup_ht9_legacy for LIFE-METAL adopted parameters or method burnup_ht9_opt for recent Dakota-optimized parameters based on improved EBR-II operating parameters available in FIPD. Originally, the time-varying sodium flow flux information from FIPD was used with BISON's CoolantChannelAction to set up the temperature boundary condition. In that case, Qb_ht9_opt and k0_ht9_opt are 1.7959350973E5 and 3.3786089381E1 , respectively. This temperature boundary condition approach was later replaced by directly adopting the time-varying cladding outer surface temperature available in FIPD by FIPDAxialProfileFunction as the Dirichlet boundary condition. The default parameters are based on the latter temperature boundary condition approach.

Derived from the burnup dependent mode, a new correlation (burnup_ht9_gap) that takes gap closure into account was developed, which has the following expression, where and are the atomic percent burnup and time corresponding to the initiation of FCCI formation. has two components: (1) gap closure burnup, , which is tracked by the local penetration depth of the contact module. A minimum gap width (min_gap) is used as the criterion for gap closure. (2) rind consumption burnup delay, , which can be defined by an empirical function rind_burnup_delay. The model will be further optimized once the advanced swelling and contact models are developed.

In the flux-burnup dependent mode, the growth rate of wastage is calculated as This is a combination of the previous two models. Currently this mode is only available for HT9 cladding through the method flux_burnup_ht9.

The parameters were fitted by using PIE data from EBR-II.

The flux mode was calibrated using PIE data from X419, X420, X421, X425, X429, X441, and X447. The burnup mode and flux-burnup mode were calibrated using the PIE data from X447.

Example Input Syntax

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [wastage]
    type = MetallicFuelWastage<<<{"description": "Computes wastage thickness on the fuel-cladding interface.", "href": "MetallicFuelWastage.html"}>>>
    method<<<{"description": "Wastage thickness method choice. Options: flux_ht9 flux_ss316 flux_d9 burnup_ht9_legacy burnup_ht9_opt burnup_ht9_gap flux_burnup_ht9"}>>> = flux_ht9
    temperature<<<{"description": "The coupled temperature (K)"}>>> = temp
    burnup<<<{"description": "Burnup material name"}>>> = 0.0
    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/MetallicFuelWastage/flux_ht9.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

  • D_0bf_ht90.000930106Coefficient term of the diffusion expression in burnup and flux related expression (m^2/s)

    Default:0.000930106

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of the diffusion expression in burnup and flux related expression (m^2/s)

  • D_d9_07.88529Coefficient term of the diffusion expression (m^2/s)

    Default:7.88529

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of the diffusion expression (m^2/s)

  • D_d9_io6.39784e-38Coefficient term of radiation enhanced diffusion expression (m^4)

    Default:6.39784e-38

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of radiation enhanced diffusion expression (m^4)

  • D_ht9_00.000112186Coefficient term of the diffusion expression (m^2/s)

    Default:0.000112186

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of the diffusion expression (m^2/s)

  • D_ht9_io1.79211e-39Coefficient term of radiation enhanced diffusion expression (m^4)

    Default:1.79211e-39

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of radiation enhanced diffusion expression (m^4)

  • D_iobf_ht97.86737e-43Coefficient term of radiation enhanced diffusion expression in burnup and flux related expression (m^4)

    Default:7.86737e-43

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of radiation enhanced diffusion expression in burnup and flux related expression (m^4)

  • D_ss316_02.41935Coefficient term of the diffusion expression (m^2/s)

    Default:2.41935

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of the diffusion expression (m^2/s)

  • D_ss316_io1.9534e-38Coefficient term of radiation enhanced diffusion expression (m^4)

    Default:1.9534e-38

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Coefficient term of radiation enhanced diffusion expression (m^4)

  • Q_d9266102Activation energy term of the diffusion expression in J/mol

    Default:266102

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term of the diffusion expression in J/mol

  • Q_ht9201782Activation energy term of the diffusion expression in J/mol

    Default:201782

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term of the diffusion expression in J/mol

  • Q_ss316266102Activation energy term of the diffusion expression in J/mol

    Default:266102

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term of the diffusion expression in J/mol

  • Qb_ht9_gap195398Activation energy term used in burnup related expression (J/mol)

    Default:195398

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term used in burnup related expression (J/mol)

  • Qb_ht9_lg179592Activation energy term used in burnup related expression (J/mol)

    Default:179592

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term used in burnup related expression (J/mol)

  • Qb_ht9_opt179069Activation energy term used in burnup related expression (J/mol)

    Default:179069

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term used in burnup related expression (J/mol)

  • Qbf_ht9252253Activation energy term used in burnup and flux related expression (J/mol)

    Default:252253

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Activation energy term used in burnup and flux related expression (J/mol)

  • R8.31446Gas constant in J/mol/K

    Default:8.31446

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Gas constant in J/mol/K

  • 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

  • burnupburnupBurnup material name

    Default:burnup

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:Burnup material name

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

  • fast_neutron_fluxfast_neutron_fluxThe fast neutron flux (neutrons/m^2-sec)

    Default:fast_neutron_flux

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:The fast neutron flux (neutrons/m^2-sec)

  • k0_ht9_gap298.835Preexponential coefficient used in burnup related expression (m/s^0.5)

    Default:298.835

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Preexponential coefficient used in burnup related expression (m/s^0.5)

  • k0_ht9_lg32.973Preexponential coefficient used in burnup related expression (m/s^0.5)

    Default:32.973

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Preexponential coefficient used in burnup related expression (m/s^0.5)

  • k0_ht9_opt39.1273Preexponential coefficient used in burnup related expression (m/s^0.5)

    Default:39.1273

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Preexponential coefficient used in burnup related expression (m/s^0.5)

  • methodburnup_ht9_legacyWastage thickness method choice. Options: flux_ht9 flux_ss316 flux_d9 burnup_ht9_legacy burnup_ht9_opt burnup_ht9_gap flux_burnup_ht9

    Default:burnup_ht9_legacy

    C++ Type:MooseEnum

    Options:flux_ht9, flux_ss316, flux_d9, burnup_ht9_legacy, burnup_ht9_opt, burnup_ht9_gap, flux_burnup_ht9

    Controllable:No

    Description:Wastage thickness method choice. Options: flux_ht9 flux_ss316 flux_d9 burnup_ht9_legacy burnup_ht9_opt burnup_ht9_gap flux_burnup_ht9

  • min_gap1e-05Minimum fuel-cladding gap to initiate FCCI

    Default:1e-05

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Minimum fuel-cladding gap to initiate FCCI

  • penetration_depth0negative value of gap width (m).

    Default:0

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:negative value of gap width (m).

  • rind_burnup_delay0delayed burnup due to rind after gap closure.

    Default:0

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:delayed burnup due to rind after gap closure.

  • 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_initial_timeFalseFlag to use an initial time corresponding to initial burnup that is determined by gap width

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Flag to use an initial time corresponding to initial burnup that is determined by gap width

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