WeibullFailureProbabilityUsingCorrelation

Computes failure probability of a TRISO layer using a correlation function.

Description

WeibullFailureProbabilityUsingCorrelation can be used to compute Weibull failure probability using a a stress correlation function obtained from a higher-fidelity model. This Postprocessor computes a Weibull failure probability value that is between 0 and 1. If the new value computed from current time step exceeds the old value from previous time step, the new value will be returned. Otherwise, the old value will be returned. In general, the failure probability mode will be used in direct integration TRISOFailureProbabilityIntegration approach.

The Weibull failure probability is computed as: $var element = document.getElementById("moose-equation-323ca72f-0611-4872-8b8c-4c716f4ecc12");katex.render("P_f = 1-\\exp{\\left(-\\left[ \\frac{\\sigma_c}{\\sigma_{ms}}\\right]^m\\right)},", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-c5c7caab-f756-4bab-a0d8-48a8f6367dd6");katex.render("\\sigma_c", element, {displayMode:false,throwOnError:false});$ is the maximum stress value, $var element = document.getElementById("moose-equation-599a7ad3-0049-4e49-8125-caa3faa8bf52");katex.render("m", element, {displayMode:false,throwOnError:false});$ is the Weibull modulus and $var element = document.getElementById("moose-equation-26af43db-3e3f-47b1-b513-ccf98d8a4d35");katex.render("\\sigma_{ms}", element, {displayMode:false,throwOnError:false});$ is the effective mean strength WeibullEffectiveMeanStrength. To account for high dimensional failure behaviors, the correlation function of $var element = document.getElementById("moose-equation-a77c9c14-bf88-4b5b-8a2b-2e0c8ec07098");katex.render("\\sigma_c", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-a153e345-b9c6-4fdf-9bac-8696d2e76827");katex.render("\\sigma_{ms}", element, {displayMode:false,throwOnError:false});$ will be obtained from a higher-fidelity model, see WeibullFailureOutputUsingCorrelation for the details.

This postprocessor can be used in conjunction with AD variables through using ADWeibullFailureProbabilityUsingCorrelation.

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]
    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
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability<<<{"description": "Computes Weibull failure probability using principle of independent action (PIA) model.", "href": "WeibullFailureProbability.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = IPyC
    weibull_modulus<<<{"description": ",Weibull modulus"}>>> = 9.5
    characteristic_strength<<<{"description": "The material name of characteristic strength."}>>> = characteristic_strength
  []
  [weibull_failure_probability_SiC_crackedIPyC]
    type = WeibullFailureProbabilityUsingCorrelation<<<{"description": "Computes failure probability of a TRISO layer using a correlation function.", "href": "WeibullFailureProbabilityUsingCorrelation.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."}>>> = stress_yy
    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'
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbabilityUsingCorrelation<<<{"description": "Computes failure probability of a TRISO layer using a correlation function.", "href": "WeibullFailureProbabilityUsingCorrelation.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."}>>> = stress_yy
    high_fidelity_analysis_strength<<<{"description": "The effective mean strength obtained from a high-fidelity analysis."}>>> = 'high_fidelity_strength_asphericity'
    stress_correlation_function<<<{"description": "The stress correlation function obtained from a high-fidelity analysis."}>>> = 'stress_correlation_asphericity'
    stress_change_correlation_function<<<{"description": "Correlation function of stress change from minumum to maximum during the irradiation history obtained from a high-fidelity analysis."}>>> = 'stress_change_correlation_asphericity'
  []
[]

Input Parameters

stress_nameName of the stress property.
C++ Type:std::string
Controllable:No
Description:Name of the stress property.
weibull_modulus,Weibull modulus
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:,Weibull modulus

Required Parameters

blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
effective_mean_strengthThe name of the postprocessor calculating the effective mean strength.
C++ Type:PostprocessorName
Unit:(no unit assumed)
Controllable:No
Description:The name of the postprocessor calculating the effective mean strength.
high_fidelity_analysis_strengthThe effective mean strength obtained from a high-fidelity analysis.
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:The effective mean strength obtained from a high-fidelity analysis.
stress_change_correlation_functionCorrelation function of stress change from minumum to maximum during the irradiation history obtained from a high-fidelity analysis.
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Correlation function of stress change from minumum to maximum during the irradiation history obtained from a high-fidelity analysis.
stress_correlation_functionThe stress correlation function obtained from a high-fidelity analysis.
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:The stress correlation function obtained from a high-fidelity analysis.

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

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/triso_failure/triso_1d_ipyc_weibull_probability.i)
(examples/TRISO/failure_probability_direct_integration/triso_1d.i)
(test/tests/triso_failure/ad_triso_1d_ipyc_weibull_probability.i)

(test/tests/triso_failure/triso_1d_ipyc_weibull_probability.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 = 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 = TRISO1DMeshGenerator
    elem_type = EDGE3
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} ${coordinates5}'
    mesh_density = '5 3 0 5 3 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
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
  []
[]

[Variables]
  [temperature]
    initial_condition = 900
  []
[]

[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 = ConstantFunction
    value = '1363350801.3058'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 7017 -2.368e8'
    polynomial_coefficients_SiC = '1 1.492e4 -3.802e7'
    polynomial_coefficients_OPyC = '1 -1.273e4 1.849e8'
    correlation_factor = -1.1824630660785265
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1086690814.283'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -2070 3.458e7'
    polynomial_coefficients_SiC = '1 -868.9 -1.368e7'
    polynomial_coefficients_OPyC = '1 1734 -1.988e7'
    correlation_factor = 1.0626986695756293
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -856 1.593e7'
    polynomial_coefficients_SiC = '1 1774 -5.253e7'
    polynomial_coefficients_OPyC = '1 456.4 -1.459e7'
    correlation_factor = 1.0113764663823708
  []
[]

[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 = 900
    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
  []
  [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
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC_crackedIPyC]
    type = WeibullFailureProbabilityUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbabilityUsingCorrelation
    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'
  []
[]

[Outputs]
  show = 'weibull_failure_probability_IPyC weibull_failure_probability_SiC weibull_failure_probability_SiC_crackedIPyC'
  print_linear_residuals = false
  time_step_interval =  1
  csv = true
  exodus = false
  perf_graph = false
  print_linear_converged_reason = false
  print_nonlinear_converged_reason = false
[]

(test/tests/triso_failure/triso_1d_ipyc_weibull_probability.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 = 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 = TRISO1DMeshGenerator
    elem_type = EDGE3
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} ${coordinates5}'
    mesh_density = '5 3 0 5 3 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
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
  []
[]

[Variables]
  [temperature]
    initial_condition = 900
  []
[]

[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 = ConstantFunction
    value = '1363350801.3058'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 7017 -2.368e8'
    polynomial_coefficients_SiC = '1 1.492e4 -3.802e7'
    polynomial_coefficients_OPyC = '1 -1.273e4 1.849e8'
    correlation_factor = -1.1824630660785265
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1086690814.283'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -2070 3.458e7'
    polynomial_coefficients_SiC = '1 -868.9 -1.368e7'
    polynomial_coefficients_OPyC = '1 1734 -1.988e7'
    correlation_factor = 1.0626986695756293
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -856 1.593e7'
    polynomial_coefficients_SiC = '1 1774 -5.253e7'
    polynomial_coefficients_OPyC = '1 456.4 -1.459e7'
    correlation_factor = 1.0113764663823708
  []
[]

[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 = 900
    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
  []
  [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
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC_crackedIPyC]
    type = WeibullFailureProbabilityUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbabilityUsingCorrelation
    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'
  []
[]

[Outputs]
  show = 'weibull_failure_probability_IPyC weibull_failure_probability_SiC weibull_failure_probability_SiC_crackedIPyC'
  print_linear_residuals = false
  time_step_interval =  1
  csv = true
  exodus = false
  perf_graph = false
  print_linear_converged_reason = false
  print_nonlinear_converged_reason = false
[]

(examples/TRISO/failure_probability_direct_integration/triso_1d.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}'

initial_fuel_density = 10966

[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 = TRISO1DMeshGenerator
    elem_type = EDGE3
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} ${coordinates5}'
    mesh_density = '5 3 0 5 3 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
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
  []
[]

[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 = '1363350801.3058'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 7017 -2.368e8'
    polynomial_coefficients_SiC = '1 1.492e4 -3.802e7'
    polynomial_coefficients_OPyC = '1 -1.273e4 1.849e8'
    correlation_factor = -1.1824630660785265
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1086690814.283'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -2070 3.458e7'
    polynomial_coefficients_SiC = '1 -868.9 -1.368e7'
    polynomial_coefficients_OPyC = '1 1734 -1.988e7'
    correlation_factor = 1.0626986695756293
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -856 1.593e7'
    polynomial_coefficients_SiC = '1 1774 -5.253e7'
    polynomial_coefficients_OPyC = '1 456.4 -1.459e7'
    correlation_factor = 1.0113764663823708
  []
[]

[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]
  [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 = ${initial_fuel_density}
    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 = ${initial_fuel_density}
  []
  [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

  #num_steps = 1

  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
  []
  [weibull_failure_probability_IPyC]
    type = WeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC_crackedIPyC]
    type = WeibullFailureProbabilityUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [weibull_failure_probability_SiC]
    type = WeibullFailureProbabilityUsingCorrelation
    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'
  []
[]

[Outputs]
  show = 'weibull_failure_probability_IPyC weibull_failure_probability_SiC weibull_failure_probability_SiC_crackedIPyC'
  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/ad_triso_1d_ipyc_weibull_probability.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 = 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 = TRISO1DMeshGenerator
    elem_type = EDGE3
    coordinates = '0 ${coordinates1} ${coordinates2} ${coordinates2} ${coordinates3} ${coordinates4} ${coordinates5}'
    mesh_density = '5 3 0 5 3 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
    IPyC_thickness_mean = 40.4e-6
    SiC_thickness_mean = 35.2e-6
    OPyC_thickness_mean = 43.4e-6
  []
[]

[Variables]
  [temperature]
    initial_condition = 900
  []
[]

[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 = ConstantFunction
    value = '1363350801.3058'
  []
  [stress_correlation_crackedIPyC]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 7017 -2.368e8'
    polynomial_coefficients_SiC = '1 1.492e4 -3.802e7'
    polynomial_coefficients_OPyC = '1 -1.273e4 1.849e8'
    correlation_factor = -1.1824630660785265
  []
  [high_fidelity_strength_asphericity]
    type = ConstantFunction
    value = '1086690814.283'
  []
  [stress_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -2070 3.458e7'
    polynomial_coefficients_SiC = '1 -868.9 -1.368e7'
    polynomial_coefficients_OPyC = '1 1734 -1.988e7'
    correlation_factor = 1.0626986695756293
  []
  [stress_change_correlation_asphericity]
    type = TRISOStressCorrelationFunction
    triso_geometry = particle_geometry
    polynomial_coefficients_IPyC = '1 -856 1.593e7'
    polynomial_coefficients_SiC = '1 1774 -5.253e7'
    polynomial_coefficients_OPyC = '1 456.4 -1.459e7'
    correlation_factor = 1.0113764663823708
  []
[]

[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'
    use_automatic_differentiation = true
  []
  [buffer]
    block = buffer
    eigenstrain_names = 'Buffer_IIDC_strain Buffer_TE_strain'
    extra_vector_tags = 'ref'
    use_automatic_differentiation = true
  []
  [IPyC]
    block = IPyC
    eigenstrain_names = 'IPyC_IIDC_strain IPyC_TE_strain'
    extra_vector_tags = 'ref'
    use_automatic_differentiation = true
  []
  [SiC]
    block = SiC
    eigenstrain_names = 'SiC_thermal_eigenstrain'
    extra_vector_tags = 'ref'
    use_automatic_differentiation = true
  []
  [OPyC]
    block = OPyC
    eigenstrain_names = 'OPyC_IIDC_strain OPyC_TE_strain'
    extra_vector_tags = 'ref'
    use_automatic_differentiation = true
  []
[]

[Kernels]
  [heat_ie]
    type = ADHeatConductionTimeDerivative
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat]
    type = ADHeatConduction
    variable = temperature
    extra_vector_tags = 'ref'
  []
  [heat_source]
    type = ADNeutronHeatSource
    variable = temperature
    block = fuel
    fission_rate = fission_rate
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [fissionrate]
    type = ADMaterialRealAux
    variable = fission_rate
    property = fission_rate
    block = fuel
    execute_on = timestep_begin
  []
  [burnup]
    type = ADMaterialRealAux
    variable = burnup
    property = burnup
    block = fuel
    execute_on = timestep_begin
  []
  [fast_neutron_flux]
    type = ADMaterialRealAux
    variable = fast_neutron_flux
    property = fast_neutron_flux
    execute_on = timestep_begin
  []
  [fast_neutron_fluence]
    type = ADMaterialRealAux
    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
    use_automatic_differentiation = true
  []
[]

[BCs]
  [no_disp_x]
    type = ADDirichletBC
    variable = disp_x
    boundary = xzero
    value = 0.0
  []
  [freesurf_temp]
    type = ADDirichletBC
    variable = temperature
    value = 900
    boundary = exterior
  []
  [exterior_pressure_x]
    type = ADPressure
    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
      use_automatic_differentiation = true
    []
  []
[]

[Materials]
  [fission_rate]
    type = ADGenericFunctionMaterial
    prop_names = fission_rate
    prop_values = fission_rate
    block = fuel
  []
  [fast_neutron_flux]
    type = ADFastNeutronFlux
    calculate_fluence = true
    factor = 6.2425e+17
  []
  [UCO_burnup]
    type = ADTRISOBurnup
    initial_density = 10966
    block = fuel
  []
  [UCO_thermal]
    type = ADUCOThermal
    block = fuel
    temperature = temperature
  []
  [UCO_elasticity_tensor]
    type = ADUCOElasticityTensor
    block = fuel
    temperature = temperature
  []
  [UCO_stress]
    type = ADComputeFiniteStrainElasticStress
    block = fuel
  []
  [UCO_VolumetricSwellingEigenstrain]
    type = ADUCOVolumetricSwellingEigenstrain
    block = fuel
    eigenstrain_name = UCO_swelling_eigenstrain
  []
  [fuel_thermal_expansion]
    type = ADComputeThermalExpansionEigenstrain
    block = fuel
    thermal_expansion_coeff = 10.0e-6
    temperature = temperature
    eigenstrain_name = UCO_TE_strain
  []
  [UCO_density]
    type = ADStrainAdjustedDensity
    block = fuel
    strain_free_density = 10966
  []
  [fission_gas_release]
    type = ADUCOFGR
    block = fuel
    average_grain_radius = 10e-6
    temperature = temperature
    triso_geometry = particle_geometry
    cutoff_neutron_flux = 0.0
  []
  [BAF_IPyC]
    type = ADBaconAnisotropyFactor
    initial_BAF = 1.0465
    block = IPyC
  []
  [BAF_OPyC]
    type = ADBaconAnisotropyFactor
    initial_BAF = 1.0429
    block = OPyC
  []
  [buffer_elasticity_tensor]
    type = ADBufferElasticityTensor
    block = buffer
    temperature = temperature
  []
  [buffer_stress]
    type = ADBufferCEGACreep
    block = buffer
    temperature = temperature
  []
  [buffer_thermal]
    type = ADBufferThermal
    block = buffer
    initial_density = 1050.0
  []
  [buffer_density]
    type = ADStrainAdjustedDensity
    block = buffer
    strain_free_density = 1050.0
  []
  [buffer_TE]
    type = ADBufferThermalExpansionEigenstrain
    block = buffer
    eigenstrain_name = Buffer_TE_strain
    temperature = temperature
  []
  [buffer_IIDC]
    type = ADBufferCEGAIrradiationEigenstrain
    block = buffer
    eigenstrain_name = Buffer_IIDC_strain
    temperature = temperature
  []
  [IPyC_elasticity_tensor]
    type = ADPyCElasticityTensor
    block = IPyC
    temperature = temperature
  []
  [IPyC_stress]
    type = ADPyCCEGACreep
    block = IPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []

  [IPyC_thermal]
    type = ADHeatConductionMaterial
    block = IPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [IPyC_density]
    type = ADGenericConstantMaterial
    block = IPyC
    prop_names = 'density'
    prop_values = 1890
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = 'buffer IPyC OPyC'
  []
  [IPyC_IIDC]
    type = ADPyCCEGAIrradiationEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [IPyC_TE]
    type = ADPyCThermalExpansionEigenstrain
    block = IPyC
    eigenstrain_name = IPyC_TE_strain
    temperature = temperature
  []
  [SiC_elasticity_tensor]
    type = ADMonolithicSiCElasticityTensor
    block = SiC
    temperature = temperature
    elastic_modulus_model = miller
  []
  [SiC_stress]
    type = ADComputeFiniteStrainElasticStress
    block = SiC
  []
  [SiC_thermal]
    type = ADMonolithicSiCThermal
    block = SiC
    temperature = temperature
    thermal_conductivity_model = miller
  []
  [SiC_density]
    type = ADStrainAdjustedDensity
    block = SiC
    strain_free_density = 3200.0
  []
  [SiC_thermal_expansion]
    type = ADComputeThermalExpansionEigenstrain
    block = SiC
    thermal_expansion_coeff = 4.9e-6
    temperature = temperature
    eigenstrain_name = SiC_thermal_eigenstrain
  []
  [OPyC_elasticity_tensor]
    type = ADPyCElasticityTensor
    block = OPyC
    temperature = temperature
    initial_BAF = 1.0
  []
  [OPyC_stress]
    type = ADPyCCEGACreep
    block = OPyC
    creep_rate_scale_factor = 1
    temperature = temperature
  []
  [OPyC_thermal_conductivity]
    type = ADHeatConductionMaterial
    block = OPyC
    thermal_conductivity = 4.0
    specific_heat = 720.0
  []
  [OPyC_density]
    type = ADGenericConstantMaterial
    block = OPyC
    prop_names = 'density'
    prop_values = 1900
  []
  [OPyC_IIDC]
    type = ADPyCCEGAIrradiationEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_IIDC_strain
    temperature = temperature
    irradiation_eigenstrain_scale_factor = 1
  []
  [OPyC_TE]
    type = ADPyCThermalExpansionEigenstrain
    block = OPyC
    eigenstrain_name = OPyC_TE_strain
    temperature = temperature
  []
  [characteristic_strength_SiC]
    type = ADGenericConstantMaterial
    prop_values = '9640000'
    block = SiC
    prop_names = 'characteristic_strength'
  []
  [characteristic_strength_PyC]
    type = ADPyCCharacteristicStrength
    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 = ADElementIntegralMaterialProperty
    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 = ADVoidVolume
    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 = ADVoidVolume
    block = fuel
    theoretical_density = kernel_th_density
    execute_on = 'initial timestep_end'
    use_displaced_mesh = true
  []
  [particle_power]
    type = ADElementIntegralPower
    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 = ADElementExtremeMaterialProperty
    block = SiC
    value_type = min
    mat_prop = stress_yy
  []
  [weibull_failure_probability_IPyC]
    type = ADWeibullFailureProbability
    block = IPyC
    weibull_modulus = 9.5
    characteristic_strength = characteristic_strength
  []
  [weibull_failure_probability_SiC_crackedIPyC]
    type = ADWeibullFailureProbabilityUsingCorrelation
    block = SiC
    weibull_modulus = 6
    stress_name = stress_yy
    high_fidelity_analysis_strength = 'high_fidelity_strength_crackedIPyC'
    stress_correlation_function = 'stress_correlation_crackedIPyC'
  []
  [weibull_failure_probability_SiC]
    type = ADWeibullFailureProbabilityUsingCorrelation
    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'
  []
[]

[Outputs]
  show = 'weibull_failure_probability_IPyC weibull_failure_probability_SiC weibull_failure_probability_SiC_crackedIPyC'
  print_linear_residuals = false
  time_step_interval =  1
  csv = true
  exodus = false
  perf_graph = false
  print_linear_converged_reason = false
  print_nonlinear_converged_reason = false
[]

Description
Example Input Syntax
Input Parameters
Input Files