MonolithicSiCVolumetricSwellingEigenstrain

Computes an eigenstrain due to volumetric swelling for Monolithic SiC based on temperature and fluence.

Description

The MonolithicSiCVolumetricSwellingEigenstrain material model computes an eigenstrain tensor that accounts for volumetric swelling in monolithic silicon carbide (SiC). The details about this model can be found in Katoh et al. (2018). The model computes the volumetric swelling strain (unitless) $var element = document.getElementById("moose-equation-df87875f-38cc-4776-a79d-275080b52480");katex.render("\\Delta V/V", element, {displayMode:false,throwOnError:false});$ : $var element = document.getElementById("moose-equation-8b811a76-8169-4628-86af-a6274a9c5e03");katex.render(" \\frac{\\Delta V}{V} = S_s \\left[ 1 - \\exp\\left( -\\frac{\\gamma}{\\gamma_c} \\right) \\right]^{2/3}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-99ea1198-fcde-402d-a7df-7c4c14990f38");katex.render("S_s", element, {displayMode:false,throwOnError:false});$ is the saturated swelling strain (volumetric), $var element = document.getElementById("moose-equation-60a1fcdc-26bc-46e1-8d7b-022888bfb92c");katex.render("\\gamma", element, {displayMode:false,throwOnError:false});$ is the dose in dpa, $var element = document.getElementById("moose-equation-2acfdf70-44df-47e4-a20f-c062f311a665");katex.render("\\gamma_c", element, {displayMode:false,throwOnError:false});$ is the characteristic dose in dpa for swelling saturation. The saturated swelling and the characteristic dose are polynomial functions of temperature $var element = document.getElementById("moose-equation-c7cae89b-f807-48a8-a97d-22df39091927");katex.render("T", element, {displayMode:false,throwOnError:false});$ as: $var element = document.getElementById("moose-equation-8a761ecf-6f43-486a-ba12-90074e020e44");katex.render(" S_s = b_0 + b_1 T^1 + b_2 T^2 + b_3 T^3", element, {displayMode:true,throwOnError:false});$

$var element = document.getElementById("moose-equation-4637b43f-89c9-4d05-8925-51804a371aea");katex.render(" \\gamma_c = a_0 + a_1 T^1 + a_2 T^2 + a_3 T^3", element, {displayMode:true,throwOnError:false});$ The model is valid for the temperature range 493K < $var element = document.getElementById("moose-equation-ecb8b27e-d4e4-4eea-a7a6-8592a97b7bab");katex.render("T", element, {displayMode:false,throwOnError:false});$ < 713K.

Table 1: Constant values

Constant	Value
$var element = document.getElementById("moose-equation-65a4e5f0-af95-4748-ab66-174e63a17bb1");katex.render("a_0", element, {displayMode:false,throwOnError:false});$	-0.57533
$var element = document.getElementById("moose-equation-1d4da715-caae-4b70-89bd-5da73adae7d4");katex.render("a_1", element, {displayMode:false,throwOnError:false});$	3.3342e-3
$var element = document.getElementById("moose-equation-6d629c3c-1000-4bdb-b80c-7d28fcc642f9");katex.render("a_2", element, {displayMode:false,throwOnError:false});$	-5.3970e-6
$var element = document.getElementById("moose-equation-aaeb7cd4-eb60-4795-ba92-6127fa27f530");katex.render("a_3", element, {displayMode:false,throwOnError:false});$	2.9754e-9
$var element = document.getElementById("moose-equation-5cb60210-3212-41d7-95c6-99f8d9e5fce3");katex.render("b_0", element, {displayMode:false,throwOnError:false});$	5.8366e-2
$var element = document.getElementById("moose-equation-7af27a09-237c-412e-9c1e-1c867a8e029f");katex.render("b_1", element, {displayMode:false,throwOnError:false});$	-1.0089e-4
$var element = document.getElementById("moose-equation-fb05f211-1a37-4dbd-89a5-8764f3cc8fc1");katex.render("b_2", element, {displayMode:false,throwOnError:false});$	6.9368e-8
$var element = document.getElementById("moose-equation-22c4d7e4-f96a-4ae1-bd50-861136915059");katex.render("b_3", element, {displayMode:false,throwOnError:false});$	-1.8152e-11

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [swelling]
    type = MonolithicSiCVolumetricSwellingEigenstrain<<<{"description": "Computes an eigenstrain due to volumetric swelling for Monolithic SiC based on temperature and fluence.", "href": "MonolithicSiCVolumetricSwellingEigenstrain.html"}>>>
    eigenstrain_name<<<{"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."}>>> = swelling_strain
    temperature<<<{"description": "Coupled temperature (K)"}>>> = temperature
    fast_neutron_fluence<<<{"description": "fast neutron fluence"}>>> = fluence
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
[]

The eigenstrain name must also be passed to the strain calculator as shown below:

[Physics<<<{"href": "../../../syntax/Physics/index.html"}>>>]
  [SolidMechanics<<<{"href": "../../../syntax/Physics/SolidMechanics/index.html"}>>>]
    [QuasiStatic<<<{"href": "../../../syntax/Physics/SolidMechanics/QuasiStatic/index.html"}>>>]
      [all]
        strain<<<{"description": "Strain formulation"}>>> = SMALL
        decomposition_method<<<{"description": "Methods to calculate the finite strain and rotation increments"}>>> = EigenSolution
        add_variables<<<{"description": "Add the displacement variables"}>>> = true
        eigenstrain_names<<<{"description": "List of eigenstrains to be applied in this strain calculation"}>>> = swelling_strain
        use_automatic_differentiation<<<{"description": "Flag to use automatic differentiation (AD) objects when possible"}>>> = false
      []
    []
  []
[]

Input Parameters

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
Controllable:No
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.
fast_neutron_fluencefast neutron fluence
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:fast neutron fluence
temperatureCoupled temperature (K)
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled temperature (K)

Required Parameters

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
Controllable:No
Description:Optional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases
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.

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

(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/exact.i)
(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/ad_exact.i)

References

Y. Katoh, T. Koyanagi, J. L. McDuffee, L.L. Snead, and K. Yueh. Dimensional stability and anisotropy of sic and sic-based composites in transition swelling regime. Journal of Nuclear Materials, 499:471–479, 2018. doi:10.1016/j.jnucmat.2017.12.009.

@article{KATOH2018471,
    author = "Katoh, Y. and Koyanagi, T. and McDuffee, J. L. and Snead, L.L. and Yueh, K.",
    title = "Dimensional stability and anisotropy of SiC and SiC-based composites in transition swelling regime",
    journal = "Journal of Nuclear Materials",
    volume = "499",
    pages = "471--479",
    year = "2018",
    doi = "10.1016/j.jnucmat.2017.12.009",
    publisher = "Elsevier"
}

(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/exact.i)

# This tests MonolithicSiCVolumetricSwellingEigenstrain as a function of temperature and fluence.
# Exact solution is calculated as a function in the input file.

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [gen]
    type = ExamplePatchMeshGenerator
    dim = 3
  []
[]

[AuxVariables]
  [temperature]
  []
  [fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [volumetric_strain]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temperature_ramp]
    type = PiecewiseLinear
    x = '0 1'
    y = '300 1400'
  []
  [fluence_ramp]
    type = ConstantFunction
    value = 2e25
  []
[]

[AuxKernels]
  [temperature_aux]
    type = FunctionAux
    variable = temperature
    function = temperature_ramp
  []
  [fluence_aux]
    type = FunctionAux
    variable = fluence
    function = fluence_ramp
  []
  [volumetric_strain]
    type = RankTwoScalarAux
    rank_two_tensor = total_strain
    variable = volumetric_strain
    scalar_type = VolumetricStrain
  []
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = SMALL
  decomposition_method = EigenSolution
  add_variables = true
  eigenstrain_names = swelling_strain
  use_automatic_differentiation = false
[]

[BCs]
  [left]
    type = DirichletBC
    variable = disp_x
    value = 0
    boundary = left
  []
  [bottom]
    type = DirichletBC
    variable = disp_y
    value = 0
    boundary = bottom
  []
  [back]
    type = DirichletBC
    variable = disp_z
    value = 0
    boundary = back
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1
    poissons_ratio = 0.3
  []
  [swelling]
    type = MonolithicSiCVolumetricSwellingEigenstrain
    eigenstrain_name = swelling_strain
    temperature = temperature
    fast_neutron_fluence = fluence
    outputs = all
  []
  [stress]
    type = ComputeLinearElasticStress
  []
[]

[Executioner]
  type = Transient
  num_steps = 10
  dt = 0.1
[]

[Postprocessors]
  [temperature_pp]
    type = ElementAverageValue
    variable = temperature
  []
  [fluence_pp]
    type = ElementAverageValue
    variable = fluence
  []
  [volumetric_strain]
    type = ElementAverageValue
    variable = volumetric_strain
  []
  [swelling_exact]
    type = ParsedPostprocessor
    pp_names = 'temperature_pp fluence_pp'
    expression = 'characteristic_fluence := -0.57533 + 3.3342e-3*temperature_pp - 5.3970e-6*temperature_pp^2 + 2.9754e-9*temperature_pp^3;
                saturated_swelling := 5.8366e-2 - 1.0089e-4*temperature_pp + 6.9368e-8*temperature_pp^2 - 1.8152e-11*temperature_pp^3;
                saturated_swelling * (1 - exp(-1 * fluence_pp / (1e25 * characteristic_fluence)))^(2.0/3.0)'
  []
  [swelling_diff]
    type = ParsedPostprocessor
    pp_names = 'volumetric_strain swelling_exact'
    expression = '(volumetric_strain - swelling_exact) / swelling_exact'
    outputs = none
  []
  [swelling_max_diff]
    type = TimeExtremeValue
    postprocessor = swelling_diff
    value_type = abs_max
  []
[]

[Outputs]
  csv = true
[]

(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/exact.i)

# This tests MonolithicSiCVolumetricSwellingEigenstrain as a function of temperature and fluence.
# Exact solution is calculated as a function in the input file.

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [gen]
    type = ExamplePatchMeshGenerator
    dim = 3
  []
[]

[AuxVariables]
  [temperature]
  []
  [fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [volumetric_strain]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temperature_ramp]
    type = PiecewiseLinear
    x = '0 1'
    y = '300 1400'
  []
  [fluence_ramp]
    type = ConstantFunction
    value = 2e25
  []
[]

[AuxKernels]
  [temperature_aux]
    type = FunctionAux
    variable = temperature
    function = temperature_ramp
  []
  [fluence_aux]
    type = FunctionAux
    variable = fluence
    function = fluence_ramp
  []
  [volumetric_strain]
    type = RankTwoScalarAux
    rank_two_tensor = total_strain
    variable = volumetric_strain
    scalar_type = VolumetricStrain
  []
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = SMALL
  decomposition_method = EigenSolution
  add_variables = true
  eigenstrain_names = swelling_strain
  use_automatic_differentiation = false
[]

[BCs]
  [left]
    type = DirichletBC
    variable = disp_x
    value = 0
    boundary = left
  []
  [bottom]
    type = DirichletBC
    variable = disp_y
    value = 0
    boundary = bottom
  []
  [back]
    type = DirichletBC
    variable = disp_z
    value = 0
    boundary = back
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1
    poissons_ratio = 0.3
  []
  [swelling]
    type = MonolithicSiCVolumetricSwellingEigenstrain
    eigenstrain_name = swelling_strain
    temperature = temperature
    fast_neutron_fluence = fluence
    outputs = all
  []
  [stress]
    type = ComputeLinearElasticStress
  []
[]

[Executioner]
  type = Transient
  num_steps = 10
  dt = 0.1
[]

[Postprocessors]
  [temperature_pp]
    type = ElementAverageValue
    variable = temperature
  []
  [fluence_pp]
    type = ElementAverageValue
    variable = fluence
  []
  [volumetric_strain]
    type = ElementAverageValue
    variable = volumetric_strain
  []
  [swelling_exact]
    type = ParsedPostprocessor
    pp_names = 'temperature_pp fluence_pp'
    expression = 'characteristic_fluence := -0.57533 + 3.3342e-3*temperature_pp - 5.3970e-6*temperature_pp^2 + 2.9754e-9*temperature_pp^3;
                saturated_swelling := 5.8366e-2 - 1.0089e-4*temperature_pp + 6.9368e-8*temperature_pp^2 - 1.8152e-11*temperature_pp^3;
                saturated_swelling * (1 - exp(-1 * fluence_pp / (1e25 * characteristic_fluence)))^(2.0/3.0)'
  []
  [swelling_diff]
    type = ParsedPostprocessor
    pp_names = 'volumetric_strain swelling_exact'
    expression = '(volumetric_strain - swelling_exact) / swelling_exact'
    outputs = none
  []
  [swelling_max_diff]
    type = TimeExtremeValue
    postprocessor = swelling_diff
    value_type = abs_max
  []
[]

[Outputs]
  csv = true
[]

(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/exact.i)

# This tests MonolithicSiCVolumetricSwellingEigenstrain as a function of temperature and fluence.
# Exact solution is calculated as a function in the input file.

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [gen]
    type = ExamplePatchMeshGenerator
    dim = 3
  []
[]

[AuxVariables]
  [temperature]
  []
  [fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [volumetric_strain]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temperature_ramp]
    type = PiecewiseLinear
    x = '0 1'
    y = '300 1400'
  []
  [fluence_ramp]
    type = ConstantFunction
    value = 2e25
  []
[]

[AuxKernels]
  [temperature_aux]
    type = FunctionAux
    variable = temperature
    function = temperature_ramp
  []
  [fluence_aux]
    type = FunctionAux
    variable = fluence
    function = fluence_ramp
  []
  [volumetric_strain]
    type = RankTwoScalarAux
    rank_two_tensor = total_strain
    variable = volumetric_strain
    scalar_type = VolumetricStrain
  []
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = SMALL
  decomposition_method = EigenSolution
  add_variables = true
  eigenstrain_names = swelling_strain
  use_automatic_differentiation = false
[]

[BCs]
  [left]
    type = DirichletBC
    variable = disp_x
    value = 0
    boundary = left
  []
  [bottom]
    type = DirichletBC
    variable = disp_y
    value = 0
    boundary = bottom
  []
  [back]
    type = DirichletBC
    variable = disp_z
    value = 0
    boundary = back
  []
[]

[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1
    poissons_ratio = 0.3
  []
  [swelling]
    type = MonolithicSiCVolumetricSwellingEigenstrain
    eigenstrain_name = swelling_strain
    temperature = temperature
    fast_neutron_fluence = fluence
    outputs = all
  []
  [stress]
    type = ComputeLinearElasticStress
  []
[]

[Executioner]
  type = Transient
  num_steps = 10
  dt = 0.1
[]

[Postprocessors]
  [temperature_pp]
    type = ElementAverageValue
    variable = temperature
  []
  [fluence_pp]
    type = ElementAverageValue
    variable = fluence
  []
  [volumetric_strain]
    type = ElementAverageValue
    variable = volumetric_strain
  []
  [swelling_exact]
    type = ParsedPostprocessor
    pp_names = 'temperature_pp fluence_pp'
    expression = 'characteristic_fluence := -0.57533 + 3.3342e-3*temperature_pp - 5.3970e-6*temperature_pp^2 + 2.9754e-9*temperature_pp^3;
                saturated_swelling := 5.8366e-2 - 1.0089e-4*temperature_pp + 6.9368e-8*temperature_pp^2 - 1.8152e-11*temperature_pp^3;
                saturated_swelling * (1 - exp(-1 * fluence_pp / (1e25 * characteristic_fluence)))^(2.0/3.0)'
  []
  [swelling_diff]
    type = ParsedPostprocessor
    pp_names = 'volumetric_strain swelling_exact'
    expression = '(volumetric_strain - swelling_exact) / swelling_exact'
    outputs = none
  []
  [swelling_max_diff]
    type = TimeExtremeValue
    postprocessor = swelling_diff
    value_type = abs_max
  []
[]

[Outputs]
  csv = true
[]

(test/tests/solid_mechanics/monolithicSiC_volumetric_swelling/ad_exact.i)

# This tests ADMonolithicSiCVolumetricSwellingEigenstrain as a function of temperature and fluence.
# Exact solution is calculated as a function in the input file.

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [gen]
    type = ExamplePatchMeshGenerator
    dim = 3
  []
[]

[AuxVariables]
  [temperature]
  []
  [fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [volumetric_strain]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temperature_ramp]
    type = PiecewiseLinear
    x = '0 1'
    y = '300 1400'
  []
  [fluence_ramp]
    type = ConstantFunction
    value = 2e25
  []
[]

[AuxKernels]
  [temperature_aux]
    type = FunctionAux
    variable = temperature
    function = temperature_ramp
  []
  [fluence_aux]
    type = FunctionAux
    variable = fluence
    function = fluence_ramp
  []
  [volumetric_strain]
    type = ADRankTwoScalarAux
    rank_two_tensor = total_strain
    variable = volumetric_strain
    scalar_type = VolumetricStrain
  []
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = SMALL
  decomposition_method = EigenSolution
  add_variables = true
  eigenstrain_names = swelling_strain
  use_automatic_differentiation = true
[]

[BCs]
  [left]
    type = ADDirichletBC
    variable = disp_x
    value = 0
    boundary = left
  []
  [bottom]
    type = ADDirichletBC
    variable = disp_y
    value = 0
    boundary = bottom
  []
  [back]
    type = ADDirichletBC
    variable = disp_z
    value = 0
    boundary = back
  []
[]

[Materials]
  [elasticity_tensor]
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 1
    poissons_ratio = 0.3
  []
  [swelling]
    type = ADMonolithicSiCVolumetricSwellingEigenstrain
    eigenstrain_name = swelling_strain
    temperature = temperature
    fast_neutron_fluence = fluence
    outputs = all
  []
  [stress]
    type = ADComputeLinearElasticStress
  []
[]

[Executioner]
  type = Transient
  num_steps = 10
  dt = 0.1
[]

[Postprocessors]
  [temperature_pp]
    type = ElementAverageValue
    variable = temperature
  []
  [fluence_pp]
    type = ElementAverageValue
    variable = fluence
  []
  [volumetric_strain]
    type = ElementAverageValue
    variable = volumetric_strain
  []
  [swelling_exact]
    type = ParsedPostprocessor
    pp_names = 'temperature_pp fluence_pp'
    expression = 'characteristic_fluence := -0.57533 + 3.3342e-3*temperature_pp - 5.3970e-6*temperature_pp^2 + 2.9754e-9*temperature_pp^3;
                saturated_swelling := 5.8366e-2 - 1.0089e-4*temperature_pp + 6.9368e-8*temperature_pp^2 - 1.8152e-11*temperature_pp^3;
                saturated_swelling * (1 - exp(-1 * fluence_pp / (1e25 * characteristic_fluence)))^(2.0/3.0)'
  []
  [swelling_diff]
    type = ParsedPostprocessor
    pp_names = 'volumetric_strain swelling_exact'
    expression = '(volumetric_strain - swelling_exact) / swelling_exact'
    outputs = none
  []
  [swelling_max_diff]
    type = TimeExtremeValue
    postprocessor = swelling_diff
    value_type = abs_max
  []
[]

[Outputs]
  csv = true
[]

Description
Example Input Syntax
Input Parameters
Input Files
References