SELF_Mesh_3D_t Module

Components

Type-Bound Procedures

Evaluates the bilinear map of face s of a trilinear (nGeo=1) element at the face coordinates (u,v) in [-1,1]^2, using the face trace-coordinate convention of BoundaryInterp: faces 1,6 -> (xi1,xi2); faces 2,4 -> (xi1,xi3); faces 3,5 -> (xi2,xi3).

Arguments

Return Value real(kind=prec)(1:3)

Identifies which of the six domain boundary planes of the SimpleMortarMesh contains face s of the element, and returns the corresponding boundary condition id. The domain is [0,3dx] x [0,2dx] x [0,2dx] with the region x > 2dx, above/below the small elements, outside the mesh; the x = 3dx plane is the domain east.

Arguments

Return Value integer

Identifies which of the six planes of the box [0,extent(1)] x [0,extent(2)] x [0,extent(3)] contains face s of the element, and returns the corresponding boundary condition id (StructuredMesh side ordering).

Arguments

Return Value integer

Determines whether face sSmall of the small element coincides with half h of face sBig of the big element, split along big-face trace coordinate split (1 or 2), under the given flip: face_small(F_flip(u,v)) == face_big(half_h(u,v)) at the sampled points.

Arguments

Return Value logical

Determines whether face sSmall of the small element coincides with quadrant q of face sBig of the big element under the given flip: face_small(F_flip(u,v)) == face_big(quadrant_q(u,v)) at the sampled points.

Arguments

Return Value logical

Arguments

Return Value integer

Map a global element position (gx,gy,gz), with gx in [1,nX] etc., to the tile-ordered element id used throughout the structured mesh. This is the inverse of the (i,ti) decomposition: gx = i + nxpertile*(ti-1).

Arguments

Return Value integer

Fills the level heights for ExtrudedMesh with nzFlat = size(depths) flat interfaces near the surface over nzS = size(sigma)-1 terrain-following layers. With d_0 = 0 and d_F = depths(nzFlat) (d_F = 0 when nzFlat = 0), the nz = nzS+nzFlat levels at every geometry node are

Arguments

Builds a layered hexahedral mesh by extruding a surface mesh downward through nz = size(zLevels,3)-1 layers bounded by the level surfaces z_0 < z_1 < ... < z_nz.

Arguments

Fills the level heights for ExtrudedMesh with terrain-following (sigma) levels,

Arguments

Coordinate form of MortarFaceMap: maps receiver-face coordinates (u,v) to donor-face coordinates (u2,v2) for each of the eight SELF face flips.

Arguments

Arguments

Returns the nGeo+1 one-dimensional geometry node positions, in [-1,1], of a mesh whose node coordinates are sampled on quadrature points. These are the points GenerateFromMesh interpolates the node coordinates from (see Lagrange_t%Init), so a coordinate sampled here is reproduced exactly by the geometry interpolant.

Arguments

Arguments

Determines whether face sA of element A coincides with face sB of element B, and if so with which flip: face_B(F_flip(u,v)) == face_A(u,v) at the sampled points. Bilinear faces are pinned by their four corners plus one interior point.

Arguments

Maps receiver-face node indices (i,j) to donor-face node indices (i2,j2) for each of the eight SELF face flips (see the CGNStoSELFflip table above; indices here are 1-based, so the 0-based rule i2 = N-i1 reads i2 = N+2-i). Mortar staging uses this to reorient a small face's trace into the big face's coordinates and to scatter restricted traces back.

Arguments

Reader for HOHQMesh 3-D (hexahedral) text mesh files in the ISM and ISM-MM formats, as written by HOHQMesh's WriteISMHexMeshFile (Source/3DSource/Mesh3DOutputMethods.f90). Unlike the 2-D writer, the 3-D writer emits NO format header line: the first line is always the count line "nNodes nElems polyOrder". The two variants differ only in the per-element corner-node line: * ISM : 8 corner-node ids * ISM-MM : 8 corner-node ids followed by a material-name string The variant is auto-detected from the presence of the 9th token.

Arguments

Arguments

Arguments

This method can be used to reset all of the boundary elements boundary condition type to the desired value.

Arguments

Create the smallest 3D anisotropic (2-subface) nonconforming mesh: one "big" element whose east face is shared with the faces of two "small" dx x dx x dx elements that split it along one face coordinate only. The mesh is conforming everywhere except at the single mortar interface.

Arguments

Create the smallest 3D 2:1 nonconforming (mortar) mesh: one "big" element of size 2dx x 2dx x 2*dx whose east face is shared with the faces of four "small" dx x dx x dx elements stacked 2x2 in (y,z). The mesh is conforming everywhere except at the single mortar interface.

Arguments

Create a fully triply-periodic structured hexahedral mesh and store it in SELF's unstructured mesh format. Element geometry and ordering are identical to UniformStructuredMesh; the only difference is connectivity. The faces on the six domain boundaries are wired as interior faces whose neighbor is the element on the opposite side of the domain. This realises the triply periodic box T^3 = [0,Lx] x [0,Ly] x [0,Lz] required by, e.g., the Arnold-Beltrami-Childress (ABC) flow benchmark.

Arguments

Create a structured mesh and store it in SELF's unstructured mesh format. The mesh is created in tiles of size (tnx,tny,tnz). Tiling is used to determine the element ordering.

Arguments

Arguments

Arguments

Fill mesh%nodeCoords(:,:,:,:,e) for one hexahedral element by transfinite (Coons) interpolation of its six face grids. Faces flagged in the mesh file use the file's face-point grids; the remaining faces are bilinear patches of their four corner nodes evaluated at Chebyshev-Gauss-Lobatto parametric coordinates. Edge curves are extracted from the face grids (preferring a flagged face when an edge borders one flagged and one bilinear face) and the standard Boolean-sum formula x = Pxi + Peta + Pzeta - PxiPeta - PxiPzeta - PetaPzeta + PxiPeta*Pzeta combines face, edge, and corner contributions. For an element with all-straight faces this reduces to trilinear interpolation of the eight corners.

Arguments