ErrorFractionLookAheadMarker

Description

ErrorFractionLookAheadMarker extends the ErrorFractionMarker from MOOSE to include the statistical relative error from fields in stochastic solvers (i.e. OpenMC). If the spatial error from indicator is within the error fraction provided in refine and the post-refinement statistical relative error is below the threshold provided in rel_error_refine, the element is marked for refinement. If the spatial error is within the error fraction coarsen or the statistical relative error is above the threshold provided in rel_error_refine, the element is marked for coarsening. The post-refinement error is estimated as: $var element = document.getElementById("moose-equation-7c3b9d25-66d1-4776-9359-cd3123829265");katex.render("e_{i+1} = \\sqrt{N_{c}}e_{i}", element, {displayMode:false,throwOnError:false});$ ; where $var element = document.getElementById("moose-equation-c93eb922-f68a-49ea-8ea6-04097043b153");katex.render("e_{i+1}", element, {displayMode:false,throwOnError:false});$ is the post-refinement statistical relative error, $var element = document.getElementById("moose-equation-184dbaa3-a665-4bbd-9e81-8fa574110af2");katex.render("N_{c}", element, {displayMode:false,throwOnError:false});$ is the number of child elements generated when refining, and $var element = document.getElementById("moose-equation-9606466a-dbc9-4631-93d6-ccd117631a33");katex.render("e_{i}", element, {displayMode:false,throwOnError:false});$ is the current statistical relative error. The statistical error is taken from stat_error_indicator.

Example Input File Syntax

[Adaptivity<<<{"href": "../../syntax/Adaptivity/index.html"}>>>]
  [Markers<<<{"href": "../../syntax/Adaptivity/Markers/index.html"}>>>]
    [lh]
      type = ErrorFractionLookAheadMarker<<<{"description": "Marks elements for refinement or coarsening based on the fraction of the min/max error from the supplied indicator and the relative error of a tally value.", "href": "ErrorFractionLookAheadMarker.html"}>>>
      # Statistical error
      rel_error_refine<<<{"description": "The relative error refinement threshold."}>>> = 0.2
      stat_error_indicator<<<{"description": "The name of the statistical relative error Indicator that this Marker uses."}>>> = 'stat_err'
      # Spatial error
      indicator<<<{"description": "The name of the Indicator that this Marker uses."}>>> = 'od'
      refine<<<{"description": "Elements within this percentage of the max error will be refined.  Must be between 0 and 1!"}>>> = 0.3
      coarsen<<<{"description": "Elements within this percentage of the min error will be coarsened.  Must be between 0 and 1!"}>>> = 0.0
    []
  []
[]

Input Parameters

indicatorThe name of the Indicator that this Marker uses.
C++ Type:IndicatorName
Controllable:No
Description:The name of the Indicator that this Marker uses.
rel_error_refineThe relative error refinement threshold.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The relative error refinement threshold.
stat_error_indicatorThe name of the statistical relative error Indicator that this Marker uses.
C++ Type:IndicatorName
Controllable:No
Description:The name of the statistical relative error Indicator that this Marker uses.

Required Parameters

blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
clear_extremesTrueWhether or not to clear the extremes during each error calculation. Changing this to `false` will result in the global extremes ever encountered during the run to be used as the min and max error.
Default:True
C++ Type:bool
Controllable:No
Description:Whether or not to clear the extremes during each error calculation. Changing this to `false` will result in the global extremes ever encountered during the run to be used as the min and max error.
coarsen0Elements within this percentage of the min error will be coarsened. Must be between 0 and 1!
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Elements within this percentage of the min error will be coarsened. Must be between 0 and 1!
refine0Elements within this percentage of the max error will be refined. Must be between 0 and 1!
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Elements within this percentage of the max error will be refined. Must be between 0 and 1!

Optional 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:No
Description:Set the enabled status of the MooseObject.
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Advanced Parameters

(cardinal/test/tests/neutronics/markers/look_ahead/openmc.i)

[Mesh]
  [sphere]
    type = FileMeshGenerator
    file = ../../meshes/sphere.e
  []
  [solid1]
    type = SubdomainIDGenerator
    input = sphere
    subdomain_id = '100'
  []

  allow_renumbering = false
  parallel_type = replicated
[]

[Adaptivity]
  [Markers/lh]
    type = ErrorFractionLookAheadMarker
    # Statistical error
    rel_error_refine = 0.2
    stat_error_indicator = 'stat_err'
    # Spatial error
    indicator = 'od'
    refine = 0.3
    coarsen = 0.0
  []
  [Indicators]
    [od]
      type = ElementOpticalDepthIndicator
      rxn_rate = 'fission'
      h_type = 'cube_root'
    []
    [stat_err]
      type = StatRelErrorIndicator
      score = 'kappa_fission'
    []
  []
[]

[Problem]
  type = OpenMCCellAverageProblem
  verbose = true
  normalize_by_global_tally = false

  power = 100.0
  source_rate_normalization = 'kappa_fission'

  [Tallies]
    [Mesh_1]
      type = MeshTally
      score = 'kappa_fission'
      output = 'unrelaxed_tally_rel_error'
    []
    [Mesh_2]
      type = MeshTally
      score = 'flux fission'
    []
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  execute_on = final
  exodus = true
  hide = 'kappa_fission kappa_fission_rel_error fission flux od stat_err'
[]

Description
Example Input File Syntax
Input Parameters