Skip to content

Commit

Permalink
check: generation of hexahedral meshes
Browse files Browse the repository at this point in the history
  • Loading branch information
schuettem committed Apr 22, 2024
1 parent 7ccec87 commit 5840ddc
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 3 deletions.
4 changes: 2 additions & 2 deletions docs/documentation/tutorials/generationofhexahedralmeshes.md
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
# Generation of Hexahedral Meshes

As many solvers require purely hexahedral meshes, HOPR implements a subdivision strategy to split meshes consisting of tetrahedra, prisms and hexahedra into purely hexahedral meshes. This feature is activated using the parameter `splitToHex=T`. Note, that pyramids cannot be decomposed to hexahedra in a straightforward way, thus this feature cannot be applied to meshes containing pyramids. Note also that this option is up to now restricted to linear meshes.
Since many solvers require purely hexahedral meshes, HOPR implements a subdivision strategy to split meshes consisting of tetrahedra, prisms and hexahedra into purely hexahedral meshes. This feature is activated using the parameter `splitToHex=T`. Note, that pyramids cannot be decomposed to hexahedra in a straightforward way, thus this feature cannot be applied to meshes containing pyramids. Also note that this option is currently limited to linear meshes.

<table align="center" style="width:100%">
<tr>
Expand All @@ -25,6 +25,6 @@ As many solvers require purely hexahedral meshes, HOPR implements a subdivision

<h4>Parameter File<a class="headerlink" href="#parameter-file" title="Permalink to this heading"></a></h4>

To test this feature, set elemtype to either 104 or 106 and add `splitToHex=T` to the parameter file, which is found in
To test this feature, set `elemtype` to either `104` (Tetrahedron) or `106` (Prism with triangular base) and add `splitToHex=T` to the parameter file, which is found in

tutorials/1-01-cartbox/parameter.ini
2 changes: 1 addition & 1 deletion docs/documentation/userguide/parameters.md
Original file line number Diff line number Diff line change
Expand Up @@ -53,7 +53,7 @@ width: 50px
| `R_INF` | `R_INF=20` | Real | 1 | MANDATORY<br>if Mode=11 | Outer radius of curved structured mesh |
| `SpaceQuandt` | `SpaceQuandt=1.0` | Real | 1 | 0.1 | Characteristic length of the mesh |
| `SplitElemFile` | `SplitElemFile=filename` | Str | 1 | MANDATORY<br>if curvingMethod=3 | Name of suvdivided surface mesh |
| `SplitToHex` | `SplitToHex=T` | Logical | 1 | F | If `SplitToHex=F` then tetrahedra and prisms are subdivided to yield pure hexahedral meshes meshes. <br>Hexahedra are also subdivided to guarantee a conforming connection. <br>NOTE: Pyramids cannot be split, this function cannot be used with pyramids present in the mesh. |
| `SplitToHex` | `SplitToHex=T` | Logical | 1 | F | If `SplitToHex=F` then tetrahedra and prisms are subdivided to yield pure hexahedral meshes. <br>Hexahedra are also subdivided to guarantee a conforming connection. <br>NOTE: Pyramids cannot be split, this function cannot be used with pyramids present in the mesh. |
| `stretchType` | `stretchType=(/3,1,0/)` | Int | 3 | (/0,0,0/) | This parameter manages the (de)activation of the stretching functions for all axis. For this reason the parameter <br>is a vector with three components.<br>0: Stretching in direction of the axis is deactivated<br>1: Stretching in direction of the axis is activated<br>3: Stretching in direction of the axis is activated and from the second half of the mesh distance on the stretching <br>factor is multiplied by -1. |
| `useCurveds` | `useCurveds=T` | Logical | 1 | F | T (True): If curved boundaries are defined<br>F (False): If no curved boundaries are defined |
| `vv` | `vv=(/0,0,1./)` | Real | 3 | (/0,0,0/) | The displacement vector has to specify in the three-dimensional cartesian coordinate system and has to be <br>normal to a surface the vector was assigned to. In addition the displacement vector has to show to the inside <br>of the cartesian box. In case of two parallel surface-planes, both with periodic boundary conditions, just one <br>displacement vector has to be defined. Therefore the different directions of the vectors can be compensated <br>by switching the sign of alpha, the fourth component of the `BoundaryType` vector.<br>The displacement vector has to be as long as the distance between the surfaces the vector was assigned to. <br>The index of a displacement vector is defined by the position of its definition like the parameter `BCIndex`. <br>Several definitions of boundary conditions between two definitions of displacement vectors will not affect <br>the index of the displacement vectors. |
Expand Down

0 comments on commit 5840ddc

Please sign in to comment.