AxisymmetricPointScalarFlux

construction:Undocumented Class

The AxisymmetricPointScalarFlux has not been documented. The content listed below should be used as a starting point for documenting the class, which includes the typical automatic documentation associated with a MooseObject; however, what is contained is ultimately determined by what is necessary to make the documentation clear for users.

A function which computes the analytical scalar flux for a point source surrounded by axisymmetric shielding media.

Overview

Example Input File Syntax

Input Parameters

num_dimsThe dimensionality of the solution. 2D yields a flatland, 3D yields spherical cartesian.
C++ Type:unsigned int
Controllable:No
Description:The dimensionality of the solution. 2D yields a flatland, 3D yields spherical cartesian.

Required Parameters

axisymmetric_cross_sectionsThe macroscopic total cross-sections for each radial region.
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:The macroscopic total cross-sections for each radial region.
axisymmetric_radiiThe ending radii for each axisymmetric region.
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:The ending radii for each axisymmetric region.
src_location0 0 0The location of the point source.
Default:0 0 0
C++ Type:libMesh::Point
Controllable:No
Description:The location of the point source.
src_strength1The emission rate of the point source.
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The emission rate of the point source.

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.

Advanced Parameters

Input Files

(examples/2D_subcritical_reactor_4_groups/neutronics_merge.i)
(examples/new_uncollided_method/spatial_convergence_2D/quads/spatial_convergence.i)
(examples/new_uncollided_method/spatial_convergence_2D/tris/spatial_convergence.i)

(examples/2D_subcritical_reactor_4_groups/neutronics_merge.i)

# An input deck which merges an analytically computed near-source
# region flux with the numerically computed SASF approach.

[Mesh]
  [NeutronicsDomain]
    type = FileMeshGenerator
    file = 2D_quarter_core.e
  []
  uniform_refine = 2
[]

[AuxVariables]
  [Flux_Moment_1_0_0_uncollided_NSR]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
  []
  [Flux_Moment_1_0_0_uncollided_FSR]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
  []

  [uncollided_flux_moment_1_0_0]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
    initial_condition = '0.0'
  []
  [uncollided_flux_moment_2_0_0]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
    initial_condition = '0.0'
  []
  [uncollided_flux_moment_3_0_0]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
    initial_condition = '0.0'
  []
  [uncollided_flux_moment_4_0_0]
    type = MooseVariable
    order = FIRST
    family = L2_LAGRANGE
    initial_condition = '0.0'
  []
[]

[AuxKernels]
  [NSR_Analytical_Flux]
    type = FunctionAux
    variable = Flux_Moment_1_0_0_uncollided_NSR
    function = AnalyticalFlux
    block = src_air
    execute_on = INITIAL
  []

  [FSR_Sum]
    type = SumAux
    variable = uncollided_flux_moment_1_0_0
    values = 'Flux_Moment_1_0_0_uncollided_NSR Flux_Moment_1_0_0_uncollided_FSR'
    execute_on = TIMESTEP_END
  []
[]

[Functions]
  [AnalyticalFlux]
    type = AxisymmetricPointScalarFlux
    num_dims = 2
    src_strength = 1.5
    src_location = '0 0 0'
    axisymmetric_radii = '15.0'
    axisymmetric_cross_sections = '9.683916522303286e-05'
  []
[]

[Transfers]
  [Group_1]
    type = MultiAppGeneralFieldShapeEvaluationTransfer
    variable = Flux_Moment_1_0_0_uncollided_FSR
    source_variable = uncollided_flux_moment_1_0_0
    from_multi_app = Uncollided
  []
[]

[Problem]
  type = FEProblem
  solve = false
  kernel_coverage_check = false
[]

[Executioner]
  type = Steady
[]

[MultiApps]
  [Uncollided]
    type = FullSolveMultiApp
    app_type = GnatApp
    input_files = 'neutronics_uncollided.i'
    execute_on = INITIAL
  []
[]

(examples/new_uncollided_method/spatial_convergence_2D/quads/spatial_convergence.i)

# A simple test case with a purely absorbing medium and a point source in the
# middle of the domain.

src_x = 0.0
src_y = 0.0
src_z = 0.0

shield_xs = 1e0
void_xs = 0.0

[GlobalParams]
  source_location = '${src_x} ${src_y} ${src_z}'
[]

[Mesh]
  [domain]
    type = FileMeshGenerator
    file = 'spatial_convergence_2_quad.e'
  []
[]

[AuxVariables]
  [AnalyticalFluxField]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
  [UncollidedError]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
[]

[AuxKernels]
  [CompAnalytical]
    type = FunctionAux
    variable = AnalyticalFluxField
    function = AnalyticalFlux
  []
  [SpatialError]
    type = ParsedAux
    variable = UncollidedError
    coupled_variables = 'UncollidedFlux AnalyticalFluxField'
    expression = 'abs(UncollidedFlux - AnalyticalFluxField) / AnalyticalFluxField'
  []
[]

[Variables]
  [UncollidedFlux]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
[]

[Kernels]
  [UncollidedFluxAdvection]
    type = SASFAdvection
    variable = UncollidedFlux
    group_index = 0
    transport_system = ''
  []
  [UncollidedFluxRemoval]
    type = SASFRemoval
    variable = UncollidedFlux
    group_index = 0
    transport_system = ''
  []
[]

[BCs]
  [UncollidedFluxVacuumBC]
    type = SASFVacuumBC
    variable = UncollidedFlux
    boundary = 'vacuum'
  []
  [UncollidedFluxSourceBC]
    type = SASFAnalyticalFluxBC
    variable = UncollidedFlux
    boundary = 'analytical_flux'
    group_index = 0
    group_source = 1.0
    group_total = 0.0
  []
[]

[Materials]
  [Vacuum]
    type = AbsorbingNeutronicsMaterial
    transport_system = ''
    num_groups = 1
    group_total = '${void_xs}'
    block = 'empty'
    saaf_eta = 1.0
  []
  [Shield]
    type = AbsorbingNeutronicsMaterial
    transport_system = ''
    num_groups = 1
    group_total = '${shield_xs}'
    block = 'shield'
    saaf_eta = 1.0
  []
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
  petsc_options_iname = '-pc_type -pc_factor_shift_type'
  petsc_options_value = ' lu       NONZERO'

  l_max_its = 50
  nl_abs_tol = 1e-12

  automatic_scaling = true
  off_diagonals_in_auto_scaling = true
  compute_scaling_once = false
[]

[Functions]
  [AnalyticalFlux]
    type = AxisymmetricPointScalarFlux
    num_dims = 2
    src_strength = 1.0
    src_location = '0 0 0'
    axisymmetric_radii = '2.0 4.0 6.0 8.0 20.0'
    axisymmetric_cross_sections = '${void_xs} ${shield_xs} ${void_xs} ${shield_xs} ${void_xs}'
  []
[]

[Postprocessors]
  [Error]
    type = ElementL2Error
    variable = UncollidedFlux
    function = AnalyticalFlux
  []
  [H]
    type = AverageElementMinSize
  []
[]

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

(examples/new_uncollided_method/spatial_convergence_2D/tris/spatial_convergence.i)

# A simple test case with a purely absorbing medium and a point source in the
# middle of the domain.

src_x = 0.0
src_y = 0.0
src_z = 0.0

shield_xs = 1.0
void_xs = 0.0

[GlobalParams]
  source_location = '${src_x} ${src_y} ${src_z}'
[]

[Mesh]
  [domain]
    type = FileMeshGenerator
    file = 'spatial_convergence_0_tris.e'
  []
[]

[AuxVariables]
  [AnalyticalFluxField]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
  [UncollidedError]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
[]

[AuxKernels]
  [CompAnalytical]
    type = FunctionAux
    variable = AnalyticalFluxField
    function = AnalyticalFlux
  []
  [SpatialError]
    type = ParsedAux
    variable = UncollidedError
    coupled_variables = 'UncollidedFlux AnalyticalFluxField'
    expression = 'abs(UncollidedFlux - AnalyticalFluxField) / AnalyticalFluxField'
  []
[]

[Variables]
  [UncollidedFlux]
    type = MooseVariable
    order = FIRST
    family = LAGRANGE
  []
[]

[Kernels]
  [UncollidedFluxAdvection]
    type = SASFAdvection
    variable = UncollidedFlux
    group_index = 0
    transport_system = ''
  []
  [UncollidedFluxRemoval]
    type = SASFRemoval
    variable = UncollidedFlux
    group_index = 0
    transport_system = ''
  []
[]

[BCs]
  [UncollidedFluxVacuumBC]
    type = SASFVacuumBC
    variable = UncollidedFlux
    boundary = 'vacuum'
  []
  [UncollidedFluxSourceBC]
    type = SASFAnalyticalFluxBC
    variable = UncollidedFlux
    boundary = 'analytical_flux'
    group_index = 0
    group_source = 1.0
    group_total = 0.0
  []
[]

[Materials]
  [Vacuum]
    type = AbsorbingNeutronicsMaterial
    transport_system = ''
    num_groups = 1
    group_total = '${void_xs}'
    block = 'empty'
    saaf_eta = 1.0
  []
  [Shield]
    type = AbsorbingNeutronicsMaterial
    transport_system = ''
    num_groups = 1
    group_total = '${shield_xs}'
    block = 'shield'
    saaf_eta = 1.0
  []
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
  petsc_options_iname = '-pc_type -pc_factor_shift_type'
  petsc_options_value = ' lu       NONZERO'

  l_max_its = 50
  nl_abs_tol = 1e-12

  automatic_scaling = true
  off_diagonals_in_auto_scaling = true
  compute_scaling_once = false
[]

[Functions]
  [AnalyticalFlux]
    type = AxisymmetricPointScalarFlux
    num_dims = 2
    src_strength = 1.0
    src_location = '0 0 0'
    axisymmetric_radii = '2.0 4.0 6.0 8.0 20.0'
    axisymmetric_cross_sections = '${void_xs} ${shield_xs} ${void_xs} ${shield_xs} ${void_xs}'
  []
[]

[Postprocessors]
  [Error]
    type = ElementL2Error
    variable = UncollidedFlux
    function = AnalyticalFlux
  []
  [H]
    type = AverageElementMinSize
  []
[]

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

Overview
Example Input File Syntax
Input Parameters
Input Files