B4CSwellingEigenstrain

Computes an eigenstrain from an incremental swelling fraction in B4C due to neutron capture.

Description

This computes the total irradiation swelling $var element = document.getElementById("moose-equation-56a43fbf-b09c-43d7-a36a-52299c0a0a0b");katex.render("\\left(\\Delta V / V_0 \\right)", element, {displayMode:false,throwOnError:false});$ for boron carbide (B4C) due to neutron capture. The porosity is then computed from this swelling, though this is an approximation. Two models are provided to calculate the swelling from neutron capture labeled as "Maruyama" (Maruyama et al., 1997) and "covolume" (Qu et al., 2022). The automatic differentiation (AD) version is also available by adding the prefix AD to the name for the "type" parameter.

Maruyama Swelling Data Fit

The fit to swelling data is divided into two regions suggested by Maruyama et al. (1997). SmootherStep is used to smoothly transition from one region to the next from $var element = document.getElementById("moose-equation-8e6cd025-09fb-4207-a1dd-56e612a38f44");katex.render("80 \\times 10^{26}", element, {displayMode:false,throwOnError:false});$ to $var element = document.getElementById("moose-equation-76d36bd1-8f25-458a-8b33-734088d12bbe");katex.render("100 \\times 10^{26}", element, {displayMode:false,throwOnError:false});$ neutron captures per m $var element = document.getElementById("moose-equation-9834ebee-6283-4ac5-a439-a22375457840");katex.render("^3", element, {displayMode:false,throwOnError:false});$ (cap/m $var element = document.getElementById("moose-equation-00c927ec-7b8c-4cd5-b179-97c5537611e5");katex.render("^3", element, {displayMode:false,throwOnError:false});$ ). The swelling fraction for both regions is

$var element = document.getElementById("moose-equation-69d701ba-7d87-49c4-8c64-5a94e0079b58");katex.render(" \\frac{\\Delta V}{V_0} = \\begin{cases} 1.1 \\times 10^{-29} \\Phi, & \\Phi \\le 80 \\times 10^{26} (\\text{cap} / \\text{m}^3) \\\\ 7.432 \\times 10^{-30} \\Phi, & \\Phi \\ge 100 \\times 10^{26} (\\text{cap} / \\text{m}^3) \\end{cases},", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-37c0cc97-fa87-479e-a3df-4aac2b9cb4be");katex.render("\\Phi", element, {displayMode:false,throwOnError:false});$ is the total neutron captures per m $var element = document.getElementById("moose-equation-fd4ace30-6c76-4559-8c0b-3c4e87c8d145");katex.render("^3", element, {displayMode:false,throwOnError:false});$ .

Covolume Approach

The covolume approach is an approximation from Qu et al. (2022). The swelling fraction is calculated as

$var element = document.getElementById("moose-equation-f6e2e73c-7a67-4d29-b6a8-d2aa9682d3ba");katex.render(" \\frac{\\Delta V}{V_0} = \\begin{cases} 0, & \\Phi \\le 10 \\times 10^{26} (\\text{cap} / \\text{m}^3) \\\\ 1.5 \\times 10^{-29} \\left( \\Phi - 10 \\times 10^{26} \\right), & \\Phi > 10 \\times 10^{26} (\\text{cap} / \\text{m}^3) \\end{cases}.", element, {displayMode:true,throwOnError:false});$

Porosity

Porosity, $var element = document.getElementById("moose-equation-222b2a41-260b-4237-91c6-a920afa0d3b2");katex.render("\\phi", element, {displayMode:false,throwOnError:false});$ , is approximated as:

$var element = document.getElementById("moose-equation-386144f0-1520-46e4-8238-813d47530bd5");katex.render(" \\phi = \\frac{\\phi_0 + \\frac{\\Delta V}{V_0}}{1 + \\frac{\\Delta V}{V_0}},", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-30415129-19f3-4c46-895c-b4ee08a1df27");katex.render("\\phi_0", element, {displayMode:false,throwOnError:false});$ is the initial porosity.

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [b4c_swell_mar]
    type = B4CSwellingEigenstrain<<<{"description": "Computes an eigenstrain from an incremental swelling fraction in B4C due to neutron capture.", "href": "B4CSwellingEigenstrain.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."}>>> = irradiation_eigenstrain
    capture_density_rate<<<{"description": "Name of neutron capture density rate material property"}>>> = capture_density_rate
    porosity_name<<<{"description": "Name of porosity material property calculated by neutron capture strain"}>>> = porosity_mar
    capture_density<<<{"description": "Name of material property tracking neutron capture density in cap/m^3."}>>> = capture_density_mar
    swelling_name<<<{"description": "Name of material property set here to the neutron capture strain"}>>> = irr_swell_mar
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
[]

Input Parameters

capture_density_rateName of neutron capture density rate material property
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Name of neutron capture density rate material property
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.

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
capture_densitycapture_densityName of material property tracking neutron capture density in cap/m^3.
Default:capture_density
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Name of material property tracking neutron capture density in cap/m^3.
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.
initial_density2269.55Initial density of B4C in kg/m^3.
Default:2269.55
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Initial density of B4C in kg/m^3.
porosity_nameporosityName of porosity material property calculated by neutron capture strain
Default:porosity
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Name of porosity material property calculated by neutron capture strain
swelling_modelmaruyamaVolumetric swelling model. Choices are: maruyama covolume
Default:maruyama
C++ Type:MooseEnum
Options:maruyama, covolume
Controllable:No
Description:Volumetric swelling model. Choices are: maruyama covolume
swelling_namevolumetric_swelling_strainName of material property set here to the neutron capture strain
Default:volumetric_swelling_strain
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Name of material property set here to the neutron capture strain
theoretical_density2389Theorectical maximum density of B4C in kg/m^3.
Default:2389
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Theorectical maximum density of B4C in kg/m^3.

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

total_swelling_scaling_factor1Scaling factor to be applied to the swelling strain before porosity calcuation. Used for sensitivity and calibration studies
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor to be applied to the swelling strain before porosity calcuation. Used for sensitivity and calibration studies

Advanced: Scaling Factors Parameters

Input Files

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

References

Tadashi Maruyama, Shoji Onose, Takeji Kaito, and Hiroto Horiuchi. Effect of fast neutron irradiation on the properties of boron carbide pellet. Journal of Nuclear Science and Technology, 34(10):1006–1014, 1997. doi:10.1080/18811248.1997.9733777.

@article{MARUYAMA19971006,
    author = "Maruyama, Tadashi and Onose, Shoji and Kaito, Takeji and Horiuchi, Hiroto",
    title = "Effect of Fast Neutron Irradiation on the Properties of Boron Carbide Pellet",
    journal = "Journal of Nuclear Science and Technology",
    volume = "34",
    number = "10",
    pages = "1006-1014",
    year = "1997",
    publisher = "Taylor \& Francis",
    doi = "10.1080/18811248.1997.9733777"
}

Zhixue Qu, Chuanjin Yu, Yitong Wei, Xiping Su, and Aibing Du. Thermal conductivity of boron carbide under fast neutron irradiation. Journal of Advanced Ceramics, 11:482–494, 3 2022. doi:10.1007/s40145-022-0572-8.

@article{Qu2022482,
    author = "Qu, Zhixue and Yu, Chuanjin and Wei, Yitong and Su, Xiping and Du, Aibing",
    title = "Thermal conductivity of boron carbide under fast neutron irradiation",
    doi = "10.1007/s40145-022-0572-8",
    journal = "Journal of Advanced Ceramics",
    month = "3",
    year = "2022",
    volume = "11",
    pages = "482-494"
}

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

# Compares B4CSwellingEigenstrain to exact calculations of each swelling component.

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

[Mesh]
  use_displaced_mesh = false
  [mesh]
    type = ExamplePatchMeshGenerator
    dim = 3
    x_length = 1
    y_length = 1
    z_length = 1
  []
[]

[Functions]
  [func_n_cap]
    type = PiecewiseLinear
    x = '0.0 1.0  2.0  8.0 10.0'
    y = '0   1e27 8e26 1e27 1e28'
  []
[]

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

[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]
  [set_captures]
    type = GenericFunctionMaterial
    prop_names = 'capture_density_rate'
    prop_values = func_n_cap
    outputs = all
  []
  [elasticity_tensor] # Fake values
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 304.4e9
    poissons_ratio = 0.3
  []
  [b4c_swell_mar]
    type = B4CSwellingEigenstrain
    eigenstrain_name = irradiation_eigenstrain
    capture_density_rate = capture_density_rate
    porosity_name = porosity_mar
    capture_density = capture_density_mar
    swelling_name = irr_swell_mar
    outputs = all
  []
  [b4c_swell_vol]
    type = B4CSwellingEigenstrain
    swelling_model = covolume
    eigenstrain_name = irradiation_eigenstrain_vol
    capture_density_rate = capture_density_rate
    porosity_name = porosity_vol
    capture_density = capture_density_vol
    total_swelling_scaling_factor = 1.01
    swelling_name = irr_swell_vol
    outputs = all
  []
  [stress]
    type = ComputeLinearElasticStress
  []
  [calc_solid_swell_mar]
    type = ParsedMaterial
    property_name = 'calc_solid_swelling_mar'
    material_property_names = 'capture_density_mar'
    expression = 'below := 1.1e-29 * capture_density_mar; above := 7.432e-30 * capture_density_mar; x := (capture_density_mar - 8e27)/(1e28 - 8e27); y := if(capture_density_mar <= 8e27, 0, if(capture_density_mar >= 1e28, 1, 6.0 * x^5 - 15.0 * x^4 +10.0 * x^3)); y * above + (1.0 - y) * below'
    outputs = all
  []
  [calc_porosity_mar]
    type = ParsedMaterial
    property_name = 'calc_porosity_mar'
    material_property_names = 'calc_solid_swelling_mar'
    expression = '(0.05 + calc_solid_swelling_mar) / (1.0 + calc_solid_swelling_mar)'
    outputs = all
  []
  [calc_solid_swell_vol]
    type = ParsedMaterial
    property_name = 'calc_solid_swelling_vol'
    material_property_names = 'capture_density_vol'
    expression = 'if(capture_density_vol < 1e27, 0, 1.01 * 1.5e-29 * (capture_density_vol - 1e27))'
    outputs = all
  []
  [calc_porosity_vol]
    type = ParsedMaterial
    property_name = 'calc_porosity_vol'
    material_property_names = 'calc_solid_swelling_vol'
    expression = '(0.05 + calc_solid_swelling_vol) / (1.0 + calc_solid_swelling_vol)'
    outputs = all
  []
[]

[Postprocessors]
  [cap_density_both]
    type = ElementAverageMaterialProperty
    mat_prop = capture_density_mar
  []
  [disp_x_right]
    type = SideAverageValue
    boundary = right
    variable = disp_x
    execute_on = 'initial timestep_end'
  []
  [solid_swell_post_mar]
    type = ElementAverageMaterialProperty
    mat_prop = irr_swell_mar
  []
  [solid_swell_post_vol]
    type = ElementAverageMaterialProperty
    mat_prop = irr_swell_vol
  []
  [porosity_post_mar]
    type = ElementAverageMaterialProperty
    mat_prop = porosity_mar
  []
  [porosity_post_vol]
    type = ElementAverageMaterialProperty
    mat_prop = porosity_vol
  []
  [solid_swell_calc_mar]
    type = ElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_mar
  []
  [solid_swell_calc_vol]
    type = ElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_vol
  []
  [porosity_calc_mar]
    type = ElementAverageMaterialProperty
    mat_prop = calc_porosity_mar
  []
  [porosity_calc_vol]
    type = ElementAverageMaterialProperty
    mat_prop = calc_porosity_vol
  []
[]

[Executioner]
  type = Transient
  num_steps = 20
  dt = 0.5
  nl_abs_tol = 1e-19
[]

[Outputs]
  csv = true
[]

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

# Compares B4CSwellingEigenstrain to exact calculations of each swelling component.

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

[Mesh]
  use_displaced_mesh = false
  [mesh]
    type = ExamplePatchMeshGenerator
    dim = 3
    x_length = 1
    y_length = 1
    z_length = 1
  []
[]

[Functions]
  [func_n_cap]
    type = PiecewiseLinear
    x = '0.0 1.0  2.0  8.0 10.0'
    y = '0   1e27 8e26 1e27 1e28'
  []
[]

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

[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]
  [set_captures]
    type = GenericFunctionMaterial
    prop_names = 'capture_density_rate'
    prop_values = func_n_cap
    outputs = all
  []
  [elasticity_tensor] # Fake values
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 304.4e9
    poissons_ratio = 0.3
  []
  [b4c_swell_mar]
    type = B4CSwellingEigenstrain
    eigenstrain_name = irradiation_eigenstrain
    capture_density_rate = capture_density_rate
    porosity_name = porosity_mar
    capture_density = capture_density_mar
    swelling_name = irr_swell_mar
    outputs = all
  []
  [b4c_swell_vol]
    type = B4CSwellingEigenstrain
    swelling_model = covolume
    eigenstrain_name = irradiation_eigenstrain_vol
    capture_density_rate = capture_density_rate
    porosity_name = porosity_vol
    capture_density = capture_density_vol
    total_swelling_scaling_factor = 1.01
    swelling_name = irr_swell_vol
    outputs = all
  []
  [stress]
    type = ComputeLinearElasticStress
  []
  [calc_solid_swell_mar]
    type = ParsedMaterial
    property_name = 'calc_solid_swelling_mar'
    material_property_names = 'capture_density_mar'
    expression = 'below := 1.1e-29 * capture_density_mar; above := 7.432e-30 * capture_density_mar; x := (capture_density_mar - 8e27)/(1e28 - 8e27); y := if(capture_density_mar <= 8e27, 0, if(capture_density_mar >= 1e28, 1, 6.0 * x^5 - 15.0 * x^4 +10.0 * x^3)); y * above + (1.0 - y) * below'
    outputs = all
  []
  [calc_porosity_mar]
    type = ParsedMaterial
    property_name = 'calc_porosity_mar'
    material_property_names = 'calc_solid_swelling_mar'
    expression = '(0.05 + calc_solid_swelling_mar) / (1.0 + calc_solid_swelling_mar)'
    outputs = all
  []
  [calc_solid_swell_vol]
    type = ParsedMaterial
    property_name = 'calc_solid_swelling_vol'
    material_property_names = 'capture_density_vol'
    expression = 'if(capture_density_vol < 1e27, 0, 1.01 * 1.5e-29 * (capture_density_vol - 1e27))'
    outputs = all
  []
  [calc_porosity_vol]
    type = ParsedMaterial
    property_name = 'calc_porosity_vol'
    material_property_names = 'calc_solid_swelling_vol'
    expression = '(0.05 + calc_solid_swelling_vol) / (1.0 + calc_solid_swelling_vol)'
    outputs = all
  []
[]

[Postprocessors]
  [cap_density_both]
    type = ElementAverageMaterialProperty
    mat_prop = capture_density_mar
  []
  [disp_x_right]
    type = SideAverageValue
    boundary = right
    variable = disp_x
    execute_on = 'initial timestep_end'
  []
  [solid_swell_post_mar]
    type = ElementAverageMaterialProperty
    mat_prop = irr_swell_mar
  []
  [solid_swell_post_vol]
    type = ElementAverageMaterialProperty
    mat_prop = irr_swell_vol
  []
  [porosity_post_mar]
    type = ElementAverageMaterialProperty
    mat_prop = porosity_mar
  []
  [porosity_post_vol]
    type = ElementAverageMaterialProperty
    mat_prop = porosity_vol
  []
  [solid_swell_calc_mar]
    type = ElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_mar
  []
  [solid_swell_calc_vol]
    type = ElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_vol
  []
  [porosity_calc_mar]
    type = ElementAverageMaterialProperty
    mat_prop = calc_porosity_mar
  []
  [porosity_calc_vol]
    type = ElementAverageMaterialProperty
    mat_prop = calc_porosity_vol
  []
[]

[Executioner]
  type = Transient
  num_steps = 20
  dt = 0.5
  nl_abs_tol = 1e-19
[]

[Outputs]
  csv = true
[]

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

# Compares B4CSwellingEigenstrain to exact calculations of each swelling component.

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

[Mesh]
  use_displaced_mesh = false
  [mesh]
    type = ExamplePatchMeshGenerator
    dim = 3
    x_length = 1
    y_length = 1
    z_length = 1
  []
[]

[Functions]
  [func_n_cap]
    type = PiecewiseLinear
    x = '0.0 1.0  2.0  8.0 10.0'
    y = '0   1e27 8e26 1e27 1e28'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    strain = SMALL
    decomposition_method = EigenSolution
    add_variables = true
    eigenstrain_names = irradiation_eigenstrain
    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]
  [set_captures]
    type = ADGenericFunctionMaterial
    prop_names = 'capture_density_rate'
    prop_values = func_n_cap
    outputs = all
  []
  [elasticity_tensor] # Fake values
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 304.4e9
    poissons_ratio = 0.3
  []
  [b4c_swell_mar]
    type = ADB4CSwellingEigenstrain
    eigenstrain_name = irradiation_eigenstrain
    capture_density_rate = capture_density_rate
    porosity_name = porosity_mar
    capture_density = capture_density_mar
    swelling_name = irr_swell_mar
    outputs = all
  []
  [b4c_swell_vol]
    type = ADB4CSwellingEigenstrain
    swelling_model = covolume
    eigenstrain_name = irradiation_eigenstrain_vol
    capture_density_rate = capture_density_rate
    porosity_name = porosity_vol
    capture_density = capture_density_vol
    total_swelling_scaling_factor = 1.01
    swelling_name = irr_swell_vol
    outputs = all
  []
  [stress]
    type = ADComputeLinearElasticStress
  []
  [calc_solid_swell_mar]
    type = ADParsedMaterial
    property_name = 'calc_solid_swelling_mar'
    material_property_names = 'capture_density_mar'
    expression = 'below := 1.1e-29 * capture_density_mar; above := 7.432e-30 * capture_density_mar; x := (capture_density_mar - 8e27)/(1e28 - 8e27); y := if(capture_density_mar <= 8e27, 0, if(capture_density_mar >= 1e28, 1, 6.0 * x^5 - 15.0 * x^4 +10.0 * x^3)); y * above + (1.0 - y) * below'
    outputs = all
  []
  [calc_porosity_mar]
    type = ADParsedMaterial
    property_name = 'calc_porosity_mar'
    material_property_names = 'calc_solid_swelling_mar'
    expression = '(0.05 + calc_solid_swelling_mar) / (1.0 + calc_solid_swelling_mar)'
    outputs = all
  []
  [calc_solid_swell_vol]
    type = ADParsedMaterial
    property_name = 'calc_solid_swelling_vol'
    material_property_names = 'capture_density_vol'
    expression = 'if(capture_density_vol < 1e27, 0, 1.01 * 1.5e-29 * (capture_density_vol - 1e27))'
    outputs = all
  []
  [calc_porosity_vol]
    type = ADParsedMaterial
    property_name = 'calc_porosity_vol'
    material_property_names = 'calc_solid_swelling_vol'
    expression = '(0.05 + calc_solid_swelling_vol) / (1.0 + calc_solid_swelling_vol)'
    outputs = all
  []
[]

[Postprocessors]
  [cap_density_both]
    type = ADElementAverageMaterialProperty
    mat_prop = capture_density_mar
  []
  [disp_x_right]
    type = SideAverageValue
    boundary = right
    variable = disp_x
    execute_on = 'initial timestep_end'
  []
  [solid_swell_post_mar]
    type = ADElementAverageMaterialProperty
    mat_prop = irr_swell_mar
  []
  [solid_swell_post_vol]
    type = ADElementAverageMaterialProperty
    mat_prop = irr_swell_vol
  []
  [porosity_post_mar]
    type = ADElementAverageMaterialProperty
    mat_prop = porosity_mar
  []
  [porosity_post_vol]
    type = ADElementAverageMaterialProperty
    mat_prop = porosity_vol
  []
  [solid_swell_calc_mar]
    type = ADElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_mar
  []
  [solid_swell_calc_vol]
    type = ADElementAverageMaterialProperty
    mat_prop = calc_solid_swelling_vol
  []
  [porosity_calc_mar]
    type = ADElementAverageMaterialProperty
    mat_prop = calc_porosity_mar
  []
  [porosity_calc_vol]
    type = ADElementAverageMaterialProperty
    mat_prop = calc_porosity_vol
  []
[]

[Executioner]
  type = Transient
  num_steps = 20
  dt = 0.5
  nl_abs_tol = 1e-19
[]

[Outputs]
  csv = true
[]

Description
Example Input Syntax
Input Parameters
Input Files
References