AnnularMesh

For rmin>0: creates an annular mesh of QUAD4 elements. For rmin=0: creates a disc mesh of QUAD4 and TRI3 elements. Boundary sidesets are created at rmax and rmin, and given these names. If dmin!0 and dmax!360, a sector of an annulus or disc is created. In this case boundary sidesets are also created a dmin and dmax, and given these names

Description

The AnnularMesh mesh generator builds simple 2D annular and disc meshes. They are created by drawing radial lines and concentric circles, and the mesh consists of the quadrilaterals thus formed. Therefore, no sophisticated paving is used to construct the mesh.

The inner radius and the outer radius must be specified. If the inner radius is zero a disc mesh is created, while if it is positive an annulus is created. The annulus has just one subdomain (block number = 0), whereas the disc has two subdomains: subdomain zero consists of the outer quadrilaterals, while the other (block number = 1) consists of the triangular elements that emanate from the origin.

The minimum and maximum angle may also be specified. These default to zero and 360, respectively. If other values are chosen, a sector of an annulus, or a sector of a disc will be created. Both angles are measured anti-clockwise from the $var element = document.getElementById("moose-equation-9a790d84-453e-4e63-9394-ac3deb36a703");katex.render("x", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ axis.

The number of elements in the radial direction and the angular direction may be specified. In addition, a growth factor on the element size in the radial direction may be chosen. The element-size (in the radial direction) is multiplied by this factor for each concentric ring of elements, moving from the inner to the outer radius.

Sidesets are also created:

* Sideset 0 is called "rmin" and is the set of sides at the minimum radius (which is zero for the disc). * Sideset 1 is called "rmax" and is the set of sides at the maximum radius. * Sideset 2 is called "dmin" and is the set of sides at the minimum angle, which is created only in the case of a sector of an annulus (or disc) * Sideset 3 is called "dmax" and is the set of sides at the maximum angle, which is created only in the case of a sector of an annulus (or disc)

Example Syntax

A full annulus with minimum radius 1 and maximum radius 5, with smaller elements near the inside of the annulus. (A disc would be created by setting rmin to zero.) !listing test/tests/mesh/mesh_generation/annulus.i block=Mesh

A sector of an annulus, sitting between 45 and 135 degrees. (A sector of a disc would be created by setting rmin to zero.) !listing test/tests/mesh/mesh_generation/annulus_sector.i block=Mesh

An example of using sidesets !listing test/tests/mesh/mesh_generation/annulus_sector.i block=BCs

Input Parameters

ntNumber of elements in the angular direction
C++ Type:unsigned int
Options:
Description:Number of elements in the angular direction
rmaxOuter radius
C++ Type:double
Options:
Description:Outer radius
rminInner radius. If rmin=0 then a disc mesh (with no central hole) will be created.
C++ Type:double
Options:
Description:Inner radius. If rmin=0 then a disc mesh (with no central hole) will be created.

Required Parameters

allow_renumberingTrueIf allow_renumbering=false, node and element numbers are kept fixed until deletion
Default:True
C++ Type:bool
Options:
Description:If allow_renumbering=false, node and element numbers are kept fixed until deletion
dmax360Maximum angle, measured in degrees anticlockwise from x axis. If dmin=0 and dmax=360 an annular mesh is created. Otherwise, only a sector of an annulus is created
Default:360
C++ Type:double
Options:
Description:Maximum angle, measured in degrees anticlockwise from x axis. If dmin=0 and dmax=360 an annular mesh is created. Otherwise, only a sector of an annulus is created
dmin0Minimum degree, measured in degrees anticlockwise from x axis
Default:0
C++ Type:double
Options:
Description:Minimum degree, measured in degrees anticlockwise from x axis
ghosting_patch_sizeThe number of nearest neighbors considered for ghosting purposes when 'iteration' patch update strategy is used. Default is 5 * patch_size.
C++ Type:unsigned int
Options:
Description:The number of nearest neighbors considered for ghosting purposes when 'iteration' patch update strategy is used. Default is 5 * patch_size.
growth_r1The ratio of radial sizes of successive rings of elements
Default:1
C++ Type:double
Options:
Description:The ratio of radial sizes of successive rings of elements
max_leaf_size10The maximum number of points in each leaf of the KDTree used in the nearest neighbor search. As the leaf size becomes larger,KDTree construction becomes faster but the nearest neighbor searchbecomes slower.
Default:10
C++ Type:unsigned int
Options:
Description:The maximum number of points in each leaf of the KDTree used in the nearest neighbor search. As the leaf size becomes larger,KDTree construction becomes faster but the nearest neighbor searchbecomes slower.
nr1Number of elements in the radial direction
Default:1
C++ Type:unsigned int
Options:
Description:Number of elements in the radial direction
parallel_typeDEFAULTDEFAULT: Use libMesh::ReplicatedMesh unless --distributed-mesh is specified on the command line REPLICATED: Always use libMesh::ReplicatedMesh DISTRIBUTED: Always use libMesh::DistributedMesh
Default:DEFAULT
C++ Type:MooseEnum
Options:DEFAULT REPLICATED DISTRIBUTED
Description:DEFAULT: Use libMesh::ReplicatedMesh unless --distributed-mesh is specified on the command line REPLICATED: Always use libMesh::ReplicatedMesh DISTRIBUTED: Always use libMesh::DistributedMesh
quad_subdomain_id0The subdomain ID given to the QUAD4 elements
Default:0
C++ Type:unsigned short
Options:
Description:The subdomain ID given to the QUAD4 elements
tri_subdomain_id1The subdomain ID given to the TRI3 elements (these exist only if rmin=0, and they exist at the center of the disc
Default:1
C++ Type:unsigned short
Options:
Description:The subdomain ID given to the TRI3 elements (these exist only if rmin=0, and they exist at the center of the disc

Optional Parameters

centroid_partitioner_directionSpecifies the sort direction if using the centroid partitioner. Available options: x, y, z, radial
C++ Type:MooseEnum
Options:x y z radial
Description:Specifies the sort direction if using the centroid partitioner. Available options: x, y, z, radial
partitionerdefaultSpecifies a mesh partitioner to use when splitting the mesh for a parallel computation.
Default:default
C++ Type:MooseEnum
Options:default metis parmetis linear centroid hilbert_sfc morton_sfc
Description:Specifies a mesh partitioner to use when splitting the mesh for a parallel computation.

Partitioning Parameters

construct_node_list_from_side_listTrueWhether or not to generate nodesets from the sidesets (usually a good idea).
Default:True
C++ Type:bool
Options:
Description:Whether or not to generate nodesets from the sidesets (usually a good idea).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
dim1This is only required for certain mesh formats where the dimension of the mesh cannot be autodetected. In particular you must supply this for GMSH meshes. Note: This is completely ignored for ExodusII meshes!
Default:1
C++ Type:MooseEnum
Options:1 2 3
Description:This is only required for certain mesh formats where the dimension of the mesh cannot be autodetected. In particular you must supply this for GMSH meshes. Note: This is completely ignored for ExodusII meshes!
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
nemesisFalseIf nemesis=true and file=foo.e, actually reads foo.e.N.0, foo.e.N.1, ... foo.e.N.N-1, where N = # CPUs, with NemesisIO.
Default:False
C++ Type:bool
Options:
Description:If nemesis=true and file=foo.e, actually reads foo.e.N.0, foo.e.N.1, ... foo.e.N.N-1, where N = # CPUs, with NemesisIO.
patch_size40The number of nodes to consider in the NearestNode neighborhood.
Default:40
C++ Type:unsigned int
Options:
Description:The number of nodes to consider in the NearestNode neighborhood.
patch_update_strategyneverHow often to update the geometric search 'patch'. The default is to never update it (which is the most efficient but could be a problem with lots of relative motion). 'always' will update the patch for all slave nodes at the beginning of every timestep which might be time consuming. 'auto' will attempt to determine at the start of which timesteps the patch for all slave nodes needs to be updated automatically.'iteration' updates the patch at every nonlinear iteration for a subset of slave nodes for which penetration is not detected. If there can be substantial relative motion between the master and slave surfaces during the nonlinear iterations within a timestep, it is advisable to use 'iteration' option to ensure accurate contact detection.
Default:never
C++ Type:MooseEnum
Options:never always auto iteration
Description:How often to update the geometric search 'patch'. The default is to never update it (which is the most efficient but could be a problem with lots of relative motion). 'always' will update the patch for all slave nodes at the beginning of every timestep which might be time consuming. 'auto' will attempt to determine at the start of which timesteps the patch for all slave nodes needs to be updated automatically.'iteration' updates the patch at every nonlinear iteration for a subset of slave nodes for which penetration is not detected. If there can be substantial relative motion between the master and slave surfaces during the nonlinear iterations within a timestep, it is advisable to use 'iteration' option to ensure accurate contact detection.

Advanced Parameters

Input Files

test/tests/mesh/mesh_generation/disc.i
test/tests/mesh/mesh_generation/annular_except2.i
test/tests/mesh/mesh_generation/annular_except5_deprecated.i
test/tests/mesh/mesh_generation/annular_except6.i
test/tests/mesh/mesh_generation/annular_except5.i
test/tests/mesh/mesh_generation/annular_except2_deprecated.i
test/tests/mesh/mesh_generation/annular_except3_deprecated.i
test/tests/mesh/mesh_generation/annulus.i
test/tests/mesh/mesh_generation/annulus_sector.i
test/tests/mesh/mesh_generation/annular_except1.i
test/tests/mesh/mesh_generation/annular_except3.i
test/tests/mesh/mesh_generation/annulus_sector_deprecated.i
test/tests/mesh/mesh_generation/disc_sector_deprecated.i
test/tests/mesh/mesh_generation/annular_except1_deprecated.i
test/tests/mesh/mesh_generation/annular_except4.i
test/tests/mesh/mesh_generation/disc_sector.i

test/tests/mesh/mesh_generation/disc.i

# Generates a Disc Mesh
# Radius of outside circle=5
# Solves the diffusion equation with u=-5 at origin, and u=0 on outside
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 0
  rmax = 5
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = -5.0
    boundary = rmin
  [../]
  [./outer]
    type = DirichletBC
    variable = u
    value = 0.0
    boundary = rmax
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

test/tests/mesh/mesh_generation/annular_except2.i

# Exception testing of AnnularMesh.
# In this file dmax < dmin
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  dmin = 135
  dmax = 45
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except5_deprecated.i

# Exception testing of AnnularMesh.
# In this file nt < (tmax - tmin) / Pi
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 1
  rmin = 1
  rmax = 5
  tmax = 4
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except6.i

# Exception testing of AnnularMesh.
# In this file tri_subdomain_id = quad_subdomain_id
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 3
  rmin = 0
  rmax = 5
  quad_subdomain_id = 1
  tri_subdomain_id = 1
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except5.i

# Exception testing of AnnularMesh.
# In this file nt < (dmax - dmin) / Pi
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 1
  rmin = 1
  rmax = 5
  dmax = 230
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except2_deprecated.i

# Exception testing of AnnularMesh.
# In this file tmax < tmin
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  tmin = 2.356194490
  tmax = 0.785398163
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except3_deprecated.i

# Exception testing of AnnularMesh.
# In this file tmax - tmin > 2 Pi
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  tmin = 0.785398163
  tmax = 433
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annulus.i

# Generates an Annular Mesh
# Radius of inside circle=1
# Radius of outside circle=5
# Solves the diffusion equation with
# u=0 on inside
# u=log(5) on outside

[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = 0.0
    boundary = rmin
  [../]
  [./outer]
    type = FunctionDirichletBC
    variable = u
    function = log(5)
    boundary = rmax
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

test/tests/mesh/mesh_generation/annulus_sector.i

# Generates a sector of an Annular Mesh between angle=Pi/4 and angle=3Pi/4
# Radius of inside circle=1
# Radius of outside circle=5
# Solves the diffusion equation with
# u=0 on inside
# u=log(5) on outside
# u=log(r) at angle=Pi/4 and angle=3Pi/4

[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  dmin = 45
  dmax = 135
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = 0.0
    boundary = rmin
  [../]
  [./outer]
    type = FunctionDirichletBC
    variable = u
    function = log(5)
    boundary = rmax
  [../]
  [./min_angle]
    type = FunctionDirichletBC
    variable = u
    function = 'log(sqrt(x*x + y*y))'
    boundary = 'dmin dmax'
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

test/tests/mesh/mesh_generation/annular_except1.i

# Exception testing of AnnularMesh.
# In this file rmax < rmin
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 5
  rmax = 1
  dmin = 45
  dmax = 135
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except3.i

# Exception testing of AnnularMesh.
# In this file dmax - dmin > 360
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  dmin = 45
  dmax = 422
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annulus_sector_deprecated.i

# Generates a sector of an Annular Mesh between angle=Pi/4 and angle=3Pi/4
# Radius of inside circle=1
# Radius of outside circle=5
# Solves the diffusion equation with
# u=0 on inside
# u=log(5) on outside
# u=log(r) at angle=Pi/4 and angle=3Pi/4

[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 1
  rmax = 5
  tmin = 0.785398163
  tmax = 2.356194490
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = 0.0
    boundary = rmin
  [../]
  [./outer]
    type = FunctionDirichletBC
    variable = u
    function = log(5)
    boundary = rmax
  [../]
  [./min_angle]
    type = FunctionDirichletBC
    variable = u
    function = 'log(sqrt(x*x + y*y))'
    boundary = 'tmin tmax'
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

test/tests/mesh/mesh_generation/disc_sector_deprecated.i

# Generates a sector of a Disc Mesh between angle=Pi/4 and angle=3Pi/4
# Radius of outside circle=5
# Solves the diffusion equation with u=-5 at origin, and u=0 on outside
# as well as u=-5+r at angle=Pi/4 and u=-5+r^4/125 at angle=3Pi/4

[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 0
  rmax = 5
  tmin = 0.785398163
  tmax = 2.356194490
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = -5.0
    boundary = rmin
  [../]
  [./outer]
    type = FunctionDirichletBC
    variable = u
    function = 0
    boundary = rmax
  [../]
  [./tmin]
    type = FunctionDirichletBC
    variable = u
    function = '-5.0+sqrt(x*x + y*y)'
    boundary = tmin
  [../]
  [./tmax]
    type = FunctionDirichletBC
    variable = u
    function = '-5.0+pow(x*x + y*y, 2)/125'
    boundary = tmax
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

test/tests/mesh/mesh_generation/annular_except1_deprecated.i

# Exception testing of AnnularMesh.
# In this file rmax < rmin
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 5
  rmax = 1
  tmin = 45
  tmax = 135
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/annular_except4.i

# Exception testing of AnnularMesh.
# In this file nt < (dmax - dmin) / Pi
[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 2
  rmin = 1
  rmax = 5
  growth_r = 1.3
[]

test/tests/mesh/mesh_generation/disc_sector.i

# Generates a sector of a Disc Mesh between angle=Pi/4 and angle=3Pi/4
# Radius of outside circle=5
# Solves the diffusion equation with u=-5 at origin, and u=0 on outside
# as well as u=-5+r at angle=Pi/4 and u=-5+r^4/125 at angle=3Pi/4

[Mesh]
  type = AnnularMesh
  nr = 10
  nt = 12
  rmin = 0
  rmax = 5
  dmin = 45
  dmax = 135
  growth_r = 1.3
[]

[Variables]
  [./u]
  [../]
[]

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

[BCs]
  [./inner]
    type = DirichletBC
    variable = u
    value = -5.0
    boundary = rmin
  [../]
  [./outer]
    type = FunctionDirichletBC
    variable = u
    function = 0
    boundary = rmax
  [../]
  [./tmin]
    type = FunctionDirichletBC
    variable = u
    function = '-5.0+sqrt(x*x + y*y)'
    boundary = dmin
  [../]
  [./tmax]
    type = FunctionDirichletBC
    variable = u
    function = '-5.0+pow(x*x + y*y, 2)/125'
    boundary = dmax
  [../]
[]

[Executioner]
  type = Steady
  solve_type = PJFNK
[]

[Outputs]
  exodus = true
[]

Description
Example Syntax
Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools