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

#include <SBMBndEdge2.h>

Inheritance diagram for SBMBndEdge2:

[legend]

Public Member Functions
	SBMBndEdge2 (const Elem *elem)
	Constructor.

Ball	computeBoundingBall () const override

bool	intersect (const libMesh::Plane &pl, Point &intersect_p) const

bool	intersect (const LineSegment &l1, Point &intersect_p) const

bool	intersect (const LineSegment &line_segment) const override

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

Point	closest_point (const Point &p) const

bool	closest_normal_point (const Point &p, Point &closest_p) const

bool	contains_point (const Point &p) const

const Point &	start () const

const Point &	end () const

void	setStart (const Point &p0)

void	setEnd (const Point &p1)

void	set (const Point &p0, const Point &p1)

Real	length () const

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

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.

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

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.

Private Member Functions
bool	closest_point (const Point &p, bool clamp_to_segment, Point &closest_p) const

Private Attributes
Point	_p0

Point	_p1

const Elem *	_elem
	< Pointer to the libMesh element representing this boundary face.

const Point	_normal

Detailed Description

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

Definition at line 19 of file SBMBndEdge2.h.

Constructor & Destructor Documentation

◆ SBMBndEdge2()

SBMBndEdge2::SBMBndEdge2 ( const Elem * elem )

explicit

Constructor.

Definition at line 13 of file SBMBndEdge2.C.

  : LineSegment(elem->point(0), elem->point(1)), SBMBndElementBase(elem, LineSegment::normal())
{
  mooseAssert(elem->type() == EDGE2, "Element must be of type EDGE2");
  mooseAssert(MooseUtils::absoluteFuzzyEqual(elem->point(0)(2), elem->point(1)(2)),
              "Currently SBMBndEdge2 must be parallel to the x-y plane, i.e., z-coordinates of "
              "points must be equal");
}

Member Function Documentation

◆ computeBoundingBall()

Ball LineSegment::computeBoundingBall ( ) const

virtual

Reimplemented from LineSegment.

◆ 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.

{
  // 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;
}

Referenced by UnsignedDistanceToSurfaceMesh::distanceVectorToSurface(), TEST(), TEST(), TEST(), TEST(), and TEST().

◆ elem()

const Elem & SBMBndElementBase::elem ( ) const

inlineinherited

Getter for the underlying element.

Definition at line 32 of file SBMBndElementBase.h.

32{ return *_elem; }

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

◆ 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; }

◆ 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();
}

Referenced by TEST().

◆ intersect() [1/3]

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

Referenced by TEST(), and TEST().

◆ intersect() [2/3]

bool LineSegment::intersect	(	const LineSegment &	l1,
		Point &	intersect_p
	)		const

◆ intersect() [3/3]

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

virtual

Reimplemented from LineSegment.

◆ normal()

const Point & SBMBndElementBase::normal ( ) const

inline

Getter for the normal vector.

Definition at line 35 of file SBMBndElementBase.h.

35{ return _normal; }

Referenced by TEST(), and TEST().

Member Data Documentation

◆ _elem

const Elem* SBMBndElementBase::_elem

privateinherited

< Pointer to the libMesh element representing this boundary face.

Derived classes access it through elem() rather than this field directly. Unit normal vector of the boundary element, computed eagerly in the derived constructor and passed through the base constructor.

Definition at line 80 of file SBMBndElementBase.h.

Referenced by SBMBndElementBase::distanceFrom(), SBMBndElementBase::elem(), SBMBndElementBase::expectedEmbeddingMeshDim(), and SBMBndElementBase::getProjectedBoundingBoxDiagonal().

◆ _normal

const Point SBMBndElementBase::_normal