KernelMigrationFailureIndicator

Determines whether the SiC layer fails based on kernel migration distance.

Description

KernelMigrationFailureIndicator determines whether or not the SiC layer failure occurs due to kernel migration. The SiC is assumed to fail once the kernel has migrated to make contact with the SiC layer, i.e., $var element = document.getElementById("moose-equation-08753b39-f03b-467b-aad3-a0c8269cebff");katex.render("kernel_{radius}+kernel\\_migration\\_distance \\ge SiC_{inner\\_radius}", element, {displayMode:false,throwOnError:false});$ . To calculate the kernel and SiC inner radius at deformed configuration, execute_on = 'INITIAL TIMESTEP_END' must be added in TRISOGeometry UserObject.

Example Input Syntax

[Postprocessors<<<{"href": "../../syntax/Postprocessors/index.html"}>>>]
  [ave_gas_temp]
    type = ElementAverageValue<<<{"description": "Computes the volumetric average of a variable", "href": "ElementAverageValue.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = buffer
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temperature
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [fis_gas_released]
    type = ElementIntegralMaterialProperty<<<{"description": "Compute the integral of the material property over the domain", "href": "ElementIntegralMaterialProperty.html"}>>>
    mat_prop<<<{"description": "The name of the material property"}>>> = fis_gas_released
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    use_displaced_mesh<<<{"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."}>>> = false
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [gap_volume]
    type = InternalVolume<<<{"description": "Computes the volume of an enclosed area by performing an integral over a user-supplied boundary.", "href": "InternalVolume.html"}>>>
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = buffer_IPyC_boundary
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial linear'
    use_displaced_mesh<<<{"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."}>>> = true
  []
  [buffer_void_volume]
    type = VoidVolume<<<{"description": "Computes void volume based on actual and theoretical density.", "href": "VoidVolume.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = buffer
    theoretical_density<<<{"description": "Theoretical density of the material (Postprocessor)."}>>> = 2250
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
    use_displaced_mesh<<<{"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."}>>> = true
  []
  [kernel_th_density]
    type = UCOTheoreticalDensity<<<{"description": "Computes UCO theoretical density.", "href": "UCOTheoreticalDensity.html"}>>>
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = initial
  []
  [kernel_void_volume]
    type = VoidVolume<<<{"description": "Computes void volume based on actual and theoretical density.", "href": "VoidVolume.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    theoretical_density<<<{"description": "Theoretical density of the material (Postprocessor)."}>>> = kernel_th_density
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
    use_displaced_mesh<<<{"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."}>>> = true
  []
  [particle_power]
    type = ElementIntegralPower<<<{"description": "Computes the power given the fission rate and energy per fission.", "href": "ElementIntegralPower.html"}>>>
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temperature
    use_material_fission_rate<<<{"description": "Flag to use the material 'fission_rate_material' instead of variable fission rate"}>>> = true
    fission_rate_material<<<{"description": "Fission rate material property name"}>>> = fission_rate
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [max_fluence]
    type = ElementExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "ElementExtremeValue.html"}>>>
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = fast_neutron_fluence
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = 'max'
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [max_burnup]
    type = ElementExtremeValue<<<{"description": "Finds either the min or max elemental value of a variable over the domain.", "href": "ElementExtremeValue.html"}>>>
    variable<<<{"description": "The name of the variable that this postprocessor operates on"}>>> = burnup
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = fuel
    value_type<<<{"description": "Type of extreme value to return. 'max' returns the maximum value. 'min' returns the minimum value. 'max_abs' returns the maximum of the absolute value."}>>> = 'max'
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [SiC_stress_max]
    type = ElementExtremeMaterialProperty<<<{"description": "Determines the minimum or maximum of a material property over a volume.", "href": "ElementExtremeMaterialProperty.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    value_type<<<{"description": "Type of extreme value to return: 'max' returns the maximum value and 'min' returns the minimum value."}>>> = max
    mat_prop<<<{"description": "Material property for which to find the extreme"}>>> = stress_yy
  []
  [SiC_stress_min]
    type = ElementExtremeMaterialProperty<<<{"description": "Determines the minimum or maximum of a material property over a volume.", "href": "ElementExtremeMaterialProperty.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    value_type<<<{"description": "Type of extreme value to return: 'max' returns the maximum value and 'min' returns the minimum value."}>>> = min
    mat_prop<<<{"description": "Material property for which to find the extreme"}>>> = stress_yy
  []
  [strength_SiC]
    type = WeibullEffectiveMeanStrength<<<{"description": "Computes Weibull effective mean strength, which is used by the WeibullFailureOutputUsingCorrelation Postprocessor.", "href": "WeibullEffectiveMeanStrength.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
  []
  [kernel_migration_distance]
    type = KernelMigrationDistance<<<{"description": "Calculates the kernel migration distance.", "href": "KernelMigrationDistance.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 'fuel buffer IPyC SiC OPyC'
    variable<<<{"description": "The name of the variable that this object operates on"}>>> = temperature
    temperature_gradient<<<{"description": "The Postprocessor that will give the temperature gradient across the particle."}>>> = 15000
    kernel_type<<<{"description": "Kernel type"}>>> = UO2
  []
  [failure_indicator_kernel_migration]
    type = KernelMigrationFailureIndicator<<<{"description": "Determines whether the SiC layer fails based on kernel migration distance.", "href": "KernelMigrationFailureIndicator.html"}>>>
    kernel_migration_distance<<<{"description": "The name of the postprocessor that computes kernel migration distance."}>>> = kernel_migration_distance
    triso_geometry<<<{"description": "User object containing the TrisoGeometry"}>>> = particle_geometry
  []
  [failure_indicator_SiC]
    type = WeibullFailureOutputUsingCorrelation<<<{"description": "Computes failure indicator of a TRISO layer using a correlation function.", "href": "WeibullFailureOutputUsingCorrelation.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
    stress_name<<<{"description": "Name of the stress property."}>>> = max_principal_stress
    effective_mean_strength<<<{"description": "The name of the postprocessor calculating the effective mean strength."}>>> = strength_SiC
  []
  [strength_IPyC]
    type = WeibullEffectiveMeanStrength<<<{"description": "Computes Weibull effective mean strength, which is used by the WeibullFailureOutputUsingCorrelation Postprocessor.", "href": "WeibullEffectiveMeanStrength.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = IPyC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
  []
  [failure_indicator_IPyC]
    type = WeibullFailureOutputUsingCorrelation<<<{"description": "Computes failure indicator of a TRISO layer using a correlation function.", "href": "WeibullFailureOutputUsingCorrelation.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = IPyC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
    stress_name<<<{"description": "Name of the stress property."}>>> = max_principal_stress
    effective_mean_strength<<<{"description": "The name of the postprocessor calculating the effective mean strength."}>>> = strength_IPyC
  []
  [strength_OPyC]
    type = WeibullEffectiveMeanStrength<<<{"description": "Computes Weibull effective mean strength, which is used by the WeibullFailureOutputUsingCorrelation Postprocessor.", "href": "WeibullEffectiveMeanStrength.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = OPyC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
  []
  [failure_indicator_OPyC]
    type = WeibullFailureOutputUsingCorrelation<<<{"description": "Computes failure indicator of a TRISO layer using a correlation function.", "href": "WeibullFailureOutputUsingCorrelation.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = OPyC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
    stress_name<<<{"description": "Name of the stress property."}>>> = max_principal_stress
    effective_mean_strength<<<{"description": "The name of the postprocessor calculating the effective mean strength."}>>> = strength_OPyC
  []
  [pd_penetration]
    type = PdPenetration<<<{"description": "Penetration depth of palladium into the SiC layer of TRISO fuel particles", "href": "PdPenetration.html"}>>>
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = SiC_inner_boundary
    variable<<<{"description": "The name of the variable which this postprocessor integrates"}>>> = temperature
    execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = 'initial timestep_end'
  []
  [failure_indicator_pd_penetration]
    type = PdPenetrationFailureIndicator<<<{"description": "Determines whether the SiC layer fails based on palladium penetration depth.", "href": "PdPenetrationFailureIndicator.html"}>>>
    triso_geometry<<<{"description": "User object containing the TrisoGeometry"}>>> = particle_geometry
    pd_penetration<<<{"description": "The name of the postprocessor that computes palladium penetration depth."}>>> = pd_penetration
  []
  [failure_indicator_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation<<<{"description": "Computes failure indicator of a TRISO layer using a correlation function.", "href": "WeibullFailureOutputUsingCorrelation.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
    stress_name<<<{"description": "Name of the stress property."}>>> = max_principal_stress
    high_fidelity_analysis_strength<<<{"description": "The effective mean strength obtained from a high-fidelity analysis."}>>> = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function<<<{"description": "The stress correlation function obtained from a high-fidelity analysis."}>>> = 'stress_correlation_crackedIPyC'
  []
  [failure_indicator_SiC_crackedOPyC]
    type = WeibullFailureOutputUsingCorrelation<<<{"description": "Computes failure indicator of a TRISO layer using a correlation function.", "href": "WeibullFailureOutputUsingCorrelation.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = SiC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 6
    stress_name<<<{"description": "Name of the stress property."}>>> = max_principal_stress
    high_fidelity_analysis_strength<<<{"description": "The effective mean strength obtained from a high-fidelity analysis."}>>> = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function<<<{"description": "The stress correlation function obtained from a high-fidelity analysis."}>>> = 'stress_correlation_crackedOPyC'
  []
  [triso_failure]
    type = TRISOFailureEvaluation<<<{"description": "Computes particle failure based on Weibull statistical theory. It returns either 1 or 0, indicating fails or not fail.", "href": "TRISOFailureEvaluation.html"}>>>
    IPyC_failure<<<{"description": "The name of the postprocessor calculating the IPyC failure."}>>> = failure_indicator_IPyC
    OPyC_failure<<<{"description": "The name of the postprocessor calculating the OPyC failure."}>>> = failure_indicator_OPyC
    SiC_failure<<<{"description": "The name of the postprocessor calculating the SiC failure."}>>> = failure_indicator_SiC
    SiC_failure_crackedIPyC<<<{"description": "The name of the postprocessor calculating the SiC failure due to cracked IPyC."}>>> = failure_indicator_SiC_crackedIPyC
    SiC_failure_crackedOPyC<<<{"description": "The name of the postprocessor calculating the SiC failure due to cracked OPyC."}>>> = failure_indicator_SiC_crackedOPyC
    SiC_failure_pd_penetration<<<{"description": "The name of the postprocessor calculating the SiC failure due to palladium penetration."}>>> = failure_indicator_pd_penetration
    SiC_failure_kernel_migration<<<{"description": "The name of the postprocessor calculating the SiC failure due to kernel migration."}>>> = failure_indicator_kernel_migration
  []
[]

Input Parameters

kernel_migration_distanceThe name of the postprocessor that computes kernel migration distance.
C++ Type:PostprocessorName
Unit:(no unit assumed)
Controllable:No
Description:The name of the postprocessor that computes kernel migration distance.
triso_geometryUser object containing the TrisoGeometry
C++ Type:UserObjectName
Controllable:No
Description:User object containing the TrisoGeometry

Required Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:XFEM_MARK, NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, TRANSFER
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup

Execution Scheduling 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.
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
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
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

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

(examples/TRISO/failure_probability_monte_carlo/triso_1d_constant.i)
(examples/TRISO/failure_probability_monte_carlo/triso_1d_function.i)
(test/tests/triso_failure/triso_1d_kernel_migration.i)

(test/tests/triso_failure/triso_1d_kernel_migration.i)

kernel_radius = 213.35e-6
buffer_thickness = 98.9e-6
IPyC_thickness = 40.4e-6
SiC_thickness = 35.2e-6
OPyC_thickness = 43.4e-6

coordinates1 = '${fparse kernel_radius}'
coordinates2 = '${fparse coordinates1+buffer_thickness}'
coordinates3 = '${fparse coordinates2+IPyC_thickness}'
coordinates4 = '${fparse coordinates3+SiC_thickness}'
coordinates5 = '${fparse coordinates4+OPyC_thickness}'

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x'
  initial_enrichment = 0.14029 # [wt-]
  flux_conversion_factor = 1.0 # convert E>0.10 to E>0.18 MeV
  stress_free_temperature = 1573 # used for thermal expansion
  energy_per_fission = 3.204e-11 # [J/fission]
  O_U = 1.428 # Initial Oxygen to Uranium atom ratio
  C_U = 0.392 # Initial Carbon to Uranium atom ratio
[]

[Mesh]
  coord_type = RSPHERICAL
  [mesh]
    type = TRISO1DMeshGenerator
    elem_type = EDGE2
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} '
                  '${coordinates5}'
    mesh_density = '20 8 0 4 4 4'
    block_names = 'fuel buffer IPyC SiC OPyC'
  []
[]

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
[]

[UserObjects]
  [particle_geometry]
    type = TRISOGeometry
    outer_OPyC = OPyC_outer_boundary
    outer_SiC = SiC_outer_boundary
    outer_IPyC = IPyC_outer_boundary
    inner_IPyC = IPyC_inner_boundary
    outer_buffer = buffer_outer_boundary
    outer_kernel = fuel_outer_boundary
    include_particle = true
    include_pebble = false
    execute_on = 'INITIAL TIMESTEP_END'
  []
[]

[Variables]
  [temperature]
    initial_condition = 1573
  []
[]

[AuxVariables]
  [fission_rate]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_flux]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [fission_rate]
    type = ConstantFunction
    value = 5.75e19
  []
  [high_fidelity_strength_crackedIPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1000000 1000000'
  []
  [high_fidelity_strength_crackedOPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1000000 1000000'
  []
  [stress_correlation_crackedOPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1 1'
  []
  [stress_correlation_crackedIPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '100 100'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  generate_output = 'stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz '
                    'max_principal_stress'
  add_variables = true
  strain = FINITE
  incremental = true
  [fuel]
    block = fuel
    eigenstrain_names = 'UCO_swelling_eigenstrain UCO_TE_strain'
    extra_vector_tags = 'ref'
  []
  [buffer]
    block = buffer
    eigenstrain_names = 'Buffer_IIDC_strain Buffer_TE_strain'
    extra_vector_tags = 'ref'
  []
  [IPyC]
    block = IPyC
    eigenstrain_names = 'IPyC_IIDC_strain IPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
  [SiC]
    block = SiC
    eigenstrain_names = 'SiC_thermal_eigenstrain'
    extra_vector_tags = 'ref'
  []
  [OPyC]
    block = OPyC
    eigenstrain_names = 'OPyC_IIDC_strain OPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
[]

[Kernels]
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat]
    type = HeatConduction
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temperature
    block = fuel
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [fissionrate]
    type = MaterialRealAux
    variable = fission_rate
    property = fission_rate
    block = fuel
    execute_on = timestep_begin
  []
  [burnup]
    type = MaterialRealAux
    variable = burnup
    property = burnup
    block = fuel
    execute_on = timestep_begin
  []
  [fast_neutron_flux]
    type = MaterialRealAux
    variable = fast_neutron_flux
    property = fast_neutron_flux
    execute_on = timestep_begin
  []
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_begin
  []
[]

[ThermalContact]
  [thermal_contact]
    type = GasGapHeatTransfer
    variable = temperature
    primary = IPyC_inner_boundary
    secondary = buffer_outer_boundary
    initial_moles = initial_moles
    gas_released = 'fis_gas_released'
    released_gas_types = 'Kr Xe'
    released_fractions = '0.185 0.815'
    tangential_tolerance = 1e-6
    quadrature = false
    min_gap = 1e-7
    max_gap = 50e-6
    gap_geometry_type = sphere
  []
[]

[BCs]
  [no_disp_x]
    type = DirichletBC
    variable = disp_x
    boundary = xzero
    value = 0.0
  []
  [freesurf_temp]
    type = DirichletBC
    variable = temperature
    value = 1573
    boundary = exterior
  []
  [exterior_pressure_x]
    type = Pressure
    variable = disp_x
    boundary = exterior
    factor = 0.1e6
  []
  [PlenumPressure]
    [plenumPressure]
      boundary = buffer_IPyC_boundary
      startup_time = 1e4
      initial_pressure = 0
      output_initial_moles = initial_moles
      temperature = ave_gas_temp
      volume = 'gap_volume buffer_void_volume kernel_void_volume'
      material_input = 'fis_gas_released'
      output = gas_pressure
    []
  []
[]

[Materials]
  [fission_rate]
    type = GenericFunctionMaterial
    prop_names = fission_rate
    prop_values = fission_rate
    block = fuel
  []
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    factor = 6.2425e+17
  []
  [UCO_burnup]
    type = TRISOBurnup
    initial_density = 10966
    block = fuel
  []
  [UCO_thermal]
    type = UCOThermal
    block = fuel
    temperature = temperature
  []
  [UCO_elasticity_tensor]
    type = UCOElasticityTensor
    block = fuel
    temperature = temperature
  []
  [UCO_stress]
    type = ComputeFiniteStrainElasticStress
    block = fuel
  []
  [UCO_VolumetricSwellingEigenstrain]
    type = UCOVolumetricSwellingEigenstrain
    block = fuel
    eigenstrain_name = UCO_swelling_eigenstrain
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = fuel
    thermal_expansion_coeff = 10.0e-6
    temperature = temperature
    eigenstrain_name = UCO_TE_strain
  []
  [UCO_density]
    type = StrainAdjustedDensity
    block = fuel
    strain_free_density = 10966
  []
  [fission_gas_release]
    type = UCOFGR
    block = fuel
    average_grain_radius = 10e-6
    temperature = temperature
    triso_geometry = particle_geometry
    cutoff_neutron_flux = 0.0
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = 'buffer IPyC OPyC'
  []
  [BAF_IPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0465
    block = IPyC
  []
  [BAF_OPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0429
    block = OPyC
  []
  [buffer_elasticity_tensor]
    type = BufferElasticityTensor
    block = buffer
    temperature = temperature
  []
  [buffer_stress]
    type = BufferCEGACreep
    block = buffer
    temperature = temperature
  []
  [buffer_thermal]
    type = BufferThermal
    block = buffer
    initial_density = 1050.0
  []
  [buffer_density]
    type = StrainAdjustedDensity
    block = buffer
    strain_free_density = 1050.0
  []
  [buffer_TE]
    type = BufferThermalExpansionEigenstrain
    block = buffer
    eigenstrain_name = Buffer_TE_strain
    temperature = temperature
  []
  [buffer_IIDC]
    type = BufferCEGAIrradiationEigenstrain
    block = buffer
    eigenstrain_name = Buffer_IIDC_strain
    temperature = temperature
  []
  [IPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = IPyC
    temperature = temperature
  []
  [IPyC_stress]
    type = PyCCEGACreep
    block = IPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []

  [IPyC_thermal]
    type = HeatConductionMaterial
    block = IPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [IPyC_density]
    type = GenericConstantMaterial
    block = IPyC
    prop_names = 'density'
    prop_values = 1890
  []
  [IPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [IPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_TE_strain
    temperature = temperature
  []
  [SiC_elasticity_tensor]
    type = MonolithicSiCElasticityTensor
    block = SiC
    temperature = temperature
    elastic_modulus_model = miller
  []
  [SiC_stress]
    type = ComputeFiniteStrainElasticStress
    block = SiC
  []
  [SiC_thermal]
    type = MonolithicSiCThermal
    block = SiC
    temperature = temperature
    thermal_conductivity_model = miller
  []
  [SiC_density]
    type = StrainAdjustedDensity
    block = SiC
    strain_free_density = 3200.0
  []
  [SiC_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = SiC
    thermal_expansion_coeff = 4.9e-6
    temperature = temperature
    eigenstrain_name = SiC_thermal_eigenstrain
  []
  [OPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = OPyC
    temperature = temperature
    initial_BAF = 1.0
  []
  [OPyC_stress]
    type = PyCCEGACreep
    block = OPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []
  [OPyC_thermal_conductivity]
    type = HeatConductionMaterial
    block = OPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [OPyC_density]
    type = GenericConstantMaterial
    block = OPyC
    prop_names = 'density'
    prop_values = 1900
  []
  [OPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [OPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_TE_strain
    temperature = temperature
  []
  [characteristic_strength_SiC]
    type = GenericConstantMaterial
    prop_values = '9640000'
    block = SiC
    prop_names = 'characteristic_strength'
  []
  [characteristic_strength_PyC]
    type = PyCCharacteristicStrength
    temperature = temperature
    X = 1.02
    block = 'IPyC OPyC'
  []
[]

[Dampers]
  [temp]
    type = MaxIncrement
    variable = temperature
    max_increment = 100
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  line_search = 'none'

  nl_rel_tol = 5e-6
  nl_abs_tol = 1e-8
  nl_max_its = 20

  l_tol = 1e-4
  l_max_its = 50

  start_time = 0.0
  num_steps = 10

  dtmin = 1e-4
  dt = 1e4
[]

[Postprocessors]
  [ave_gas_temp]
    type = ElementAverageValue
    block = buffer
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fis_gas_released]
    type = ElementIntegralMaterialProperty
    mat_prop = fis_gas_released
    block = fuel
    use_displaced_mesh = false
    execute_on = 'initial timestep_end'
  []
  [gap_volume]
    type = InternalVolume
    boundary = buffer_IPyC_boundary
    execute_on = 'initial linear'
    use_displaced_mesh = true
  []
  [buffer_void_volume]
    type = VoidVolume
    block = buffer
    theoretical_density = 2250
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [kernel_th_density]
    type = UCOTheoreticalDensity
    execute_on = initial
  []
  [kernel_void_volume]
    type = VoidVolume
    block = fuel
    theoretical_density = kernel_th_density
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [particle_power]
    type = ElementIntegralPower
    variable = temperature
    use_material_fission_rate = true
    fission_rate_material = fission_rate
    block = fuel
    execute_on = 'initial timestep_end'
  []
  [max_fluence]
    type = ElementExtremeValue
    variable = fast_neutron_fluence
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [max_burnup]
    type = ElementExtremeValue
    variable = burnup
    block = fuel
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [SiC_stress_max]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = max
    mat_prop = stress_yy
  []
  [SiC_stress_min]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = min
    mat_prop = stress_yy
  []
  [strength_SiC]
    type = WeibullEffectiveMeanStrength
    block = SiC
    weibull_modulus = 6
  []
  [kernel_migration_distance]
    type = KernelMigrationDistance
    block = 'fuel buffer IPyC SiC OPyC'
    variable = temperature
    temperature_gradient = 15000
    kernel_type = UO2
  []
  [failure_indicator_kernel_migration]
    type = KernelMigrationFailureIndicator
    kernel_migration_distance = kernel_migration_distance
    triso_geometry = particle_geometry
  []
  [failure_indicator_SiC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_SiC
  []
  [strength_IPyC]
    type = WeibullEffectiveMeanStrength
    block = IPyC
    weibull_modulus = 6
  []
  [failure_indicator_IPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = IPyC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_IPyC
  []
  [strength_OPyC]
    type = WeibullEffectiveMeanStrength
    block = OPyC
    weibull_modulus = 6
  []
  [failure_indicator_OPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = OPyC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_OPyC
  []
  [pd_penetration]
    type = PdPenetration
    boundary = SiC_inner_boundary
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [failure_indicator_pd_penetration]
    type = PdPenetrationFailureIndicator
    triso_geometry = particle_geometry
    pd_penetration = pd_penetration
  []
  [failure_indicator_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [failure_indicator_SiC_crackedOPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedOPyC'
  []
  [triso_failure]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    OPyC_failure = failure_indicator_OPyC
    SiC_failure = failure_indicator_SiC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure_crackedOPyC = failure_indicator_SiC_crackedOPyC
    SiC_failure_pd_penetration = failure_indicator_pd_penetration
    SiC_failure_kernel_migration = failure_indicator_kernel_migration
  []
[]

[Outputs]
  show = 'kernel_migration_distance failure_indicator_kernel_migration triso_failure '
  print_linear_residuals = false
  time_step_interval =  1
  csv = true
  perf_graph = true
[]

(examples/TRISO/failure_probability_monte_carlo/triso_1d_constant.i)


initial_fuel_density = 5

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x'
  initial_enrichment = 0.14029 # [wt-]
  flux_conversion_factor = 1.0 # convert E>0.10 to E>0.18 MeV
  stress_free_temperature = 481 # used for thermal expansion
  energy_per_fission = 3.204e-11 # [J/fission]
  O_U = 1.428 # Initial Oxygen to Uranium atom ratio
  C_U = 0.392 # Initial Carbon to Uranium atom ratio
[]

[Mesh]
  coord_type = RSPHERICAL
  [gen]
    type = TRISO1DFiveLayerMeshGenerator
    elem_type = EDGE3
    kernel_radius = 213.35e-6
    buffer_thickness = 98.9e-6
    IPyC_thickness = 40.4e-6
    SiC_thickness = 35.2e-6
    OPyC_thickness = 43.4e-6
    kernel_mesh_density = ${initial_fuel_density}
    buffer_mesh_density = 3
    IPyC_mesh_density = 5
    SiC_mesh_density = 3
    OPyC_mesh_density = 4
  []
[]

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
[]

[UserObjects]
  [particle_geometry]
    type = TRISOGeometry
    outer_OPyC = OPyC_outer_boundary
    outer_SiC = SiC_outer_boundary
    outer_IPyC = IPyC_outer_boundary
    inner_IPyC = IPyC_inner_boundary
    outer_buffer = buffer_outer_boundary
    outer_kernel = fuel_outer_boundary
    include_particle = true
    include_pebble = false
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
    execute_on = 'INITIAL TIMESTEP_END'
  []
  [sic_failure_terminator]
    type = Terminator
    expression = 'sic_failure_overall > 0'
  []
[]

[Variables]
  [temperature]
    initial_condition = 481
  []
[]

[AuxVariables]
  [fission_rate]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_flux]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temp_bc]
    type = PiecewiseLinear
    data_file = outer_temp.csv
    x_index_in_file = 0
    y_index_in_file = 1
    format = columns
  []
  [fission_rate]
    type = ConstantFunction
    value = 5.75e19
  []
  [high_fidelity_strength_crackedIPyC]
    type = ConstantFunction
    value = '1403604095.5707'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 0 0'
    polynomial_coefficients_SiC = '1 0 0'
    polynomial_coefficients_OPyC = '1 0 0'
    correlation_factor = -1.2447543103484047
  []
  [high_fidelity_strength_debonding]
    type = ConstantFunction
    value = '1705800293.3578'
  []
  [stress_correlation_debonding]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 0 0'
    polynomial_coefficients_SiC = '1 0 0'
    polynomial_coefficients_OPyC = '1 0 0'
    correlation_factor = -0.14916368684964607
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1371700766.8875'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 0 0'
    polynomial_coefficients_SiC = '1 0 0'
    polynomial_coefficients_OPyC = '1 0 0'
    correlation_factor = 1.5191967987843993
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 0 0'
    polynomial_coefficients_SiC = '1 0 0'
    polynomial_coefficients_OPyC = '1 0 0'
    correlation_factor = 1.391516859626456
  []
  [sic_crackedipyc_stress_strength]
    type = ParsedFunction
    expression = 'a-b'
    symbol_names = 'a b'
    symbol_values = 'stress_SiC_crackedIPyC actual_strength_SiC_crackedIPyC'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  generate_output = 'stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz max_principal_stress'
  add_variables = true
  strain = FINITE
  incremental = true
  [fuel]
    block = fuel
    eigenstrain_names = 'UCO_swelling_eigenstrain UCO_TE_strain'
    extra_vector_tags = 'ref'
  []
  [buffer]
    block = buffer
    eigenstrain_names = 'Buffer_IIDC_strain Buffer_TE_strain'
    extra_vector_tags = 'ref'
  []
  [IPyC]
    block = IPyC
    eigenstrain_names = 'IPyC_IIDC_strain IPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
  [SiC]
    block = SiC
    eigenstrain_names = 'SiC_thermal_eigenstrain'
    extra_vector_tags = 'ref'
  []
  [OPyC]
    block = OPyC
    eigenstrain_names = 'OPyC_IIDC_strain OPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
[]

[Kernels]
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat]
    type = HeatConduction
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temperature
    block = fuel
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [fissionrate]
    type = MaterialRealAux
    variable = fission_rate
    property = fission_rate
    block = fuel
    execute_on = timestep_begin
  []
  [burnup]
    type = MaterialRealAux
    variable = burnup
    property = burnup
    block = fuel
    execute_on = timestep_begin
  []
  [fast_neutron_flux]
    type = MaterialRealAux
    variable = fast_neutron_flux
    property = fast_neutron_flux
    execute_on = timestep_begin
  []
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_begin
  []
[]

[ThermalContact]
  [thermal_contact]
    type = GasGapHeatTransfer
    variable = temperature
    primary = IPyC_inner_boundary
    secondary = buffer_outer_boundary
    initial_moles = initial_moles
    gas_released = 'fis_gas_released'
    released_gas_types = 'Kr Xe'
    released_fractions = '0.185 0.815'
    tangential_tolerance = 1e-6
    quadrature = false
    min_gap = 1e-7
    max_gap = 50e-6
    gap_geometry_type = sphere
  []
[]

[BCs]
  [no_disp_x]
    type = DirichletBC
    variable = disp_x
    boundary = xzero
    value = 0.0
  []
  [freesurf_temp]
    type = FunctionDirichletBC
    variable = temperature
    function = temp_bc
    boundary = exterior
  []
  [exterior_pressure_x]
    type = Pressure
    variable = disp_x
    boundary = exterior
    factor = 0.1e6
  []
  [PlenumPressure]
    [plenumPressure]
      boundary = buffer_IPyC_boundary
      startup_time = 1e4
      initial_pressure = 0
      R = 8.3145
      output_initial_moles = initial_moles
      temperature = ave_gas_temp
      volume = 'gap_volume buffer_void_volume kernel_void_volume'
      material_input = 'fis_gas_released'
      output = gas_pressure
    []
  []
[]

[Materials]
  [radial_stress]
    type = RankTwoCylindricalComponent
    rank_two_tensor = stress
    cylindrical_axis_point1 = '0 0 0'
    cylindrical_axis_point2 = '0 0 1'
    cylindrical_component = RadialStress
    property_name = radial_stress
    outputs = all
  []
  [fission_rate]
    type = GenericFunctionMaterial
    prop_names = fission_rate
    prop_values = fission_rate
    block = fuel
  []
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    factor = 6.2425e+17
  []
  [UCO_burnup]
    type = TRISOBurnup
    initial_density = 10966
    block = fuel
  []
  [UCO_thermal]
    type = UCOThermal
    block = fuel
    temperature = temperature
  []
  [UCO_elasticity_tensor]
    type = UCOElasticityTensor
    block = fuel
    temperature = temperature
  []
  [UCO_stress]
    type = ComputeFiniteStrainElasticStress
    block = fuel
  []
  [UCO_VolumetricSwellingEigenstrain]
    type = UCOVolumetricSwellingEigenstrain
    block = fuel
    eigenstrain_name = UCO_swelling_eigenstrain
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = fuel
    thermal_expansion_coeff = 10.0e-6
    temperature = temperature
    eigenstrain_name = UCO_TE_strain
  []
  [UCO_density]
    type = StrainAdjustedDensity
    block = fuel
    strain_free_density = 10966
  []
  [fission_gas_release]
    type = UCOFGR
    block = fuel
    average_grain_radius = 10e-6
    temperature = temperature
    triso_geometry = particle_geometry
    cutoff_neutron_flux = 0.0
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = 'buffer IPyC OPyC'
  []
  [BAF_IPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0465
    block = IPyC
  []
  [BAF_OPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0429
    block = OPyC
  []
  [buffer_elasticity_tensor]
    type = BufferElasticityTensor
    block = buffer
    temperature = temperature
  []
  [buffer_stress]
    type = BufferCEGACreep
    block = buffer
    temperature = temperature
  []
  [buffer_thermal]
    type = BufferThermal
    block = buffer
    initial_density = 1050.0
  []
  [buffer_density]
    type = StrainAdjustedDensity
    block = buffer
    strain_free_density = 1050.0
  []
  [buffer_TE]
    type = BufferThermalExpansionEigenstrain
    block = buffer
    eigenstrain_name = Buffer_TE_strain
    temperature = temperature
  []
  [buffer_IIDC]
    type = BufferCEGAIrradiationEigenstrain
    block = buffer
    eigenstrain_name = Buffer_IIDC_strain
    temperature = temperature
  []
  [IPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = IPyC
    temperature = temperature
  []
  [IPyC_stress]
    type = PyCCEGACreep
    block = IPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []

  [IPyC_thermal]
    type = HeatConductionMaterial
    block = IPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [IPyC_density]
    type = GenericConstantMaterial
    block = IPyC
    prop_names = 'density'
    prop_values = 1890
  []
  [IPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [IPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_TE_strain
    temperature = temperature
  []
  [SiC_elasticity_tensor]
    type = MonolithicSiCElasticityTensor
    block = SiC
    temperature = temperature
    elastic_modulus_model = miller
  []
  [SiC_stress]
    type = ComputeFiniteStrainElasticStress
    block = SiC
  []
  [SiC_thermal]
    type = MonolithicSiCThermal
    block = SiC
    temperature = temperature
    thermal_conductivity_model = miller
  []
  [SiC_density]
    type = StrainAdjustedDensity
    block = SiC
    strain_free_density = 3200.0
  []
  [SiC_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = SiC
    thermal_expansion_coeff = 4.9e-6
    temperature = temperature
    eigenstrain_name = SiC_thermal_eigenstrain
  []
  [OPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = OPyC
    temperature = temperature
    initial_BAF = 1.0
  []
  [OPyC_stress]
    type = PyCCEGACreep
    block = OPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []
  [OPyC_thermal_conductivity]
    type = HeatConductionMaterial
    block = OPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [OPyC_density]
    type = GenericConstantMaterial
    block = OPyC
    prop_names = 'density'
    prop_values = 1900
  []
  [OPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [OPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_TE_strain
    temperature = temperature
  []
  [characteristic_strength_SiC]
    type = GenericConstantMaterial
    prop_values = '9640000'
    block = SiC
    prop_names = 'characteristic_strength'
  []
  [characteristic_strength_PyC]
    type = PyCCharacteristicStrength
    temperature = temperature
    X = 1.02
    block = 'IPyC OPyC'
  []
[]

[Dampers]
  [temp]
    type = MaxIncrement
    variable = temperature
    max_increment = 100
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  line_search = 'none'

  nl_rel_tol = 5e-6
  nl_abs_tol = 1e-8
  nl_max_its = 20

  l_tol = 1e-4
  l_max_its = 50

  start_time = 0.0
  end_time = 4.831315e7

  dtmin = 1e-4
  dt = 5e5
[]

[Postprocessors]
  [ave_gas_temp]
    type = ElementAverageValue
    block = buffer
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fis_gas_released]
    type = ElementIntegralMaterialProperty
    mat_prop = fis_gas_released
    block = fuel
    use_displaced_mesh = false
    execute_on = 'initial timestep_end'
  []
  [gap_volume]
    type = InternalVolume
    boundary = buffer_IPyC_boundary
    execute_on = 'initial linear'
    use_displaced_mesh = true
  []
  [buffer_void_volume]
    type = VoidVolume
    block = buffer
    theoretical_density = 2250
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [kernel_th_density]
    type = UCOTheoreticalDensity
    execute_on = initial
  []
  [kernel_void_volume]
    type = VoidVolume
    block = fuel
    theoretical_density = kernel_th_density
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [particle_power]
    type = ElementIntegralPower
    variable = temperature
    use_material_fission_rate = true
    fission_rate_material = fission_rate
    block = fuel
    execute_on = 'initial timestep_end'
  []
  [max_fluence]
    type = ElementExtremeValue
    variable = fast_neutron_fluence
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [max_burnup]
    type = ElementExtremeValue
    variable = burnup
    block = fuel
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [SiC_stress]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = min
    mat_prop = stress_yy
  []
  [strength_SiC]
    type = WeibullEffectiveMeanStrength
    block = SiC
    weibull_modulus = 6
  []
  [failure_indicator_SiC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_asphericity'
    stress_correlation_function = 'stress_correlation_asphericity'
    stress_change_correlation_function = 'stress_change_correlation_asphericity'
  []
  [strength_IPyC]
    type = WeibullEffectiveMeanStrength
    block = IPyC
    weibull_modulus = 9.5
  []
  [failure_indicator_IPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = IPyC
    weibull_modulus = 9.5
    stress_name = max_principal_stress
    effective_mean_strength = strength_IPyC
  []
  [failure_indicator_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [failure_indicator_debonding]
    type = TRISODebondingFailureIndicator
    boundary = IPyC_outer_boundary
    bond_strength = 10e6
    stress_name = radial_stress
  []
  [failure_indicator_SiC_debonding]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_debonding'
    stress_correlation_function = 'stress_correlation_debonding'
  []
  [sic_failure_overall]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    SiC_failure_pd_penetration = failure_indicator_pd_penetration
    SiC_failure_kernel_migration = failure_indicator_kernel_migration
    failure_type = SIC_FAILURE_OVERALL
  []
  [ipyc_cracking]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = IPYC_CRACKING
  []
  [sic_failure_due_to_pressure]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = SIC_FAILURE_DUE_TO_PRESSURE
  []
  [sic_failure_due_to_ipyc_cracking]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = SIC_FAILURE_DUE_TO_IPYC_CRACKING
  []
  [stress_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
    output_type = 'stress'
  []
  [actual_strength_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
    output_type = 'strength'
  []
  [SiC_crackedIPyC_stressminusstrength]
    type = FunctionValuePostprocessor
    function = 'sic_crackedipyc_stress_strength'
  []
  [debonding]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    IPyC_SiC_debonding = failure_indicator_debonding
    SiC_failure_debonding = failure_indicator_SiC_debonding
    failure_type = IPYC_SIC_DEBONDING
  []
  [fluence_at_failure]
    type = TRISOFailureOccurrenceStatus
    failure_evaluation = ipyc_cracking
    failure_information = max_fluence
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbability
    block = SiC
    weibull_modulus = 6
    characteristic_strength = characteristic_strength
  []
  [pd_penetration]
    type = PdPenetration
    boundary = SiC_inner_boundary
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [failure_indicator_pd_penetration]
    type = PdPenetrationFailureIndicator
    triso_geometry = particle_geometry
    pd_penetration = pd_penetration
  []
  [kernel_migration_distance]
    type = KernelMigrationDistance
    block = 'fuel buffer IPyC SiC OPyC'
    variable = temperature
    temperature_gradient = 15000
    kernel_type = UCO
  []
  [failure_indicator_kernel_migration]
    type = KernelMigrationFailureIndicator
    kernel_migration_distance = kernel_migration_distance
    triso_geometry = particle_geometry
  []
[]

[Outputs]
  print_linear_residuals = false
  time_step_interval =  1
  csv = false
  exodus = false
  perf_graph = true
  print_linear_converged_reason = false
  print_nonlinear_converged_reason = false
[]

(examples/TRISO/failure_probability_monte_carlo/triso_1d_function.i)


initial_fuel_density = 5

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x'
  initial_enrichment = 0.14029 # [wt-]
  flux_conversion_factor = 1.0 # convert E>0.10 to E>0.18 MeV
  stress_free_temperature = 481 # used for thermal expansion
  energy_per_fission = 3.204e-11 # [J/fission]
  O_U = 1.428 # Initial Oxygen to Uranium atom ratio
  C_U = 0.392 # Initial Carbon to Uranium atom ratio
[]

[Mesh]
  coord_type = RSPHERICAL
  [gen]
    type = TRISO1DFiveLayerMeshGenerator
    elem_type = EDGE3
    kernel_radius = 213.35e-6
    buffer_thickness = 98.9e-6
    IPyC_thickness = 40.4e-6
    SiC_thickness = 35.2e-6
    OPyC_thickness = 43.4e-6
    kernel_mesh_density = ${initial_fuel_density}
    buffer_mesh_density = 3
    IPyC_mesh_density = 5
    SiC_mesh_density = 3
    OPyC_mesh_density = 4
  []
[]

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
[]

[UserObjects]
  [particle_geometry]
    type = TRISOGeometry
    outer_OPyC = OPyC_outer_boundary
    outer_SiC = SiC_outer_boundary
    outer_IPyC = IPyC_outer_boundary
    inner_IPyC = IPyC_inner_boundary
    outer_buffer = buffer_outer_boundary
    outer_kernel = fuel_outer_boundary
    include_particle = true
    include_pebble = false
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
    execute_on = 'INITIAL TIMESTEP_END'
  []
  [sic_failure_terminator]
    type = Terminator
    expression = 'sic_failure_overall > 0'
  []
[]

[Variables]
  [temperature]
    initial_condition = 481
  []
[]

[AuxVariables]
  [fission_rate]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_flux]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temp_bc]
    type = PiecewiseLinear
    data_file = outer_temp.csv
    x_index_in_file = 0
    y_index_in_file = 1
    format = columns
  []
  [fission_rate]
    type = ConstantFunction
    value = 5.75e19
  []
  [high_fidelity_strength_crackedIPyC]
    type = ConstantFunction
    value = '1403604095.0794'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 5.95176524e3 -2.25337303e8'
    polynomial_coefficients_SiC = '1 1.43220859e4 -5.17689523e7'
    polynomial_coefficients_OPyC = '1 -1.25870267e4 1.81620484e8'
    correlation_factor = -1.2447543093270736
  []
  [high_fidelity_strength_debonding]
    type = ConstantFunction
    value = '1705800293.3578'
  []
  [stress_correlation_debonding]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 0 0'
    polynomial_coefficients_SiC = '1 0 0'
    polynomial_coefficients_OPyC = '1 0 0'
    correlation_factor = -0.14916368684964607
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1371700806.9481'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 1.00595402e3 1.43530004e7'
    polynomial_coefficients_SiC = '1 3.27925856e3 -2.02308753e8'
    polynomial_coefficients_OPyC = '1 2.07404580e3 -6.12612615e6'
    correlation_factor = 1.5191967993808713
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -5.81891553e3 -2.81628655e7'
    polynomial_coefficients_SiC = '1 1.00990700e4 -5.55290343e8'
    polynomial_coefficients_OPyC = '1 -3.59151050e3 -2.65952373e7'
    correlation_factor = 1.3915168526633837
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  generate_output = 'stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz max_principal_stress'
  add_variables = true
  strain = FINITE
  incremental = true
  [fuel]
    block = fuel
    eigenstrain_names = 'UCO_swelling_eigenstrain UCO_TE_strain'
    extra_vector_tags = 'ref'
  []
  [buffer]
    block = buffer
    eigenstrain_names = 'Buffer_IIDC_strain Buffer_TE_strain'
    extra_vector_tags = 'ref'
  []
  [IPyC]
    block = IPyC
    eigenstrain_names = 'IPyC_IIDC_strain IPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
  [SiC]
    block = SiC
    eigenstrain_names = 'SiC_thermal_eigenstrain'
    extra_vector_tags = 'ref'
  []
  [OPyC]
    block = OPyC
    eigenstrain_names = 'OPyC_IIDC_strain OPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
[]

[Kernels]
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat]
    type = HeatConduction
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temperature
    block = fuel
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [fissionrate]
    type = MaterialRealAux
    variable = fission_rate
    property = fission_rate
    block = fuel
    execute_on = timestep_begin
  []
  [burnup]
    type = MaterialRealAux
    variable = burnup
    property = burnup
    block = fuel
    execute_on = timestep_begin
  []
  [fast_neutron_flux]
    type = MaterialRealAux
    variable = fast_neutron_flux
    property = fast_neutron_flux
    execute_on = timestep_begin
  []
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_begin
  []
[]

[ThermalContact]
  [thermal_contact]
    type = GasGapHeatTransfer
    variable = temperature
    primary = IPyC_inner_boundary
    secondary = buffer_outer_boundary
    initial_moles = initial_moles
    gas_released = 'fis_gas_released'
    released_gas_types = 'Kr Xe'
    released_fractions = '0.185 0.815'
    tangential_tolerance = 1e-6
    quadrature = false
    min_gap = 1e-7
    max_gap = 50e-6
    gap_geometry_type = sphere
  []
[]

[BCs]
  [no_disp_x]
    type = DirichletBC
    variable = disp_x
    boundary = xzero
    value = 0.0
  []
  [freesurf_temp]
    type = FunctionDirichletBC
    variable = temperature
    function = temp_bc
    boundary = exterior
  []
  [exterior_pressure_x]
    type = Pressure
    variable = disp_x
    boundary = exterior
    factor = 0.1e6
  []
  [PlenumPressure]
    [plenumPressure]
      boundary = buffer_IPyC_boundary
      startup_time = 1e4
      initial_pressure = 0
      R = 8.3145
      output_initial_moles = initial_moles
      temperature = ave_gas_temp
      volume = 'gap_volume buffer_void_volume kernel_void_volume'
      material_input = 'fis_gas_released'
      output = gas_pressure
    []
  []
[]

[Materials]
  [radial_stress]
    type = RankTwoCylindricalComponent
    rank_two_tensor = stress
    cylindrical_axis_point1 = '0 0 0'
    cylindrical_axis_point2 = '0 0 1'
    cylindrical_component = RadialStress
    property_name = radial_stress
    outputs = all
  []
  [fission_rate]
    type = GenericFunctionMaterial
    prop_names = fission_rate
    prop_values = fission_rate
    block = fuel
  []
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    factor = 6.2425e+17
  []
  [UCO_burnup]
    type = TRISOBurnup
    initial_density = 10966
    block = fuel
  []
  [UCO_thermal]
    type = UCOThermal
    block = fuel
    temperature = temperature
  []
  [UCO_elasticity_tensor]
    type = UCOElasticityTensor
    block = fuel
    temperature = temperature
  []
  [UCO_stress]
    type = ComputeFiniteStrainElasticStress
    block = fuel
  []
  [UCO_VolumetricSwellingEigenstrain]
    type = UCOVolumetricSwellingEigenstrain
    block = fuel
    eigenstrain_name = UCO_swelling_eigenstrain
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = fuel
    thermal_expansion_coeff = 10.0e-6
    temperature = temperature
    eigenstrain_name = UCO_TE_strain
  []
  [UCO_density]
    type = StrainAdjustedDensity
    block = fuel
    strain_free_density = 10966
  []
  [fission_gas_release]
    type = UCOFGR
    block = fuel
    average_grain_radius = 10e-6
    temperature = temperature
    triso_geometry = particle_geometry
    cutoff_neutron_flux = 0.0
  []
  [BAF_IPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0465
    block = IPyC
  []
  [BAF_OPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0429
    block = OPyC
  []
  [buffer_elasticity_tensor]
    type = BufferElasticityTensor
    block = buffer
    temperature = temperature
  []
  [buffer_stress]
    type = BufferCEGACreep
    block = buffer
    temperature = temperature
  []
  [buffer_thermal]
    type = BufferThermal
    block = buffer
    initial_density = 1050.0
  []
  [buffer_density]
    type = StrainAdjustedDensity
    block = buffer
    strain_free_density = 1050.0
  []
  [buffer_TE]
    type = BufferThermalExpansionEigenstrain
    block = buffer
    eigenstrain_name = Buffer_TE_strain
    temperature = temperature
  []
  [buffer_IIDC]
    type = BufferCEGAIrradiationEigenstrain
    block = buffer
    eigenstrain_name = Buffer_IIDC_strain
    temperature = temperature
  []
  [IPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = IPyC
    temperature = temperature
  []
  [IPyC_stress]
    type = PyCCEGACreep
    block = IPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []

  [IPyC_thermal]
    type = HeatConductionMaterial
    block = IPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [IPyC_density]
    type = GenericConstantMaterial
    block = IPyC
    prop_names = 'density'
    prop_values = 1890
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = 'buffer IPyC OPyC'
  []
  [IPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [IPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_TE_strain
    temperature = temperature
  []
  [SiC_elasticity_tensor]
    type = MonolithicSiCElasticityTensor
    block = SiC
    temperature = temperature
    elastic_modulus_model = miller
  []
  [SiC_stress]
    type = ComputeFiniteStrainElasticStress
    block = SiC
  []
  [SiC_thermal]
    type = MonolithicSiCThermal
    block = SiC
    temperature = temperature
    thermal_conductivity_model = miller
  []
  [SiC_density]
    type = StrainAdjustedDensity
    block = SiC
    strain_free_density = 3200.0
  []
  [SiC_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = SiC
    thermal_expansion_coeff = 4.9e-6
    temperature = temperature
    eigenstrain_name = SiC_thermal_eigenstrain
  []
  [OPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = OPyC
    temperature = temperature
    initial_BAF = 1.0
  []
  [OPyC_stress]
    type = PyCCEGACreep
    block = OPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []
  [OPyC_thermal_conductivity]
    type = HeatConductionMaterial
    block = OPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [OPyC_density]
    type = GenericConstantMaterial
    block = OPyC
    prop_names = 'density'
    prop_values = 1900
  []
  [OPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [OPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_TE_strain
    temperature = temperature
  []
  [characteristic_strength_SiC]
    type = GenericConstantMaterial
    prop_values = '9640000'
    block = SiC
    prop_names = 'characteristic_strength'
  []
  [characteristic_strength_PyC]
    type = PyCCharacteristicStrength
    temperature = temperature
    X = 1.02
    block = 'IPyC OPyC'
  []
[]

[Dampers]
  [temp]
    type = MaxIncrement
    variable = temperature
    max_increment = 100
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  line_search = 'none'

  nl_rel_tol = 5e-6
  nl_abs_tol = 1e-8
  nl_max_its = 20

  l_tol = 1e-4
  l_max_its = 50

  start_time = 0.0
  end_time = 4.831315e7

  dtmin = 1e-4
  dt = 5e5
[]

[Postprocessors]
  [ave_gas_temp]
    type = ElementAverageValue
    block = buffer
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fis_gas_released]
    type = ElementIntegralMaterialProperty
    mat_prop = fis_gas_released
    block = fuel
    use_displaced_mesh = false
    execute_on = 'initial timestep_end'
  []
  [gap_volume]
    type = InternalVolume
    boundary = buffer_IPyC_boundary
    execute_on = 'initial linear'
    use_displaced_mesh = true
  []
  [buffer_void_volume]
    type = VoidVolume
    block = buffer
    theoretical_density = 2250
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [kernel_th_density]
    type = UCOTheoreticalDensity
    execute_on = initial
  []
  [kernel_void_volume]
    type = VoidVolume
    block = fuel
    theoretical_density = kernel_th_density
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [particle_power]
    type = ElementIntegralPower
    variable = temperature
    use_material_fission_rate = true
    fission_rate_material = fission_rate
    block = fuel
    execute_on = 'initial timestep_end'
  []
  [max_fluence]
    type = ElementExtremeValue
    variable = fast_neutron_fluence
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [max_burnup]
    type = ElementExtremeValue
    variable = burnup
    block = fuel
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [SiC_stress]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = min
    mat_prop = stress_yy
  []
  [strength_SiC]
    type = WeibullEffectiveMeanStrength
    block = SiC
    weibull_modulus = 6
  []
  [failure_indicator_SiC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_asphericity'
    stress_correlation_function = 'stress_correlation_asphericity'
    stress_change_correlation_function = 'stress_change_correlation_asphericity'
  []
  [strength_IPyC]
    type = WeibullEffectiveMeanStrength
    block = IPyC
    weibull_modulus = 9.5
  []
  [failure_indicator_IPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = IPyC
    weibull_modulus = 9.5
    stress_name = max_principal_stress
    effective_mean_strength = strength_IPyC
  []
  [failure_indicator_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [failure_indicator_debonding]
    type = TRISODebondingFailureIndicator
    boundary = IPyC_outer_boundary
    bond_strength = 10e6
    stress_name = radial_stress
  []
  [failure_indicator_SiC_debonding]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_debonding'
    stress_correlation_function = 'stress_correlation_debonding'
  []
  [sic_failure_overall]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    SiC_failure_pd_penetration = failure_indicator_pd_penetration
    SiC_failure_kernel_migration = failure_indicator_kernel_migration
    failure_type = SIC_FAILURE_OVERALL
  []
  [ipyc_cracking]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = IPYC_CRACKING
  []
  [sic_failure_due_to_pressure]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = SIC_FAILURE_DUE_TO_PRESSURE
  []
  [sic_failure_due_to_ipyc_cracking]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    failure_type = SIC_FAILURE_DUE_TO_IPYC_CRACKING
  []
  [debonding]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure = failure_indicator_SiC
    IPyC_SiC_debonding = failure_indicator_debonding
    SiC_failure_debonding = failure_indicator_SiC_debonding
    failure_type = IPYC_SIC_DEBONDING
  []
  [fluence_at_failure]
    type = TRISOFailureOccurrenceStatus
    failure_evaluation = ipyc_cracking
    failure_information = max_fluence
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbability
    block = SiC
    weibull_modulus = 6
    characteristic_strength = characteristic_strength
  []
  [pd_penetration]
    type = PdPenetration
    boundary = SiC_inner_boundary
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [failure_indicator_pd_penetration]
    type = PdPenetrationFailureIndicator
    triso_geometry = particle_geometry
    pd_penetration = pd_penetration
  []
  [kernel_migration_distance]
    type = KernelMigrationDistance
    block = 'fuel buffer IPyC SiC OPyC'
    variable = temperature
    temperature_gradient = 15000
    kernel_type = UCO
  []
  [failure_indicator_kernel_migration]
    type = KernelMigrationFailureIndicator
    kernel_migration_distance = kernel_migration_distance
    triso_geometry = particle_geometry
  []
[]

[Outputs]
  print_linear_residuals = false
  time_step_interval =  1
  csv = false
  exodus = false
  perf_graph = true
  print_linear_converged_reason = false
  print_nonlinear_converged_reason = false
[]

(test/tests/triso_failure/triso_1d_kernel_migration.i)

kernel_radius = 213.35e-6
buffer_thickness = 98.9e-6
IPyC_thickness = 40.4e-6
SiC_thickness = 35.2e-6
OPyC_thickness = 43.4e-6

coordinates1 = '${fparse kernel_radius}'
coordinates2 = '${fparse coordinates1+buffer_thickness}'
coordinates3 = '${fparse coordinates2+IPyC_thickness}'
coordinates4 = '${fparse coordinates3+SiC_thickness}'
coordinates5 = '${fparse coordinates4+OPyC_thickness}'

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x'
  initial_enrichment = 0.14029 # [wt-]
  flux_conversion_factor = 1.0 # convert E>0.10 to E>0.18 MeV
  stress_free_temperature = 1573 # used for thermal expansion
  energy_per_fission = 3.204e-11 # [J/fission]
  O_U = 1.428 # Initial Oxygen to Uranium atom ratio
  C_U = 0.392 # Initial Carbon to Uranium atom ratio
[]

[Mesh]
  coord_type = RSPHERICAL
  [mesh]
    type = TRISO1DMeshGenerator
    elem_type = EDGE2
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} '
                  '${coordinates5}'
    mesh_density = '20 8 0 4 4 4'
    block_names = 'fuel buffer IPyC SiC OPyC'
  []
[]

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
[]

[UserObjects]
  [particle_geometry]
    type = TRISOGeometry
    outer_OPyC = OPyC_outer_boundary
    outer_SiC = SiC_outer_boundary
    outer_IPyC = IPyC_outer_boundary
    inner_IPyC = IPyC_inner_boundary
    outer_buffer = buffer_outer_boundary
    outer_kernel = fuel_outer_boundary
    include_particle = true
    include_pebble = false
    execute_on = 'INITIAL TIMESTEP_END'
  []
[]

[Variables]
  [temperature]
    initial_condition = 1573
  []
[]

[AuxVariables]
  [fission_rate]
    order = CONSTANT
    family = MONOMIAL
  []
  [burnup]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_flux]
    order = CONSTANT
    family = MONOMIAL
  []
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [fission_rate]
    type = ConstantFunction
    value = 5.75e19
  []
  [high_fidelity_strength_crackedIPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1000000 1000000'
  []
  [high_fidelity_strength_crackedOPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1000000 1000000'
  []
  [stress_correlation_crackedOPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '1 1'
  []
  [stress_correlation_crackedIPyC]
    type = PiecewiseLinear
    x = '0 1.0e11'
    y = '100 100'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  generate_output = 'stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz '
                    'max_principal_stress'
  add_variables = true
  strain = FINITE
  incremental = true
  [fuel]
    block = fuel
    eigenstrain_names = 'UCO_swelling_eigenstrain UCO_TE_strain'
    extra_vector_tags = 'ref'
  []
  [buffer]
    block = buffer
    eigenstrain_names = 'Buffer_IIDC_strain Buffer_TE_strain'
    extra_vector_tags = 'ref'
  []
  [IPyC]
    block = IPyC
    eigenstrain_names = 'IPyC_IIDC_strain IPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
  [SiC]
    block = SiC
    eigenstrain_names = 'SiC_thermal_eigenstrain'
    extra_vector_tags = 'ref'
  []
  [OPyC]
    block = OPyC
    eigenstrain_names = 'OPyC_IIDC_strain OPyC_TE_strain'
    extra_vector_tags = 'ref'
  []
[]

[Kernels]
  [heat_ie]
    type = HeatConductionTimeDerivative
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat]
    type = HeatConduction
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = NeutronHeatSource
    variable = temperature
    block = fuel
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [fissionrate]
    type = MaterialRealAux
    variable = fission_rate
    property = fission_rate
    block = fuel
    execute_on = timestep_begin
  []
  [burnup]
    type = MaterialRealAux
    variable = burnup
    property = burnup
    block = fuel
    execute_on = timestep_begin
  []
  [fast_neutron_flux]
    type = MaterialRealAux
    variable = fast_neutron_flux
    property = fast_neutron_flux
    execute_on = timestep_begin
  []
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_begin
  []
[]

[ThermalContact]
  [thermal_contact]
    type = GasGapHeatTransfer
    variable = temperature
    primary = IPyC_inner_boundary
    secondary = buffer_outer_boundary
    initial_moles = initial_moles
    gas_released = 'fis_gas_released'
    released_gas_types = 'Kr Xe'
    released_fractions = '0.185 0.815'
    tangential_tolerance = 1e-6
    quadrature = false
    min_gap = 1e-7
    max_gap = 50e-6
    gap_geometry_type = sphere
  []
[]

[BCs]
  [no_disp_x]
    type = DirichletBC
    variable = disp_x
    boundary = xzero
    value = 0.0
  []
  [freesurf_temp]
    type = DirichletBC
    variable = temperature
    value = 1573
    boundary = exterior
  []
  [exterior_pressure_x]
    type = Pressure
    variable = disp_x
    boundary = exterior
    factor = 0.1e6
  []
  [PlenumPressure]
    [plenumPressure]
      boundary = buffer_IPyC_boundary
      startup_time = 1e4
      initial_pressure = 0
      output_initial_moles = initial_moles
      temperature = ave_gas_temp
      volume = 'gap_volume buffer_void_volume kernel_void_volume'
      material_input = 'fis_gas_released'
      output = gas_pressure
    []
  []
[]

[Materials]
  [fission_rate]
    type = GenericFunctionMaterial
    prop_names = fission_rate
    prop_values = fission_rate
    block = fuel
  []
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    factor = 6.2425e+17
  []
  [UCO_burnup]
    type = TRISOBurnup
    initial_density = 10966
    block = fuel
  []
  [UCO_thermal]
    type = UCOThermal
    block = fuel
    temperature = temperature
  []
  [UCO_elasticity_tensor]
    type = UCOElasticityTensor
    block = fuel
    temperature = temperature
  []
  [UCO_stress]
    type = ComputeFiniteStrainElasticStress
    block = fuel
  []
  [UCO_VolumetricSwellingEigenstrain]
    type = UCOVolumetricSwellingEigenstrain
    block = fuel
    eigenstrain_name = UCO_swelling_eigenstrain
  []
  [fuel_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = fuel
    thermal_expansion_coeff = 10.0e-6
    temperature = temperature
    eigenstrain_name = UCO_TE_strain
  []
  [UCO_density]
    type = StrainAdjustedDensity
    block = fuel
    strain_free_density = 10966
  []
  [fission_gas_release]
    type = UCOFGR
    block = fuel
    average_grain_radius = 10e-6
    temperature = temperature
    triso_geometry = particle_geometry
    cutoff_neutron_flux = 0.0
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = 'buffer IPyC OPyC'
  []
  [BAF_IPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0465
    block = IPyC
  []
  [BAF_OPyC]
    type = BaconAnisotropyFactor
    initial_BAF = 1.0429
    block = OPyC
  []
  [buffer_elasticity_tensor]
    type = BufferElasticityTensor
    block = buffer
    temperature = temperature
  []
  [buffer_stress]
    type = BufferCEGACreep
    block = buffer
    temperature = temperature
  []
  [buffer_thermal]
    type = BufferThermal
    block = buffer
    initial_density = 1050.0
  []
  [buffer_density]
    type = StrainAdjustedDensity
    block = buffer
    strain_free_density = 1050.0
  []
  [buffer_TE]
    type = BufferThermalExpansionEigenstrain
    block = buffer
    eigenstrain_name = Buffer_TE_strain
    temperature = temperature
  []
  [buffer_IIDC]
    type = BufferCEGAIrradiationEigenstrain
    block = buffer
    eigenstrain_name = Buffer_IIDC_strain
    temperature = temperature
  []
  [IPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = IPyC
    temperature = temperature
  []
  [IPyC_stress]
    type = PyCCEGACreep
    block = IPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []

  [IPyC_thermal]
    type = HeatConductionMaterial
    block = IPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [IPyC_density]
    type = GenericConstantMaterial
    block = IPyC
    prop_names = 'density'
    prop_values = 1890
  []
  [IPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [IPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_TE_strain
    temperature = temperature
  []
  [SiC_elasticity_tensor]
    type = MonolithicSiCElasticityTensor
    block = SiC
    temperature = temperature
    elastic_modulus_model = miller
  []
  [SiC_stress]
    type = ComputeFiniteStrainElasticStress
    block = SiC
  []
  [SiC_thermal]
    type = MonolithicSiCThermal
    block = SiC
    temperature = temperature
    thermal_conductivity_model = miller
  []
  [SiC_density]
    type = StrainAdjustedDensity
    block = SiC
    strain_free_density = 3200.0
  []
  [SiC_thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = SiC
    thermal_expansion_coeff = 4.9e-6
    temperature = temperature
    eigenstrain_name = SiC_thermal_eigenstrain
  []
  [OPyC_elasticity_tensor]
    type = PyCElasticityTensor
    block = OPyC
    temperature = temperature
    initial_BAF = 1.0
  []
  [OPyC_stress]
    type = PyCCEGACreep
    block = OPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []
  [OPyC_thermal_conductivity]
    type = HeatConductionMaterial
    block = OPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [OPyC_density]
    type = GenericConstantMaterial
    block = OPyC
    prop_names = 'density'
    prop_values = 1900
  []
  [OPyC_IIDC]
    type = PyCCEGAIrradiationEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [OPyC_TE]
    type = PyCThermalExpansionEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_TE_strain
    temperature = temperature
  []
  [characteristic_strength_SiC]
    type = GenericConstantMaterial
    prop_values = '9640000'
    block = SiC
    prop_names = 'characteristic_strength'
  []
  [characteristic_strength_PyC]
    type = PyCCharacteristicStrength
    temperature = temperature
    X = 1.02
    block = 'IPyC OPyC'
  []
[]

[Dampers]
  [temp]
    type = MaxIncrement
    variable = temperature
    max_increment = 100
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  line_search = 'none'

  nl_rel_tol = 5e-6
  nl_abs_tol = 1e-8
  nl_max_its = 20

  l_tol = 1e-4
  l_max_its = 50

  start_time = 0.0
  num_steps = 10

  dtmin = 1e-4
  dt = 1e4
[]

[Postprocessors]
  [ave_gas_temp]
    type = ElementAverageValue
    block = buffer
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fis_gas_released]
    type = ElementIntegralMaterialProperty
    mat_prop = fis_gas_released
    block = fuel
    use_displaced_mesh = false
    execute_on = 'initial timestep_end'
  []
  [gap_volume]
    type = InternalVolume
    boundary = buffer_IPyC_boundary
    execute_on = 'initial linear'
    use_displaced_mesh = true
  []
  [buffer_void_volume]
    type = VoidVolume
    block = buffer
    theoretical_density = 2250
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [kernel_th_density]
    type = UCOTheoreticalDensity
    execute_on = initial
  []
  [kernel_void_volume]
    type = VoidVolume
    block = fuel
    theoretical_density = kernel_th_density
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [particle_power]
    type = ElementIntegralPower
    variable = temperature
    use_material_fission_rate = true
    fission_rate_material = fission_rate
    block = fuel
    execute_on = 'initial timestep_end'
  []
  [max_fluence]
    type = ElementExtremeValue
    variable = fast_neutron_fluence
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [max_burnup]
    type = ElementExtremeValue
    variable = burnup
    block = fuel
    value_type = 'max'
    execute_on = 'initial timestep_end'
  []
  [SiC_stress_max]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = max
    mat_prop = stress_yy
  []
  [SiC_stress_min]
    type = ElementExtremeMaterialProperty
    block = SiC
    value_type = min
    mat_prop = stress_yy
  []
  [strength_SiC]
    type = WeibullEffectiveMeanStrength
    block = SiC
    weibull_modulus = 6
  []
  [kernel_migration_distance]
    type = KernelMigrationDistance
    block = 'fuel buffer IPyC SiC OPyC'
    variable = temperature
    temperature_gradient = 15000
    kernel_type = UO2
  []
  [failure_indicator_kernel_migration]
    type = KernelMigrationFailureIndicator
    kernel_migration_distance = kernel_migration_distance
    triso_geometry = particle_geometry
  []
  [failure_indicator_SiC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_SiC
  []
  [strength_IPyC]
    type = WeibullEffectiveMeanStrength
    block = IPyC
    weibull_modulus = 6
  []
  [failure_indicator_IPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = IPyC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_IPyC
  []
  [strength_OPyC]
    type = WeibullEffectiveMeanStrength
    block = OPyC
    weibull_modulus = 6
  []
  [failure_indicator_OPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = OPyC
    weibull_modulus = 6
    stress_name = max_principal_stress
    effective_mean_strength = strength_OPyC
  []
  [pd_penetration]
    type = PdPenetration
    boundary = SiC_inner_boundary
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [failure_indicator_pd_penetration]
    type = PdPenetrationFailureIndicator
    triso_geometry = particle_geometry
    pd_penetration = pd_penetration
  []
  [failure_indicator_SiC_crackedIPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [failure_indicator_SiC_crackedOPyC]
    type = WeibullFailureOutputUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = max_principal_stress
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedOPyC'
  []
  [triso_failure]
    type = TRISOFailureEvaluation
    IPyC_failure = failure_indicator_IPyC
    OPyC_failure = failure_indicator_OPyC
    SiC_failure = failure_indicator_SiC
    SiC_failure_crackedIPyC = failure_indicator_SiC_crackedIPyC
    SiC_failure_crackedOPyC = failure_indicator_SiC_crackedOPyC
    SiC_failure_pd_penetration = failure_indicator_pd_penetration
    SiC_failure_kernel_migration = failure_indicator_kernel_migration
  []
[]

[Outputs]
  show = 'kernel_migration_distance failure_indicator_kernel_migration triso_failure '
  print_linear_residuals = false
  time_step_interval =  1
  csv = true
  perf_graph = true
[]

Description
Example Input Syntax
Input Parameters
Input Files