GraphiteGradeIrradiationEigenstrain

Computes irradiation-induced dimensional changes (IIDC) for various graphite grades such as H-451 (anisotropic) and IG-110 (isotropic).

Description

Two modeling options are available in the code for computation of the irradiation-induced dimensional change for various graphite types (e.g., H-451 and IG-110). Following sections detail the material models of each graphite type.

note:Constant temperature behavior

These models do not include any correction for changing temperature. The response is as though the material experienced the current temperature for its entire history.

Grade H-451 graphite

Grade H-451 graphite is manufactured in cylindrical shapes, and it is recommended that the material be used only with cylindrical geometries. Expressions for radial and axial irradiation strains are provided.

The irradiation-induced dimensional change (IIDC), $var element = document.getElementById("moose-equation-30d631f4-d0ec-4f97-99db-6427392f34b3");katex.render("\\varepsilon_{i}", element, {displayMode:false,throwOnError:false});$ (%) ( $var element = document.getElementById("moose-equation-2a17765b-197e-4c8d-8033-8e8c1da95b5e");katex.render("i=", element, {displayMode:false,throwOnError:false});$ radial, axial) of grade H-451 graphite is expressed by Ho (1991) as: $var element = document.getElementById("moose-equation-abea97e9-616c-4643-a58e-2ff26600deab");katex.render(" \\begin{split} \\varepsilon_i & = \\left( c_1 + c_2 T_C+ c_3 T_C^2 + c_4 T_C^3 + c_5 T_C^4 \\right) \\phi \\\\ & + \\left( c_6 + c_7 T_C + c_8 T_C^2 + c_9 T_C^3 + c_{10} T_C^4 + c_{16} T_C^5 \\right) \\phi^2 \\\\ & + \\left( c_{11} + c_{12} T_C + c_{13} T_C^2 + c_{14} T_C^3 + c_{15} T_C^4 + c_{17} T_C^5 + c_{18} T_C^6\\right) \\phi^3 \\end{split}", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-a50cad75-f93a-44b3-a1bd-779afc4624a8");katex.render("T_C", element, {displayMode:false,throwOnError:false});$ ( $var element = document.getElementById("moose-equation-78c76074-d8b5-4b92-a58d-f94ed970b227");katex.render("\\degree{C}", element, {displayMode:false,throwOnError:false});$ ) is the irradiation temperature, and $var element = document.getElementById("moose-equation-4b065679-ed5b-4396-bfb2-cbfeb09e3f4d");katex.render("\\phi", element, {displayMode:false,throwOnError:false});$ (10 $var element = document.getElementById("moose-equation-a049a550-2117-4a58-99b5-1489f61c8a0f");katex.render("^{25}", element, {displayMode:false,throwOnError:false});$ n/m $var element = document.getElementById("moose-equation-7dd94f6f-5919-4230-9754-24d928c38953");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-d25ac60c-ba1f-4768-80b5-a664dcf992b3");katex.render("E>", element, {displayMode:false,throwOnError:false});$ 0.18 MeV) is the fast neutron fluence. The empirical coefficients ( $var element = document.getElementById("moose-equation-c030d267-9819-4680-a2b7-2cec0a25f2ec");katex.render("c_i", element, {displayMode:false,throwOnError:false});$ for $var element = document.getElementById("moose-equation-6fd04ea3-8793-4eba-82dd-dbe148a4a678");katex.render("i=1,2,\\ldots,18", element, {displayMode:false,throwOnError:false});$ ) are tabulated in Table 1.

This model is considered valid over the temperature range of 623-1573 K and to fast neutron fluence up to 10 $var element = document.getElementById("moose-equation-1b119c17-9684-4565-92bf-2068505b246e");katex.render("\\times 10^{25}", element, {displayMode:false,throwOnError:false});$ n/m $var element = document.getElementById("moose-equation-f4e3a8ea-cfbd-4dda-a542-3a6ce3d92af8");katex.render("^2", element, {displayMode:false,throwOnError:false});$ . The correlation uses fast fluence with neutron energy threshold $var element = document.getElementById("moose-equation-a0c3a456-284c-4b8b-a3b0-7f1d9d0c5ccd");katex.render("E>", element, {displayMode:false,throwOnError:false});$ 0.18 MeV. The model performs the fast fluence conversion from $var element = document.getElementById("moose-equation-4e822db4-5359-4399-ae61-6d67fb6a6256");katex.render("E>", element, {displayMode:false,throwOnError:false});$ 0.10 MeV to $var element = document.getElementById("moose-equation-b3347dca-527e-46aa-a861-bf2f8d3a8ddd");katex.render("E>", element, {displayMode:false,throwOnError:false});$ 0.18 MeV using flux conversion factor.

Table 1: Empirical coefficients, $var element = document.getElementById("moose-equation-1788d514-ce5d-4baa-ad73-6ee9d788d3fc");katex.render("c_i", element, {displayMode:false,throwOnError:false});$ for dimensional change of H-451 graphite (Ho, 1991).

$var element = document.getElementById("moose-equation-bb4d7d86-de78-4ac0-9e5e-e52cca917336");katex.render("\\mathbf{i}", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-d9725d6f-2f08-426b-9a68-7e4afe570c45");katex.render("\\mathbf{c_i}", element, {displayMode:false,throwOnError:false});$ (axial)	$var element = document.getElementById("moose-equation-bca84881-d476-46da-8466-de3ffa3f7777");katex.render("\\mathbf{c_i}", element, {displayMode:false,throwOnError:false});$ (radial)	$var element = document.getElementById("moose-equation-25eee83b-406b-4715-bafe-344cfa7f45dc");katex.render("\\mathbf{i}", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-a6094b22-e118-4499-ad53-6a39fb34adb1");katex.render("\\mathbf{c_i}", element, {displayMode:false,throwOnError:false});$ (axial)	$var element = document.getElementById("moose-equation-d118f80b-f723-4c56-9e88-39b2062a9cd2");katex.render("\\mathbf{c_i}", element, {displayMode:false,throwOnError:false});$ (radial)
1	1.11617	1.15132	10	0.20435 $var element = document.getElementById("moose-equation-223c6a1f-612c-46ff-b2f1-65feabf31ad4");katex.render("\\times 10^{-10}", element, {displayMode:false,throwOnError:false});$	0.14233 $var element = document.getElementById("moose-equation-2665b677-abe0-47aa-a858-85986931d543");katex.render("\\times 10^{-10}", element, {displayMode:false,throwOnError:false});$
2	-0.92197 $var element = document.getElementById("moose-equation-02a9eccf-57e9-44ee-b5dc-806596425b8c");katex.render("\\times 10^{-2}", element, {displayMode:false,throwOnError:false});$	-0.82968 $var element = document.getElementById("moose-equation-e7b36d50-2547-43f9-8413-aabdefede26d");katex.render("\\times 10^{-2}", element, {displayMode:false,throwOnError:false});$	11	-0.59274 $var element = document.getElementById("moose-equation-df685478-b243-4c3a-b00b-fd5c509d0ee4");katex.render("\\times 10^{-2}", element, {displayMode:false,throwOnError:false});$	0.13110 $var element = document.getElementById("moose-equation-06cb9639-b5d4-422f-a380-cc6df7059f7e");katex.render("\\times 10^{-1}", element, {displayMode:false,throwOnError:false});$
3	0.20463 $var element = document.getElementById("moose-equation-08958dd5-8ea0-4865-bfed-d06f1aa1e02b");katex.render("\\times 10^{-4}", element, {displayMode:false,throwOnError:false});$	0.17060 $var element = document.getElementById("moose-equation-dd02da37-0abf-47cd-a161-35b08f9b9267");katex.render("\\times 10^{-4}", element, {displayMode:false,throwOnError:false});$	12	-0.65404 $var element = document.getElementById("moose-equation-1de71220-8060-496e-8ada-453cfde5044d");katex.render("\\times 10^{-5}", element, {displayMode:false,throwOnError:false});$	-0.18768 $var element = document.getElementById("moose-equation-a168ac4f-7942-41ed-9edd-a8704f7f36f5");katex.render("\\times 10^{-3}", element, {displayMode:false,throwOnError:false});$
4	-0.16458 $var element = document.getElementById("moose-equation-32fbe47d-f580-4ad9-8b0b-1e202dd654ee");katex.render("\\times 10^{-7}", element, {displayMode:false,throwOnError:false});$	-0.12645 $var element = document.getElementById("moose-equation-c6c76c1e-7e14-4254-8850-4d3fa7f611ab");katex.render("\\times 10^{-7}", element, {displayMode:false,throwOnError:false});$	13	0.32751 $var element = document.getElementById("moose-equation-1f451379-c021-43a3-8f90-267e5399d6ce");katex.render("\\times 10^{-6}", element, {displayMode:false,throwOnError:false});$	0.10199 $var element = document.getElementById("moose-equation-b5a85784-6ceb-4f0b-8627-8dac98ac9ec7");katex.render("\\times 10^{-5}", element, {displayMode:false,throwOnError:false});$
5	0.40809 $var element = document.getElementById("moose-equation-c92c797a-b80c-47fa-b0cc-28734c7fb060");katex.render("\\times 10^{-11}", element, {displayMode:false,throwOnError:false});$	0.27657 $var element = document.getElementById("moose-equation-8e221c07-0e0d-49ae-8314-067e837d22dc");katex.render("\\times 10^{-11}", element, {displayMode:false,throwOnError:false});$	14	-0.13449 $var element = document.getElementById("moose-equation-144feeaa-cfbe-46f9-9673-24a4e84fd883");katex.render("\\times 10^{-8}", element, {displayMode:false,throwOnError:false});$	-0.27004 $var element = document.getElementById("moose-equation-d716fade-1571-40fa-9bc1-bf4fcb00daa8");katex.render("\\times 10^{-8}", element, {displayMode:false,throwOnError:false});$
6	0.64947	0.39177	15	0.22245 $var element = document.getElementById("moose-equation-c1ba1069-250e-4386-85a6-674dcc68d71d");katex.render("\\times 10^{-11}", element, {displayMode:false,throwOnError:false});$	0.36498 $var element = document.getElementById("moose-equation-6a5180b7-c1f8-42b3-9e7f-9718ff80cb60");katex.render("\\times 10^{-11}", element, {displayMode:false,throwOnError:false});$
7	-0.56929 $var element = document.getElementById("moose-equation-50e7b810-437e-4e0b-92cb-852877c934ae");katex.render("\\times 10^{-2}", element, {displayMode:false,throwOnError:false});$	-0.36540 $var element = document.getElementById("moose-equation-5c4d7630-f36c-4df3-b85d-38b6f35d1b40");katex.render("\\times 10^{-2}", element, {displayMode:false,throwOnError:false});$	16	-0.51973 $var element = document.getElementById("moose-equation-f61b5417-4756-4ff5-871e-19a2b87fac32");katex.render("\\times 10^{-14}", element, {displayMode:false,throwOnError:false});$	-0.35768 $var element = document.getElementById("moose-equation-2fa0f766-9d3a-4041-a1d0-da812eea5f85");katex.render("\\times 10^{-14}", element, {displayMode:false,throwOnError:false});$
8	0.18972 $var element = document.getElementById("moose-equation-c7827b5e-45ba-4131-92dd-3c48fbb3ba87");katex.render("\\times 10^{-4}", element, {displayMode:false,throwOnError:false});$	0.12750 $var element = document.getElementById("moose-equation-be18792e-be2c-4cf2-abe3-b6fd404b9d47");katex.render("\\times 10^{-4}", element, {displayMode:false,throwOnError:false});$	17	-0.16038 $var element = document.getElementById("moose-equation-a4c1630c-b2e1-4636-b03d-0246082968b7");katex.render("\\times 10^{-14}", element, {displayMode:false,throwOnError:false});$	-0.23579 $var element = document.getElementById("moose-equation-62bb32ab-82e2-4162-8f9b-4a20e3d82957");katex.render("\\times 10^{-14}", element, {displayMode:false,throwOnError:false});$
9	-0.29277 $var element = document.getElementById("moose-equation-38257267-0d96-445c-adaa-d0c1633dc6da");katex.render("\\times 10^{-7}", element, {displayMode:false,throwOnError:false});$	-0.20230 $var element = document.getElementById("moose-equation-5879df34-6e43-469c-8993-f1afc9f804e7");katex.render("\\times 10^{-7}", element, {displayMode:false,throwOnError:false});$	18	0.41756 $var element = document.getElementById("moose-equation-de97f12b-34b4-4170-864f-6c25c54acb4d");katex.render("\\times 10^{-18}", element, {displayMode:false,throwOnError:false});$	0.57329 $var element = document.getElementById("moose-equation-df229dbd-a76c-4b05-80ea-aa78880f316d");katex.render("\\times 10^{-18}", element, {displayMode:false,throwOnError:false});$

Grade IG-110 graphite

The IIDC, $var element = document.getElementById("moose-equation-f96fab6c-ab57-4d5e-a970-3485053967c3");katex.render("\\varepsilon_{i}", element, {displayMode:false,throwOnError:false});$ (%) of the grade IG-110 graphite is expressed in Mohanty and Majumdar (2011) as: $var element = document.getElementById("moose-equation-9087116c-bb32-41de-835e-3157757ee3bb");katex.render(" \\varepsilon_i = a_1 \\phi^2 + a_2\\phi", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-72cc1d1b-ff26-4e64-88c3-89daebccb69a");katex.render("\\phi", element, {displayMode:false,throwOnError:false});$ (10 $var element = document.getElementById("moose-equation-701b3719-a91b-44d2-83c4-bee479a326d3");katex.render("^{26}", element, {displayMode:false,throwOnError:false});$ n/m $var element = document.getElementById("moose-equation-3a2f712d-1348-4594-9207-828986e338e6");katex.render("^2", element, {displayMode:false,throwOnError:false});$ ) is the fast neutron fluence and $var element = document.getElementById("moose-equation-4e581855-efc3-4179-aa7f-40c76ce16e3f");katex.render("a_1", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-dc2479f7-65df-4b80-aa9f-54a5e6e676ab");katex.render("a_2", element, {displayMode:false,throwOnError:false});$ are the empirical coefficients (see Table 2).

Table 2: Empirical coefficients, $var element = document.getElementById("moose-equation-652510ad-4226-4366-99f7-9eaf888db4ca");katex.render("a_1", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-64a6de0a-9f88-415b-8c98-a33988486ca3");katex.render("a_2", element, {displayMode:false,throwOnError:false});$ for dimensional change of IG-110 graphite (Mohanty and Majumdar, 2011).

	400 $var element = document.getElementById("moose-equation-a5416b23-3514-4da5-b08c-ced5e3d52ae4");katex.render("\\degree{C}", element, {displayMode:false,throwOnError:false});$	600 $var element = document.getElementById("moose-equation-d428714b-a891-488b-84d1-b5b683872e44");katex.render("\\degree{C}", element, {displayMode:false,throwOnError:false});$	800 $var element = document.getElementById("moose-equation-daa5392c-db09-457d-ba28-323c0958a15e");katex.render("\\degree{C}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-aa08b0ea-5010-4b70-ad50-a89620392f3c");katex.render("a_1", element, {displayMode:false,throwOnError:false});$	0.279	0.450	0.821
$var element = document.getElementById("moose-equation-e041b933-e04f-4f7d-b244-34fbd77c43a3");katex.render("a_2", element, {displayMode:false,throwOnError:false});$	-1.64	-1.86	-2.19

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [graphite_irrad_strain]
    type = GraphiteGradeIrradiationEigenstrain<<<{"description": "Computes irradiation-induced dimensional changes (IIDC) for various graphite grades such as H-451 (anisotropic) and IG-110 (isotropic).", "href": "GraphiteGradeIrradiationEigenstrain.html"}>>>
    temperature<<<{"description": "Coupled temperature"}>>> = temperature
    graphite_grade<<<{"description": "The grade of graphite H_451 IG_110"}>>> = H_451
    flux_conversion_factor<<<{"description": "Convert fast neutron flux E>0.10 to E>0.18 MeV"}>>> = 1.0
    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."}>>> = irrad_strain
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = exodus
  []
[]

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.
flux_conversion_factorConvert fast neutron flux E>0.10 to E>0.18 MeV
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Convert fast neutron flux E>0.10 to E>0.18 MeV
temperatureCoupled temperature
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled temperature

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.
fast_neutron_fluencefast_neutron_fluenceCoupled fast (E>0.18 MeV) neutron fluence (n/m^2)
Default:fast_neutron_fluence
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Coupled fast (E>0.18 MeV) neutron fluence (n/m^2)
graphite_gradeH_451The grade of graphite H_451 IG_110
Default:H_451
C++ Type:MooseEnum
Options:H_451, IG_110
Controllable:No
Description:The grade of graphite H_451 IG_110

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

irradiation_eigenstrain_scale_factor1Scale factor for graphite grade IIDC
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scale factor for graphite grade IIDC

Advanced: Scaling Factors 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/graphite_eigenstrains/irradiation_strain/ig110_test.i)
(test/tests/solid_mechanics/graphite_eigenstrains/irradiation_strain/h451_test.i)

References

F.H. Ho. Graphite design handbook. Technical Report DOE-HTGR-88111, General Atomics, 1991. URL: https://www.osti.gov/servlets/purl/714896, doi:10.2172/714896.

@techreport{GraphiteDesignHandbook,
    author = "Ho, F.H.",
    title = "Graphite Design Handbook",
    institution = "General Atomics",
    number = "DOE-HTGR-88111",
    doi = "10.2172/714896",
    url = "https://www.osti.gov/servlets/purl/714896",
    year = "1991"
}

S. Mohanty and S. Majumdar. HTGR graphite core component stress analysis research program - task 1. Technical Letter Report ANL-11/04, Argonne National Laboratory, 9 2011. URL: https://www.nrc.gov/docs/ML1127/ML11276A009.pdf.

@techreport{nrcHTGR2011,
    author = "Mohanty, S. and Majumdar, S.",
    title = "{HTGR} Graphite Core Component Stress Analysis Research Program - Task 1",
    institution = "Argonne National Laboratory",
    type = "Technical Letter Report",
    number = "ANL-11/04",
    month = "9",
    url = "https://www.nrc.gov/docs/ML1127/ML11276A009.pdf",
    year = "2011"
}

(test/tests/solid_mechanics/graphite_eigenstrains/irradiation_strain/h451_test.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  coord_type = RZ
  use_displaced_mesh = false
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    elem_type = quad8
    xmin = 0.039
    xmax = 0.064
    ymin = 0
    ymax = 0.025
  []
[]

[Variables]
  [disp_x]
    order = SECOND
    family = LAGRANGE
  []
  [disp_y]
    order = SECOND
    family = LAGRANGE
  []
[]

[AuxVariables]
  [temperature]
    initial_condition = 1073.15
  []
[]

[Functions]
  [fluence]
    type = PiecewiseLinear
    x = '0 10'
    y = '0 10e25'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    strain = small
    incremental = true
    eigenstrain_names = irrad_strain
  []
[]

[BCs]
  [x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_values = fluence
    prop_names = fast_neutron_fluence
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    normal_x = 1
    normal_y = 0
    outputs = exodus
  []
  [graphite_irrad_strain]
    type = GraphiteGradeIrradiationEigenstrain
    temperature = temperature
    graphite_grade = H_451
    flux_conversion_factor = 1.0
    eigenstrain_name = irrad_strain
    outputs = exodus
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1e6
    poissons_ratio = 0.3
  []
[]

[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton

  petsc_options = '-snes_ksp '
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'
  line_search = none

  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6

  start_time = 0.0
  dt = 0.1
  end_time = 10
[]

[Outputs]
  file_base = 'H_451_out'
  exodus = true
[]

(test/tests/solid_mechanics/graphite_eigenstrains/irradiation_strain/ig110_test.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  coord_type = RZ
  use_displaced_mesh = false
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    elem_type = quad8
    xmin = 0.039
    xmax = 0.064
    ymin = 0
    ymax = 0.025
  []
[]

[Variables]
  [disp_x]
    order = SECOND
    family = LAGRANGE
  []
  [disp_y]
    order = SECOND
    family = LAGRANGE
  []
[]

[AuxVariables]
  [temperature]
    initial_condition = 673.15
  []
[]

[Functions]
  [fluence]
    type = PiecewiseLinear
    x = '0 10'
    y = '0 10e26'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    strain = small
    incremental = true
    eigenstrain_names = irrad_strain
  []
[]

[BCs]
  [x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_values = fluence
    prop_names = fast_neutron_fluence
  []
  [graphite_irrad_strain]
    type = GraphiteGradeIrradiationEigenstrain
    temperature = temperature
    graphite_grade = IG_110
    flux_conversion_factor = 1.0
    eigenstrain_name = irrad_strain
    outputs = exodus
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1e6
    poissons_ratio = 0.3
  []
[]

[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton

  petsc_options = '-snes_ksp '
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'
  line_search = none

  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6

  start_time = 0.0
  dt = 0.1
  end_time = 6
[]

[Outputs]
  file_base = 'IG_110_out'
  exodus = true
[]

(test/tests/solid_mechanics/graphite_eigenstrains/irradiation_strain/h451_test.i)

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  coord_type = RZ
  use_displaced_mesh = false
  [mesh]
    type = GeneratedMeshGenerator
    dim = 2
    elem_type = quad8
    xmin = 0.039
    xmax = 0.064
    ymin = 0
    ymax = 0.025
  []
[]

[Variables]
  [disp_x]
    order = SECOND
    family = LAGRANGE
  []
  [disp_y]
    order = SECOND
    family = LAGRANGE
  []
[]

[AuxVariables]
  [temperature]
    initial_condition = 1073.15
  []
[]

[Functions]
  [fluence]
    type = PiecewiseLinear
    x = '0 10'
    y = '0 10e25'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    strain = small
    incremental = true
    eigenstrain_names = irrad_strain
  []
[]

[BCs]
  [x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_values = fluence
    prop_names = fast_neutron_fluence
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    normal_x = 1
    normal_y = 0
    outputs = exodus
  []
  [graphite_irrad_strain]
    type = GraphiteGradeIrradiationEigenstrain
    temperature = temperature
    graphite_grade = H_451
    flux_conversion_factor = 1.0
    eigenstrain_name = irrad_strain
    outputs = exodus
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 1e6
    poissons_ratio = 0.3
  []
[]

[Preconditioning]
  [SMP]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton

  petsc_options = '-snes_ksp '
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'
  line_search = none

  l_max_its = 60
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6

  start_time = 0.0
  dt = 0.1
  end_time = 10
[]

[Outputs]
  file_base = 'H_451_out'
  exodus = true
[]

Description
Example Input Syntax
Input Parameters
Input Files
References