CircularCrossSectionMeshGenerator

Generates a circular cross section mesh of an axisymmetric fuel rod.

Description

This mesh generator class generates a sliced, circular cross section mesh of fuel and cladding. It allows the user to specify radii, mesh density, and gap sizes. A coating layer may be added to the exterior of the cladding, or multiple claddings may be defined. The fuel or cladding portion of the mesh can be disabled if desired.

The bias parameter affects mesh density moving from the inner edge to the outer edge. For example, a bias value of 2 for a given block doubles element size from the inner to the outer edge, whereas a bias value of 0.5 cuts element size in half going from inner to outer edge. The dual_bias parameter may be used to affect mesh density on both ends of the block. It is also applied from inner to outer edge. When dual_bias is given, bias affects the first half of the elements, and dual_bias affects the right half of the elements. If the two biases are the same, the result will be the same as the original, single-value bias approach.

The resulting mesh will have the following sidesets.

Sideset	Description
`2`	outer surface (named `outer`)
`5`	For a typical fuel and cladding mesh, this is the interior surface of the cladding.
`7`	For a typical fuel and cladding mesh, this is the interior surface of the cladding.
`8`	For a typical fuel and cladding mesh, this is the exterior surface of the fuel.
`9`	For a typical fuel and cladding mesh, this is the union of `7` and `8`.
`10`	For a typical fuel and cladding mesh, this is the exterior surface of the fuel.
`101`	bottom surface (named `bottom`)
`102`	left surface (named `left`)
`1001`, `1002`, `1003`, $var element = document.getElementById("moose-equation-adaf94a2-632f-4ace-8ea2-726cbd731fb2");katex.render("\\ldots", element, {displayMode:false,throwOnError:false});$	every outside and inside radial surface of mesh blocks will receive a sideset, starting with `1001`.
`2001`, `2002`, `2003`, $var element = document.getElementById("moose-equation-76640679-0fa5-4805-9afd-6e1aa308c039");katex.render("\\ldots", element, {displayMode:false,throwOnError:false});$	if `all_bottom_left=true`, every bottom surface of mesh blocks will receive a sideset, starting with `2001`.
`3001`, `3002`, `3003`, $var element = document.getElementById("moose-equation-e77d3ae9-aeb0-4cf7-8ceb-a4adc8253ff6");katex.render("\\ldots", element, {displayMode:false,throwOnError:false});$	if `all_bottom_left=true`, every left surface of mesh blocks will receive a sideset, starting with `3001`.

Meshes of more than one quadrant retain the bottom and left sidesets as defined for the top-right quadrant case.

Figure 1: Typical mesh from CircularCrossSectionMeshGenerator.

The resulting mesh will have the following nodesets:

Nodeset	Description
`10000`, `10001`, `10002`, `10003`, $var element = document.getElementById("moose-equation-58bdfc71-12ae-4212-8f94-a00bea600390");katex.render("\\ldots", element, {displayMode:false,throwOnError:false});$	corresponds to specific nodes at two locations at the interior and exterior of every ring as well as the central node if the first ring is solid. Nodeset `10000` is reserved for the central node of the first ring if it is solid. See Figure 2 for an illustration of the nodeset locations. The nodeset id continues to increase in the same fashion if additional rings are added.

Figure 2: Nodeset ids and locations for a typical circular cross section mesh.

[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
  [ccsmg]
    type = CircularCrossSectionMeshGenerator<<<{"description": "Generates a circular cross section mesh of an axisymmetric fuel rod.", "href": "CircularCrossSectionMeshGenerator.html"}>>>
    num_sectors<<<{"description": "Number of azimuthal sectors in each quadrant. 'num_sectors' must be an even number."}>>> = 20
    elements_per_ring<<<{"description": "A list giving the number of elements in each ring (will be zero for a gap). The first entry is the number of radial elements outside the inner quadrilateral portion of the fuel mesh."}>>> = '5 0 3'
    block_names<<<{"description": "A list giving the block names.  Must be the same length as elements_per_ring.  If an entry in 'elements_per_ring' is zero, the corresponding entry in block_names will be ignored."}>>> = 'fuel null clad'
    coordinates<<<{"description": "Radial coordinates of mesh block boundaries."}>>> = '0.0041 0.00418 0.00474'
  []
[]

Input Parameters

block_namesA list giving the block names. Must be the same length as elements_per_ring. If an entry in 'elements_per_ring' is zero, the corresponding entry in block_names will be ignored.
C++ Type:std::vector<std::string>
Controllable:No
Description:A list giving the block names. Must be the same length as elements_per_ring. If an entry in 'elements_per_ring' is zero, the corresponding entry in block_names will be ignored.
coordinatesRadial coordinates of mesh block boundaries.
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Radial coordinates of mesh block boundaries.
elements_per_ringA list giving the number of elements in each ring (will be zero for a gap). The first entry is the number of radial elements outside the inner quadrilateral portion of the fuel mesh.
C++ Type:std::vector<unsigned int>
Controllable:No
Description:A list giving the number of elements in each ring (will be zero for a gap). The first entry is the number of radial elements outside the inner quadrilateral portion of the fuel mesh.
num_sectorsNumber of azimuthal sectors in each quadrant. 'num_sectors' must be an even number.
C++ Type:unsigned int
Controllable:No
Description:Number of azimuthal sectors in each quadrant. 'num_sectors' must be an even number.

Required Parameters

all_bottom_leftFalseWhether to generate sidesets for every bottom and left surface of mesh blocks.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to generate sidesets for every bottom and left surface of mesh blocks.
axial_directionzThe direction perpendicular to the circular cross section.
Default:z
C++ Type:MooseEnum
Options:x, y, z
Controllable:No
Description:The direction perpendicular to the circular cross section.
axial_offset0The offset of the generated mesh in the axial direction of the fuel rod from the origin.
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The offset of the generated mesh in the axial direction of the fuel rod from the origin.
biasMesh bias. Same length as elements_per_ring, each entry b: 0.5 <= b <= 2.0.
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Mesh bias. Same length as elements_per_ring, each entry b: 0.5 <= b <= 2.0.
compress_boundary_idsFalseWhether to compress the boundary IDs into a consecutive list.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to compress the boundary IDs into a consecutive list.
dual_biasMesh dual bias. If given, bias must be given. Same length as bias, each entry b: 0.5 <= b <= 2.0.
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Mesh dual bias. If given, bias must be given. Same length as bias, each entry b: 0.5 <= b <= 2.0.
elem_typequad8Whether to generate quad8 or quad4 elements.
Default:quad8
C++ Type:MooseEnum
Options:quad4, quad8
Controllable:No
Description:Whether to generate quad8 or quad4 elements.
name_prefixIf provided, prefix the built in boundary names with this string
C++ Type:std::string
Controllable:No
Description:If provided, prefix the built in boundary names with this string
offset0 0 The in-plane offset of the fuel (to generate eccentric fuel.
Default:0 0
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:The in-plane offset of the fuel (to generate eccentric fuel.
portiontop_rightControl which part of the cross section is meshed.
Default:top_right
C++ Type:MooseEnum
Options:full, top_right, right_half, top_half
Controllable:No
Description:Control which part of the cross section is meshed.

Optional Parameters

enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.
save_with_nameKeep the mesh from this mesh generator in memory with the name specified
C++ Type:std::string
Controllable:No
Description:Keep the mesh from this mesh generator in memory with the name specified

Advanced Parameters

nemesisFalseWhether or not to output the mesh file in the nemesisformat (only if output = true)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file in the nemesisformat (only if output = true)
outputFalseWhether or not to output the mesh file after generating the mesh
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file after generating the mesh
show_infoFalseWhether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)

Debugging Parameters

Input Files

(test/tests/triso/pyc_elasticity_tensor/anisotropic_exact.i)
(test/tests/triso_pebble/sphere_2d_rz.i)
(test/tests/solid_mechanics/graphite_grade_elasticity_tensor/test.i)
(test/tests/circular_cross_section_mesh/circular_cross_section_mesh_generator.i)

Child Objects

(include/meshgenerators/MPSCircularCrossSectionMeshGenerator.h)

References

No citations exist within this document.

(test/tests/circular_cross_section_mesh/circular_cross_section_mesh_generator.i)

#
# Test of capability to generate cross-section meshes.
#
# Volume check with r_fuel = 0.0041, r_clad_inner = 0.00418,
#   r_clad_outer = 0.00474 for a quarter of a fuel rod:
#
#                Fuel         Gap          Total
#        BISON:  1.320254e-5  5.202477e-7  1.764601e-5
#   Analytical:  1.320254e-5  5.202477e-7  1.764601e-5
#

[Mesh]
  [ccsmg]
    type = CircularCrossSectionMeshGenerator
    num_sectors = 20
    elements_per_ring = '5 0 3'
    block_names = 'fuel null clad'
    coordinates = '0.0041 0.00418 0.00474'
  []
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = 9
    value = 100
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = 2
    value = 200
  []
[]

[Executioner]
  type = Steady
[]

[Postprocessors]
  [fuelVolume]
    type = InternalVolume
    boundary = 8
  []
  [gapVolume]
    type = InternalVolume
    boundary = 9
  []
  [totalVolume]
    type = InternalVolume
    boundary = 2
  []
[]

[Outputs]
  [out]
    type = Exodus
  []
[]

(test/tests/triso/pyc_elasticity_tensor/anisotropic_exact.i)

# This test is to verify the implementation of PyCElasticityTensor material.

# It focuses on the full anisotropic elasticity tensor, including rotations.

# It is also used to verify exceptions and the isotropic and anisotropic Jacobians calculated by AD.

[GlobalParams]
  displacements = 'disp_x disp_y'
  flux_conversion_factor = 0.85
[]

[Mesh]
  coord_type = RZ
  use_displaced_mesh = false
  [mesh]
    type = CircularCrossSectionMeshGenerator
    num_sectors = 30
    offset = '0.0 0.0'
    elements_per_ring = '0 4 0'
    block_names = 'null PyC_quad null2'
    coordinates = '0.037 0.038 0.039'
  []
[]

[AuxVariables]
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [density]
    order = CONSTANT
    family = MONOMIAL
  []
  [temperature]
    initial_condition = 673.15
  []
[]

[Functions]
  [flux_history]
    type = PiecewiseLinear
    x = '0    1e4    1e8'
    y = '0 7.5e18 7.5e18'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [perm_PyC]
    block = PyC_quad
    add_variables = true
    strain = SMALL
    generate_output = 'vonmises_stress stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz'
  []
[]

[AuxKernels]
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_end
  []
  [density]
    type = MaterialRealAux
    variable = density
    property = density
    block = PyC_quad
    execute_on = 'initial linear'
  []
[]

[BCs]
  [no_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0
  []
  [no_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0
  []
  [Pressure]
    [inside_pressure]
      boundary = 1001
      factor = 1e4
    []
  []
[]

[Materials]
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    flux_function = flux_history
  []
  [stress]
    type = ComputeLinearElasticStress
    block = PyC_quad
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
    block = PyC_quad
  []
  [elasticity_tensor]
    type = PyCElasticityTensor
    block = PyC_quad
    initial_BAF = 1.05
    temperature = temperature
  []
  [density]
    type = StrainAdjustedDensity
    block = PyC_quad
    strain_free_density = 1900.0
  []
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
  petsc_options_value = 'lu       superlu_dist                  51'
  line_search = 'none'

  l_tol = 1e-5
  nl_max_its = 10
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-10
  start_time = 0.0
  num_steps = 1
  dt = 1e3
[]

[Postprocessors]
  [temperature]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fluence]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = fast_neutron_fluence
    execute_on = 'initial timestep_end'
  []
  [sigma_x]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = stress_xx
    execute_on = 'initial timestep_end'
  []
  [sigma_y]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = stress_yy
    execute_on = 'initial timestep_end'
  []
  [sigma_z]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = stress_zz
    execute_on = 'initial timestep_end'
  []
  [strain_x]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = strain_xx
    execute_on = 'initial timestep_end'
  []
  [strain_y]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = strain_yy
    execute_on = 'initial timestep_end'
  []
  [strain_z]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = strain_zz
    execute_on = 'initial timestep_end'
  []
  [density]
    type = ElementExtremeValue
    value_type = 'max'
    block = PyC_quad
    variable = density
    execute_on = 'initial timestep_end'
  []
[]

[Outputs]
  csv = true
  exodus = true
[]

(test/tests/triso_pebble/sphere_2d_rz.i)

[Mesh]
  coord_type = RZ
  [gen]
    type = CircularCrossSectionMeshGenerator
    num_sectors = 20
    elements_per_ring = '5 3'
    block_names = 'core fuel'
    coordinates = '0.75 1.0'
    portion = right_half
  []
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
  [time]
    type = TimeDerivative
    variable = u
  []
[]

[BCs]
  [outer]
    type = DirichletBC
    variable = u
    boundary = outer
    value = 0
  []
[]

[DiracKernels]
  [vpp_point_source]
    type = TRISOMonteCarloPointSource
    variable = u
    point_source_values = constant
    value_name = source
    point_source_location = point_source_location
  []
[]

[UserObjects]
  [point_source_location]
    type = TRISOMonteCarloPointSourceLocation
    min_radius = 0.75
    max_radius = 1.0
    geometry = SPHERE
    sampler = sample
    fuel_element_blocks = 2
    particle_space = 0.05
    execute_on = 'INITIAL'
  []
[]

[Distributions]
  [random]
    type = TruncatedNormal
    mean = 0
    standard_deviation = 0.5
  []
[]

[Samplers]
  [sample]
    type = MonteCarlo
    num_rows = 10
    distributions = 'random'
    execute_on = 'PRE_MULTIAPP_SETUP'
  []
[]

[VectorPostprocessors]
  [constant]
    type = ConstantVectorPostprocessor
    value = '1 2 3 4 5 6 7 8 9 10'
    vector_names = source
    execute_on = initial
  []
  [point_source_output]
    type = TRISOMonteCarloPointSourceOutput
    vector_postprocessor = constant
    value_name = source
    point_source_location = point_source_location
  []
[]

[Executioner]
  type = Transient
  solve_type = 'NEWTON'

  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'

  line_search = 'none'

  nl_rel_tol = 1e-6
  nl_abs_tol = 1e-6
  nl_max_its = 20

  l_tol = 1e-4
  l_max_its = 50

  start_time = 0.0
  end_time = 4.831315e7

  num_steps = 1

  dt = 400000
[]

[Outputs]
  csv = true
  exodus = true
[]

(test/tests/solid_mechanics/graphite_grade_elasticity_tensor/test.i)

# This test case is prepared to test the elastic properties of
# grade H-451 grade material model.

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  coord_type = RZ
  use_displaced_mesh = false
  [mesh]
    type = CircularCrossSectionMeshGenerator
    num_sectors = 30
    offset = '0.0 0.0'
    elements_per_ring = '0 4 0'
    block_names = 'null shell null2'
    coordinates = '0.039 0.064 0.065'
  []
[]

[AuxVariables]
  [temperature]
    initial_condition = 673.15
  []
  [fast_neutron_fluence]
    order = CONSTANT
    family = MONOMIAL
  []
  [density]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Functions]
  [temp_fcn]
    type = ParsedFunction
    expression = '600 + 0.075 * t'
  []
  [flux_history]
    type = PiecewiseLinear
    x = '0   1e4'
    y = '0 16e21'
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [perm]
    add_variables = true
    strain = SMALL
    generate_output = 'vonmises_stress stress_xx stress_yy stress_zz strain_xx strain_yy strain_zz'
  []
[]

[AuxKernels]
  [temp_ramp]
    type = FunctionAux
    variable = temperature
    function = temp_fcn
  []
  [fast_neutron_fluence]
    type = MaterialRealAux
    variable = fast_neutron_fluence
    property = fast_neutron_fluence
    execute_on = timestep_end
  []
  [density]
    type = MaterialRealAux
    variable = density
    property = density
    execute_on = 'initial linear'
  []
[]

[BCs]
  [no_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0
  []
  [no_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0
  []
  [Pressure]
    [inside_pressure]
      boundary = 1001
      factor = 1e4
    []
  []
[]

[Materials]
  [fast_neutron_flux]
    type = FastNeutronFlux
    calculate_fluence = true
    flux_function = flux_history
  []
  [stress]
    type = ComputeLinearElasticStress
  []
  [normal_vectors_triso]
    type = NormalVectorsTRISO
  []
  [elasticity_tensor]
    type = GraphiteGradeElasticityTensor
    temperature = temperature
    graphite_grade = H_451
    fast_neutron_fluence = fast_neutron_fluence
    flux_conversion_factor = 0.85
  []
  [thermal_props]
    type = GraphiteMatrixThermal
    graphite_grade = H_451
    temperature = temperature
    packing_fraction = 0.0
    specific_heat_scale_factor = 1.0
    thermal_conductivity_scale_factor = 1.0
    flux_conversion_factor = 0.85
  []
  [density]
    type = StrainAdjustedDensity
    strain_free_density = 1900.0
  []
[]

[Postprocessors]
  [temperature]
    type = ElementExtremeValue
    value_type = 'max'
    variable = temperature
    execute_on = 'initial timestep_end'
  []
  [fluence]
    type = ElementExtremeValue
    value_type = 'max'
    variable = fast_neutron_fluence
    execute_on = 'initial timestep_end'
  []
  [sigma_x]
    type = ElementExtremeValue
    value_type = 'max'
    variable = stress_xx
    execute_on = 'initial timestep_end'
  []
  [sigma_y]
    type = ElementExtremeValue
    value_type = 'max'
    variable = stress_yy
    execute_on = 'initial timestep_end'
  []
  [sigma_z]
    type = ElementExtremeValue
    value_type = 'max'
    variable = stress_zz
    execute_on = 'initial timestep_end'
  []
  [strain_x]
    type = ElementExtremeValue
    value_type = 'max'
    variable = strain_xx
    execute_on = 'initial timestep_end'
  []
  [strain_y]
    type = ElementExtremeValue
    value_type = 'max'
    variable = strain_yy
    execute_on = 'initial timestep_end'
  []
  [strain_z]
    type = ElementExtremeValue
    value_type = 'max'
    variable = strain_zz
    execute_on = 'initial timestep_end'
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -ksp_gmres_restart'
  petsc_options_value = 'lu       superlu_dist                  51'
  line_search = 'none'

  l_tol = 1e-5
  nl_max_its = 10
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-10
  start_time = 0.0
  num_steps = 10
  dt = 1e3
[]

[Outputs]
  csv = true
[]

(test/tests/circular_cross_section_mesh/circular_cross_section_mesh_generator.i)

#
# Test of capability to generate cross-section meshes.
#
# Volume check with r_fuel = 0.0041, r_clad_inner = 0.00418,
#   r_clad_outer = 0.00474 for a quarter of a fuel rod:
#
#                Fuel         Gap          Total
#        BISON:  1.320254e-5  5.202477e-7  1.764601e-5
#   Analytical:  1.320254e-5  5.202477e-7  1.764601e-5
#

[Mesh]
  [ccsmg]
    type = CircularCrossSectionMeshGenerator
    num_sectors = 20
    elements_per_ring = '5 0 3'
    block_names = 'fuel null clad'
    coordinates = '0.0041 0.00418 0.00474'
  []
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = 9
    value = 100
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = 2
    value = 200
  []
[]

[Executioner]
  type = Steady
[]

[Postprocessors]
  [fuelVolume]
    type = InternalVolume
    boundary = 8
  []
  [gapVolume]
    type = InternalVolume
    boundary = 9
  []
  [totalVolume]
    type = InternalVolume
    boundary = 2
  []
[]

[Outputs]
  [out]
    type = Exodus
  []
[]

(include/meshgenerators/MPSCircularCrossSectionMeshGenerator.h)

/*************************************************/
/*           DO NOT MODIFY THIS HEADER           */
/*                                               */
/*                     BISON                     */
/*                                               */
/*    (c) 2015 Battelle Energy Alliance, LLC     */
/*            ALL RIGHTS RESERVED                */
/*                                               */
/*   Prepared by Battelle Energy Alliance, LLC   */
/*     Under Contract No. DE-AC07-05ID14517      */
/*     With the U. S. Department of Energy       */
/*                                               */
/*     See COPYRIGHT for full restrictions       */
/*************************************************/

#pragma once

#include "CircularCrossSectionMeshGenerator.h"

/**
 * Generates a circular cross section mesh of an axisymmetric fuel rod with an MPS.
 */
class MPSCircularCrossSectionMeshGenerator : public CircularCrossSectionMeshGenerator
{
public:
  static InputParameters validParams();
  MPSCircularCrossSectionMeshGenerator(const InputParameters & parameters);

  std::unique_ptr<MeshBase> generate() override;
};

Description
Input Parameters
Input Files
Child Objects
References