TMAP8 Requirements Traceability Matrix

[equilibriumBC_field_variable]
  type = 'CSVDiff'
  input = 'equilibriumBC.i'
  csvdiff = 'equilibriumBC_out.csv'
  issues = '#134'
  design = 'EquilibriumBC.md'
  requirement = 'The system shall compute the equilibrium flux condition between an enclosure and a structure when a field variable is used for the enclosure variable.'
[]

[Shimada2024_input_check]
  type = RunApp
  input = 'divertor_monoblock.i'
  requirement = 'The system shall maintain a working input file to model heat and tritium transport in a divertor monoblock during pulsed operation.'
  check_input = True
  method = opt
[]

[Shimada2024_run]
  type = Exodiff
  input = 'divertor_monoblock.i'
  exodiff = 'divertor_monoblock_exodus.e'
  cli_args = 'Outputs/exodus/sync_only=true Executioner/end_time=1600 Mesh/ccmg/rings="1 10 8 40" Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-11'
  requirement = 'The system shall model heat and tritium transport in a divertor monoblock during pulsed operation.'
  rel_err = 8e-4 # increasing the relative error tolerance because the test is sensitive to the environment. This error is still small for physical problems.
  heavy = true
  min_parallel = 2
  max_time = 1000
[]

[Abdou2021]
  type = CSVDiff
  input = 'fuel_cycle.i'
  csvdiff = 'fuel_cycle_out.csv'
  requirement = 'The system shall reproduce a consistent solution to an ODE system of equations modeling the tritium fuel cycle.'
[]

[fuel_cycle_Abdou_comparison]
  type = RunCommand
  command = 'python3 plot_comparison.py'
  requirement = 'The system shall be able to generate comparison plots between the simulated solution from TMAP8 and Abdou et al. (2020), modeling tritium fuel cycle.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[gui]
  type = RunCommand
  command = 'python3 fuel_cycle_gui.py --test'
  installation_type = 'in_tree' # See #207
  requirement = 'The system shall be able to open a graphical interface for the tritium fuel cycle example for user training.'
  required_python_packages = 'tempfile tkinter re subprocess matplotlib numpy scipy atexit shutil atexit os'
[]

[fuel_cycle_benchmarking_csvdiff]
  type = CSVDiff
  input = 'fuel_cycle.i'
  cli_args = 'accuracy_time="${units 17000 s}" Outputs/file_base=fuel_cycle_limited_out'
  csvdiff = fuel_cycle_limited_out.csv
  requirement = 'The system shall be able to model the tritium fuel cycle from Meschini et al. (2023).'
[]

[fuel_cycle_benchmarking_heavy_csvdiff]
  type = CSVDiff
  input = 'fuel_cycle.i'
  csvdiff = fuel_cycle_out.csv
  heavy = true
  requirement = 'The system shall be able to the model tritium fuel cycle from Meschini et al. (2023) with fine time step.'
[]

[fuel_cycle_benchmarking_comparison]
  type = RunCommand
  command = 'python3 comparison_fuel_cycle_benchmark.py'
  requirement = 'The system shall be able to generate comparison plots between the simulated solution from TMAP8 and Meschini et al. (2023), modeling tritium fuel cycle.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[Sievert_non_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'Outputs/file_base=interface_sorption_Sievert_non_ad_out'
  exodiff = 'interface_sorption_Sievert_non_ad_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Sievert_non_ad_penalty]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using a penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e5
                Outputs/file_base=interface_sorption_Sievert_non_ad_penalty_out'
  exodiff = 'interface_sorption_Sievert_non_ad_penalty_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Sievert_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using automatic differentiation.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'Kernels/u1/type=ADMatDiffusion
                Kernels/u2/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_Sievert_ad_out'
  exodiff = 'interface_sorption_Sievert_ad_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Sievert_ad_penalty]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using automatic differentiation and a penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'Kernels/u1/type=ADMatDiffusion
                Kernels/u2/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e5
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_Sievert_ad_penalty_out'
  exodiff = 'interface_sorption_Sievert_ad_penalty_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Henry_non_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'InterfaceKernels/interface/n_sorption=1 Outputs/file_base=interface_sorption_Henry_non_ad_out'
  exodiff = 'interface_sorption_Henry_non_ad_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Henry_non_ad_penalty]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using a penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e5
                InterfaceKernels/interface/n_sorption=1
                Outputs/file_base=interface_sorption_Henry_non_ad_penalty_out'
  exodiff = 'interface_sorption_Henry_non_ad_penalty_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Henry_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using automatic differentiation.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'Kernels/u1/type=ADMatDiffusion
                Kernels/u2/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                InterfaceKernels/interface/n_sorption=1
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_Henry_ad_out'
  exodiff = 'interface_sorption_Henry_ad_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Henry_ad_penalty]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using automatic differentiation and a penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption.i'
  cli_args = 'Kernels/u1/type=ADMatDiffusion
                Kernels/u2/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                InterfaceKernels/interface/n_sorption=1
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e5
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_Henry_ad_penalty_out'
  exodiff = 'interface_sorption_Henry_ad_penalty_out.e'
  abs_zero = 1e-7
  mesh_mode = REPLICATED
[]

[Sievert_transient_non_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions during transient simulations.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'Outputs/file_base=interface_sorption_transient_Sievert_non_ad_out'
  exodiff = 'interface_sorption_transient_Sievert_non_ad_out.e'
  mesh_mode = REPLICATED
[]

[Sievert_transient_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using automatic differentiation during transient simulations.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_transient_Sievert_ad_out
                Executioner/nl_rel_tol=1e-7
                Executioner/nl_abs_tol=1e-12'
  exodiff = 'interface_sorption_transient_Sievert_ad_out.e'
  mesh_mode = REPLICATED
[]

[Sievert_transient_non_ad_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'unit_scale=4
                unit_scale_neighbor=6
                Outputs/file_base=interface_sorption_transient_Sievert_non_ad_scaling_out'
  exodiff = 'interface_sorption_transient_Sievert_non_ad_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Sievert_transient_ad_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using automatic differentiation during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'unit_scale=4
                unit_scale_neighbor=6
                Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_transient_Sievert_ad_scaling_out
                Executioner/nl_abs_tol=5e-12'
  exodiff = 'interface_sorption_transient_Sievert_ad_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Sievert_transient_non_ad_penalty_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using a penalty-enforced flux balance during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'unit_scale=4
                unit_scale_neighbor=6
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Sievert_non_ad_penalty_scaling_out'
  exodiff = 'interface_sorption_transient_Sievert_non_ad_penalty_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Sievert_transient_non_ad_penalty_scaling_comparison]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using a penalty-enforced flux balance during transient simulations with unit scaling on both variables and provide similar results to the approach without the penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'unit_scale=4
                unit_scale_neighbor=6
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Sievert_non_ad_penalty_scaling_comparison_out'
  exodiff = 'interface_sorption_transient_Sievert_non_ad_penalty_scaling_comparison_out.e'
  # compared to Sievert_transient_non_ad_scaling because it is supposed to be equivalent,
  # but the residuals are expected to be different when the penalty is used, so the relative error tolerance is increased here
  rel_err = 8.5e-2
  mesh_mode = REPLICATED
[]

[Sievert_transient_ad_penalty_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Sievert law in isothermal conditions using automatic differentiation and a penalty-enforced flux balance during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'unit_scale=4
                unit_scale_neighbor=6
                Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Sievert_ad_penalty_scaling_out
                Executioner/nl_abs_tol=1e-12'
  exodiff = 'interface_sorption_transient_Sievert_ad_penalty_scaling_out.e'
  # compared to Sievert_transient_non_ad_scaling because it is supposed to be equivalent,
  # but the residuals are expected to be different when the penalty is used, so the relative error tolerance is increased here
  rel_err = 8.5e-2
  mesh_mode = REPLICATED
[]

[Henry_transient_non_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions during transient simulations.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                Outputs/file_base=interface_sorption_transient_Henry_non_ad_out'
  exodiff = 'interface_sorption_transient_Henry_non_ad_out.e'
  mesh_mode = REPLICATED
[]

[Henry_transient_ad]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using automatic differentiation during transient simulations.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_transient_Henry_ad_out
                Executioner/nl_rel_tol=1e-7
                Executioner/nl_abs_tol=1e-12'
  exodiff = 'interface_sorption_transient_Henry_ad_out.e'
  mesh_mode = REPLICATED
[]

[Henry_transient_non_ad_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                unit_scale=4
                unit_scale_neighbor=6
                Outputs/file_base=interface_sorption_transient_Henry_non_ad_scaling_out'
  exodiff = 'interface_sorption_transient_Henry_non_ad_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Henry_transient_ad_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using automatic differentiation during transient simulation with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                unit_scale=4
                unit_scale_neighbor=6
                Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                Outputs/file_base=interface_sorption_transient_Henry_ad_scaling_out
                Executioner/nl_abs_tol=5e-12'
  exodiff = 'interface_sorption_transient_Henry_ad_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Henry_transient_non_ad_penalty_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using a penalty-enforced flux balance during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                unit_scale=4
                unit_scale_neighbor=6
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Henry_non_ad_penalty_scaling_out'
  exodiff = 'interface_sorption_transient_Henry_non_ad_penalty_scaling_out.e'
  mesh_mode = REPLICATED
[]

[Henry_transient_non_ad_penalty_scaling_comparison]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using a penalty-enforced flux balance during transient simulations with unit scaling on both variables and provide similar results to the approach without the penalty-enforced flux balance.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                unit_scale=4
                unit_scale_neighbor=6
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Henry_non_ad_penalty_scaling_comparison_out'
  exodiff = 'interface_sorption_transient_Henry_non_ad_penalty_scaling_comparison_out.e'
  # compared to Henry_transient_non_ad_scaling because it is supposed to be equivalent,
  # but the residuals are expected to be different when the penalty is used, so the relative error tolerance is increased here
  rel_err = 5.4e-2
  mesh_mode = REPLICATED
[]

[Henry_transient_ad_penalty_scaling]
  requirement = 'The system shall have the capability to enforce interfacial conditions based on the Henry law in isothermal conditions using automatic differentiation and a penalty-enforced flux balance during transient simulations with unit scaling on both variables.'
  type = 'Exodiff'
  input = 'interface_sorption_transient.i'
  cli_args = 'n_sorption=1
                unit_scale=4
                unit_scale_neighbor=6
                Kernels/u1_time_derivative/type=ADTimeDerivative
                Kernels/u2_time_derivative/type=ADTimeDerivative
                Kernels/temperature/type=ADTimeDerivative
                Kernels/u1_diffusion/type=ADMatDiffusion
                Kernels/u2_diffusion/type=ADMatDiffusion
                Kernels/temperature/type=ADHeatConduction
                InterfaceKernels/interface/type=ADInterfaceSorption
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_1/type=ADGenericConstantMaterial
                Materials/properties_2/type=ADGenericConstantMaterial
                Postprocessors/flux_inner/type=ADSideDiffusiveFluxIntegral
                Postprocessors/flux_outer/type=ADSideDiffusiveFluxIntegral
                InterfaceKernels/interface/use_flux_penalty=true
                InterfaceKernels/interface/flux_penalty=1e1
                Outputs/file_base=interface_sorption_transient_Henry_ad_penalty_scaling_out
                Executioner/nl_abs_tol=1e-12'
  exodiff = 'interface_sorption_transient_Henry_ad_penalty_scaling_out.e'
  # compared to Henry_transient_non_ad_scaling because it is supposed to be equivalent,
  # but the residuals are expected to be different when the penalty is used, so the relative error tolerance is increased here
  rel_err = 5.4e-2
  mesh_mode = REPLICATED
[]

[ad_mat_coupled_annihilation]
  type = 'Exodiff'
  input = 'ad_mat_coupled_annihilation_test.i'
  exodiff = 'ad_mat_coupled_annihilation_test_out.e'
  issues = '#35'
  design = '/ADMatCoupledDefectAnnihilation.md'
  requirement = 'The system shall compute the annihilation of a concentration variable towards its equilibrium in the material.'
[]

[./ad_mat_reaction_flexible]
  type = 'Exodiff'
  input = 'ad_mat_reaction_flexible_test.i'
  exodiff = 'ad_mat_reaction_flexible_test_out.e'
  issues = '#23'
  design = '/ADMatReactionFlexible.md'
  requirement = 'The system shall compute the reaction contribution in the material.'
[../]

[pore_scale_microstructure_formation_slice]
  type = Exodiff
  input = '2D_microstructure_reader_smoothing_base.i pore_structure_closed.params'
  cli_args = "output_name=pore_structure_closed_slice/pore_structure_closed_slice
                num_nodes_x=60
                num_nodes_y=2
                domain_start_x=4000
                domain_end_x=5000
                domain_end_y=1
                Outputs/exodus/execute_on='INITIAL TIMESTEP_END FINAL'"
  exodiff = pore_structure_closed_slice/pore_structure_closed_slice.e
  requirement = 'The system shall be able to read a microstructure image and import it to perform a simulation to smoothen the interfaces.'
  capabilities = 'vtk'
[]

[pore_scale_transport_slice]
  type = Exodiff
  input = '2D_absorption_base.i pore_structure_closed_absorption.params'
  cli_args = "input_name=gold/pore_structure_closed_slice/pore_structure_closed_slice.e
                output_name=pore_structure_closed_slice/pore_structure_closed_absorption_slice
                BCs/tritium_2g_sides_d/boundary=right
                Executioner/num_steps=20"
  exodiff = pore_structure_closed_slice/pore_structure_closed_absorption_slice.e
  custom_cmp = '2D_absorption_closed_slice.exodiff' # neeeded for -p2 tests
  requirement = 'The system shall be able to import an existing microstructure and perform a simulation of tritium transport during
                   absorption at the pore scale.'
[]

[val-2a_csvdiff]
  type = CSVDiff
  input = 'val-2a.i'
  csvdiff = val-2a_out.csv
  requirement = 'The system shall be able to model deuterium ion implantation in a steel alloy for comparison with experimental results, particularly focusing on the permeation flux.'
[]

[val-2a_exodiff]
  type = Exodiff
  input = 'val-2a.i'
  exodiff = val-2a_out.e
  prereq = val-2a_csvdiff
  should_execute = false # this test relies on the output files from val-2a_csvdiff, so it shouldn't be run twice
  requirement = 'The system shall be able to model deuterium ion implantation in a steel alloy for comparison with experimental results, focused on the full set of simulation output including deuterium concentration, recombination coefficient, and dissociation coefficient.'
[]

[val-2a_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2a.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of validation case 2a, modeling deuterium ion implantation in a steel alloy.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[val-2b_heavy_csv]
  type = CSVDiff
  input = val-2b.i
  cli_args = 'Outputs/file_base=val-2b_heavy_out'
  csvdiff = val-2b_heavy_out.csv
  heavy = true
  requirement = 'The system shall be able to model diffusion of deuterium in a beryllium sample and generate CSV data output for comparison to experimental results.'
[]

[val-2b_heavy_exodus]
  type = Exodiff
  input = val-2b.i
  cli_args = 'Outputs/file_base=val-2b_heavy_out'
  exodiff = val-2b_heavy_out.e
  prereq = val-2b_heavy_csv
  should_execute = false # this test relies on the output files from val-2b_heavy_csv, so it shouldn't be run twice
  heavy = true
  abs_zero = 1e-8
  requirement = 'The system shall be able to model diffusion of deuterium in a beryllium sample and generate field and material property data output in the Exodus format for comparison to experimental results.'
[]

[val-2b_exodus]
  type = Exodiff
  input = val-2b.i
  cli_args = "Outputs/file_base=val-2b_out
                Executioner/nl_abs_tol=1e-11
                charge_time='${units 10 h -> s}'
                cooldown_duration=$'{units 1 h -> s}'
                dt_max_large=$'{units 500 s}'
                dt_max_small=$'{units 200 s}'
                dt_start_charging=$'{units 300 s}'
                dt_start_cooldown=$'{units 100 s}'
                dt_start_desorption=$'{units 10 s}'"
  exodiff = val-2b_out.e
  abs_zero = 1e-8
  requirement = 'The system shall be able to model diffusion of deuterium in beryllium sample with a short runtime suitable for regular regression testing.'
[]

[val-2b_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2b.py'
  requirement = 'The system shall be able to generate comparison plots between simulated solutions and experimental data of validation case val-2b, modeling diffusion and release of deuterium in a beryllium sample.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[val-2c_immediate_injection_csv]
  type = CSVDiff
  input = val-2c_immediate_injection.i
  csvdiff = val-2c_immediate_injection_csv.csv
  requirement = 'The system shall be able to model the Test Cell Release Experiment (val-2c) with immediate T2 injection.'
  max_parallel = 1 # see #200
  recover = false # see #196
[]

[val-2c_immediate_injection_exodus]
  type = Exodiff
  input = val-2c_immediate_injection.i
  prereq = val-2c_immediate_injection_csv
  should_execute = false # this test relies on the output files from val-2c_immediate_injection_csv, so it shouldn't be run twice
  exodiff = val-2c_immediate_injection_out.e
  custom_cmp = 'val-2c_exodus.exodiff'
  requirement = 'The system shall be able to model the Test Cell Release Experiment (val-2c) with immediate T2 injection and properly compute the exodus file.'
  max_parallel = 1 # see #200
  recover = false # see #196
[]

[val-2c_delay_csv]
  type = CSVDiff
  input = val-2c_delay.i
  csvdiff = val-2c_delay_csv.csv
  requirement = 'The system shall be able to model the Test Cell Release Experiment (val-2c) with delayed T2 injection.'
  abs_zero = 1e-8
  max_parallel = 1 # see #200
  recover = false # see #196
[]

[val-2c_delay_exodus]
  type = Exodiff
  input = val-2c_delay.i
  prereq = val-2c_delay_csv
  should_execute = false # this test relies on the output files from val-2c_delay_csv, so it shouldn't be run twice
  exodiff = val-2c_delay_out.e
  custom_cmp = 'val-2c_exodus.exodiff'
  requirement = 'The system shall be able to model the Test Cell Release Experiment (val-2c) with delayed T2 injection and properly compute the exodus file.'
  max_parallel = 1 # see #200
  recover = false # see #196
[]

[val-2c_delay_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2c.py'
  requirement = 'The system shall be able to generate comparison plots between simulated solutions and experimental data of validation cases val-2c, modeling a Test Cell Release Experiment.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[val-2d_csvdiff]
  type = CSVDiff
  input = 'val-2d.i'
  cli_args = 'Mesh/active=cartesian_mesh_coarse Postprocessors/max_time_step_size/function=max_dt_size_function_coarse Executioner/nl_rel_tol=8e-9 Outputs/file_base=val-2d_limited_out'
  csvdiff = val-2d_limited_out.csv
  requirement = 'The system shall be able to model thermal desorption spectroscopy on Tungsten.'
[]

[val-2d_exodiff]
  type = Exodiff
  input = 'val-2d.i'
  cli_args = 'Mesh/active=cartesian_mesh_coarse Postprocessors/max_time_step_size/function=max_dt_size_function_coarse Executioner/nl_rel_tol=8e-9 Outputs/file_base=val-2d_limited_out'
  exodiff = val-2d_limited_out.e
  prereq = val-2d_csvdiff
  should_execute = false # this test relies on the output files from val-2d_csvdiff, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal desorption spectroscopy on Tungsten to include full set of simulation outputs, including tritium concentration, diffusion flux, and trapping properties.'
[]

[val-2d_heavy_csvdiff]
  type = CSVDiff
  input = 'val-2d.i'
  csvdiff = val-2d_out.csv
  heavy = true
  requirement = 'The system shall be able to model thermal desorption spectroscopy on Tungsten with fine mesh and time step to compare with the desorption flux from experiment results.'
[]

[val-2d_heavy_exodiff]
  type = Exodiff
  input = 'val-2d.i'
  exodiff = val-2d_out.e
  prereq = val-2d_csvdiff
  heavy = true
  should_execute = false # this test relies on the output files from val-2d_csvdiff, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal desorption spectroscopy on Tungsten with fine mesh and time step to include full set of simulation outputs, including tritium concentration, diffusion flux, and trapping properties.'
[]

[val-2d_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2d.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of validation case 2d, modeling thermal desorption spectroscopy on Tungsten.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[val-2ea_csvdiff]
  type = CSVDiff
  input = val-2ea.i
  csvdiff = val-2ea_out.csv
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.05 mm thick membrane at 825 K to generate CSV data for use in comparisons with the experimental data.'
[]

[val-2ea]
  type = Exodiff
  input = val-2ea.i
  should_execute = False # this test relies on the output files from val-2ea_csvdiff, so it shouldn't be run twice
  exodiff = val-2ea_out.e
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.05 mm thick membrane at 825 K.'
[]

[val-2eb_csvdiff]
  type = CSVDiff
  input = val-2ea.i
  csvdiff = val-2eb_out.csv
  cli_args = 'slab_thickness="${units 2.5e-5 m -> mum}" file_name="val-2eb_out"'
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.025 mm thin membrane at 825 K to generate CSV data for use in comparisons with the experimental data.'
[]

[val-2eb]
  type = Exodiff
  input = val-2ea.i
  should_execute = False # this test relies on the output files from val-2eb_csvdiff, so it shouldn't be run twice
  exodiff = val-2eb_out.e
  cli_args = 'slab_thickness="${units 2.5e-5 m -> mum}" file_name="val-2eb_out"'
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.025 mm thin membrane at 825 K.'
[]

[val-2ec_csvdiff]
  type = CSVDiff
  input = val-2ea.i
  csvdiff = val-2ec_out.csv
  cli_args = 'temperature="${units 865 K}" slab_thickness="${units 2.5e-5 m -> mum}" file_name="val-2ec_out"'
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.025 mm thin membrane at 865 K to generate CSV data for use in comparisons with the experimental data.'
[]

[val-2ec]
  type = Exodiff
  input = val-2ea.i
  should_execute = False # this test relies on the output files from val-2ec_csvdiff, so it shouldn't be run twice
  exodiff = val-2ec_out.e
  cli_args = 'temperature="${units 865 K}" slab_thickness="${units 2.5e-5 m -> mum}" file_name="val-2ec_out"'
  requirement = 'The system shall be able to model permeation of Deuterium from a 0.025 mm thin membrane at 865 K and generate an exodus file.'
[]

[val-2ed_csvdiff]
  type = CSVDiff
  input = val-2ed.i
  csvdiff = val-2ed_out.csv
  requirement = 'The system shall be able to model permeation of mixture gas from a 0.025 mm thin membrane at 870 K using lawdep boundary conditions to generate CSV data for use in comparisons with the experimental data.'
[]

[val-2ed]
  type = Exodiff
  input = val-2ed.i
  should_execute = False # this test relies on the output files from val-2ed_csvdiff, so it shouldn't be run twice
  exodiff = val-2ed_out.e
  requirement = 'The system shall be able to model permeation of mixture gas with chemical reaction from a  0.025 mm thin membrane at 870 K using lawdep boundary conditions and generate an exodus file.'
[]

[val-2ee_csvdiff]
  type = CSVDiff
  input = val-2ee.i
  csvdiff = val-2ee_out.csv
  requirement = 'The system shall be able to model permeation of mixture gas from a 0.025 mm thin membrane at 870 K using ratedep boundary conditions to generate CSV data for use in comparisons with the experimental data.'
[]

[val-2ee]
  type = Exodiff
  input = val-2ee.i
  should_execute = False # this test relies on the output files from val-2ee_csvdiff, so it shouldn't be run twice
  exodiff = val-2ee_out.e
  requirement = 'The system shall be able to model permeation of mixture gas with chemical reaction from a 0.025 mm thin membrane at 870 K using ratedep boundary conditions and generate an exodus file.'
[]

[ver-2e_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2e.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and experimental data of validation case 2e, modeling the permeation of Deuterium from a membrane.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[val-2f_light_csv]
  type = CSVDiff
  input = val-2f.i
  cli_args = "sample_thickness='${units 0.9e-4 m -> mum}' ix1=1 ix2=1 ix3=1 ix4=1 ix5=1 charge_time='${units 1e-1 s}' cooldown_duration='${units 1e-1 s}' Executioner/TimeStepper/timestep_limiting_postprocessor=max_time_step_size_coarse Outputs/file_base=val-2f_light_out"
  csvdiff = val-2f_light_out.csv
  requirement = 'The system shall be able to model self-damaged tungsten effects on deuterium transport and generate CSV data output with a short runtime and coarse mesh testing.'
  max_parallel = 1 # see #200
[]

[val-2f_light_exodus]
  type = Exodiff
  input = val-2f.i
  cli_args = "sample_thickness='${units 0.9e-4 m -> mum}' ix1=1 ix2=1 ix3=1 ix4=1 ix5=1 charge_time='${units 1e-1 s}' cooldown_duration='${units 1e-1 s}' Executioner/TimeStepper/timestep_limiting_postprocessor=max_time_step_size_coarse Outputs/file_base=val-2f_light_out"
  exodiff = val-2f_light_out.e
  prereq = val-2f_light_csv
  should_execute = false # this test relies on the output files from val-2f_light_csv, so it shouldn't be run twice
  requirement = 'The system shall be able to model self-damaged tungsten effects on deuterium transport with a short runtime and coarse mesh testing.'
  max_parallel = 1 # see #200
[]

[val-2f_heavy_csv]
  type = CSVDiff
  input = val-2f.i
  csvdiff = val-2f_out.csv
  heavy = true
  requirement = 'The system shall be able to model self-damaged tungsten effects on deuterium transport and generate CSV data output.'
  max_time = 2000
[]

[val-2f_heavy_exodus]
  type = Exodiff
  input = val-2f.i
  exodiff = val-2f_out.e
  heavy = true
  prereq = val-2f_heavy_csv
  should_execute = false # this test relies on the output files from val-2f_heavy_csv, so it shouldn't be run twice
  requirement = 'The system shall be able to model self-damaged tungsten effects on deuterium transport.'
[]

[val-2f_heavy_csv_inf_recombination]
  type = CSVDiff
  input = val-2f.i
  cli_args = "ix1='${fparse 500}' ix4='${fparse 500}'
                AuxVariables/active='bounds_dummy temperature'
                BCs/active='left_concentration_sieverts right_concentration_sieverts'
                NodalKernels/trapping_1/Ct0='trap_distribution_function_1_inf'
                NodalKernels/trapping_2/Ct0='trap_distribution_function_2_inf'
                NodalKernels/trapping_3/Ct0='trap_distribution_function_3_inf'
                NodalKernels/trapping_4/Ct0='trap_distribution_function_4_inf'
                NodalKernels/trapping_5/Ct0='trap_distribution_function_5_inf'
                Materials/active='diffusivity_W_func diffusivity_nonAD'
                Postprocessors/active='integral_source_deuterium scaled_implanted_deuterium integral_deuterium_concentration scaled_mobile_deuterium flux_surface_left_sieverts scaled_flux_surface_left_sieverts flux_surface_right_sieverts scaled_flux_surface_right_sieverts temperature diffusion_W max_time_step_size max_time_step_size_coarse integral_trapped_concentration_1 scaled_trapped_deuterium_1 integral_trapped_concentration_2 scaled_trapped_deuterium_2 integral_trapped_concentration_3 scaled_trapped_deuterium_3 integral_trapped_concentration_4 scaled_trapped_deuterium_4 integral_trapped_concentration_5 scaled_trapped_deuterium_5 integral_trapped_concentration_intrinsic scaled_trapped_deuterium_intrinsic'
                Postprocessors/max_time_step_size/function=max_dt_size_function_inf
                Outputs/file_base='val-2f_out_inf_recombination'"
  csvdiff = val-2f_out_inf_recombination.csv
  heavy = true
  requirement = 'The system shall be able to model self-damaged tungsten effects on deuterium transport and generate CSV data output, for the infinite recombination case.'
  max_time = 3600
[]

[val-2f_comparison]
  type = RunCommand
  command = 'python3 comparison_val-2f.py'
  requirement = 'The system shall be able to generate comparison plots between simulated solutions and experimental data of validation case val-2f, modeling self-damaged tungsten effects on deuterium transport.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1a]
  type = Exodiff
  input = ver-1a.i
  exodiff = ver-1a_test_out.e
  cli_args = 'Executioner/dt=5 Mesh/nx=10 Outputs/file_base=ver-1a_test_out'
  requirement = 'The system shall be able to model species diffusion through a structure, originating from a depleting source enclosure.'
[]

[ver-1a_heavy]
  type = Exodiff
  input = ver-1a.i
  exodiff = ver-1a_out.e
  requirement = 'The system shall be able to model species diffusion through a structure, originating from a depleting source enclosure, with the fine mesh and timestep required to match the analytical solution.'
  heavy = true
[]

[ver-1a_heavy_csvdiff]
  type = CSVDiff
  input = ver-1a.i
  should_execute = False # this test relies on the output files from ver-1a_heavy, so it shouldn't be run twice
  csvdiff = ver-1a_csv.csv
  requirement = 'The system shall be able to model species diffusion through a structure, originating from a depleting source enclosure, with the fine mesh and timestep required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution.'
  heavy = true
  prereq = ver-1a_heavy
[]

[ver-1a_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1a.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1a, modeling species diffusion through a structure, originating from a depleting source enclosure.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1b]
  type = Exodiff
  input = ver-1b.i
  exodiff = ver-1b_test_out.e
  cli_args = 'Executioner/dt=10 Mesh/nx=100 Outputs/exodus/file_base=ver-1b_test_out'
  requirement = 'The system shall be able to model transient diffusion through a slab with a constant concentration boundary condition as the species source.'
[]

[ver-1b_heavy]
  type = Exodiff
  input = ver-1b.i
  exodiff = ver-1b_out.e
  heavy = true
  requirement = 'The system shall be able to model transient diffusion through a slab with a constant concentration boundary condition as the species source with the fine mesh and time step required to match the analytical solution.'
[]

[ver-1b_heavy_csvdiff]
  type = CSVDiff
  input = ver-1b.i
  should_execute = False # this test relies on the output files from ver-1b_heavy, so it shouldn't be run twice
  csvdiff = ver-1b_csv.csv
  requirement = 'The system shall be able to model transient diffusion through a slab with a constant concentration boundary condition as the species source, with the fine mesh and time step required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution over time.'
  heavy = true
  prereq = ver-1b_heavy
[]

[ver-1b_heavy_lineplot]
  type = CSVDiff
  input = ver-1b.i
  should_execute = False # this test relies on the output files from ver-1b_heavy, so it shouldn't be run twice
  csvdiff = ver-1b_vector_postproc_line_0250.csv
  requirement = 'The system shall be able to model transient diffusion through a slab with a constant concentration boundary condition as the species source, with the fine mesh and timestep required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution for the profile concentration.'
  heavy = true
  prereq = ver-1b_heavy
[]

[ver-1b_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1b.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1b, modeling transient diffusion through a slab with a constant concentration boundary condition as the species source.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1c_TMAP4]
  type = Exodiff
  input = 'ver-1c.i'
  exodiff = ver-1c_tmap4.e
  cli_args = 'BCs/lhs/type=NeumannBC Outputs/file_base=ver-1c_tmap4'
  requirement = 'The system shall be able to model species permeation into an unloaded portion of a slab from a pre-loaded portion with boundary conditions consistent with TMAP4.'
[]

[ver-1c_TMAP7]
  type = Exodiff
  input = 'ver-1c.i'
  exodiff = ver-1c_tmap7.e
  cli_args = 'BCs/lhs/type=DirichletBC Outputs/file_base=ver-1c_tmap7'
  requirement = 'The system shall be able to model species permeation into an unloaded portion of a slab from a pre-loaded portion with boundary conditions consistent with TMAP7'
[]

[ver-1c_TMAP4_csvdiff]
  type = CSVDiff
  input = ver-1c.i
  should_execute = False # this test relies on the output files from ver-1c, so it shouldn't be run twice
  csvdiff = ver-1c_tmap4.csv
  requirement = 'The system shall be able to model species permeation into an unloaded portion of a slab from a pre-loaded portion to generate CSV data for use in comparisons with the analytic solution over time for the TMAP4 verification case.'
  prereq = ver-1c_TMAP4
[]

[ver-1c_TMAP7_csvdiff]
  type = CSVDiff
  input = ver-1c.i
  should_execute = False # this test relies on the output files from ver-1c, so it shouldn't be run twice
  csvdiff = ver-1c_tmap7.csv
  requirement = 'The system shall be able to model species permeation into an unloaded portion of a slab from a pre-loaded portion to generate CSV data for use in comparisons with the analytic solution over time for the TMAP7 verification case.'
  prereq = ver-1c_TMAP7
[]

[ver-1c_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1c.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1c, modeling species permeation into an unloaded portion of a slab from a pre-loaded portion for both the TMAP4 and TMAP7 verification cases.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1d_diffusion_limited]
  type = Exodiff
  input = ver-1d-diffusion.i
  exodiff = ver-1d_diffusion_limited_test_exodus.e
  cli_args = "Mesh/nx=20 Executioner/num_steps=300 Outputs/file_base=ver-1d_diffusion_limited_test_exodus"
  requirement = 'The system shall be able to model a breakthrough problem where diffusion is the rate limiting process.'
[]

[ver-1d_diffusion_limited_heavy]
  type = Exodiff
  heavy = true
  input = ver-1d-diffusion.i
  exodiff = ver-1d-diffusion_out.e
  requirement = 'The system shall be able to model a breakthrough problem where diffusion is the rate limiting process, with the fine mesh and time step required to match the analytical solution for the verification case.'
[]

[ver-1d_diffusion_limited_heavy_csvdiff]
  type = CSVDiff
  heavy = true
  input = ver-1d-diffusion.i
  should_execute = false # this test relies on the output files from ver-1d_diffusion_limited_heavy, so it shouldn't be run twice
  csvdiff = ver-1d-diffusion_out.csv
  prereq = ver-1d_diffusion_limited_heavy
  requirement = 'The system shall be able to model a breakthrough problem where diffusion is the rate limiting process,
    with the fine mesh and time step required to match the analytical solution for the verification case and generate CSV data for use in comparisons with the analytic solution.'
[]

[ver-1d_trapping_limited]
  type = Exodiff
  input = ver-1d-trapping.i
  exodiff = ver-1d_trapping_limited_test_exodus.e
  requirement = 'The system shall be able to model a breakthrough problem where trapping is the rate limiting process.'
  cli_args = "Mesh/nx=20 Executioner/num_steps=300 Outputs/file_base=ver-1d_trapping_limited_test_exodus"
[]

[ver-1d_trapping_limited_heavy]
  type = Exodiff
  heavy = true
  input = ver-1d-trapping.i
  exodiff = ver-1d-trapping_out.e
  requirement = 'The system shall be able to model a breakthrough problem where trapping is the rate limiting process with the fine mesh and time step required to match the analytical solution for the verification case.'
[]

[ver-1d_trapping_limited_heavy_csvdiff]
  type = CSVDiff
  heavy = true
  input = ver-1d-trapping.i
  should_execute = false # this test relies on the output files from ver-1d_trapping_limited_heavy, so it shouldn't be run twice
  csvdiff = ver-1d-trapping_out.csv
  prereq = ver-1d_trapping_limited_heavy
  requirement = 'The system shall be able to model a breakthrough problem where trapping is the rate limiting process with the fine mesh and time step required to match the analytical solution for the verification case and generate CSV data for use in comparisons with the analytic solution.'
[]

[ver-1d_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1d.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solutions of verification cases 1d, modeling a breakthrough problem where diffusion and trapping are the rate limiting processes.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1dc_limited]
  type = Exodiff
  input = ver-1dc.i
  exodiff = ver-1dc_limited_test_exodus.e
  cli_args = "nx_num=20 Executioner/num_steps=300 Outputs/file_base=ver-1dc_limited_test_exodus"
  requirement = 'The system shall be able to model a breakthrough problem of multiple traps.'
[]

[ver-1dc_limited_heavy]
  type = Exodiff
  heavy = true
  input = ver-1dc.i
  exodiff = ver-1dc_out.e
  requirement = 'The system shall be able to model a breakthrough problem of multiple traps, with the fine mesh and time step required to match the analytical solution for the verification case.'
[]

[ver-1dc_limited_heavy_csvdiff]
  type = CSVDiff
  heavy = true
  input = ver-1dc.i
  should_execute = false # this test relies on the output files from ver-1dc_limited_heavy, so it shouldn't be run twice
  csvdiff = ver-1dc_out.csv
  prereq = ver-1dc_limited_heavy
  requirement = 'The system shall be able to model a breakthrough problem of multiple traps,
    with the fine mesh and time step required to match the analytical solution for the verification case and generate CSV data for use in comparisons with the analytic solution.'
[]

[ver-1d_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1dc.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solutions of verification cases 1dc, modeling a breakthrough problem of multiple traps.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1dc-mms]
  type = PythonUnitTest
  input = test.py
  test_case = TestMMS
  heavy = true
  requirement = 'The system shall show second order spatial convergence for a diffusion-trapping-release test case.'
  method = '!dbg'
  required_python_packages = 'pandas matplotlib'
  # skip test if test is being run out-of-tree. Issue Ref: idaholab/moose#25279
  installation_type = in_tree
[]

[ver-1dd]
  type = Exodiff
  input = ver-1dd.i
  exodiff = ver-1dd_out.e
  requirement = 'The system shall be able to model a breakthrough problem without traps.'
[]

[ver-1dd_csvdiff]
  type = CSVDiff
  input = ver-1dd.i
  should_execute = false # this test relies on the output files from ver-1dd, so it shouldn't be run twice
  csvdiff = ver-1dd_out.csv
  prereq = ver-1dd
  requirement = 'The system shall be able to model a breakthrough problem without traps, and generate CSV data for use in comparisons with the analytic solution.'
[]

[ver-1d_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1dd.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solutions of verification cases 1dd, modeling a breakthrough problem without traps.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1e]
  type = Exodiff
  input = ver-1e.i
  exodiff = ver-1e_test_exodus.e
  cli_args = 'Executioner/dt=0.2 Executioner/end_time=50 Outputs/exodus/file_base=ver-1e_test_exodus'
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source.'
  exodiff_opts = '-match_ids'
  map = False
[]

[ver-1e_TMAP4_heavy]
  type = Exodiff
  input = ver-1e.i
  exodiff = TMAP4.e
  heavy = true
  cli_args = 'T_SiC=63e-6 D_ver=8e-6 Outputs/exodus/file_base=TMAP4 Outputs/csv/file_base=TMAP4 Outputs/vector_postproc/file_base=TMAP4_vector_postproc'
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and time step required to match the analytical solution for the TMAP4 verification case.'
  exodiff_opts = '-match_ids'
  map = False
[]

[ver-1e_TMAP7_heavy]
  type = Exodiff
  input = ver-1e.i
  exodiff = TMAP7.e
  heavy = true
  cli_args = 'T_SiC=66e-6 D_ver=15.75e-6 Outputs/exodus/file_base=TMAP7 Outputs/csv/file_base=TMAP7 Outputs/vector_postproc/file_base=TMAP7_vector_postproc'
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and time step required to match the analytical solution for the TMAP7 verification case.'
  exodiff_opts = '-match_ids'
  map = False
[]

[ver-1e_TMAP4_heavy_csvdiff]
  type = CSVDiff
  input = ver-1e.i
  should_execute = False # this test relies on the output files from ver-1e_TMAP4_heavy, so it shouldn't be run twice
  csvdiff = TMAP4.csv
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and time step required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution over time for the TMAP4 verification case.'
  heavy = true
  prereq = ver-1e_TMAP4_heavy
[]

[ver-1e_TMAP7_heavy_csvdiff]
  type = CSVDiff
  input = ver-1e.i
  should_execute = False # this test relies on the output files from ver-1e_TMAP7_heavy, so it shouldn't be run twice
  csvdiff = TMAP7.csv
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and time step required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution over time for the TMAP7 verification case.'
  heavy = true
  prereq = ver-1e_TMAP7_heavy
[]

[ver-1e_TMAP4_heavy_lineplot]
  type = CSVDiff
  input = ver-1e.i
  should_execute = False # this test relies on the output files from ver-1e_TMAP4_heavy, so it shouldn't be run twice
  csvdiff = TMAP4_vector_postproc_line_0548.csv
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and timestep required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution for the profile concentration for the TMAP4 verification case.'
  heavy = true
  prereq = ver-1e_TMAP4_heavy
[]

[ver-1e_TMAP7_heavy_lineplot]
  type = CSVDiff
  input = ver-1e.i
  should_execute = False # this test relies on the output files from ver-1e_TMAP7_heavy, so it shouldn't be run twice
  csvdiff = TMAP7_vector_postproc_line_0548.csv
  requirement = 'The system shall be able to model transient diffusion through a composite slab with a constant concentration boundary condition as the species source, with the fine mesh and timestep required to match the analytical solution to generate CSV data for use in comparisons with the analytic solution for the profile concentration for the TMAP7 verification case.'
  heavy = true
  prereq = ver-1e_TMAP7_heavy
[]

[ver-1e_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1e.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1e, modeling transient diffusion through a composite slab with a constant concentration boundary condition as the species source for both the TMAP4 and TMAP7 verification cases.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1fa]
  type = Exodiff
  input = ver-1fa.i
  exodiff = ver-1fa_out.e
  requirement = 'The system shall be able to model heat conduction in a slab that has heat generation'
[]

[ver-1fa_lineplot]
  type = CSVDiff
  input = ver-1fa.i
  should_execute = False # this test relies on the output files from ver-1fa, so it shouldn't be run twice
  csvdiff = ver-1fa_csv_line_0011.csv
  requirement = 'The system shall be able to model heat conduction in a slab that has heat generation to generate CSV data for use in comparisons with the analytic solution for the profile concentration.'
  prereq = ver-1fa
[]

[ver-1fa_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1fa.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1fa, to model heat conduction in a slab that has heat generation.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1fb]
  type = Exodiff
  input = ver-1fb.i
  exodiff = ver-1fb_out.e
  requirement = 'The system shall be able to model thermal transient in a slab that has temperatures fixed at both the ends'
[]

[ver-1fb_csvdiff_0pt1sec]
  type = CSVDiff
  input = ver-1fb.i
  csvdiff = ver-1fb_u_vs_x_line_0010.csv
  should_execute = False # this test relies on the output files from ver-1fb, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal transient in a slab that has temperatures fixed at both the ends to generate CSV data at time of 0.1 s for use in comparison with analytical solution.'
  prereq = ver-1fb
[]

[ver-1fb_csvdiff_0pt5sec]
  type = CSVDiff
  input = ver-1fb.i
  csvdiff = ver-1fb_u_vs_x_line_0050.csv
  should_execute = False # this test relies on the output files from ver-1fb, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal transient in a slab that has temperatures fixed at both the ends to generate CSV data at time of 0.5 s for use in comparison with analytical solution.'
  prereq = ver-1fb
[]

[ver-1fb_csvdiff_1pt0sec]
  type = CSVDiff
  input = ver-1fb.i
  csvdiff = ver-1fb_u_vs_x_line_0100.csv
  should_execute = False # this test relies on the output files from ver-1fb, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal transient in a slab that has temperatures fixed at both the ends to generate CSV data at time of 1.0 s for use in comparison with analytical solution.'
  prereq = ver-1fb
[]

[ver-1fb_csvdiff_5pt0sec]
  type = CSVDiff
  input = ver-1fb.i
  csvdiff = ver-1fb_u_vs_x_line_0500.csv
  should_execute = False # this test relies on the output files from ver-1fb, so it shouldn't be run twice
  requirement = 'The system shall be able to model thermal transient in a slab that has temperatures fixed at both the ends to generate CSV data at time of 5.0 s for use in comparison with analytical solution.'
  prereq = ver-1fb
[]

[ver-1fb_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1fb.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1fb, modeling thermal transient in a slab with fixed temperatures at both the ends.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1fc]
  type = Exodiff
  input = ver-1fc.i
  exodiff = ver-1fc_out.e
  requirement = 'The system shall be able to model conduction in a composite structure with constant surface temperatures.'
[]

[ver-1fc_csvdiff_transient]
  type = CSVDiff
  input = ver-1fc.i
  should_execute = False # this test relies on the output files from ver-1fc, so it shouldn't be run twice
  csvdiff = ver-1fc_temperature_at_x0.09.csv
  requirement = 'The system shall be able to model conduction in a composite structure with constant surface temperatures to generate CSV data for use in comparisons with ABAQUS during transient at x=0.09 m.'
  prereq = ver-1fc
[]

[ver-1fc_csvdiff_transient_profile]
  type = CSVDiff
  input = ver-1fc.i
  should_execute = False # this test relies on the output files from ver-1fc, so it shouldn't be run twice
  csvdiff = ver-1fc_vector_postproc_line_0032.csv
  requirement = 'The system shall be able to model conduction in a composite structure with constant surface temperatures to generate CSV data for use in comparisons with ABAQUS during transient at t=150 s.'
  prereq = ver-1fc
[]

[ver-1fc_csvdiff_steady_state]
  type = CSVDiff
  input = ver-1fc.i
  should_execute = False # this test relies on the output files from ver-1fc, so it shouldn't be run twice
  csvdiff = ver-1fc_vector_postproc_line_0063.csv
  requirement = 'The system shall be able to model conduction in a composite structure with constant surface temperatures to generate CSV data for use in comparisons with ABAQUS and an analytical solution at steady state (t=10000 s).'
  prereq = ver-1fc
[]

[ver-1fc_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1fc.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution, ABAQUS data, and simulated solution of verification case 1fc, modeling conduction in a composite structure with constant surface temperatures.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1fd]
  type = Exodiff
  input = ver-1fd.i
  exodiff = ver-1fd_out.e
  requirement = 'The system shall be able to model convective heating.'
[]

[ver-1fd_csv]
  type = CSVDiff
  input = ver-1fd.i
  csvdiff = ver-1fd_out.csv
  should_execute = False # this test relies on the output files from ver-1fd, so it shouldn't be run twice
  requirement = 'The system shall be able to model convective heating to generate CSV data for use in comparisons with the analytic solution.'
  prereq = ver-1fd
[]

[ver-1fd_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1fd.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1fd, modeling convective heating.'
  required_python_packages = 'matplotlib scipy numpy pandas os'
[]

[binary_reaction_equal_concentrations]
  type = Exodiff
  input = 'ver-1g.i equal_conc.i'
  exodiff = equal_conc_out.e
  requirement = 'The system shall be able to model a chemical reaction between two species with the same concentrations and calculate the concentrations of reactants and product as a function of time'
[]

[binary_reaction_diff_concentrations_TMAP4]
  type = Exodiff
  input = 'ver-1g.i diff_conc_TMAP4.i'
  exodiff = diff_conc_TMAP4_out.e
  requirement = 'The system shall be able to model a chemical reaction between two species with different concentrations and calculate the concentrations of reactants and product as a function of time using the initial conditions from the TMAP4 case'
[]

[binary_reaction_diff_concentrations_TMAP7]
  type = Exodiff
  input = 'ver-1g.i diff_conc_TMAP7.i'
  exodiff = diff_conc_TMAP7_out.e
  requirement = 'The system shall be able to model a chemical reaction between two species with different concentrations and calculate the concentrations of reactants and product as a function of time using the initial conditions from the TMAP7 case'
[]

[binary_reaction_equal_concentrations_csv_diff]
  type = CSVDiff
  input = 'ver-1g.i equal_conc.i'
  should_execute = False # this test relies on the output files from binary_reaction_equal_concentrations, so it shouldn't be run twice
  csvdiff = equal_conc_out.csv
  requirement = 'The system shall be able to model a chemical reaction between two species with the same concentrations and calculate the concentrations of reactants and product as a function of time, to match the analytical solution to generate CSV data for use in comparisons with the analytic solution over time.'
  prereq = binary_reaction_equal_concentrations
[]

[binary_reaction_diff_concentrations_csv_diff_TMAP4]
  type = CSVDiff
  input = 'ver-1g.i diff_conc_TMAP4.i'
  should_execute = False # this test relies on the output files from binary_reaction_equal_concentrations, so it shouldn't be run twice
  csvdiff = diff_conc_TMAP4_out.csv
  requirement = 'The system shall be able to model a chemical reaction between two species with different concentrations using the initial conditions from the TMAP4 case and calculate the concentrations of reactants and product as a function of time to generate CSV data for use in comparisons with the analytic solution over time.'
  prereq = binary_reaction_diff_concentrations_TMAP4
[]

[binary_reaction_diff_concentrations_csv_diff_TMAP7]
  type = CSVDiff
  input = 'ver-1g.i diff_conc_TMAP7.i'
  should_execute = False # this test relies on the output files from binary_reaction_equal_concentrations, so it shouldn't be run twice
  csvdiff = diff_conc_TMAP7_out.csv
  requirement = 'The system shall be able to model a chemical reaction between two species with different concentrations using the initial conditions from the TMAP7 case and calculate the concentrations of reactants and product as a function of time to generate CSV data for use in comparisons with the analytic solution over time.'
  prereq = binary_reaction_diff_concentrations_TMAP7
[]

[ver-1g_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1g.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a chemical reaction between two species with same or different concentrations, using the initial conditions from both TMAP4 and TMAP7 cases.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1gc]
  type = Exodiff
  input = 'ver-1gc.i'
  exodiff = ver-1gc_out.e
  requirement = 'The system shall be able to model a series of chemical reactions involving three species and calculate the concentrations of each species as a function of time.'
[]

[ver-1gc_csv]
  type = CSVDiff
  input = 'ver-1gc.i'
  should_execute = False # this test relies on the output files from ver-1gc, so it shouldn't be run twice
  csvdiff = ver-1gc_out.csv
  requirement = 'The system shall be able to model a series of chemical reactions involving three species and calculate the concentrations of each species as a function of time and to generate CSV data for use in comparisons with the analytic solution over time.'
  prereq = ver-1gc
[]

[ver-1gc_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1gc.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a series of chemical reactions involving three species and calculate the concentrations of each species as a function of time and to generate CSV data for use in comparisons with the analytic solution over time.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1ha_csv]
  type = CSVDiff
  input = 'ver-1ha.i'
  csvdiff = ver-1ha_out.csv
  requirement = 'The system shall be able to model a convective outflow problem and calculate the pressure and concentration of the gas in the second and third enclosure and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1ha_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ha.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a convective outflow problem involving three enclosures and calculate the pressure and concentration of the gas in the second and third enclosure.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1hb_csv]
  type = CSVDiff
  input = 'ver-1hb.i'
  csvdiff = ver-1hb_out.csv
  requirement = 'The system shall be able to model a convective outflow problem and calculate the pressure and concentration of tritium and deuterium gas in the first and second enclosure and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1hb_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1hb.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a convective outflow problem involving two enclosures and two different gases and calculate the pressure and concentration of the gases in the enclosures.'
  prereq = 'ver-1hb_csv'
  required_python_packages = 'csv matplotlib numpy pandas scipy os math'
[]

[ver-1ia_csv]
  type = CSVDiff
  input = 'ver-1ia.i'
  csvdiff = ver-1ia_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface with equal starting pressures in ratedep condition and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1ia_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ia.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface with equal starting pressures in ratedep condition'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1ib_csv]
  type = CSVDiff
  input = '../ver-1ia/ver-1ia.i'
  cli_args = 'p0_B2="${units 1e5 Pa}" time_interval="${units 0.05 s}" Outputs/file_base="ver-1ib_out"' # these are the only things that need to be updated from ver-1ia to ver-1ib.
  csvdiff = ver-1ib_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface with unequal starting pressures in ratedep condition and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1ib_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ib.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface in ratedep condition with unequal starting pressures.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1ic_csv]
  type = CSVDiff
  input = 'ver-1ic.i'
  csvdiff = ver-1ic_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface in surfdep conditions with low barrier energy and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1ic_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ic.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface in surfdep condition with low barrier energy.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1id_csv]
  type = CSVDiff
  input = '../ver-1ic/ver-1ic.i'
  cli_args = "E_x='${units 0.2 eV -> J}' Outputs/file_base='ver-1id_out'" # these are the only things that need to be updated from ver-1ic to ver-1id.
  csvdiff = ver-1id_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface in surfdep conditions with high barrier energy and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1id_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1id.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface in surfdep condition with high barrier energy.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1ie_csv]
  type = CSVDiff
  input = 'ver-1ie.i'
  csvdiff = ver-1ie_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface in lawdep condition with equal starting pressures and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1ie_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ie.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface in lawdep condition with equal starting pressures.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1if_csv]
  type = CSVDiff
  input = '../ver-1ie/ver-1ie.i'
  cli_args = 'p0_B2="${units 1e5 Pa}" Outputs/file_base="ver-1if_out"' # these are the only things that need to be updated from ver-1ie to ver-1if.
  csvdiff = ver-1if_out.csv
  requirement = 'The system shall be able to model a equilibration problem on a reactive surface in lawdep condition with unequal starting pressures and to generate CSV data for use in comparisons with the analytic solution over time.'
[]

[ver-1if_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1if.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of a equilibration on a reactive surface in lawdep condition with unequal starting pressures.'
  required_python_packages = 'matplotlib numpy pandas scipy os'
[]

[ver-1ja_csvdiff]
  type = CSVDiff
  input = ver-1ja.i
  csvdiff = ver-1ja_out.csv
  requirement = 'The system shall be able to model decay of tritium and associated growth of helium in a diffusion segment and generate CSV data for use in comparisons with the analytic solution.'
[]

[ver-1ja_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ja.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1ja, which models decay of tritium and associated growth of helium in a diffusion segment.'
  required_python_packages = 'matplotlib numpy pandas os'
[]

[ver-1jb_csvdiff]
  type = CSVDiff
  input = ver-1jb.i
  csvdiff = ver-1jb_time_dependent_out.csv
  requirement = 'The system shall be able to model decay of tritium and associated growth of He in a diffusion segment with
                   distributed traps, with the fine mesh and timestep required to match the analytical solution to generate CSV
                   data for use in comparisons with the analytic solution.'
[]

[ver-1jb_csvdiff_profile]
  type = CSVDiff
  input = ver-1jb.i
  csvdiff = ver-1jb_profile_out_line_0048.csv
  should_execute = false # uses ver-1jb_csvdiff
  prereq = ver-1jb_csvdiff
  requirement = 'The system shall be able to model decay of tritium and associated growth of He in a diffusion segment with
                   distributed traps and output the profiles of concentrations.'
[]

[ver-1jb_csvdiff_equivalent_concentrations]
  type = CSVDiff
  input = ver-1jb.i
  cli_args = 'tritium_mobile_concentration_initial=1e25
                trap_per_free=1
                Outputs/file_base=ver-1jb_equivalent_concentrations_out
                Outputs/time_dependent_out/file_base=ver-1jb_equivalent_concentrations_time_dependent_out
                Outputs/profile_out/file_base=ver-1jb_equivalent_concentrations_profile_out'
  csvdiff = ver-1jb_equivalent_concentrations_time_dependent_out.csv
  requirement = 'The system shall be able to model decay of tritium and associated growth of He in a diffusion segment with
                   distributed traps with equivalent initial mobile and trapped tritium concentration,
                   with the fine mesh and timestep required to match the analytical solution to generate CSV data for use in
                   comparisons with the analytic solution.'
[]

[ver-1jb_csvdiff_profile_equivalent_concentrations]
  type = CSVDiff
  input = ver-1jb.i
  cli_args = 'tritium_mobile_concentration_initial=1e25
                trap_per_free=1
                Outputs/file_base=ver-1jb_equivalent_concentrations_out
                Outputs/time_dependent_out/file_base=ver-1jb_equivalent_concentrations_time_dependent_out
                Outputs/profile_out/file_base=ver-1jb_equivalent_concentrations_profile_out'
  csvdiff = ver-1jb_equivalent_concentrations_profile_out_line_0048.csv
  should_execute = false # uses ver-1jb_csvdiff_equivalent_concentrations
  prereq = ver-1jb_csvdiff_equivalent_concentrations
  requirement = 'The system shall be able to model decay of tritium and associated growth of He in a diffusion segment with
                   distributed traps with equivalent initial mobile and trapped tritium concentration and output the
                   profiles of concentrations.'
[]

[ver-1jb_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1jb.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution
                   when modeling decay of tritium and associated growth of He in a diffusion segment with distributed traps.'
  required_python_packages = 'csv matplotlib numpy pandas os'
[]

[ver-1ka_csv]
  type = CSVDiff
  input = ver-1ka.i
  csvdiff = ver-1ka_out.csv
  requirement = 'The system shall be able to model a tritium volumetric source in one enclosure'
[]

[ver-1ka_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1ka.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1ka, modeling a tritium volumetric source in one enclosure.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[ver-1kb_csv_without_concentration_jump]
  type = CSVDiff
  input = ver-1kb.i
  csvdiff = ver-1kb_out_k1.csv
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Henry’s law without any concentration jump at the interface.'
[]

[ver-1kb_csv_concentration_jump]
  type = CSVDiff
  input = ver-1kb.i
  cli_args = "solubility='${fparse 10/(R*temperature)}'
                Outputs/file_base=ver-1kb_out_k10"
  csvdiff = ver-1kb_out_k10.csv
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Henry’s law with a concentration jump at the interface.'
[]

[ver-1kb_exodus_concentration_jump]
  type = Exodiff
  input = ver-1kb.i
  exodiff = ver-1kb_out_k10.e
  cli_args = "solubility='${fparse 10/(R*temperature)}'
                Outputs/file_base=ver-1kb_out_k10"
  prereq = ver-1kb_csv_concentration_jump
  should_execute = false # this test relies on the output files from ver-1kb_csv_concentration_jump, so it shouldn't be run twice
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Henry’s law with a concentration jump at the interface'
[]

[ver-1kb_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1kb.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1kb, modeling a diffusion across a membrane separating two enclosures in accordance with Henry’s law.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[ver-1kc-1_csv]
  type = CSVDiff
  input = ver-1kc-1.i
  cli_args = "nb_segments_TMAP8=1e1
                simulation_time=4
                Executioner/nl_abs_tol=1e-5
                Executioner/nl_rel_tol=1e-4
                Outputs/exodus=false
                Outputs/file_base=ver-1kc-1_out_k10_light"
  csvdiff = ver-1kc-1_out_k10_light.csv
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface.'
[]

[ver-1kc-1_csv_heavy]
  type = CSVDiff
  heavy = true
  input = ver-1kc-1.i
  csvdiff = ver-1kc-1_out_k10.csv
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface with a fine mesh and tight tolerances for higher accuracy.'
[]

[ver-1kc-1_exodus_heavy]
  type = Exodiff
  heavy = true
  input = ver-1kc-1.i
  exodiff = ver-1kc-1_out_k10.e
  prereq = ver-1kc-1_csv_heavy
  should_execute = false # this test relies on the output files from ver-1kc-1_csv_concentration_jump, so it shouldn't be run twice
  requirement = 'The system shall be able to model the diffusion of T2 across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface with a fine mesh and tight tolerances for higher accuracy and generate an exodus file.'
[]

[ver-1kc-1_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1kc-1.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1kc-1, modeling a diffusion across a membrane separating two enclosures in accordance with Sieverts’ law.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[ver-1kc-2_csv]
  type = CSVDiff
  input = ver-1kc-2.i
  cli_args = "simulation_time=0.05
                Executioner/nl_rel_tol=1e-6
                Outputs/exodus=false
                Outputs/file_base=ver-1kc-2_out_k10_light"
  csvdiff = ver-1kc-2_out_k10_light.csv
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface.'
[]

[ver-1kc-2_csv_heavy]
  type = CSVDiff
  heavy = true
  input = ver-1kc-2.i
  csvdiff = ver-1kc-2_out_k10.csv
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface with tight tolerances for higher accuracy.'
[]

[ver-1kc-2_exodus_heavy]
  type = Exodiff
  heavy = true
  input = ver-1kc-2.i
  exodiff = ver-1kc-2_out_k10.e
  prereq = ver-1kc-2_csv_heavy
  should_execute = false # this test relies on the output files from ver-1kc-2_csv_heavy, so it shouldn't be run twice
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface and generate an exodus file with tight tolerances for higher accuracy.'
[]

[ver-1kc-2_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1kc-2.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1kc-2, modeling a diffusion across a membrane separating two enclosures in accordance with Sieverts’ law.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[ver-1kd_csv]
  type = CSVDiff
  input = ver-1kd.i
  cli_args = "simulation_time=0.05
                Executioner/nl_rel_tol=1e-6
                Outputs/exodus=false
                Outputs/file_base=ver-1kd_out_k10_light"
  csvdiff = ver-1kd_out_k10_light.csv
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface and a T2 volumetric source term.'
[]

[ver-1kd_csv_heavy]
  type = CSVDiff
  heavy = true
  input = ver-1kd.i
  csvdiff = ver-1kd_out_k10.csv
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface and a T2 volumetric source term with tight tolerances for higher accuracy.'
[]

[ver-1kd_exodus_heavy]
  type = Exodiff
  heavy = true
  input = ver-1kd.i
  exodiff = ver-1kd_out_k10.e
  prereq = ver-1kd_csv_heavy
  should_execute = false # this test relies on the output files from ver-1kd_csv_heavy, so it shouldn't be run twice
  requirement = 'The system shall be able to model the diffusion of T2, H2 and HT across a membrane separating two enclosures in accordance with Sieverts’ law with a concentration jump at the interface and a T2 volumetric source term and generate an exodus file with tight tolerances for higher accuracy.'
[]

[ver-1kd_comparison]
  type = RunCommand
  command = 'python3 comparison_ver-1kd.py'
  requirement = 'The system shall be able to generate comparison plots between the analytical solution and simulated solution of verification case 1kd, modeling a diffusion across a membrane separating two enclosures in accordance with Sieverts’ law  and a T2 volumetric source term.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[YHx_PCT_csv]
  type = CSVDiff
  input = YHx_PCT.i
  cli_args = 'InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  csvdiff = YHx_PCT_out.csv
  max_parallel = 1 # see #200
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature.'
[]

[YHx_PCT_exodus]
  type = Exodiff
  input = YHx_PCT.i
  cli_args = 'InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  exodiff = YHx_PCT_out.e
  prereq = YHx_PCT_csv
  should_execute = false # this test relies on the output files from YHx_PCT_csv, so it shouldn't be run twice
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature and generate an exodus file.'
[]

[YHx_PCT_T1173_P1e3_csv]
  type = CSVDiff
  input = YHx_PCT.i
  cli_args = 'output_file_base=YHx_PCT_T1173_P1e3_out
                dt_init=1e1
                temperature=1173.15
                initial_pressure_H2_enclosure_1=1e3
                initial_atomic_fraction=1.55
                Executioner/num_steps=150
                InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  csvdiff = YHx_PCT_T1173_P1e3_out.csv
  max_parallel = 1 # see #200
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature for P=1e3 Pa and T=1173.15 K.'
[]

[YHx_PCT_T1173_P1e4_csv]
  type = CSVDiff
  input = YHx_PCT.i
  cli_args = 'output_file_base=YHx_PCT_T1173_P1e4_out
                temperature=1173.15
                initial_pressure_H2_enclosure_1=1e4
                initial_atomic_fraction=1.65
                Executioner/num_steps=200
                InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  csvdiff = YHx_PCT_T1173_P1e4_out.csv
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature for P=1e4 Pa and T=1173.15 K.'
[]

[YHx_PCT_T1173_P5e4_csv]
  type = CSVDiff
  input = YHx_PCT.i
  cli_args = 'output_file_base=YHx_PCT_T1173_P5e4_out
                temperature=1173.15
                initial_pressure_H2_enclosure_1=5e4
                initial_atomic_fraction=1.85
                Executioner/num_steps=200
                InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  csvdiff = YHx_PCT_T1173_P5e4_out.csv
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature for P=5e4 Pa and T=1173.15 K.'
[]

[YHx_PCT_T1273_P3e3_csv]
  type = CSVDiff
  input = YHx_PCT.i
  cli_args = 'output_file_base=YHx_PCT_T1273_P3e3_out
                temperature=1273.15
                initial_pressure_H2_enclosure_1=3e3
                initial_atomic_fraction=1.41
                Executioner/num_steps=200
                InterfaceKernels/interface_reaction_YHx_PCT/silence_warnings=true' # should be removed once PCT curve implementation for low pressure is complete
  csvdiff = YHx_PCT_T1273_P3e3_out.csv
  max_parallel = 1 # see #200
  requirement = 'The system shall be able to model the PCT curves of YHx to determine the surface atomic fraction as a function of pressure and temperature for P=3e3 Pa and T=1273.15 K.'
[]

[YHx_PCT_comparison]
  type = RunCommand
  command = 'python3 comparison_YHx_PCT.py'
  requirement = 'The system shall be able to generate comparison plots between experimental PCT curves, the model used in TMAP8, and TMAP8 predictions.'
  required_python_packages = 'matplotlib numpy pandas os git'
[]

[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).'
[]