MooseMeshXYCuttingUtils.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "libmesh/mesh_base.h"
#include "libmesh/boundary_info.h"
#include "libmesh/replicated_mesh.h"

#include "MooseUtils.h"
#include "MooseTypes.h"
#include "FaceInfo.h"

namespace MooseMeshXYCuttingUtils
{
void lineRemoverMoveNode(libMesh::ReplicatedMesh & mesh,
                         const std::vector<Real> & bdry_pars,
                         const subdomain_id_type block_id_to_remove,
                         const std::set<subdomain_id_type> & subdomain_ids_set,
                         const boundary_id_type trimming_section_boundary_id,
                         const boundary_id_type external_boundary_id,
                         const std::vector<boundary_id_type> & other_boundaries_to_conform =
                             std::vector<boundary_id_type>(),
                         const bool assign_ext_to_new = false,
                         const bool side_to_remove = true);

bool lineSideDeterminator(const Real px,
                          const Real py,
                          const Real param_1,
                          const Real param_2,
                          const Real param_3,
                          const bool direction_param,
                          const Real dis_tol = libMesh::TOLERANCE);

Point twoLineIntersection(const Real param_11,
                          const Real param_12,
                          const Real param_13,
                          const Real param_21,
                          const Real param_22,
                          const Real param_23);

Point twoPointandLineIntersection(const Point & pt1,
                                  const Point & pt2,
                                  const Real param_1,
                                  const Real param_2,
                                  const Real param_3);

bool
quasiTriElementsFixer(libMesh::ReplicatedMesh & mesh,
                      const std::set<subdomain_id_type> & subdomain_ids_set,
                      const subdomain_id_type tri_elem_subdomain_shift = Moose::INVALID_BLOCK_ID,
                      const SubdomainName tri_elem_subdomain_name_suffix = "tri");

std::vector<std::pair<Real, unsigned int>> vertex_angles(const Elem & elem);

std::vector<std::pair<Real, unsigned int>> vertex_distances(const Elem & elem);

void triElemSplitter(libMesh::ReplicatedMesh & mesh,
                     const dof_id_type elem_id,
                     const unsigned short node_shift,
                     const dof_id_type nid_3,
                     const dof_id_type nid_4,
                     const subdomain_id_type single_elem_side_id,
                     const subdomain_id_type double_elem_side_id);

void triElemSplitter(libMesh::ReplicatedMesh & mesh,
                     const dof_id_type elem_id,
                     const unsigned short node_shift,
                     const dof_id_type nid_m,
                     const subdomain_id_type first_elem_side_id,
                     const subdomain_id_type second_elem_side_id);

void quadElemSplitter(libMesh::ReplicatedMesh & mesh,
                      const dof_id_type elem_id,
                      const subdomain_id_type tri_elem_subdomain_shift);

void quadToTriOnLine(libMesh::ReplicatedMesh & mesh,
                     const std::vector<Real> & cut_line_params,
                     const dof_id_type tri_subdomain_id_shift,
                     const SubdomainName tri_elem_subdomain_name_suffix);

void lineRemoverCutElemTri(libMesh::ReplicatedMesh & mesh,
                           const std::vector<Real> & cut_line_params,
                           const subdomain_id_type block_id_to_remove,
                           const boundary_id_type new_boundary_id);

void lineRemoverCutElem(libMesh::ReplicatedMesh & mesh,
                        const std::vector<Real> & cut_line_params,
                        const dof_id_type tri_subdomain_id_shift,
                        const SubdomainName tri_elem_subdomain_name_suffix,
                        const subdomain_id_type block_id_to_remove,
                        const boundary_id_type new_boundary_id,
                        const bool improve_boundary_tri_elems = false);

void boundaryTriElemImprover(libMesh::ReplicatedMesh & mesh,
                             const boundary_id_type boundary_to_improve);

void makeImprovedTriElement(
    libMesh::ReplicatedMesh & mesh,
    const dof_id_type node_id_0,
    const dof_id_type node_id_1,
    const dof_id_type node_id_2,
    const subdomain_id_type subdomain_id,
    const std::vector<dof_id_type> & extra_elem_ids,
    const std::vector<boundary_id_type> & boundary_ids_for_side_1 = std::vector<boundary_id_type>(),
    const std::vector<boundary_id_type> & boundary_ids_for_side_0 = std::vector<boundary_id_type>(),
    const std::vector<boundary_id_type> & boundary_ids_for_side_2 =
        std::vector<boundary_id_type>());

bool elemSideLocator(libMesh::ReplicatedMesh & mesh,
                     const dof_id_type elem_id,
                     const dof_id_type node_id_0,
                     const dof_id_type node_id_1,
                     unsigned short & side_id,
                     bool & is_inverse);
}