MoabSkinner

Re-generate the OpenMC geometry on-the-fly according to changes in the mesh geometry and/or contours in temperature and density

Description

This user object bins the spatial domain according to two characteristics:

any number of user-specified scalar fields (e.g., temperature, density), by grouping elements into uniform bins with equal spacing between each field's fields_min and fields_max bounds, as set in n_field_bins
subdomain

Fields are binned in the order they are listed in fields. Each field increases the total number of bins; the combined bin index of an element accounts for its subdomain bin and its bin in every field.

For example, suppose the domain consists of two subdomains, with temperature and density distributions as shown in the top row of Figure 1. The elements are then grouped into individualized bins, according to unique combinations of the temperature bin, the density bin, and the subdomain bin. In the lower right of Figure 1 is then shown the surface skins produced internally, which can be used to define the boundaries of new regions on which to track particles.

Figure 1: Mesh skinning operation

tip

You can visualize this user object (i.e., the bin indices) using a SpatialUserObjectAux.

Supported Mesh Element Types

MoabSkinner natively supports meshes consisting of TET4 and TET10 elements, since DAGMC geometry requires triangulated surfaces (i.e. skins of tetrahedra). Meshes containing HEX, PYRAMID, PRISM elements are also supported: when these elements are detected, an internal all-TET4 copy of the mesh is built automatically, and all skinning and geometry operations are performed on that copy. The original MOOSE mesh is never modified; field and subdomain binning are always evaluated on the original mesh where the auxiliary variables live. Higher order versions of these element types are also supported.

Example Input Syntax

Below is an example input file that skins a mesh, generating the bin distributions shown in Figure 1.

[UserObjects<<<{"href": "../../syntax/UserObjects/index.html"}>>>]
  [moab]
    type = MoabSkinner<<<{"description": "Re-generate the OpenMC geometry on-the-fly according to changes in the mesh geometry and/or contours in temperature and density", "href": "MoabSkinner.html"}>>>
    fields<<<{"description": "Variable(s) to bin elements by, e.g. temperature, density, etc."}>>> = 'temp rho'
    fields_min<<<{"description": "Lower bounds of the bins for each field"}>>> = '445 445'
    fields_max<<<{"description": "Upper bounds of the bins for each field"}>>> = '655 655'
    n_field_bins<<<{"description": "Number of bins for each field"}>>> = '3 5'
    verbose<<<{"description": "Whether to print diagnostic information"}>>> = true
    material_blocks<<<{"description": "List of mesh subdomain names (or IDs) for which to assign material names in the generated DAGMC geometry. Must be provided together with 'material_names' where both have the same length and are listed in the same order.Any subdomain not listed in 'material_blocks' will have no material assignment and OpenMC willdefault to void for that region. This parameter is optional when used in combination with OpenMCCellAverageProblembecause in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning."}>>> = "1 3"
    material_names<<<{"description": "Material names (or IDs) to assign to subdomains in the generated DAGMC geometry. Must be provided together with 'material_blocks' where both have the same lengthand listed in the same order. This parameter is optional when used in combination with OpenMCCellAverageProblem because in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning."}>>> = "mat mat"

    output_skins<<<{"description": "Whether the skinned MOAB mesh (skins generated from the libMesh [Mesh]) should be written to a file. The files will be named moab_skins_<n>.h5m, where <n> is the time step index. You can then visualize these files by running 'mbconvert'."}>>> = true
  []
[]

The binning, colored on a per-element basis as shown in Figure 1 can be visualized with the SkinnedBins auxiliary kernel.

Material Names and Blocks

For the vast majority of use cases, you will be using the MoabSkinner in tandem with OpenMCCellAverageProblem. In this case, you can rely on OpenMC to auto-detect the mapping between materials and subdomains from the loaded .h5m model, so that when the geometry is re-skinned, OpenMC will re-apply the original materials corresponding to the regions of space. However, if you'd like to instead force the use of a _different_ material from what was present in the original .h5m file, you can provide material_blocks and material_names to have the skinner apply specific materials on each block.

The MoabSkinner can also be used independently of OpenMCCellAverageProblem to generate an initial DAGMC geometry (.h5m file) directly from a MOOSE mesh, without needing any pre-existing DAGMC model. In this mode, material_blocks and material_names are required in order to tag the generated volumes with OpenMC material assignments, and output_full should be set to true to write the generated geometry to a file, then point OpenMC's geometry.xml at the generated file for subsequent coupled runs.

Both mesh subdomains (in material_blocks) and materials (in material_names) may be specified either by name or by numeric ID. The two parameters must be the same length and are matched in order. Any mesh subdomain not listed in material_blocks receives no material assignment, and the corresponding region will be void in OpenMC; a warning listing these subdomains is printed at setup.

Below is an example of generating an initial DAGMC geometry directly from a MOOSE mesh

[UserObjects<<<{"href": "../../syntax/UserObjects/index.html"}>>>]
  [moab]
    type = MoabSkinner<<<{"description": "Re-generate the OpenMC geometry on-the-fly according to changes in the mesh geometry and/or contours in temperature and density", "href": "MoabSkinner.html"}>>>
    fields<<<{"description": "Variable(s) to bin elements by, e.g. temperature, density, etc."}>>> = 'temp'
    fields_min<<<{"description": "Lower bounds of the bins for each field"}>>> = '0'
    fields_max<<<{"description": "Upper bounds of the bins for each field"}>>> = '1000'
    n_field_bins<<<{"description": "Number of bins for each field"}>>> = '1'
    build_graveyard<<<{"description": "Whether to build a graveyard around the geometry"}>>> = true
    graveyard_scale_inner<<<{"description": "Multiplier on mesh bounding box to form inner graveyard surface"}>>> = 1.1
    graveyard_scale_outer<<<{"description": "Multiplier on mesh bounding box to form outer graveyard surface"}>>> = 1.2
    material_blocks<<<{"description": "List of mesh subdomain names (or IDs) for which to assign material names in the generated DAGMC geometry. Must be provided together with 'material_names' where both have the same length and are listed in the same order.Any subdomain not listed in 'material_blocks' will have no material assignment and OpenMC willdefault to void for that region. This parameter is optional when used in combination with OpenMCCellAverageProblembecause in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning."}>>> = 'block_1 block_2'
    material_names<<<{"description": "Material names (or IDs) to assign to subdomains in the generated DAGMC geometry. Must be provided together with 'material_blocks' where both have the same lengthand listed in the same order. This parameter is optional when used in combination with OpenMCCellAverageProblem because in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning."}>>> = '1 U35_U38'
    output_full<<<{"description": "Whether the MOAB mesh (copied from the libMesh [Mesh]) should be written to a file. The files will be named moab_full_<n>.h5m, where <n> is the time step index. You can then visualize these files by running 'mbconvert'."}>>> = true
  []
[]

Input Parameters

build_graveyardFalseWhether to build a graveyard around the geometry
Default:False
C++ Type:bool
Controllable:No
Description:Whether to build a graveyard around the geometry
faceting_tol0.0001Faceting tolerance for DagMC
Default:0.0001
C++ Type:Real
Unit:(no unit assumed)
Range:faceting_tol > 0
Controllable:No
Description:Faceting tolerance for DagMC
fieldsVariable(s) to bin elements by, e.g. temperature, density, etc.
C++ Type:std::vector<std::string>
Controllable:No
Description:Variable(s) to bin elements by, e.g. temperature, density, etc.
fields_maxUpper bounds of the bins for each field
C++ Type:std::vector<Real>
Unit:(no unit assumed)
Controllable:No
Description:Upper bounds of the bins for each field
fields_minLower bounds of the bins for each field
C++ Type:std::vector<Real>
Unit:(no unit assumed)
Controllable:No
Description:Lower bounds of the bins for each field
geom_tol1e-06Geometry tolerance for DagMC
Default:1e-06
C++ Type:Real
Unit:(no unit assumed)
Range:geom_tol > 0
Controllable:No
Description:Geometry tolerance for DagMC
graveyard_scale_inner1.01Multiplier on mesh bounding box to form inner graveyard surface
Default:1.01
C++ Type:Real
Unit:(no unit assumed)
Range:graveyard_scale_inner > 1
Controllable:No
Description:Multiplier on mesh bounding box to form inner graveyard surface
graveyard_scale_outer1.1Multiplier on mesh bounding box to form outer graveyard surface
Default:1.1
C++ Type:Real
Unit:(no unit assumed)
Controllable:No
Description:Multiplier on mesh bounding box to form outer graveyard surface
implicit_complement_materialAssigns OpenMC material name or ID to the implicit complement region. If not provided, void material is assigned by default.
C++ Type:std::string
Controllable:No
Description:Assigns OpenMC material name or ID to the implicit complement region. If not provided, void material is assigned by default.
material_blocksList of mesh subdomain names (or IDs) for which to assign material names in the generated DAGMC geometry. Must be provided together with 'material_names' where both have the same length and are listed in the same order.Any subdomain not listed in 'material_blocks' will have no material assignment and OpenMC willdefault to void for that region. This parameter is optional when used in combination with OpenMCCellAverageProblembecause in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning.
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:List of mesh subdomain names (or IDs) for which to assign material names in the generated DAGMC geometry. Must be provided together with 'material_names' where both have the same length and are listed in the same order.Any subdomain not listed in 'material_blocks' will have no material assignment and OpenMC willdefault to void for that region. This parameter is optional when used in combination with OpenMCCellAverageProblembecause in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning.
material_namesMaterial names (or IDs) to assign to subdomains in the generated DAGMC geometry. Must be provided together with 'material_blocks' where both have the same lengthand listed in the same order. This parameter is optional when used in combination with OpenMCCellAverageProblem because in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning.
C++ Type:std::vector<std::string>
Controllable:No
Description:Material names (or IDs) to assign to subdomains in the generated DAGMC geometry. Must be provided together with 'material_blocks' where both have the same lengthand listed in the same order. This parameter is optional when used in combination with OpenMCCellAverageProblem because in this case, the default is to auto-detect the materials in the original h5m file and re-apply them during skinning.
n_field_binsNumber of bins for each field
C++ Type:std::vector<unsigned int>
Controllable:No
Description:Number of bins for each field
output_fullFalseWhether the MOAB mesh (copied from the libMesh [Mesh]) should be written to a file. The files will be named moab_full_.h5m, where is the time step index. You can then visualize these files by running 'mbconvert'.
Default:False
C++ Type:bool
Controllable:No
Description:Whether the MOAB mesh (copied from the libMesh [Mesh]) should be written to a file. The files will be named moab_full_.h5m, where is the time step index. You can then visualize these files by running 'mbconvert'.
output_skinsFalseWhether the skinned MOAB mesh (skins generated from the libMesh [Mesh]) should be written to a file. The files will be named moab_skins_.h5m, where is the time step index. You can then visualize these files by running 'mbconvert'.
Default:False
C++ Type:bool
Controllable:No
Description:Whether the skinned MOAB mesh (skins generated from the libMesh [Mesh]) should be written to a file. The files will be named moab_skins_.h5m, where is the time step index. You can then visualize these files by running 'mbconvert'.
reflective_bcs_surfacesMesh sideset names or numeric sideset IDs to assign DAGMC reflective boundary conditions to. Both string names and integer IDs are accepted. If not specified, the surface defaults to transmission.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:Mesh sideset names or numeric sideset IDs to assign DAGMC reflective boundary conditions to. Both string names and integer IDs are accepted. If not specified, the surface defaults to transmission.
vacuum_bcs_surfacesMesh sideset names or numeric sideset IDs to assign DAGMC vacuum boundary conditions to. Both string names and integer IDs are accepted. If not specified, the surface defaults to transmission.
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:Mesh sideset names or numeric sideset IDs to assign DAGMC vacuum boundary conditions to. Both string names and integer IDs are accepted. If not specified, the surface defaults to transmission.
verboseFalseWhether to print diagnostic information
Default:False
C++ Type:bool
Controllable:No
Description:Whether to print diagnostic information

Optional Parameters

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

Execution Scheduling Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
use_displaced_meshFalseWhether the skinned mesh should be generated from a displaced mesh
Default:False
C++ Type:bool
Controllable:No
Description:Whether the skinned mesh should be generated from a displaced mesh

Advanced Parameters

prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

(test/tests/neutronics/dagmc/missing_graveyard.i)
(test/tests/userobjects/moab_skinner/wrong_type.i)
(test/tests/neutronics/dagmc/implicit_comp_material/without_graveyard/openmc.i)
(test/tests/neutronics/dagmc/add_bcs/openmc.i)
(test/tests/neutronics/dagmc/density_skin/openmc.i)
(test/tests/neutronics/dagmc/cell_tallies/unmapped_cells.i)
(test/tests/userobjects/moab_skinner/hex-mesh/sub_bins.i)
(test/tests/neutronics/dagmc/density_skin/only_solid.i)
(test/tests/userobjects/moab_skinner/expanding-cube/all_bins.i)
(test/tests/neutronics/dagmc/incompatible_geom/mismatch/openmc.i)
(test/tests/neutronics/dagmc/no_symmetry.i)
(test/tests/neutronics/dagmc/mesh_tallies/direct_tally.i)
(test/tests/neutronics/dagmc/with_csg/dag_cell_not_root/openmc.i)
(test/tests/neutronics/dagmc/mesh_tallies/with_relaxation.i)
(test/tests/neutronics/dagmc/cell_tallies/one_bin.i)
(tutorials/csg_to_cad/openmc.i)
(test/tests/userobjects/moab_skinner/all_bins.i)
(test/tests/userobjects/moab_skinner/graveyard/graveyard.i)
(tutorials/tokamak/openmc.i)
(test/tests/neutronics/dagmc/cell_tallies/with_relaxation.i)
(test/tests/neutronics/mg/rr_thermal_expansion/openmc.i)
(tutorials/msfr/openmc.i)
(test/tests/neutronics/kinetics/geom_reset/one_bin.i)
(test/tests/userobjects/moab_skinner/errors/skinner.i)
(test/tests/neutronics/dagmc/with_csg/dag_lattice_outer/openmc.i)
(test/tests/userobjects/moab_skinner/hex-mesh/mixed_elements.i)
(test/tests/neutronics/dagmc/properties.i)
(doc/content/vv/slab_benchmark_model/openmc.i)
(test/tests/neutronics/dagmc/with_csg/multi_cell_dag/openmc.i)
(test/tests/userobjects/moab_skinner/hex-mesh/all_bins.i)
(test/tests/neutronics/dagmc/with_csg/openmc.i)
(test/tests/neutronics/dagmc/density_skin/disjoint_bins.i)
(test/tests/neutronics/dagmc/mesh_tallies/openmc.i)
(test/tests/neutronics/moving_mesh/scale.i)
(test/tests/userobjects/moab_skinner/output/mesh.i)
(test/tests/neutronics/dagmc/mesh_tallies/disjoint_bins.i)
(test/tests/neutronics/dagmc/with_csg/dag_lattice/openmc.i)
(test/tests/userobjects/moab_skinner/three_fields.i)
(test/tests/neutronics/dagmc/cell_tallies/with_scaling.i)
(test/tests/openmc_errors/input_params/no_dag.i)
(test/tests/neutronics/dagmc/mesh_tallies/one_bin.i)
(test/tests/userobjects/moab_skinner/sub_bins.i)
(test/tests/neutronics/dagmc/assign_materials/generate_dagmc.i)
(test/tests/userobjects/moab_skinner/second_order/all_bins.i)
(test/tests/neutronics/moving_mesh/volumes.i)
(test/tests/neutronics/dagmc/assign_materials/openmc.i)
(test/tests/neutronics/dagmc/cell_tallies/null_density.i)
(test/tests/neutronics/dagmc/assign_materials/no_material_names.i)
(test/tests/neutronics/dagmc/cell_tallies/openmc.i)
(test/tests/neutronics/dagmc/with_csg/multiple_dag/openmc.i)
(test/tests/neutronics/dagmc/implicit_comp_material/openmc.i)
(test/tests/neutronics/dagmc/mesh_tallies/with_scaling.i)
(test/tests/neutronics/dagmc/add_bcs/openmc_hex.i)
(test/tests/neutronics/dagmc/with_csg/mixed_feedback/openmc.i)

(test/tests/userobjects/moab_skinner/all_bins.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '445 445'
    fields_max = '655 655'
    n_field_bins = '3 5'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"



    output_skins = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(test/tests/neutronics/dagmc/assign_materials/generate_dagmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  [block_1]
    type = ParsedSubdomainMeshGenerator
    input = slab
    combinatorial_geometry = 'x < 37.5'
    block_id = 1
    block_name = block_1
    epsilon = 1e-3
  []
  [block_2]
    type = ParsedSubdomainMeshGenerator
    input = block_1
    combinatorial_geometry = 'x >= 37.5'
    block_id = 2
    block_name = block_2
    epsilon = 1e-3
  []
  allow_renumbering = false
[]

[AuxVariables]
  [temp]
    order = CONSTANT
    family = MONOMIAL
    initial_condition = 200
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
    graveyard_scale_inner = 1.1
    graveyard_scale_outer = 1.2
    material_blocks = 'block_1 block_2'
    material_names = '1 U35_U38'
    output_full = true
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/neutronics/dagmc/missing_graveyard.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = ../meshes/tet_cube.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[Problem]
  type = OpenMCCellAverageProblem
  cell_level = 0
  temperature_blocks = '1'
  power = 1000.0
  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '900'
    n_field_bins = '1'
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/userobjects/moab_skinner/wrong_type.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '445'
    fields_max = '655'
    n_field_bins = '3'
    material_names = "mat mat"

  []
[]

[Executioner]
  type = Transient
[]

(test/tests/neutronics/dagmc/implicit_comp_material/without_graveyard/openmc.i)

[Mesh]
  type = FileMesh
  file = ../../mesh_tallies/slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = false
    implicit_complement_material="h2o"
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
[]

(test/tests/neutronics/dagmc/add_bcs/openmc.i)

[Mesh]
  [f]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -12.5
    xmax = 87.5
    ymin = -12.5
    ymax = 37.5
    zmin = -12.5
    zmax = 12.5
    nx = 8
    ny = 4
    nz = 2
  []
  [subdomains]
    type = ParsedSubdomainIDsGenerator
    input = f
    expression = 'if (y > 12.5, 1, 2)'
  []
  [name]
    type = RenameBlockGenerator
    input = subdomains
    old_block = '1 2'
    new_block = 'block_1 block_2'
  []
  # The following sidesets/nodesets are only used by the error-checking and internal-surface tests
  [left_dup] # duplicate of the 'left' sideset, to test one side belonging to two sidesets
    type = ParsedGenerateSideset
    input = name
    combinatorial_geometry = 'x < -12.49'
    normal = '-1 0 0'
    new_sideset_name = 'left_dup'
  []
  [interface_a] # the block_1/block_2 interface, attached to block_1 element sides
    type = SideSetsBetweenSubdomainsGenerator
    input = left_dup
    primary_block = 'block_1'
    paired_block = 'block_2'
    new_boundary = 'interface_a'
  []
  [interface_b] # the same interface, attached to block_2 element sides
    type = SideSetsBetweenSubdomainsGenerator
    input = interface_a
    primary_block = 'block_2'
    paired_block = 'block_1'
    new_boundary = 'interface_b'
  []
  [midplane] # a plane interior to block_2 (y = 0), which never appears in any skin
    type = ParsedGenerateSideset
    input = interface_b
    combinatorial_geometry = 'y > -0.01 & y < 0.01'
    normal = '0 1 0'
    new_sideset_name = 'midplane'
  []
  [corner_nodes] # a nodeset-only boundary, to test the nodeset error message
    type = ExtraNodesetGenerator
    input = midplane
    new_boundary = 'corner_nodes'
    coord = '-12.5 -12.5 -12.5'
  []
  [hextotet]
    type = ElementsToTetrahedronsConverter
    input = corner_nodes
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

[]

[UserObjects]
  [moab]
    # BCs are added in tests as cli_args to use this file in all different cases
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    verbose = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

[Outputs]
  exodus = false
  [file]
    type = CSV
  []
[]

(test/tests/neutronics/dagmc/density_skin/openmc.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

[AuxVariables]
  [material_id]
    family = MONOMIAL
    order = CONSTANT
  []
  [cell_density]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [material_id]
    type = CellMaterialIDAux
    variable = material_id
  []
  [cell_density]
    type = CellDensityAux
    variable = cell_density
  []
[]

rho0 = 600.0
drho = 50.0

[Functions]
  [density]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse rho0 - drho}, if (x <= ${x1}, ${rho0}, if (x <= ${x2}, ${fparse rho0 + drho}, ${fparse rho0 + 2 * drho})))'
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 500.0
  []
[]

[AuxKernels]
  [density]
    type = FunctionAux
    variable = density
    function = density
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  density_blocks = '1 2'
  power = 100.0

  skinner = moab
  cell_level = 0

  [Tallies]
    [Mesh]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '0 ${fparse rho0 - drho}'
    fields_max = '1000 ${fparse rho0 + 2 * drho}'
    n_field_bins = '1 4'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp density cell_instance cell_id'
[]

(test/tests/neutronics/dagmc/cell_tallies/unmapped_cells.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  [delete]
    type = BlockDeletionGenerator
    input = file
    block = '2'
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[Problem]
  type = OpenMCCellAverageProblem
  cell_level = 0
  temperature_blocks = '1'
  skinner = moab
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '900'
    n_field_bins = '1'
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/userobjects/moab_skinner/hex-mesh/sub_bins.i)

[Mesh]
  [cube]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -0.5
    xmax = 0.5
    ymin = -0.5
    ymax = 0.5
    zmin = -0.5
    zmax = 0.5
    nx = 5
    ny = 5
    nz = 5
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '445'
    fields_max = '655'
    n_field_bins = '3'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"

  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/neutronics/dagmc/density_skin/only_solid.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

[AuxVariables]
  [density]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  initial_properties = xml
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '0 0'
    fields_max = '1000 100'
    n_field_bins = '1 4'
    build_graveyard = true
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/userobjects/moab_skinner/expanding-cube/all_bins.i)

[Mesh]
 [cube]
    type = FileMeshGenerator
    file = '../../../neutronics/meshes/tet_cube.e'
 []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
  use_displaced_mesh = true
[]

[Problem]
  type = FEProblem
  solve = false
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []
  [disp_z]
  []
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 1000+300*x*x+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 1000+300*y*y+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '996 996'
    fields_max = '1070 1070'
    n_field_bins = '10 10'
    verbose = true
    material_blocks = "1"
    material_names = "mat"



    output_skins = true
  []
[]

[Executioner]
  type = Transient
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(test/tests/neutronics/dagmc/incompatible_geom/mismatch/openmc.i)

[Mesh]
  [load]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []
  [merge]
    type = RenameBlockGenerator
    input = load
    old_block = '2'
    new_block = '1'
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[Problem]
  type = OpenMCCellAverageProblem

  temperature_blocks = '1'
  cell_level = 0
  power = 100.0

  initial_properties = xml
  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = ../../mesh_tallies/slab.e
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '4'
    build_graveyard = true
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/neutronics/dagmc/no_symmetry.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = ../meshes/tet_cube.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[Problem]
  type = OpenMCCellAverageProblem
  cell_level = 0
  temperature_blocks = '1'
  power = 1000.0
  skinner = moab
  symmetry_mapper = sym

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = ../meshes/tet_cube.e
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '900'
    n_field_bins = '1'
    build_graveyard = true
  []
  [sym]
    type = SymmetryPointGenerator
    normal = '1.0 0.0 0.0'
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/neutronics/dagmc/mesh_tallies/direct_tally.i)

!include csg_step_2_comp.i

[Mesh]
  type = FileMesh
  file = slab.e
  allow_renumbering = false
  parallel_type = replicated
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true

  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      output = 'unrelaxed_tally_std_dev'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  execute_on = 'FINAL'
[]

(test/tests/neutronics/dagmc/with_csg/dag_cell_not_root/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/neutronics/dagmc/mesh_tallies/with_relaxation.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = slab.e
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = '500.0 + 10.0*x'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true

  temperature_blocks = '1 2'
  cell_level = 0
  power = 16.0

  skinner = moab

  relaxation = 'constant'

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = slab.e
      output = 'unrelaxed_tally'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '300'
    fields_max = '1500'
    n_field_bins = '10'
    build_graveyard = true
  []
[]

[Postprocessors]
  [heating]
    type = PointValue
    variable = kappa_fission
    point = '47 15 5'
  []
  [heating_raw]
    type = PointValue
    variable = kappa_fission_raw
    point = '47 15 5'
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  csv = true
  execute_on = 'TIMESTEP_END'
[]

(test/tests/neutronics/dagmc/cell_tallies/one_bin.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp cell_instance'
[]

(tutorials/csg_to_cad/openmc.i)

power = 3000

[Problem]
  type = OpenMCCellAverageProblem
  power = ${power}

  temperature_blocks = '1'
  cell_level = 1

  skinner = skinner

  source_rate_normalization = 'kappa_fission'

  [Tallies]
    [mesh]
      type = MeshTally
      score = 'kappa_fission flux'

      normalize_by_global_tally = false
    []
  []
[]

[Mesh]
  [fuel]
    type = FileMeshGenerator
    file = krusty_fuel.e
  []
[]

[UserObjects]
  [skinner]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '800'
    fields_max = '1000'
    n_field_bins = '4'
    build_graveyard = true
  []
[]

[MultiApps]
  [conduction]
    type = TransientMultiApp
    input_files = 'fuel.i'
    execute_on = timestep_end
  []
[]

[Transfers]
  [heat_source_from_openmc]
    type = MultiAppGeneralFieldShapeEvaluationTransfer
    to_multi_app = conduction
    source_variable = kappa_fission
    variable = power
    from_postprocessors_to_be_preserved = tally_integral
    to_postprocessors_to_be_preserved = tally_integral
  []
  [temp_to_openmc]
    type = MultiAppGeneralFieldShapeEvaluationTransfer
    from_multi_app = conduction
    variable = temp
    source_variable = T
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 800
  []
  [power]
    type = ConstantIC
    variable = kappa_fission
    value = ${fparse power/1.865e+03}
  []
[]

[Postprocessors]
  [tally_integral]
    type = ElementIntegralVariablePostprocessor
    variable = kappa_fission
    execute_on = 'initial transfer timestep_end'
  []
[]

[Executioner]
  type = Transient
  num_steps = 3
[]

[Outputs]
  exodus = true
[]

(test/tests/userobjects/moab_skinner/all_bins.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '445 445'
    fields_max = '655 655'
    n_field_bins = '3 5'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"



    output_skins = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(test/tests/userobjects/moab_skinner/graveyard/graveyard.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '445 445'
    fields_max = '655 655'
    n_field_bins = '3 5'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"



    output_skins = true
    build_graveyard = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(tutorials/tokamak/openmc.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = tokamak.e
  []
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  scaling = 100.0

  source_strength = 2e18

  cell_level = 0
  temperature_blocks = 'plasma_facing multiplier breeder divertor vacuum_vessel'

  # this is a low number of particles; you will want to increase in order to obtain
  # high-quality results
  first_iteration_particles = 1000
  relaxation = dufek_gudowski

  skinner = moab

  [Tallies]
    [tokamak]
      type = MeshTally
      mesh_template = tokamak.e
      score = 'heating_local H3_production'
      output = unrelaxed_tally_std_dev
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '2000'
    n_field_bins = '40'
    build_graveyard = true
    output_skins = true
  []
[]

[Postprocessors]
  [heating]
    type = ElementIntegralVariablePostprocessor
    variable = heating_local
  []
  [tritium_production]
    type = ElementIntegralVariablePostprocessor
    variable = H3_production
  []
  [tritium_error]
    type = TallyRelativeError
    tally_score = H3_production
    value_type = average
  []
  [heating_error]
    type = TallyRelativeError
    tally_score = heating_local
    value_type = average
  []
[]

[Executioner]
  type = Transient
[]

[Outputs]
  exodus = true
  csv = true
[]

(test/tests/neutronics/dagmc/cell_tallies/with_relaxation.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  allow_renumbering = false
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 500.0
  []
[]

[Problem]
  type = OpenMCCellAverageProblem

  temperature_blocks = '1 2'
  cell_level = 0
  power = 16.0

  skinner = moab

  relaxation = 'constant'

  [Tallies]
    [Mesh]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '300'
    fields_max = '1500'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

(test/tests/neutronics/mg/rr_thermal_expansion/openmc.i)

# Cardinal input file for slab benchmark OpenMC model to couple with thermomechanics

P = 1.0e22 # eV/s
eV_to_J = 1.602e-19 # J/eV

[Mesh]
  [file]
    type = FileMeshGenerator
    file = slab20v.e
  []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
  use_displaced_mesh = true
[]

[Variables]
  [disp_x]
  []
  [disp_y]
    initial_condition = 0.0
  []
  [disp_z]
    initial_condition = 0.0
  []
[]

[AuxVariables]
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[ICs]
  [temperature]
    type = FunctionIC
    variable = temp
    function = 293
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  source_rate_normalization = 'kappa_fission'

  inactive_batches = 100
  batches = 200
  particles = 10

  power = ${fparse P*eV_to_J}
  cell_level = 0
  temperature_blocks = '1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20'
  volume_calculation = vol
  skinner = moab
  [Tallies]
    [kappa_fission]
      type = CellTally
      name = kappa_fission
      output = unrelaxed_tally_std_dev
      normalize_by_global_tally = false
    []
    [flux]
      type = CellTally
      score = flux
      output = unrelaxed_tally_std_dev
      normalize_by_global_tally = false
    []
  []
[]

[Executioner]
  type = Transient
  start_time = 0.0
  num_steps = 5
  steady_state_detection = true
[]

[MultiApps]
  [solid]
    type = TransientMultiApp
    input_files = 'solid.i'
    execute_on = timestep_end
  []
[]

[Transfers]
  [kappa_fission_to_solid]
    type = MultiAppGeneralFieldShapeEvaluationTransfer
    to_multi_app = solid
    variable = kappa_fission
    source_variable = kappa_fission
    from_postprocessors_to_be_preserved = heat_source
    to_postprocessors_to_be_preserved = source_integral
    execute_on = timestep_end
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
  [temp_from_solid]
    type = MultiAppGeneralFieldNearestLocationTransfer
    from_multi_app = solid
    variable = temp
    source_variable = temperature_praux
    execute_on = timestep_end
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
  [dispx]
    type = MultiAppGeneralFieldNearestLocationTransfer
    source_variable = disp_praux
    variable = disp_x
    from_multi_app = solid
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
[]

[Outputs]
  csv = true
  execute_on = 'TIMESTEP_END'
[]

[Postprocessors]
  [heat_source]
    type = ElementIntegralVariablePostprocessor
    variable = kappa_fission
    execute_on = 'transfer initial timestep_end'
  []
  [k]
    type = KEigenvalue
    value_type = 'tracklength'
  []
  [k_std_dev]
    type = KEigenvalue
    value_type = 'tracklength'
    output = 'std_dev'
  []
  [k_less_than_initial]
    type = ParsedPostprocessor
    expression = 'if (k < initial_k, 1, 0)'
    pp_names = 'k'
    constant_names = 'initial_k'
    constant_expressions = '0.31361125245905'
  []
[]

[UserObjects]
  [vol]
    type = OpenMCVolumeCalculation
    n_samples = 200000
  []
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    execute_on = 'timestep_begin'
    verbose = true
    build_graveyard = true
    graveyard_scale_inner = 1.5
    graveyard_scale_outer = 1.6
  []
[]

(tutorials/msfr/openmc.i)

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = msr.e
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 948.0
  []
  [density]
    type = ConstantIC
    variable = density
    value = ${fparse -0.882*948+4983.6}
  []
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
  [cell_density]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
  [cell_density]
    type = CellDensityAux
    variable = cell_density
  []
  [density]
    type = ParsedAux
    variable = density
    expression = '-0.882*temp+4983.6'
    coupled_variables = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = false

  # this will start each Picard iteration from the fission source from the previous one
  reuse_source = true

  scaling = 100.0

  density_blocks = '1'
  temperature_blocks = '1'
  cell_level = 0

  power = 300.0e6

  #relaxation = dufek_gudowski
  #first_iteration_particles = 5000
  skip_statepoint = true

  skinner = moab

  [Tallies]
    [heat_source]
      type = MeshTally
      mesh_template = msr.e
      output = unrelaxed_tally_std_dev
    []
  []
[]

nb = 15.0
tmin = 700.0            # Moab skinner gives an error because of max density -> Minimum temperature below 800K --> This might be an NekRS error that gives negatives nondimensional temperature?
tmax = 1250.0

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '${tmin} ${fparse -0.882*tmax+4983.6}'
    fields_max = '${tmax} ${fparse -0.882*tmin+4983.6}'
    n_field_bins = '${nb} ${nb}'


    build_graveyard = true
    output_skins = true
  []
[]

[MultiApps]
  [nek]
    type = TransientMultiApp
    app_type = CardinalApp
    input_files = 'nek.i'
    sub_cycling = true
    execute_on = timestep_end
  []
[]

[Transfers]
  [temp_to_openmc]
    type = MultiAppGeneralFieldNearestLocationTransfer
    from_multi_app = nek
    variable = temp
    source_variable = temperature
  []
  [power_to_nek]
    type = MultiAppGeneralFieldNearestLocationTransfer
    to_multi_app = nek
    source_variable = kappa_fission
    variable = source
  []
  [power_integral_to_nek]
    type = MultiAppPostprocessorTransfer
    to_postprocessor = source_integral
    from_postprocessor = power
    to_multi_app = nek
  []
  [synchronize_in]
    type = MultiAppPostprocessorTransfer
    to_postprocessor = transfer_in
    from_postprocessor = synchronize
    to_multi_app = nek
  []
[]

[Postprocessors]
  [power]
    type = ElementIntegralVariablePostprocessor
    variable = kappa_fission
  []
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
  [max_tally_err]
    type = TallyRelativeError
  []
  [max_T]
    type = ElementExtremeValue
    variable = temp
  []
  [avg_T]
    type = ElementAverageValue
    variable = temp
  []
  [max_q]
    type = ElementExtremeValue
    variable = kappa_fission
  []
  [synchronize]
    type = Receiver
    default = 1.0
  []
[]

t_nek = ${fparse 2e-4 * 7.669}

[Executioner]
  type = Transient
  dt = ${fparse 2000 * t_nek}
[]

[Outputs]
  exodus = true
  csv = true
  hide = 'synchronize'
[]

(test/tests/neutronics/kinetics/geom_reset/one_bin.i)

[Mesh]
  type = FileMesh
  file = ../../dagmc/mesh_tallies/slab.e
  allow_renumbering = false
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  calc_kinetics_params = true
  ifp_generations = 5

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [lambda]
    type = LambdaEffective
  []
  [lambda_rel]
    type = LambdaEffective
    output = rel_err
  []
  [lambda_std]
    type = LambdaEffective
    output = std_dev
  []
  [beta]
    type = BetaEffective
  []
  [beta_rel]
    type = BetaEffective
    output = rel_err
  []
  [beta_std]
    type = BetaEffective
    output = std_dev
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  csv = true
  execute_on = 'TIMESTEP_END'
[]

(test/tests/userobjects/moab_skinner/errors/skinner.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../../neutronics/meshes/tet_cube.e
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
    family = MONOMIAL
    order = CONSTANT
    initial_condition = 500.0
  []
  [rho]
    family = MONOMIAL
    order = CONSTANT
    initial_condition = 1000.0
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '400 900'
    fields_max = '650 2000'
    n_field_bins = '5 3'
    material_blocks = "1"
    material_names = "mat"



    build_graveyard = true
    graveyard_scale_inner = 1.05
    graveyard_scale_outer = 1.25
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

(test/tests/neutronics/dagmc/with_csg/dag_lattice_outer/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab


  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/userobjects/moab_skinner/hex-mesh/mixed_elements.i)

[Mesh]
  [cylinder]
    type = AnnularMeshGenerator
    nr = 10
    nt = 20
    rmin = 0
    rmax = 1
    quad_subdomain_id = 2
    tri_subdomain_id = 3
    growth_r = -1.2
  []
  [extrude]
    type = AdvancedExtruderGenerator
    input = cylinder
    heights = '10'
    num_layers = '10'
    direction = '0 0 1'
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = '400 + x * 100 + 100 * t'
  []
  [rho]
    type = ParsedFunction
    expression = '400 + y * 100 + 100 * t'
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '295 400'
    fields_max = '715 715'
    n_field_bins = '3 10'
    verbose = true
    material_blocks = '2 3'
    material_names = 'mat2 mat3'



    output_skins = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
[]

(test/tests/neutronics/dagmc/properties.i)

[Mesh]
  [file]
    type = FileMeshGenerator
    file = ../meshes/tet_cube.e
  []
  parallel_type = replicated
[]

[Problem]
  type = OpenMCCellAverageProblem
  cell_level = 0
  temperature_blocks = '1'
  power = 1000.0
  skinner = moab

  initial_properties = hdf5

  [Tallies]
    [Cell]
      type = CellTally
      block = '1'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '100'
    n_field_bins = '1'
  []
[]

[Executioner]
  type = Steady
[]

(doc/content/vv/slab_benchmark_model/openmc.i)

P = ${units 1.0e22 eV/s}
eV_to_J = ${units 1.602176634e-19 J/eV}
initial_temperature = ${units 293 K}

[Mesh]
  [file]
    type = FileMeshGenerator
    file = slab20v.e
  []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
  use_displaced_mesh = true
[]

[Variables]
  [disp_x]
  []
  [disp_y]
    initial_condition = 0.0
  []
  [disp_z]
    initial_condition = 0.0
  []
[]

[AuxVariables]
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[ICs]
  [temperature]
    type = FunctionIC
    variable = temperature
    function = ${initial_temperature}
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  source_rate_normalization = 'kappa_fission'
  inactive_batches = 150
  batches = 300
  particles = 10000000
  power = '${fparse P*eV_to_J}'
  cell_level = 0
  temperature_variables = 'temperature'
  temperature_blocks = '1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20'
  volume_calculation = vol
  skinner = moab
  [Tallies]
    [kappa_fission]
      type = CellTally
      estimator = collision
      block = '1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20'
      name = kappa_fission
      output = unrelaxed_tally_std_dev
    []
    [flux]
      type = CellTally
      estimator = collision
      block = '1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20'
      score = flux
      output = unrelaxed_tally_std_dev
    []
  []
[]

[Executioner]
  type = Transient
  start_time = 0.0
  end_time = 1e5
  dt = 50
[]

[MultiApps]
  [solid]
    type = TransientMultiApp
    input_files = 'solid.i'
    execute_on = timestep_end
  []
[]

[Transfers]
  [kappa_fission_to_solid]
    type = MultiAppGeneralFieldShapeEvaluationTransfer
    to_multi_app = solid
    variable = kappa_fission
    source_variable = kappa_fission
    from_postprocessors_to_be_preserved = heat_source
    to_postprocessors_to_be_preserved = source_integral
    execute_on = timestep_end
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
  [temperature_from_solid]
    type = MultiAppGeneralFieldNearestLocationTransfer
    from_multi_app = solid
    variable = temperature
    source_variable = temperature_praux
    execute_on = timestep_end
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
  [dispx]
    type = MultiAppGeneralFieldNearestLocationTransfer
    source_variable = disp_praux
    variable = disp_x
    from_multi_app = solid
    displaced_source_mesh = true
    displaced_target_mesh = true
  []
[]

[Outputs]
  exodus = true
  csv = true
[]

[Postprocessors]
  [heat_source]
    type = ElementIntegralVariablePostprocessor
    variable = kappa_fission
    execute_on = 'transfer initial timestep_end'
  []
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[VectorPostprocessors]
  [temperature]
    type = ElementValueSampler
    variable = 'temperature'
    sort_by = x
    execute_on = 'FINAL'
  []
[]

[UserObjects]
  [vol]
    type = OpenMCVolumeCalculation
    n_samples = 200000
  []
  [moab]
    type = MoabSkinner
    fields = 'temperature'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '20'
    execute_on = 'timestep_begin'
    verbose = true
    build_graveyard = true
    graveyard_scale_inner = 1.5
    graveyard_scale_outer = 1.6
  []
[]

(test/tests/neutronics/dagmc/with_csg/multi_cell_dag/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/userobjects/moab_skinner/hex-mesh/all_bins.i)

[Mesh]
  [cube]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -0.5
    xmax = 0.5
    ymin = -0.5
    ymax = 0.5
    zmin = -0.5
    zmax = 0.5
    nx = 5
    ny = 5
    nz = 5
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '445 445'
    fields_max = '655 655'
    n_field_bins = '3 5'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"



    output_skins = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(test/tests/neutronics/dagmc/with_csg/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'

      # Since we're ignoring the CSG geometry when mapping tallies, we need to
      # disable the tally sum check.
      check_tally_sum = false
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = false
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  csv = true
[]

(test/tests/neutronics/dagmc/density_skin/disjoint_bins.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

T0 = 600.0
dT = 50.0

[Functions]
  [density]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 - dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = density
    function = density
    execute_on = timestep_begin
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 500.0
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true

  temperature_blocks = '1 2'
  density_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '0 ${fparse T0 - dT}'
    fields_max = '1000 ${fparse T0 + 2 * dT}'
    n_field_bins = '1 4'


    build_graveyard = true
    output_skins = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp density cell_instance cell_id'
[]

(test/tests/neutronics/dagmc/mesh_tallies/openmc.i)

!include csg_step_2_comp.i

[Mesh]
  type = FileMesh
  file = slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true

  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = slab.e
      output = 'unrelaxed_tally_std_dev'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  execute_on = 'FINAL'
[]

(test/tests/neutronics/moving_mesh/scale.i)

[GlobalParams]
  use_displaced_mesh = true
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [f]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -0.05
    xmax = 0.15
    ymin = -0.05
    ymax = 0.05
    zmin = -0.05
    zmax = 0.05
    nx = 20
    ny = 5
    nz = 5
  []
  [subdomains]
    type = ParsedSubdomainIDsGenerator
    input = f
    expression = 'if (x < 0.05, 1, 2)'
  []
  [name]
    type = RenameBlockGenerator
    input = subdomains
    old_block = '1 2'
    new_block = 'left right'
  []
  [tets]
    type = ElementsToTetrahedronsConverter
    input = name
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  scaling = 100.0
  temperature_blocks = 'left right'
  cell_level = 0
  skinner = skin
  verbose = true

  # settings to match the reference/ case
  #particles = 500000
  #inactive_batches = 200
  #batches = 1000
[]

[AuxVariables]
  [disp_x]
  []
  [disp_y]
  []
  [disp_z]
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 1000
  []
[]

[AuxKernels]
  [disp_x] # this executes after the [Problem]
    type = FunctionAux
    variable = disp_x
    function = disp_x
  []
[]

[Functions]
  [disp_x]
    type = ParsedFunction
    expression = 'exp((x+0.05))-1'
  []
[]

[UserObjects]
  [skin]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '2000'
    n_field_bins = '1'
    build_graveyard = true
    verbose = True
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Outputs]
  csv = true
  exodus = true
[]

(test/tests/userobjects/moab_skinner/output/mesh.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
    family = MONOMIAL
    order = CONSTANT
    initial_condition = 500.0
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '445'
    fields_max = '655'
    n_field_bins = '3'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"


    output_full = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

(test/tests/neutronics/dagmc/mesh_tallies/disjoint_bins.i)

!include csg_step_2_comp.i

[Mesh]
  type = FileMesh
  file = slab.e
  allow_renumbering = false
  parallel_type = replicated
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 - dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true

  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = slab.e
      output = 'unrelaxed_tally_std_dev'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = true
    output_skins = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  execute_on = 'FINAL'
  hide := 'kappa_fission_csg kappa_fission
           kappa_fission_std_dev_csg kappa_fission_std_dev
           k k_std_dev k_csg_pp k_std_dev_csg_pp'
[]

(test/tests/neutronics/dagmc/with_csg/dag_lattice/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab


  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/userobjects/moab_skinner/three_fields.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [zdep]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [rho_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [zdep_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [zdep]
    type = FunctionAux
    variable = zdep
    function = zdep
    execute_on = timestep_begin
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [rho_bins]
    type = SkinnedBins
    variable = rho_bins
    skinner = moab
    skin_by = rho
  []
  [zdep_bins]
    type = SkinnedBins
    variable = zdep_bins
    skinner = moab
    skin_by = zdep
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
  [zdep]
    type = ParsedFunction
    expression = 11+z*10+5*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho zdep'
    fields_min = '445 445 11'
    fields_max = '655 655 31'
    n_field_bins = '3 4 2'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"

  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp rho zdep'
[]

(test/tests/neutronics/dagmc/cell_tallies/with_scaling.i)

scale = 100.0

[Mesh]
  [file]
    type = FileMeshGenerator
    file = scale_in.e
  []
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = '500.0 + 10.0*x*${scale}'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  scaling = ${scale}

  temperature_blocks = '1 2'
  cell_level = 0
  power = 16.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '300'
    fields_max = '1500'
    n_field_bins = '10'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

[Outputs]
  exodus = true
[]

(test/tests/openmc_errors/input_params/no_dag.i)

[Mesh]
  [sphere]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/sphere.e
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  temperature_blocks = '1'
  cell_level = 0
  skinner = moab
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0.0'
    fields_max = '1000'
    n_field_bins = '5'
  []
[]

[Executioner]
  type = Transient
[]

(test/tests/neutronics/dagmc/mesh_tallies/one_bin.i)

[Mesh]
  type = FileMesh
  file = slab.e
  allow_renumbering = false
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem

  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = slab.e
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp cell_instance'
[]

(test/tests/userobjects/moab_skinner/sub_bins.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [subdomain_bins]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [temp_bins]
    type = SkinnedBins
    variable = temp_bins
    skinner = moab
    skin_by = temp
  []
  [subdomain_bins]
    type = SkinnedBins
    variable = subdomain_bins
    skinner = moab
    skin_by = subdomain
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '445'
    fields_max = '655'
    n_field_bins = '3'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"

  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/neutronics/dagmc/assign_materials/generate_dagmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  [block_1]
    type = ParsedSubdomainMeshGenerator
    input = slab
    combinatorial_geometry = 'x < 37.5'
    block_id = 1
    block_name = block_1
    epsilon = 1e-3
  []
  [block_2]
    type = ParsedSubdomainMeshGenerator
    input = block_1
    combinatorial_geometry = 'x >= 37.5'
    block_id = 2
    block_name = block_2
    epsilon = 1e-3
  []
  allow_renumbering = false
[]

[AuxVariables]
  [temp]
    order = CONSTANT
    family = MONOMIAL
    initial_condition = 200
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
    graveyard_scale_inner = 1.1
    graveyard_scale_outer = 1.2
    material_blocks = 'block_1 block_2'
    material_names = '1 U35_U38'
    output_full = true
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/userobjects/moab_skinner/second_order/all_bins.i)

[Mesh]
  [cube]
    type = FileMeshGenerator
    file = ../../../neutronics/meshes/tet_cube.e
  []
  [id1]
    type = ParsedSubdomainMeshGenerator
    input = cube
    combinatorial_geometry = 'z < 0.0'
    block_id = 1
    epsilon = 1e-8
  []
  [id2]
    type = ParsedSubdomainMeshGenerator
    input = id1
    combinatorial_geometry = 'z >= 0.0'
    block_id = 3
    epsilon = 1e-8
  []

  second_order = true
[]

[Problem]
  type = FEProblem
  solve = false
[]

[AuxVariables]
  [rho]
    family = MONOMIAL
    order = CONSTANT
  []
  [temp]
    family = MONOMIAL
    order = CONSTANT
  []
  [bins]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [rho]
    type = FunctionAux
    variable = rho
    function = rho
    execute_on = timestep_begin
  []
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [bins]
    type = SkinnedBins
    variable = bins
    skinner = moab
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = 400+x*100+100*t
  []
  [rho]
    type = ParsedFunction
    expression = 400+y*100+100*t
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp rho'
    fields_min = '445 445'
    fields_max = '655 655'
    n_field_bins = '3 5'
    verbose = true
    material_blocks = "1 3"
    material_names = "mat mat"



    output_skins = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp rho'
[]

(test/tests/neutronics/moving_mesh/volumes.i)

[GlobalParams]
  use_displaced_mesh = true
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  [f]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -0.05
    xmax = 0.15
    ymin = -0.05
    ymax = 0.05
    zmin = -0.05
    zmax = 0.05
    nx = 20
    ny = 2
    nz = 2
  []
  [subdomains]
    type = ParsedSubdomainIDsGenerator
    input = f
    expression = 'if (x < 0.05, 1, 2)'
  []
  [name]
    type = RenameBlockGenerator
    input = subdomains
    old_block = '1 2'
    new_block = 'left right'
  []
  [tets]
    type = ElementsToTetrahedronsConverter
    input = name
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  scaling = 100.0
  temperature_blocks = 'left right'
  cell_level = 0
  skinner = skin
  volume_calculation = vol
  verbose = true
[]

[AuxVariables]
  [disp_x]
  []
  [disp_y]
  []
  [disp_z]
  []

  [mapped_vol]
    family = MONOMIAL
    order = CONSTANT
  []
  [actual_openmc_vol]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[ICs]
  [temp]
    type = ConstantIC
    variable = temp
    value = 1000
  []
[]

[AuxKernels]
  [disp_x] # this executes after the [Problem]
    type = FunctionAux
    variable = disp_x
    function = disp_x
  []
  [mapped_vol]
    type = CellVolumeAux
    variable = mapped_vol
    volume_type = mapped
  []
  [actual_openmc_vol]
    type = CellVolumeAux
    variable = actual_openmc_vol
    volume_type = actual
  []
[]

[Functions]
  [disp_x]
    type = ParsedFunction
    expression = 'exp((x+0.05))-1'
  []
[]

[UserObjects]
  [skin]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '2000'
    n_field_bins = '1'
    build_graveyard = true
    verbose = True
  []
  [vol]
    type = OpenMCVolumeCalculation
    n_samples = 1000000
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Postprocessors]
  [moose_vol_left]
    type = VolumePostprocessor
    block = 'left'
  []
  [moose_vol_right]
    type = VolumePostprocessor
    block = 'right'
  []
  [mapped_vol_left]
    type = ElementExtremeValue
    variable = mapped_vol
    block = 'left'
  []
  [mapped_vol_right]
    type = ElementExtremeValue
    variable = mapped_vol
    block = 'right'
  []
  [openmc_vol_left]
    type = ElementExtremeValue
    variable = actual_openmc_vol
    block = 'left'
  []
  [openmc_vol_right]
    type = ElementExtremeValue
    variable = actual_openmc_vol
    block = 'right'
  []
[]

[Outputs]
  csv = true
[]

(test/tests/neutronics/dagmc/assign_materials/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  [block_1]
    type = ParsedSubdomainMeshGenerator
    input = slab
    combinatorial_geometry = 'x < 37.5'
    block_id = 1
    block_name = block_1
    epsilon = 1e-3
  []
  [block_2]
    type = ParsedSubdomainMeshGenerator
    input = block_1
    combinatorial_geometry = 'x >= 37.5'
    block_id = 2
    block_name = block_2
    epsilon = 1e-3
  []
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxVariables]
  [material_id]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
  [density]
    type = FunctionAux
    variable = density
    function = temp
    execute_on = timestep_begin
  []
  [material_id]
    type = CellMaterialIDAux
    variable = material_id
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  density_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '0 0'
    fields_max = '1000 1000'
    n_field_bins = '4 1'
    build_graveyard = true
    graveyard_scale_inner = 1.1
    graveyard_scale_outer = 1.2
    material_blocks = 'block_1 2'
    material_names = 'h2o 1'
  []
[]

[Postprocessors]
  [mat_in_block_1]
    type = PointValue
    variable = material_id
    point = '12.5 12.5 12.5'
  []
  [mat_in_block_2]
    type = PointValue
    variable = material_id
    point = '62.5 12.5 12.5'
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
[]

(test/tests/neutronics/dagmc/cell_tallies/null_density.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[ICs]
  [density]
    type = ConstantIC
    variable = density
    value = 11000.0
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  density_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp density'
    fields_min = '${fparse T0 - dT} 0'
    fields_max = '${fparse T0 + 2 * dT} 15000'
    n_field_bins = '4 1'
    build_graveyard = true

    # Effectively no density skinning at all, because theres just one bin.
    # Results should be the same as if there was no density skinning enabled.
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'density cell_instance'
[]

(test/tests/neutronics/dagmc/assign_materials/no_material_names.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../mesh_tallies/slab.e
  []
  allow_renumbering = false
[]

[AuxVariables]
  [temp]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Problem]
  type = FEProblem
  solve = false
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
  []
[]

[Executioner]
  type = Steady
[]

(test/tests/neutronics/dagmc/cell_tallies/openmc.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
[]

(test/tests/neutronics/dagmc/with_csg/multiple_dag/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 1
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

(test/tests/neutronics/dagmc/implicit_comp_material/openmc.i)

[Mesh]
  type = FileMesh
  file = ../mesh_tallies/slab.e
  allow_renumbering = false
[]

x0 = 12.5
x1 = 37.5
x2 = 62.5

T0 = 600.0
dT = 50.0

[Functions]
  [temp]
    type = ParsedFunction
    expression = 'if (x <= ${x0}, ${fparse T0 - dT}, if (x <= ${x1}, ${T0}, if (x <= ${x2}, ${fparse T0 + dT}, ${fparse T0 + 2 * dT})))'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab

  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '${fparse T0 - dT}'
    fields_max = '${fparse T0 + 2 * dT}'
    n_field_bins = '4'
    build_graveyard = true
    implicit_complement_material="h2o"
    graveyard_scale_inner = 1.1
    graveyard_scale_outer = 1.2
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
  [k_std_dev]
    type = KEigenvalue
    output = 'std_dev'
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
[]

(test/tests/neutronics/dagmc/mesh_tallies/with_scaling.i)

scale = 100.0

[Mesh]
  [file]
    type = FileMeshGenerator
    file = scale_in.e
  []
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[Functions]
  [temp]
    type = ParsedFunction
    expression = '500.0 + 10.0*x*${scale}'
  []
[]

[AuxKernels]
  [temp]
    type = FunctionAux
    variable = temp
    function = temp
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  scaling = ${scale}

  temperature_blocks = '1 2'
  cell_level = 0
  power = 16.0

  skinner = moab

  [Tallies]
    [Mesh]
      type = MeshTally
      mesh_template = scale_in.e
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '300'
    fields_max = '1500'
    n_field_bins = '10'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

[Outputs]
  exodus = true
[]

(test/tests/neutronics/dagmc/add_bcs/openmc_hex.i)

[Mesh]
  [f]
    type = GeneratedMeshGenerator
    dim = 3
    xmin = -12.5
    xmax = 87.5
    ymin = -12.5
    ymax = 37.5
    zmin = -12.5
    zmax = 12.5
    nx = 8
    ny = 4
    nz = 2
  []
  [subdomains]
    type = ParsedSubdomainIDsGenerator
    input = f
    expression = 'if (y > 12.5, 1, 2)'
  []
  [name]
    type = RenameBlockGenerator
    input = subdomains
    old_block = '1 2'
    new_block = 'block_1 block_2'
  []
  # The following sidesets/nodesets are only used by the error-checking and internal-surface tests
  [left_dup] # duplicate of the 'left' sideset, to test one side belonging to two sidesets
    type = ParsedGenerateSideset
    input = name
    combinatorial_geometry = 'x < -12.49'
    normal = '-1 0 0'
    new_sideset_name = 'left_dup'
  []
  [interface_a] # the block_1/block_2 interface, attached to block_1 element sides
    type = SideSetsBetweenSubdomainsGenerator
    input = left_dup
    primary_block = 'block_1'
    paired_block = 'block_2'
    new_boundary = 'interface_a'
  []
  [interface_b] # the same interface, attached to block_2 element sides
    type = SideSetsBetweenSubdomainsGenerator
    input = interface_a
    primary_block = 'block_2'
    paired_block = 'block_1'
    new_boundary = 'interface_b'
  []
  [midplane] # a plane interior to block_2 (y = 0), which never appears in any skin
    type = ParsedGenerateSideset
    input = interface_b
    combinatorial_geometry = 'y > -0.01 & y < 0.01'
    normal = '0 1 0'
    new_sideset_name = 'midplane'
  []
  [corner_nodes] # a nodeset-only boundary, to test the nodeset error message
    type = ExtraNodesetGenerator
    input = midplane
    new_boundary = 'corner_nodes'
    coord = '-12.5 -12.5 -12.5'
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2'
  cell_level = 0
  power = 100.0

  skinner = moab
[]

[UserObjects]
  [moab]
    # BCs are added in tests as cli_args to use this file in all different cases
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    verbose = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
[]

[Outputs]
  exodus = false
  [file]
    type = CSV
  []
[]

(test/tests/neutronics/dagmc/with_csg/mixed_feedback/openmc.i)

[Mesh]
  [slab]
    type = FileMeshGenerator
    file = ../../mesh_tallies/slab.e
  []
  [cube]
    type = FileMeshGenerator
    file = ../../../meshes/tet_cube.e
  []
  [rename_cube]
    type = RenameBlockGenerator
    input = cube
    old_block = '1'
    new_block = '3'
  []
  [combine]
    type = CombinerGenerator
    inputs = 'slab rename_cube'
    positions = '0 0 0 0 50 0'
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[AuxVariables]
  [cell_temperature]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [cell_temperature]
    type = CellTemperatureAux
    variable = cell_temperature
  []
[]

[AuxKernels]
  [temp]
    type = ConstantAux
    variable = temp
    value = 500.0
    execute_on = timestep_begin
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  temperature_blocks = '1 2 3'
  cell_level = 0
  power = 100.0

  skinner = moab


  [Tallies]
    [Cell]
      type = CellTally
      block = '1 2 3'
    []
  []
[]

[UserObjects]
  [moab]
    type = MoabSkinner
    fields = 'temp'
    fields_min = '0'
    fields_max = '1000'
    n_field_bins = '1'
    build_graveyard = true
  []
[]

[Postprocessors]
  [k]
    type = KEigenvalue
  []
[]

[Executioner]
  type = Transient
  num_steps = 2
[]

[Outputs]
  exodus = true
  hide = 'temp'
[]

Description
Example Input Syntax
Material Names and Blocks
Input Parameters
Input Files