LineSegmentLevelSetAux Class Reference

Calculate level set values for an interface that is defined by a set of line segments. More...

#include <LineSegmentLevelSetAux.h>

Inheritance diagram for LineSegmentLevelSetAux:

Public Member Functions
	LineSegmentLevelSetAux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue () override
virtual void	compute () override
Real	calculateSignedDistance (Point p)
	calculate the signed distance value for a given point. More...

Protected Attributes
const LineSegmentCutSetUserObject *	_linesegment_uo
	Pointer to the LineSegmentCutSetUserObject object. More...
std::vector< Real >	_cut_data
	Store the cut locations. More...

Detailed Description

Calculate level set values for an interface that is defined by a set of line segments.

Definition at line 24 of file LineSegmentLevelSetAux.h.

Constructor & Destructor Documentation

LineSegmentLevelSetAux()

LineSegmentLevelSetAux::LineSegmentLevelSetAux

(

const InputParameters &

parameters

)

Definition at line 29 of file LineSegmentLevelSetAux.C.

  : AuxKernel(parameters)
{
  if (!isNodal())
    mooseError("LineSegmentLevelSetAux: Aux variable must be nodal variable.");
 
  FEProblemBase * fe_problem = dynamic_cast<FEProblemBase *>(&_subproblem);
 
  const UserObject * uo = &(
      fe_problem->getUserObjectBase(getParam<UserObjectName>("line_segment_cut_set_user_object")));
 
  if (dynamic_cast<const LineSegmentCutSetUserObject *>(uo) == nullptr)
    mooseError(
        "Failed to cast UserObject to LineSegmentCutSetUserObject in LineSegmentLevelSetAux");
 
  _linesegment_uo = dynamic_cast<const LineSegmentCutSetUserObject *>(uo);
}

Member Function Documentation

calculateSignedDistance()

Real LineSegmentLevelSetAux::calculateSignedDistance ( Point  p )

protected

calculate the signed distance value for a given point.

Parameters

p	Coordinate of point

Returns

Signed distance

Definition at line 48 of file LineSegmentLevelSetAux.C.

{
  const int line_cut_data_len = 6;
 
  Real min_dist = std::numeric_limits<Real>::max();
 
  for (unsigned int i = 0; i < _cut_data.size() / line_cut_data_len; ++i)
  {
    Point a = Point(_cut_data[i * line_cut_data_len + 0], _cut_data[i * line_cut_data_len + 1], 0);
    Point b = Point(_cut_data[i * line_cut_data_len + 2], _cut_data[i * line_cut_data_len + 3], 0);
 
    Point c = p - a;
    Point v = (b - a) / (b - a).norm();
    Real d = (b - a).norm();
    Real t = v * c;
 
    Real dist;
    Point nearest_point;
 
    if (t < 0)
    {
      dist = (p - a).norm();
      nearest_point = a;
    }
    else if (t > d)
    {
      dist = (p - b).norm();
      nearest_point = b;
    }
    else
    {
      v *= t;
      dist = (p - a - v).norm();
      nearest_point = (a + v);
    }
 
    Point p_nearest_point = nearest_point - p;
 
    Point normal_ab = Point(-(b - a)(1), (b - a)(0), 0);
 
    if (normal_ab * p_nearest_point < 0)
      dist = -dist;
 
    if (std::abs(dist) < std::abs(min_dist))
      min_dist = dist;
  }
 
  return min_dist;
}

Referenced by computeValue().

compute()

void LineSegmentLevelSetAux::compute ( )

overrideprotectedvirtual

Definition at line 99 of file LineSegmentLevelSetAux.C.

{
  _cut_data = _linesegment_uo->getCutData();
  AuxKernel::compute();
}

computeValue()

Real LineSegmentLevelSetAux::computeValue ( )

overrideprotectedvirtual

Definition at line 106 of file LineSegmentLevelSetAux.C.

{
  return calculateSignedDistance(*_current_node);
}

Member Data Documentation

_cut_data

std::vector<Real> LineSegmentLevelSetAux::_cut_data

protected

Store the cut locations.

Definition at line 44 of file LineSegmentLevelSetAux.h.

Referenced by calculateSignedDistance(), and compute().

_linesegment_uo

const LineSegmentCutSetUserObject* LineSegmentLevelSetAux::_linesegment_uo

protected

Pointer to the LineSegmentCutSetUserObject object.

Definition at line 41 of file LineSegmentLevelSetAux.h.

Referenced by compute(), and LineSegmentLevelSetAux().

---

The documentation for this class was generated from the following files:

• xfem/include/auxkernels/LineSegmentLevelSetAux.h
• xfem/src/auxkernels/LineSegmentLevelSetAux.C