MonolithicSiCThermal

Computes thermal conductivity and specific heat of monolithic silicon carbide.

Description

The MonolithicSiCThermal model is appropriate for silicon carbide produced through chemical vapor deposition.

Thermal Conductivity

Three models exist for thermal conductivity of monolithic SiC: Snead, Stone, and Miller.

Snead Model

Snead et al. (2007) developed a model that applies to both irradiated and unirradiated monolithic SiC. The unirradiated thermal conductivity (W/m-K) is given by Koyanagi et al. (2021): where is the thermal conductivity (W/m-K) and is the temperature (K). The increase in thermal resistance with irradiation is considered by the following equation:

where is the volumetric swelling strain (unitless). The thermal conductivity of irradiated material is expressed as:

Stone Model

Stone et al. (2015) used a model that also applies to both irradiated and unirradiated monolithic SiC. The unirradiated thermal conductivity (W/m-K) is given by: where is the temperature in K. Note that the coefficient is different than the value reported in Stone et al. (2015) as it appears to be a typo as per Freeman et al. (2018). Irradiation effects are taken into account by assuming a resistance network where the total thermal conductivity is given by: where is the resistance due to irradiation effects given by: where is the total volumetric swelling strain calculated by CompositeSiCVolumetricSwellingEigenstrain.

Miller Model

The thermal conductivity used in the PARFUME code is given by Miller et al. (2018).

where (W/m-K) is the thermal conductivity and is in Kelvin.

Specific Heat Capacity

The correlation for specific heat (J/kg-K) is given by Snead et al. (2007): where is the temperature in Kelvin.

Range of Applicability and Uncertainty

Thermal Conductivity

Snead model is applicable in the point defect regime for which the temperature range is 423 - 1273 K (Koyanagi et al., 2021). The uncertainty associated with this model is not available.

The Stone model was derived from a 4th order polymonial fit from temperatures 273 - 1573 K. The uncertainty in the calculation was not published.

The Miller model is valid for all operating temperatures. The uncertainty in the model has not been published.

Specific Heat

The model used for the specific heat applies to temperatures 200 - 2400 K with uncertainty of 7% for 200 T 1000 K, and 4% for 1000 T 2400 K.

Example Input Syntax

An example of using the Snead model is given by:

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [MonolithicSiCThermal]
    type = MonolithicSiCThermal<<<{"description": "Computes thermal conductivity and specific heat of monolithic silicon carbide.", "href": "MonolithicSiCThermal.html"}>>>
    temperature<<<{"description": "Coupled Temperature"}>>> = temperature
    thermal_conductivity_model<<<{"description": "Options for the correlation used to calculate thermal conductivity"}>>> = SNEAD
  []
[]
(test/tests/monolithicSiCThermal/thermal_snead.i)

An example of using the Stone model is given by:

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [MonolithicSiCThermal]
    type = MonolithicSiCThermal<<<{"description": "Computes thermal conductivity and specific heat of monolithic silicon carbide.", "href": "MonolithicSiCThermal.html"}>>>
    temperature<<<{"description": "Coupled Temperature"}>>> = temperature
    thermal_conductivity_model<<<{"description": "Options for the correlation used to calculate thermal conductivity"}>>> = STONE # This is the default
  []
[]
(test/tests/monolithicSiCThermal/thermal_stone_irradiated.i)

An example of using the Miller model is given by:

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  [thermal]
    type = MonolithicSiCThermal<<<{"description": "Computes thermal conductivity and specific heat of monolithic silicon carbide.", "href": "MonolithicSiCThermal.html"}>>>
    temperature<<<{"description": "Coupled Temperature"}>>> = temperature
    thermal_conductivity_model<<<{"description": "Options for the correlation used to calculate thermal conductivity"}>>> = MILLER
  []
[]
(test/tests/monolithicSiCThermal/thermal_miller.i)

Input Parameters

  • temperatureCoupled Temperature

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

    Unit:(no unit assumed)

    Controllable:No

    Description:Coupled Temperature

Required Parameters

  • 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

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

  • thermal_conductivity_modelSTONEOptions for the correlation used to calculate thermal conductivity

    Default:STONE

    C++ Type:MooseEnum

    Options:SNEAD, STONE, MILLER

    Controllable:No

    Description:Options for the correlation used to calculate thermal conductivity

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

  • specific_heat_scale_factor1Scale factor to be applied to the specific heat. Used for calibration and sensitivity studies

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor to be applied to the specific heat. Used for calibration and sensitivity studies

  • thermal_conductivity_scale_factor1Scale factor to be applied to the thermal conducitivity. Used for calibration and sensitivity studies

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor to be applied to the thermal conducitivity. Used for calibration and sensitivity studies

Advanced: Scaling Factors Parameters

Input Files

References

  1. R. A. Freeman, T. Martin, E. Roberts, and T. W. Knight. Analysis of thermal creep for uranium silicide fuel using Bison. In Proceedings of the 2018 International Congress on Advances in Nuclear Power Plants (ICAPP 18). Charlotte, NC, 2018.[BibTeX]
  2. Takaaki Koyanagi, Hsin Wang, Christian M Petrie, Christian P Deck, Weon-Ju Kim, Daejong Kim, Cédric Sauder, James Braun, Yutai Katoh, and others. Thermal diffusivity and thermal conductivity of sic composite tubes: the effects of microstructure and irradiation. Journal of Nuclear Materials, 557:153217, 2021. doi:10.1016/j.jnucmat.2021.153217.[BibTeX]
  3. G.K. Miller, D.A. Petti, J.T. Maki, D.L. Knudson, and W.F. Skerjanc. PARFUME Theory and Model Basis Report. Report INL/EXT-08-14497 (Rev.1), Idaho National Laboratory, September 2018.[BibTeX]
  4. L. L. Snead, T. Nozawa, Y. Katoh, T.-S. Byun, S. Kondo, and D. A. Petti. Handbook of sic properties for fuel performance modeling. Journal of Nuclear Materials, 371:329–377, 2007.[BibTeX]
  5. J. G. Stone, R. Schleicher, C. P. Deck, G. M. Jacobsen, H. E. Khalifa, and C. A. Back. Stress analysis and probabilistic assessment of multi-layer sic-based accident tolerant nuclear fuel cladding. Journal of Nuclear Materials, 466:682–697, 2015.[BibTeX]