Difference between revisions of "Boundary conditions"

From Serpent Wiki
Jump to: navigation, search
(Boundary conditions in 2D geometries)
m
 
(17 intermediate revisions by one other user not shown)
Line 1: Line 1:
= Boundary conditions =
+
Boundary conditions are invoked when the particle crosses the boundary from the defined part of the geometry into a region labeled as "outside" in the [[Input syntax manual#cell (cell definition)|cell cards]]. The outcome then depends on the boundary conditions, as defined by the "[[Input syntax manual#set bc|set bc]]" option. Serpent provides three types of boundary conditions:
  
== Boundary conditions in 2D geometries ==
+
# Black boundary - the particle is killed
[[File:vacu2D.png|thumb|300px|Vacuum boundary conditions ("set bc 1"). See [[media:vacu2D.gif|animated particle tracks]].]]
+
# Reflective boundary - the particle is reflected back into the geometry
[[File:refl2D.png|thumb|300px|Reflective boundary conditions ("set bc 2"). See [[media:refl2D.gif|animated particle tracks]].]]
+
# Periodic boundary - the particle is moved to the opposite side of the geometry
[[File:peri2D.png|thumb|300px|Periodic boundary conditions ("set bc 3"). See [[media:peri2D.gif|animated particle tracks]].]]
+
  
=== Vacuum boundary conditions ===
+
Reflective and periodic boundary conditions can be used to construct infinite and semi-infinite lattice structures. The way these boundary conditions are handled in Serpent is somewhat different from other Monte Carlo codes. Instead of stopping the particle at the boundary surface, reflections and translations are handled by coordinate transformations. This limits the outermost boundary to a few specific surface types that form regular lattices. In 2D geometries such surfaces are:
Vakuumi jnejne...  
+
  
=== Reflective boundary conditions ===
+
*Square prism ("sqc") - forms a 2D square lattice
 +
*Hexagonal prisms ("hexxc" and "hexxyc") -- form 2D hexagonal lattices
 +
*Rectangular prism ("rect") - forms a 2D rectangular lattice
  
=== Periodic boundary conditions ===
+
And in 3D:
 +
 
 +
*Cube ("cube") - forms a 3D cubical lattice
 +
*Cuboid ("cuboid") - forms a 3D cuboidal lattice
 +
*Truncated hexagonal prisms ("hexxprism" and "hexyprism") - form a 3D hexagonal lattices
 +
 
 +
Repeated boundary conditions are invoked only on the first surface listed in the "outside" type cells. Since the geometry boundary may consist of additional surfaces, some boundaries may be treated as reflective or periodic, while others are treated as vacuum. In reactor geometries there are typically three possible configurations:
 +
 
 +
;Infinite 2D geometry: The geometry has no dependence on the z-coordinate. The outer boundary is defined by one of the 2D surfaces listed above.
 +
;Radially infinite, axially finite 3D geometry:The outer boundary is defined by one of the 2D surfaces listed above and two axial planes ("pz"). The boundary condition takes effect in the radial direction only, while the axial boundaries are vacuum.
 +
;Infinite 3D geometry: The outer boundary is defined by one of the 2D surfaces listed above. The boundary condition takes effect in all directions.
 +
 
 +
These configurations are demonstrated below.
 +
 
 +
== Infinite 2D geometries ==
 +
 
 +
In infinite 2D lattice geometries the outside part of the geometry is composed of a single cell filling the space outside the surface at which the boundary conditions are applied. For example:
 +
 
 +
<nowiki>
 +
cell 999 0 outside  b1</nowiki>
 +
 
 +
where surface b1 is an infinite square prism.
 +
 
 +
Two types of infinite square lattices created using reflective and periodic boundary conditions are illustrated below, together with a non-repetitive black boundary. The geometry plot is extended beyond the outermost surface, in which case the plot shows the vacuum region as black and applies reflections and translations accordingly. The image captions provide links to particle track animations, which illustrate the behaviour when the boundary is hit. It should be noted that the boundary surface itself is not plotted, since the material is the same on both sides.
 +
 
 +
<gallery mode="nolines" widths=320px heights=320px>
 +
Image:vacu2D.png|'''Vacuum boundary conditions''' ("set&nbsp;bc&nbsp;1", xy-plot). See [[media:vacu2D.gif|animated particle tracks]].
 +
Image:refl2D.png|'''Reflective boundary conditions''' ("set&nbsp;bc&nbsp;2", xy-plot). See [[media:refl2D.gif|animated particle tracks]].
 +
Image:peri2D.png|'''Periodic boundary conditions''' ("set&nbsp;bc&nbsp;3", xy-plot). See [[media:peri2D.gif|animated particle tracks]].
 +
</gallery>
 +
 
 +
Serpent allows using reflective and periodic boundary conditions for both square and hexagonal prisms, although it should be noted that a reflective boundary with a non-symmetrical hexagonal cell makes physically no sense.
 +
 
 +
== Radially infinite, axially finite 3D geometries ==
 +
 
 +
Three-dimensional model of a fuel assembly in an infinite lattice is an example of a geometry that is radially infinite and axially finite. The outer boundary of such geometry is formed by the surface at which the repeated boundary conditions are applied, and two axial planes that the define the bottom and top of the geometry model. The outside cells could then be defined, for example, as:
 +
 
 +
<nowiki>
 +
cell 997 0 outside  b1  b2 -b3
 +
cell 998 0 outside -b2
 +
cell 999 0 outside  b3</nowiki>
 +
 
 +
where surfaces b1, b2 and b3 are an infinite square prism and planes forming the bottom and the top of the geometry, respectively. It is important that the surface at which the boundary conditions are applied is listed first.
 +
 
 +
<gallery mode="nolines" widths=320px heights=320px>
 +
Image:refl3Df.png|'''Reflective boundary conditions''' ("set&nbsp;bc&nbsp;2", yz-plot). See [[media:refl3Df.gif|animated particle tracks]].
 +
</gallery>
 +
 
 +
== Infinite 3D geometries ==
 +
 
 +
Infinite 3D lattices are created similar to infinite 2D lattices, by defining the outside space using a single cell and surface, for example:
 +
 
 +
<nowiki>
 +
cell 999 0 outside  b1</nowiki>
 +
 
 +
The difference to the 2D case is that surface b1 is one of the 3D surfaces listed above, in this case, a cube. The difference to the radially infinite, axially finite 3D case above is that the geometry is limited to a repeated boundary in all directions. Reflective boundary conditions are illustrated below.
 +
 
 +
<gallery mode="nolines" widths=320px heights=320px>
 +
Image:refl3Di.png|'''Reflective boundary conditions''' ("set&nbsp;bc&nbsp;2", yz-plot). See [[media:refl3Di.gif|animated particle tracks]].
 +
</gallery>
 +
 
 +
[[Category:Input]]
 +
[[Category:Theory]]

Latest revision as of 14:38, 29 November 2017

Boundary conditions are invoked when the particle crosses the boundary from the defined part of the geometry into a region labeled as "outside" in the cell cards. The outcome then depends on the boundary conditions, as defined by the "set bc" option. Serpent provides three types of boundary conditions:

  1. Black boundary - the particle is killed
  2. Reflective boundary - the particle is reflected back into the geometry
  3. Periodic boundary - the particle is moved to the opposite side of the geometry

Reflective and periodic boundary conditions can be used to construct infinite and semi-infinite lattice structures. The way these boundary conditions are handled in Serpent is somewhat different from other Monte Carlo codes. Instead of stopping the particle at the boundary surface, reflections and translations are handled by coordinate transformations. This limits the outermost boundary to a few specific surface types that form regular lattices. In 2D geometries such surfaces are:

  • Square prism ("sqc") - forms a 2D square lattice
  • Hexagonal prisms ("hexxc" and "hexxyc") -- form 2D hexagonal lattices
  • Rectangular prism ("rect") - forms a 2D rectangular lattice

And in 3D:

  • Cube ("cube") - forms a 3D cubical lattice
  • Cuboid ("cuboid") - forms a 3D cuboidal lattice
  • Truncated hexagonal prisms ("hexxprism" and "hexyprism") - form a 3D hexagonal lattices

Repeated boundary conditions are invoked only on the first surface listed in the "outside" type cells. Since the geometry boundary may consist of additional surfaces, some boundaries may be treated as reflective or periodic, while others are treated as vacuum. In reactor geometries there are typically three possible configurations:

Infinite 2D geometry
The geometry has no dependence on the z-coordinate. The outer boundary is defined by one of the 2D surfaces listed above.
Radially infinite, axially finite 3D geometry
The outer boundary is defined by one of the 2D surfaces listed above and two axial planes ("pz"). The boundary condition takes effect in the radial direction only, while the axial boundaries are vacuum.
Infinite 3D geometry
The outer boundary is defined by one of the 2D surfaces listed above. The boundary condition takes effect in all directions.

These configurations are demonstrated below.

Infinite 2D geometries

In infinite 2D lattice geometries the outside part of the geometry is composed of a single cell filling the space outside the surface at which the boundary conditions are applied. For example:

cell 999 0 outside  b1

where surface b1 is an infinite square prism.

Two types of infinite square lattices created using reflective and periodic boundary conditions are illustrated below, together with a non-repetitive black boundary. The geometry plot is extended beyond the outermost surface, in which case the plot shows the vacuum region as black and applies reflections and translations accordingly. The image captions provide links to particle track animations, which illustrate the behaviour when the boundary is hit. It should be noted that the boundary surface itself is not plotted, since the material is the same on both sides.

Serpent allows using reflective and periodic boundary conditions for both square and hexagonal prisms, although it should be noted that a reflective boundary with a non-symmetrical hexagonal cell makes physically no sense.

Radially infinite, axially finite 3D geometries

Three-dimensional model of a fuel assembly in an infinite lattice is an example of a geometry that is radially infinite and axially finite. The outer boundary of such geometry is formed by the surface at which the repeated boundary conditions are applied, and two axial planes that the define the bottom and top of the geometry model. The outside cells could then be defined, for example, as:

cell 997 0 outside  b1  b2 -b3
cell 998 0 outside -b2
cell 999 0 outside  b3

where surfaces b1, b2 and b3 are an infinite square prism and planes forming the bottom and the top of the geometry, respectively. It is important that the surface at which the boundary conditions are applied is listed first.

Infinite 3D geometries

Infinite 3D lattices are created similar to infinite 2D lattices, by defining the outside space using a single cell and surface, for example:

cell 999 0 outside  b1

The difference to the 2D case is that surface b1 is one of the 3D surfaces listed above, in this case, a cube. The difference to the radially infinite, axially finite 3D case above is that the geometry is limited to a repeated boundary in all directions. Reflective boundary conditions are illustrated below.