Mesh¶

class Mesh¶

Public Functions

Mesh(std::vector<Vector> nodes, std::vector<std::vector<unsigned>> cells, std::vector<std::pair<BOUNDARY_TYPE, std::vector<unsigned>>> boundaries)¶: Constructor of mesh. Needs nodes, cells, and boundary descriptions as specified above. See LoadSU2MeshFromFile in io.cpp for setup information.

const std::vector<Vector> &GetNodes() const¶

Returns all node coordinates.

Return: dimension: numNodes x dim

const std::vector<Vector> &GetCellMidPoints() const¶

Returns the mid point coordinates of each cell.

Return: dimension: numCells x dim

const std::vector<std::vector<unsigned>> &GetCells() const¶

Returns all node IDs that construct up each cell.

Return: dimension: numCells x numNodes

const std::vector<double> &GetCellAreas() const¶

Returns the cell area of each cell.

Return: dimension: numCells

const std::vector<std::vector<unsigned>> &GetNeighbours() const¶

Returns the neighbor cell IDs for every cell.

Return: dimension: numCells x numNodes

const std::vector<std::vector<Vector>> &GetNormals() const¶

Returns the edge length scaled normal vectors of each cell.

Return: dimension: numCells x numNodes x dim

const std::vector<BOUNDARY_TYPE> &GetBoundaryTypes() const¶

Returns the boundary enum for each cell. BOUNDARY_TYPE::NONE is the default.

Return: dimension: numCells

void SetBoundaryType(int idx_cell, BOUNDARY_TYPE boundary_type)¶

Set a boundary type for a cell. Use with caution!

Return

void

Parameters

idx_cell: cell index
boundary_type: boundary_type to change the cell to

const std::vector<std::pair<double, double>> GetBounds() const¶

Returns the minimal and maximal coordinates of all nodes for each dimension.

Return: dimension: dim

const std::vector<std::vector<Vector>> GetInterfaceMidPoints() const¶

Returns the the coordinates of midpoints of all interfaces of all cells.

Return: dimension: numCells<numNeighborsPerCell<nDim>>

double GetDistanceToOrigin(unsigned idx_cell) const¶

Returns distance of a specified cells center to the coordinate systems origin.

Return: dimension: scalar

void ComputeSlopes(unsigned nq, VectorVector &psiDerX, VectorVector &psiDerY, const VectorVector &psi) const¶

ComputeSlopes calculates the slope in every cell into x and y direction using the divergence theorem.

Parameters

nq: is number of quadrature points
psiDerX: is slope in x direction (gets computed. Slope is stored here)
psiDerY: is slope in y direction (gets computed. Slope is stored here)
psi: is solution for which slope is computed

void ComputeLimiter(unsigned nSys, const VectorVector &solDx, const VectorVector &solDy, const VectorVector &sol, VectorVector &limiter) const¶

Use gauss theorem and limiters. For unstructured mesh *.

Parameters

nq: is number of quadrature points
psiDerX: is slope in x direction (gets computed. Slope is stored here)
psiDerY: is slope in y direction (gets computed. Slope is stored here)
psi: is solution for which slope is computed

void ComputeLimiter1D(unsigned nSys, const VectorVector &sol, VectorVector &limiter) const¶

Use gauss theorem and limiters. For unstructured mesh *.

Parameters

nq: is number of quadrature points
psiDerX: is slope in x direction (gets computed. Slope is stored here)
psiDerY: is slope in y direction (gets computed. Slope is stored here)
psi: is solution for which slope is computed

void ComputeSlopes1D(unsigned nq, VectorVector &psiDerX, const VectorVector &psi) const¶

ComputeSlopes calculates for 1D meshes using finite difference formula in x direction.

Parameters

nq: is number of quadrature points
psiDerX: is slope in x direction (gets computed. Slope is stored here)
psi: is solution for which slope is computed

Protected Functions

void ComputeCellAreas()¶: Computes only the areas of the mesh cells. Write to _cellAreas.

void ComputeCellMidpoints()¶: Compute only the midpoints of the cells. Write to _cellMidPoints.

void ComputeConnectivity()¶: Computes _cellNeighbors and _nodeNeighbors, i.e. neighborship relation in mesh.

Vector ComputeOutwardFacingNormal(const Vector &nodeA, const Vector &nodeB, const Vector &cellCenter)¶

Computes outward facing normal of two neighboring nodes nodeA and nodeB with common cellCellcenter. Normals are scaled with their respective edge length.

Return

outward facing normal

Parameters

nodeA: first node
nodeB: neighboring node to nodeA
cellCenter: Center of the cell that has nodeA and nodeB as nodes.

void ComputeBounds()¶: Computes the spatial bounds of a 2D domain.

Vector ComputeCellInterfaceMidpoints(const Vector &nodeA, const Vector &nodeB)¶: compute the midpoint of the edge between nodeA and nodeB

Protected Attributes

const unsigned _dim¶: spatial dimension of the mesh, i.e. 1D,2D,3D

const unsigned _numCells¶: number of cells in the mesh

const unsigned _numNodes¶: number of nodes in the mesh (for node centered view)

const unsigned _numNodesPerCell¶: number of nodes per cell

const unsigned _numBoundaries¶: number of boundary cells in the mesh

const unsigned _ghostCellID¶: Id of the ghost cell. (we use only one ghost cell). equal to _numCells and therefore has the ID of the last cell + 1.

std::vector<Vector> _nodes¶: nodes coordinates in the mesh. dimension:_numNodes<_dim>

std::vector<std::vector<unsigned>> _cells¶: node indices for each cell. dimension:_numCells<_numNodesPerCell>

std::vector<std::pair<BOUNDARY_TYPE, std::vector<unsigned>>> _boundaries¶: boundary cells in the mesh. Pair defines boundary type of the boundary nodes of the cell. numBoundaries<(1,numBoundaryNodes)>

std::vector<double> _cellAreas¶: cell areas of the mesh. dimension: numCells

std::vector<Vector> _cellMidPoints¶: cell midpoints of the mesh. dimension: numCells<dim>

std::vector<std::vector<unsigned>> _cellNeighbors¶: neighbors of each cell. dimension: numCells<numNodesPerCell>

std::vector<std::vector<Vector>> _cellInterfaceMidPoints¶: coordinates of the interface midpoints of all cells of the mesh dimension: numCells<numNeighborsPerCell<nDim>>. interfaces of each cell are in same order as _cellNeighbors

std::vector<std::vector<Vector>> _cellNormals¶: outward facing normals of each side of each cell. dimension: numCells<numNodesPerCell<dim>>, all normals are facing away from the cell center, and scaled with the edge length

std::vector<BOUNDARY_TYPE> _cellBoundaryTypes¶: Tags each cell with its boundary type. None means no boundary. dimension: numCells.

blaze::CompressedMatrix<bool> _nodeNeighbors¶: neighborshood relationship of nodes for (par-)metis