Derived class for 3-node triangular elements (Tri3) Triangle is listed first so it is initialized before SBMBndElementBase; the base class constructor then receives the unit normal from Triangle::normal() (see the .C file). More...

#include <SBMBndTri3.h>

Inheritance diagram for SBMBndTri3:

[legend]

Public Member Functions
	SBMBndTri3 (const Elem *elem)
	Constructor. More...

const Point &	normal () const
	Getter for the normal vector. More...

Ball	computeBoundingBall () const override

bool	intersect (const LineSegment &l, libMesh::Point &intersect_p) const

bool	intersect (const LineSegment &line_segment) const override

libMesh::Point	normal () const

bool	intersect (const LineSegment &l, libMesh::Point &intersect_p) const

bool	intersect (const LineSegment &line_segment) const override

Ball	computeBoundingBall () const override

const Elem &	elem () const
	Getter for the underlying element. More...

const Point &	normal () const
	Getter for the normal vector. More...

bool	intersect (const LineSegment &line_segment) const
	Check if the given line segment intersects this boundary element. More...

unsigned int	expectedEmbeddingMeshDim () const
	Getter of expected embedding solving mesh dimension Because the boundary element is a face of the embedding mesh, its dimension is one less than the dimension of the embedding mesh. More...

virtual Point	distanceFrom (const Point &pt) const
	Compute the distance vector from an arbitrary point to this boundary element. More...

Ball	computeBoundingBall () const
	Compute a bounding ball for this boundary element. More...

Real	getProjectedBoundingBoxDiagonal (const Point &normal_dir) const
	Compute the length of the element bounding-box diagonal projected onto a plane orthogonal to a given direction. More...

Detailed Description

Derived class for 3-node triangular elements (Tri3) Triangle is listed first so it is initialized before SBMBndElementBase; the base class constructor then receives the unit normal from Triangle::normal() (see the .C file).

Definition at line 19 of file SBMBndTri3.h.

Constructor & Destructor Documentation

◆ SBMBndTri3()

SBMBndTri3::SBMBndTri3 ( const Elem * elem )

explicit

Constructor.

Definition at line 13 of file SBMBndTri3.C.

   : Triangle(elem->point(0), elem->point(1), elem->point(2)),
     SBMBndElementBase(elem, Triangle::normal())
 {
   mooseAssert(elem->type() == TRI3, "Element must be of type TRI3");
 }

Member Function Documentation

◆ computeBoundingBall() [1/2]

Ball Triangle::computeBoundingBall

◆ computeBoundingBall() [2/2]

Ball SBMBndElementBase::computeBoundingBall ( ) const

inherited

Compute a bounding ball for this boundary element.

Delegates to the underlying geometry via dynamic_cast.

Definition at line 134 of file SBMBndElementBase.C.

 {
   if (const auto * geom = dynamic_cast<const GeometryBase *>(this))
     return geom->computeBoundingBall();
 
   mooseError("SBMBndElementBase::computeBoundingBall: unsupported geometry type");
 }

◆ distanceFrom()

Point SBMBndElementBase::distanceFrom ( const Point & pt ) const

virtualinherited

Compute the distance vector from an arbitrary point to this boundary element.

Default: normal-based distance vector. Override in derived class if needed. (a) First calculation of the normal-based distance and the projection point (b) If the projection point is outside the element, return the distance to the closest vertex.

Definition at line 39 of file SBMBndElementBase.C.

Referenced by UnsignedDistanceToSurfaceMesh::distanceVectorToSurface().

 {
   // Side-type precondition is validated once in the base ctor; no per-call
   // check needed here.
 
   // (a) Project pt onto the normal direction
   const auto vec_to_first = _elem->point(0) - pt;
   const auto scale = vec_to_first * _normal;
   const auto projection = _normal * scale;
 
   // Check if projection point lands inside the geometry
   if (_elem->contains_point(pt + projection))
     return projection;
 
   // (d) Point to closest edge or node
   Real min_dist = std::numeric_limits<Real>::max();
   Point closest_vec;
 
   const unsigned int n_edges = _elem->n_sides();
   for (unsigned int j = 0; j < n_edges; ++j)
   {
     std::unique_ptr<const Elem> curr_edge = _elem->build_side_ptr(j);
 
     switch (curr_edge->type())
     {
       case EDGE2:
       {
         const Point & p1 = *curr_edge->node_ptr(0);
         const Point & p2 = *curr_edge->node_ptr(1);
 
         const Point edge = p2 - p1;
         Real t = ((pt - p1) * edge) / (edge * edge);
         t = std::clamp(t, 0.0, 1.0);
         const Point proj = p1 + t * edge;
         const Real dist = (pt - proj).norm();
 
         if (dist < min_dist)
         {
           min_dist = dist;
           closest_vec = proj - pt;
         }
         break;
       }
 
       case NODEELEM:
       {
         const Point & p = *curr_edge->node_ptr(0);
         const Real dist = (pt - p).norm();
         if (dist < min_dist)
         {
           min_dist = dist;
           closest_vec = p - pt;
         }
         break;
       }
 
       default:
         mooseAssert(false, "unreachable: side type validated in SBMBndElementBase ctor");
     }
   }
 
   return closest_vec;
 }

◆ elem()

const Elem& SBMBndElementBase::elem ( ) const

inlineinherited

Getter for the underlying element.

Definition at line 32 of file SBMBndElementBase.h.

Referenced by SBMBndEdge2::SBMBndEdge2(), SBMBndElementBase::SBMBndElementBase(), and SBMBndTri3().

32 { return *_elem; }

SBMBndElementBase::_elem

const Elem * _elem

< Pointer to the libMesh element representing this boundary face.

Definition: SBMBndElementBase.h:80

◆ expectedEmbeddingMeshDim()

unsigned int SBMBndElementBase::expectedEmbeddingMeshDim ( ) const

inlineinherited

Getter of expected embedding solving mesh dimension Because the boundary element is a face of the embedding mesh, its dimension is one less than the dimension of the embedding mesh.

Definition at line 46 of file SBMBndElementBase.h.

46 { return _elem->dim() + 1; }

SBMBndElementBase::_elem

const Elem * _elem

< Pointer to the libMesh element representing this boundary face.

Definition: SBMBndElementBase.h:80

◆ getProjectedBoundingBoxDiagonal()

Real SBMBndElementBase::getProjectedBoundingBoxDiagonal ( const Point & normal_dir ) const

inherited

Compute the length of the element bounding-box diagonal projected onto a plane orthogonal to a given direction.

The projection plane is defined only by its normal direction; the specific location (offset) of the plane is irrelevant for this operation. Equivalently, this function removes the component of the bounding-box diagonal along the given normal direction and returns the norm of the remaining tangential component.

Parameters

normal_dir A direction vector defining the plane normal.

Returns: The length of the longest diagonal of the projected rectangle.

Definition at line 104 of file SBMBndElementBase.C.

 {
   BoundingBox bbox = _elem->loose_bounding_box();
 
   const Point & min_pt = bbox.first;
   const Point & max_pt = bbox.second;
 
   // Step (a): Calculate box_vec
   Point box_vec = max_pt - min_pt;
 
   // Step (b): Project box_vec onto normal_dir
   Real normal_scale = box_vec * normal_dir;
   Point normal_box_vec = normal_dir * normal_scale;
 
   // Step (c): Calculate tangent_vec and its norm
   Point tangent_vec = box_vec - normal_box_vec;
 
   return tangent_vec.norm();
 }

◆ intersect() [1/3]

bool Triangle::intersect

◆ intersect() [2/3]

bool Triangle::intersect

◆ intersect() [3/3]

bool SBMBndElementBase::intersect ( const LineSegment & line_segment ) const

inherited

Check if the given line segment intersects this boundary element.

Delegates to the underlying geometry via dynamic_cast.

Definition at line 125 of file SBMBndElementBase.C.

 {
   if (const auto * geom = dynamic_cast<const GeometryBase *>(this))
     return geom->intersect(line_segment);
 
   mooseError("SBMBndElementBase::intersect: unsupported geometry type");
 }