NekMeshGenerator.h

Go to the documentation of this file.

/********************************************************************/
/* SOFTWARE COPYRIGHT NOTIFICATION */
/* Cardinal */
/* */
/* (c) 2021 UChicago Argonne, LLC */
/* ALL RIGHTS RESERVED */
/* */
/* Prepared by UChicago Argonne, LLC */
/* Under Contract No. DE-AC02-06CH11357 */
/* With the U. S. Department of Energy */
/* */
/* Prepared by Battelle Energy Alliance, LLC */
/* Under Contract No. DE-AC07-05ID14517 */
/* With the U. S. Department of Energy */
/* */
/* See LICENSE for full restrictions */
/********************************************************************/
 
#pragma once
 
#include "MeshGenerator.h"
 
class NekMeshGenerator : public MeshGenerator
{
public:
 static InputParameters validParams();
 
 NekMeshGenerator(const InputParameters & parameters);
 
 std::unique_ptr<MeshBase> generate() override;
 
 
 void initializeElemData(std::unique_ptr<MeshBase> & mesh);
 
 void checkElementType(std::unique_ptr<MeshBase> & mesh);
 
 virtual std::pair<unsigned int, unsigned int> pairedNodesAboutMidPoint(const unsigned int & node_id) const;
 
 unsigned int getFaceNode(const unsigned int & primary_face) const;
 
 virtual bool isCornerNode(const unsigned int & node) const;
 
 void moveNodes(std::unique_ptr<MeshBase> & mesh, std::vector<Real> & polygon_layer_smoothing);
 
 const Elem * getNextLayerElem(const Elem & elem, const unsigned int & touching_face,
 unsigned int & next_touching_face) const;
 
 unsigned int midPointNodeIndex(const unsigned int & face_id, const unsigned int & face_node) const;
 
 unsigned int pairedFaceNode(const unsigned int & node_id, const unsigned int & face_id) const;
 
 std::vector<Elem *> getBoundaryLayerElems(Elem * elem, const unsigned int & n_layers,
 const unsigned int & primary_face) const;
 
 void moveElem(Elem * elem, const unsigned int & boundary_index, const unsigned int & primary_face,
 const std::vector<Real> & polygon_layer_smoothing);
 
 bool isNearCorner(const Point & pt) const;
 
 unsigned int getNodeIndex(const Elem * elem, const Point & pt) const;
 
 Point getClosestOrigin(const unsigned int & index, const Point & pt) const;
 
 void checkPointLength(const std::vector<std::vector<Real>> & points, std::string name) const;
 
 void adjustMidPointNode(const unsigned int & node_id, Elem * elem) const;
 
 BoundaryID getBoundaryID(const BoundaryName & name, const MeshBase & mesh) const;
 
 Point projectPoint(const Point & origin, const Point & pt) const;
 
 Point adjustPointToCircle(const unsigned int & node_id, Elem * elem, const Real & radius, const Point & origin) const;
 
protected:
 std::unique_ptr<MeshBase> & _input;
 
 const MooseEnum & _geometry_type;
 
 const MooseEnum & _axis;
 
 const bool & _curve_corners;
 
 const Real & _rotation_angle;
 
 const bool & _retain_original_elem_type;
 
 const bool _has_moving_boundary;
 
 unsigned int _n_start_nodes;
 
 unsigned int _n_end_nodes;
 
 unsigned int _n_start_nodes_per_side;
 
 unsigned int _n_sides;
 
 unsigned int _n_corner_nodes;
 
 std::vector<BoundaryID> _moving_boundary;
 
 std::vector<Real> _radius;
 
 std::vector<std::vector<Real>> _origin;
 
 std::vector<unsigned int> _layers;
 
 Real _corner_radius;
 
 std::vector<Point> _polygon_corners;
 
 Real _max_corner_distance;
 
 unsigned int _n_noncorner_boundaries;
 
 std::set<BoundaryID> _boundaries_to_rebuild;
 
 std::map<BoundaryID, BoundaryName> _boundary_id_to_name;
 
 std::vector<std::vector<unsigned int>> _side_ids;
 
 std::vector<std::vector<unsigned int>> _corner_nodes;
 
 std::vector<std::vector<std::pair<unsigned int, unsigned int>>> _across_pair;
 
 std::vector<unsigned int> _across_face;
 
 ElemType _etype;
 
 std::vector<std::vector<unsigned int>> _side_nodes_map;
 
 std::vector<std::vector<unsigned int>> _face_nodes_map;
};