Reconstructed Discontinuous Galerkin Failure Analysis Report

Introduction

The FAR for the Reconstructed Discontinuous Galerkin module describes the system functional and non-functional requirements that describe the expected interactions that the software shall provide.

Dependencies

The Reconstructed Discontinuous Galerkin module is developed using MOOSE and can itself be based on various MOOSE modules, as such the FAR for the Reconstructed Discontinuous Galerkin module is dependent upon the following documents.

Framework Failure Analysis Report

Failure Analysis Requirements

The following is a complete list for all the requirements related to failure analysis for the Reconstructed Discontinuous Galerkin module.

rdg: Actioncomponents
8.1.1
The system shall return an error if
1. the number of species does not match the number of species initial pressures,
2. the number of species does not match the number of scaling factors for the species equations,
3. the number of species does not match the number of equilibrium constants.
Specification(s): enclosure_0D_wrong_sizes/species_initial_p, enclosure_0D_wrong_sizes/species_scaling_factors, enclosure_0D_wrong_sizes/species_equilibrium
Design: Enclosure0D
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5
8.1.2
The system shall return an error if
1. no Physics have been defined on an enclosure,
2. more than one Physics have been defined on an enclosure, which is currently not supported,
3. an unsupported Physics has been defined on an enclosure,
Specification(s): enclosure_0D_physics_support/no_physics, enclosure_0D_physics_support/two_physics, enclosure_0D_physics_support/unsupported_physics
Design: Enclosure0D
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5
8.1.3
The system shall return an error if the parameters of the enclosure component and on the sorption exchange Physics are
1. both specified and not matching,
2. both not specified but necessary to define the equations.
Specification(s): enclosure_0D_parameter_consistency/not_matching, enclosure_0D_parameter_consistency/missing_in_both
Design: Enclosure0D Species Solubility in 0D Structures / SorptionExchangePhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5
8.1.4
The system shall return an error if
1. the number of species does not match the number of species initial pressures,
Specification(s): structure_1D_wrong_sizes/species_initial_condition
Design: Structure1D
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5
8.1.5
The system shall return an error if the parameters of the structure component and on the species trapping Physics active on it are
1. both specified and not matching for the species,
2. both specified and not matching,
3. both specified and not matching as a material property,
4. both not specified but necessary to define the equations.
Specification(s): structure_1D_parameter_consistency/not_matching_species, structure_1D_parameter_consistency/not_matching_params, structure_1D_parameter_consistency/not_matching_property, structure_1D_parameter_consistency/missing_in_both
Design: Structure1D Species Trapping Physics using a Continuous Galerkin Finite Element discretization / SpeciesTrappingPhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5 5
8.1.6
The system shall return an error if the material properties specified on the component and on the species trapping Physics active on it are
1. both specified and not matching,
2. both not specified but necessary to define the equations.
Specification(s): structure_1D_matprop_consistency/not_matching_material_prop, structure_1D_matprop_consistency/missing_material_prop
Design: Structure1D Species Trapping Physics using a Continuous Galerkin Finite Element discretization / SpeciesTrappingPhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5

rdg: Bcs
8.2.2
The system shall error if
1. the user does not specify a function or a variable for the equilibrium flux boundary condition,
2. the user does not specify the block to use for the nodal evaluation of equilibrium constant.
Specification(s): errors/missing_temp, errors/missing_block
Design: EquilibriumBC
Issue(s): #134
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5

rdg: Physics
8.8.1
The system shall return an error if
1. the number of species does not match the number of species initial pressures,
2. the number of species does not match the number of sorption equilibrium constants.
Specification(s): sorption_wrong_sizes/species_initial_p, sorption_wrong_sizes/species_equilibrium
Design: Species Solubility in 0D Structures / SorptionExchangePhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5
8.8.2
The system shall return an error if
1. the temperature specified is negative,
2. the temperature is specified as a variable which is not currently specified,
3. the species equation scaling factor is negative or null,
4. the species initial pressure is negative,
5. the sorption equilibrium constant is negative or null,
Specification(s): sorption_bad_data/temperature, sorption_bad_data/temperature_var, sorption_bad_data/scaling_factor, sorption_bad_data/initial_pressure, sorption_bad_data/eq_constant
Design: Species Solubility in 0D Structures / SorptionExchangePhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5 5 5
8.8.3
The system shall return an error if
1. a species is repeated twice in the list of species,
2. the number of species does not match the number of species initial concentrations,
3. the number of species does not match the number of species balance equation scaling factors.
4. the number of trapped species does not match the number of mobile species,
5. the number of species does not match the number of trapping coefficients,
6. the number of species does not match the number of trapping Ct0 coefficients,
7. the number of species does not match the number of release coefficients,
8. the number of species does not match the number of detrapping energies.
Specification(s): trapping_wrong_sizes/repeated_species, trapping_wrong_sizes/species_initial_concentration, trapping_wrong_sizes/scaling_factors, trapping_wrong_sizes/mobile, trapping_wrong_sizes/alpha_t, trapping_wrong_sizes/Ct0, trapping_wrong_sizes/alpha_r, trapping_wrong_sizes/detrapping
Design: Species Trapping Physics using a Continuous Galerkin Finite Element discretization / SpeciesTrappingPhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5 5 5 5 5 5
8.8.4
The system shall return an error if
1. the temperature is specified as a function when only constants and variables are supported at this time.
2. no components are specified yet the Physics is set to create variables for each component,
Specification(s): trapping_bad_data/temperature_fun, trapping_bad_data/no_component_yet_asking_for_vars_per_comp
Design: Species Trapping Physics using a Continuous Galerkin Finite Element discretization / SpeciesTrappingPhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5
8.8.5
The system shall return an error if the species trapping Physics is not used on a component and
1. the species is not specified,
2. the trapping coefficient is not specified,
3. the trapping energy is not specified,
4. the number density of the traps is not specified,
5. the trapping Ct0 coefficient is not specified,
6. the number of traps filled per free particle trapped is not specified,
7. the release coefficient is not specified,
8. the detrapping energy is not specified.
Specification(s): trapping_missing_data/species, trapping_missing_data/alpha_t, trapping_missing_data/trapping_energy, trapping_missing_data/N, trapping_missing_data/Ct0, trapping_missing_data/trap_per_free, trapping_missing_data/alpha_r, trapping_missing_data/detrapping_energy
Design: Species Trapping Physics using a Continuous Galerkin Finite Element discretization / SpeciesTrappingPhysics
Issue(s): #8
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5 5 5 5 5 5 5 5

rdg: Yttrium Hydrogen System
8.44.8The system shall be able to return a warning when the pressure and temperature are outside the range of validity of the YHxPCT model (pressure too low).
Specification(s): YHx_PCT_error_low_pressure
Design: InterfaceDiffusion ADMatInterfaceReactionYHxPCT MatDiffusion TimeDerivative
Issue(s): #261
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5
8.44.9The system shall be able to return a warning when the pressure and temperature are outside the range of validity of the YHxPCT model (pressure too high).
Specification(s): YHx_PCT_error_high_pressure
Design: InterfaceDiffusion ADMatInterfaceReactionYHxPCT MatDiffusion TimeDerivative
Issue(s): #261
Collection(s): FAILURE_ANALYSISFUNCTIONAL
Type(s): RunException
5


[species_initial_p]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/species_initial_pressures='1 2'"
  detail = "the number of species does not match the number of species initial pressures,"
  expect_err = "The number of species partial pressures must match the number of species."
[]


[species_scaling_factors]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/species_scaling_factors='1 2'"
  detail = "the number of species does not match the number of scaling factors for the species equations,"
  expect_err = "The number of species scaling factors must match the number of species."
[]


[species_equilibrium]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/equilibrium_constants='1 2'"
  detail = "the number of species does not match the number of equilibrium constants."
  expect_err = "The number of equilibrium constants must match the number of species"
[]


[no_physics]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/physics=''"
  detail = "no Physics have been defined on an enclosure,"
  expect_err = "A physics must be specified in the enclosure"
[]


[two_physics]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/physics='0d_solub multi-D'"
  detail = "more than one Physics have been defined on an enclosure, which is currently not supported,"
  expect_err = "Enclosure0D has only been implemented for a single 'SorptionExchangePhysics'"
[]


[unsupported_physics]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/physics='multi-D'"
  detail = "an unsupported Physics has been defined on an enclosure,"
  expect_err = "Physics 'multi-D' is not a 'SorptionExchangePhysics'. This component has only been implemented for 'SorptionExchangePhysics'."
[]


[not_matching]
  type = RunException
  input = enclosure0D.i
  cli_args = "ActionComponents/enc/equilibrium_constants=1 Physics/SorptionExchange/0d_solub/equilibrium_constants=2"
  detail = "both specified and not matching,"
  expect_err = " differs from 'equilibrium_constants' in SorptionExchangePhysics:"
[]


[missing_in_both]
  type = RunException
  input = enclosure0D.i
  detail = "both not specified but necessary to define the equations."
  expect_err = "This parameter should be specified, in the Physics '0d_solub' \(applying to all components\) or in component 'enc'"
[]


[species_initial_condition]
  type = RunException
  input = structure1D.i
  cli_args = "ActionComponents/structure/species_initial_concentrations='1 2'"
  detail = "the number of species does not match the number of species initial pressures,"
  expect_err = "The number of species concentrations must match the number of species."
[]


[not_matching_species]
  type = RunException
  input = structure1D.i
  cli_args = "ActionComponents/structure/species='a b c' ActionComponents/structure/species_initial_concentrations='1 1 1'"
  detail = "both specified and not matching for the species,"
  expect_err = "'species' in component 'structure' :"
[]


[not_matching_params]
  type = RunException
  input = structure1D.i
  cli_args = "ActionComponents/structure/temperature=100 Physics/SpeciesTrapping/trapped/temperature=200"
  detail = "both specified and not matching,"
  expect_err = " differs from 'temperature' in SpeciesTrappingPhysics:"
[]


[not_matching_property]
  type = RunException
  input = structure1D.i
  cli_args = "ActionComponents/structure/temperature=100 Physics/SpeciesTrapping/trapped/alpha_t=100"
  detail = "both specified and not matching as a material property,"
  expect_err = "'alpha_t' in component 'structure' :"
[]


[missing_in_both]
  type = RunException
  input = structure1D.i
  detail = "both not specified but necessary to define the equations."
  expect_err = "This parameter should be specified, in the Physics 'trapped' \(applying to all components\) or in component 'structure'"
[]


[not_matching_material_prop]
  type = RunException
  input = structure1D.i
  cli_args = "Physics/SpeciesTrapping/trapped/alpha_t=1e12 ActionComponents/structure/temperature=100"
  detail = "both specified and not matching,"
  expect_err = " differs from 'alpha_t' in SpeciesTrappingPhysics:"
[]


[missing_material_prop]
  type = RunException
  input = structure1D.i
  cli_args = "ActionComponents/structure/property_names='a b c d e f g diff' ActionComponents/structure/temperature=100"
  detail = "both not specified but necessary to define the equations."
  expect_err = "alpha_t: This parameter should be specified, in the Physics 'trapped' \(applying to all components\) or in component 'structure' using the property_name 'alpha_t'"
[]


[missing_temp]
  type = RunException
  input = 'equilibriumBC_exceptions.i'
  expect_err = 'The temperature must be specified'
  detail = 'the user does not specify a function or a variable for the equilibrium flux boundary condition,'
[]


[missing_block]
  type = RunException
  input = 'equilibriumBC_exceptions.i'
  cli_args = 'BCs/left/temperature=100 BCs/left/Ko=some_functor'
  expect_err = 'The subdomain of the enclosure should be specified if Ko and Ea are not constants'
  detail = 'the user does not specify the block to use for the nodal evaluation of equilibrium constant.'
[]


[species_initial_p]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/species_initial_pressures='1 2'"
  detail = "the number of species does not match the number of species initial pressures,"
  expect_err = "Parameter 'species' \(size 1\) and 'species_initial_pressures' \(size 2\) must be the same size if set."
[]


[species_equilibrium]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/equilibrium_constants='1 2'"
  detail = "the number of species does not match the number of sorption equilibrium constants."
  expect_err = "Parameter 'species' \(size 1\) and 'equilibrium_constants' \(size 2\) must be the same size if set"
[]


[temperature]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/temperature=-1 ActionComponents/enc/temperature=-1"
  detail = "the temperature specified is negative,"
  expect_err = "Temperature '-1' on component 'enc' less than or equal to 0 which is not physical"
[]


[temperature_var]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/temperature=temp ActionComponents/enc/temperature=temp AuxVariables/temp/family=LAGRANGE"
  detail = "the temperature is specified as a variable which is not currently specified,"
  expect_err = "Only real values are supported"
[]


[scaling_factor]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/species_scaling_factors=0 ActionComponents/enc/species_scaling_factors=0"
  detail = "the species equation scaling factor is negative or null,"
  expect_err = "Range check failed for parameter ActionComponents/enc/species_scaling_factors"
[]


[initial_pressure]
  type = RunException
  input = sorption_exchange.i
  cli_args = "initial_pressure=-1 Physics/SorptionExchange/0d_solub/species_initial_pressures='-1'"
  detail = "the species initial pressure is negative,"
  expect_err = "Initial condition '-1' for species 'v' on component 'enc' less than 0 which is not physical"
[]


[eq_constant]
  type = RunException
  input = sorption_exchange.i
  cli_args = "Physics/SorptionExchange/0d_solub/equilibrium_constants='0'"
  detail = "the sorption equilibrium constant is negative or null,"
  expect_err = "Equilibrium constant '0' for species 'v' on component 'enc' inferior or equal to 0"
[]


[repeated_species]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/species='x x'"
  detail = "a species is repeated twice in the list of species,"
  expect_err = "Item 'x' specified in vector parameter 'species' is repeated, which is not allowed."
[]


[species_initial_concentration]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/species_initial_concentrations='1 2'"
  detail = "the number of species does not match the number of species initial concentrations,"
  expect_err = "Parameter 'species' \(size 1\) and 'species_initial_concentrations' \(size 2\) must be the same size if set."
[]


[scaling_factors]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/species_scaling_factors='1 2'"
  detail = "the number of species does not match the number of species balance equation scaling factors."
  expect_err = "Parameter 'species' \(size 1\) and 'species_scaling_factors' \(size 2\) must be the same size if set"
[]


[mobile]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/mobile='a b'"
  detail = "the number of trapped species does not match the number of mobile species,"
  expect_err = "Parameter 'species' \(size 1\) and 'mobile' \(size 2\) must be the same size if set"
[]


[alpha_t]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/alpha_t='1 2'"
  detail = "the number of species does not match the number of trapping coefficients,"
  expect_err = "Parameter 'species' \(size 1\) and 'alpha_t' \(size 2\) must be the same size if set"
[]


[Ct0]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/Ct0=''"
  detail = "the number of species does not match the number of trapping Ct0 coefficients,"
  expect_err = "Parameter 'species' \(size 1\) and 'Ct0' \(size 0\) must be the same size if set"
[]


[alpha_r]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/alpha_r='1 2'"
  detail = "the number of species does not match the number of release coefficients,"
  expect_err = "Parameter 'species' \(size 1\) and 'alpha_r' \(size 2\) must be the same size if set"
[]


[detrapping]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/detrapping_energy='1 2'"
  detail = "the number of species does not match the number of detrapping energies."
  expect_err = "Parameter 'species' \(size 1\) and 'detrapping_energy' \(size 2\) must be the same size if set"
[]


[temperature_fun]
  type = RunException
  input = trapping.i
  cli_args = "Physics/SpeciesTrapping/trapped/temperature=temp_fun Functions/temp_fun/type=ConstantFunction Functions/temp_fun/value=1"
  detail = "the temperature is specified as a function when only constants and variables are supported at this time."
  expect_err = "Should be a constant or the name of a variable"
[]


[no_component_yet_asking_for_vars_per_comp]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/separate_variables_per_component=true"
  detail = "no components are specified yet the Physics is set to create variables for each component,"
  expect_err = "Physics is not defined on any Component, this parameter should be set to false"
[]


[species]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/species=''"
  detail = "the species is not specified,"
  expect_err = "Parameter 'species' should not be set to an empty vector"
[]


[alpha_t]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff"
  detail = "the trapping coefficient is not specified,"
  expect_err = "alpha_t: Should not be empty if not using Components"
[]


[trapping_energy]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1"
  detail = "the trapping energy is not specified,"
  expect_err = "trapping_energy: Should not be empty if not using Components"
[]


[N]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1 Physics/SpeciesTrapping/trapped/trapping_energy=0"
  detail = "the number density of the traps is not specified,"
  expect_err = "N: Should not be empty if not using Components"
[]


[Ct0]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1
                  Physics/SpeciesTrapping/trapped/trapping_energy=0 Physics/SpeciesTrapping/trapped/N=1"
  detail = "the trapping Ct0 coefficient is not specified,"
  expect_err = "Ct0: Should not be empty if not using Components"
[]


[trap_per_free]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1
                  Physics/SpeciesTrapping/trapped/trapping_energy=0 Physics/SpeciesTrapping/trapped/N=1 Physics/SpeciesTrapping/trapped/Ct0=1"
  detail = "the number of traps filled per free particle trapped is not specified,"
  expect_err = "trap_per_free: Should not be empty if not using Components"
[]


[alpha_r]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1
                  Physics/SpeciesTrapping/trapped/trapping_energy=0 Physics/SpeciesTrapping/trapped/N=1
                  Physics/SpeciesTrapping/trapped/Ct0=1 Physics/SpeciesTrapping/trapped/trap_per_free=1"
  detail = "the release coefficient is not specified,"
  expect_err = "alpha_r: Should not be empty if not using Components"
[]


[detrapping_energy]
  type = RunException
  input = trapping.i
  cli_args = "ActionComponents/structure/physics=mobile_diff Physics/SpeciesTrapping/trapped/alpha_t=1
                  Physics/SpeciesTrapping/trapped/trapping_energy=0 Physics/SpeciesTrapping/trapped/N=1
                  Physics/SpeciesTrapping/trapped/Ct0=1 Physics/SpeciesTrapping/trapped/trap_per_free=1
                  Physics/SpeciesTrapping/trapped/alpha_r=100"
  detail = "the detrapping energy is not specified."
  expect_err = "detrapping_energy: Should not be empty if not using Components"
[]


[YHx_PCT_error_low_pressure]
  type = RunException
  input = YHx_PCT.i
  cli_args = 'initial_pressure_H2_enclosure_1=0.01'
  expect_err = 'In YHxPCT: pressure'
  requirement = 'The system shall be able to return a warning when the pressure and temperature are outside the range of validity of the YHxPCT model (pressure too low).'
[]


[YHx_PCT_error_high_pressure]
  type = RunException
  input = YHx_PCT.i
  cli_args = 'initial_pressure_H2_enclosure_1=1e7'
  expect_err = 'In YHxPCT: pressure'
  requirement = 'The system shall be able to return a warning when the pressure and temperature are outside the range of validity of the YHxPCT model (pressure too high).'
[]

Introduction
Dependencies
Failure Analysis Requirements