Builds mapping between two aligned subdomains/boundaries. More...

#include <MeshAlignmentBase.h>

Inheritance diagram for MeshAlignmentBase:

Public Member Functions
	MeshAlignmentBase (const MooseMesh &mesh)
	Constructor. More...

const std::vector< dof_id_type > &	getPrimaryElemIDs () const
	Returns the list of element IDs on the primary boundary. More...

const std::vector< dof_id_type > &	getSecondaryElemIDs () const
	Returns the list of element IDs on the secondary boundary. More...

bool	meshesAreAligned () const
	Returns true if the primary and secondary meshes are aligned. More...

Protected Member Functions
void	extractFrom1DElements (const std::vector< dof_id_type > &elem_ids, std::vector< Point > &elem_points, std::vector< dof_id_type > &node_ids, std::vector< Point > &node_points) const
	Extracts mesh information from 1D elements. More...

void	extractFromBoundaryInfo (const std::vector< std::tuple< dof_id_type, unsigned short int >> &boundary_info, std::vector< dof_id_type > &elem_ids, std::vector< unsigned short int > &side_ids, std::vector< Point > &side_points, std::vector< dof_id_type > &node_ids, std::vector< Point > &node_points) const
	Extracts mesh information from boundary info. More...

Protected Attributes
const MooseMesh &	_mesh
	Mesh. More...

std::vector< dof_id_type >	_primary_elem_ids
	List of primary element IDs. More...

std::vector< dof_id_type >	_secondary_elem_ids
	List of secondary element IDs. More...

std::vector< Point >	_primary_elem_points
	List of primary element points. More...

std::vector< Point >	_secondary_elem_points
	List of secondary element points. More...

std::vector< unsigned short int >	_primary_side_ids
	List of primary side IDs (if any) More...

std::vector< unsigned short int >	_secondary_side_ids
	List of secondary side IDs (if any) More...

std::vector< dof_id_type >	_primary_node_ids
	List of primary node IDs. More...

std::vector< dof_id_type >	_secondary_node_ids
	List of secondary node IDs. More...

std::vector< Point >	_primary_node_points
	List of primary node points. More...

std::vector< Point >	_secondary_node_points
	List of secondary node points. More...

bool	_meshes_are_aligned
	Flag that meshes are aligned. More...

Detailed Description

Builds mapping between two aligned subdomains/boundaries.

This class handles the following cases:

1D subdomain coupled to 2D boundary
2D boundary coupled to 2D boundary

Definition at line 21 of file MeshAlignmentBase.h.

Constructor & Destructor Documentation

◆ MeshAlignmentBase()

MeshAlignmentBase::MeshAlignmentBase ( const MooseMesh & mesh )

Constructor.

Parameters

mesh[in] mesh Mesh

Definition at line 14 of file MeshAlignmentBase.C.

   : _mesh(mesh), _meshes_are_aligned(false)
 {
 }

Member Function Documentation

◆ extractFrom1DElements()

void MeshAlignmentBase::extractFrom1DElements	(	const std::vector< dof_id_type > &	elem_ids,
		std::vector< Point > &	elem_points,
		std::vector< dof_id_type > &	node_ids,
		std::vector< Point > &	node_points
	)		const

protected

Extracts mesh information from 1D elements.

Parameters

[in]	elem_ids	Vector of element IDs
[out]	elem_points	Vector of element centroids
[out]	node_ids	Vector of node IDs
[out]	node_points	Vector of node points

Definition at line 20 of file MeshAlignmentBase.C.

Referenced by MeshAlignment1D3D::initialize(), and MeshAlignment::initialize().

 {
   elem_points.clear();
   node_ids.clear();
   node_points.clear();
 
   for (const auto & elem_id : elem_ids)
   {
     const Elem * elem = _mesh.elemPtr(elem_id);
     elem_points.push_back(elem->vertex_average());
 
     for (const auto j : elem->node_index_range())
     {
       const Node & node = elem->node_ref(j);
       const auto node_id = node.id();
       if (std::find(node_ids.begin(), node_ids.end(), node_id) == node_ids.end())
       {
         node_ids.push_back(node_id);
         node_points.push_back(node);
       }
     }
   }
 }

◆ extractFromBoundaryInfo()

void MeshAlignmentBase::extractFromBoundaryInfo	(	const std::vector< std::tuple< dof_id_type, unsigned short int >> &	boundary_info,
		std::vector< dof_id_type > &	elem_ids,
		std::vector< unsigned short int > &	side_ids,
		std::vector< Point > &	side_points,
		std::vector< dof_id_type > &	node_ids,
		std::vector< Point > &	node_points
	)		const

protected

Extracts mesh information from boundary info.

Parameters

[in]	boundary_info	Vector of tuples of element ID and side ID on boundary
[out]	elem_ids	Vector of element IDs
[out]	side_ids	Vector of side IDs
[out]	side_points	Vector of side centroids
[out]	node_ids	Vector of node IDs
[out]	node_points	Vector of node points

Definition at line 48 of file MeshAlignmentBase.C.

Referenced by MeshAlignment2D2D::initialize(), MeshAlignment1D3D::initialize(), MeshAlignment2D3D::initialize(), and MeshAlignment::initialize().

 {
   elem_ids.clear();
   side_ids.clear();
   side_points.clear();
   node_ids.clear();
   node_points.clear();
 
   for (const auto & elem_id_and_side : boundary_info)
   {
     auto elem_id = std::get<0>(elem_id_and_side);
     elem_ids.push_back(elem_id);
 
     auto side = std::get<1>(elem_id_and_side);
     side_ids.push_back(side);
 
     const Elem * elem = _mesh.elemPtr(elem_id);
     const Elem * side_elem = elem->build_side_ptr(side).release();
     const Point side_center = side_elem->vertex_average();
     side_points.push_back(side_center);
 
     for (const auto j : side_elem->node_index_range())
     {
       const Node & node = side_elem->node_ref(j);
       const auto node_id = node.id();
       if (std::find(node_ids.begin(), node_ids.end(), node_id) == node_ids.end())
       {
         node_ids.push_back(node_id);
         node_points.push_back(node);
       }
     }
     delete side_elem;
   }
 }