olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE > Class Template Reference

implements a smooth particle ellipsoid in 3D with an _epsilon sector. More...

#include <smoothIndicatorF3D.h>

Inheritance diagram for olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >:

Collaboration diagram for olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >:

Public Member Functions
	SmoothIndicatorEllipsoid3D (IndicatorEllipsoid3D< S > &ind, S epsilon, Vector< S, 3 > theta=Vector< S, 3 >(0., 0., 0.))
	SmoothIndicatorEllipsoid3D (Vector< S, 3 > center, Vector< S, 3 > radius, S epsilon, Vector< S, 3 > theta=Vector< S, 3 >(0., 0., 0.))
IndicatorEllipsoid3D< S > &	getIndicator ()
Vector< S, 3 >	surfaceNormal (const Vector< S, 3 > &pos, const S meshSize) override
const S	signedDistance (const PhysR< S, 3 > input) override
S	getVolume () override
Vector< S, 4 >	calcMofiAndMass (const S density) override
bool	operator() (T output[], const S input[]) override

Detailed Description

template<typename T, typename S, bool PARTICLE = false>
class olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >

implements a smooth particle ellipsoid in 3D with an _epsilon sector.

Definition at line 59 of file smoothIndicatorF3D.h.

Constructor & Destructor Documentation

SmoothIndicatorEllipsoid3D() [1/2]

template<typename T , typename S , bool PARTICLE>

olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::SmoothIndicatorEllipsoid3D	(	IndicatorEllipsoid3D< S > &	ind,
		S	epsilon,
		Vector< S, 3 >	theta = Vector<S,3> (0.,0.,0.) )

Definition at line 134 of file smoothIndicatorF3D.hh.

  :SmoothIndicatorEllipsoid3D(ind.getCenter(), ind.getRadius(), epsilon, theta)
{ }

SmoothIndicatorEllipsoid3D() [2/2]

template<typename T , typename S , bool PARTICLE>

olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::SmoothIndicatorEllipsoid3D	(	Vector< S, 3 >	center,
		Vector< S, 3 >	radius,
		S	epsilon,
		Vector< S, 3 >	theta = Vector<S,3> (0.,0.,0.) )

Definition at line 140 of file smoothIndicatorF3D.hh.

  :_ind(center, radius)
{
  this->_epsilon = epsilon;
  if constexpr (!PARTICLE) {
    this->_pos = center;
    this->_theta = util::degreeToRadian(theta);
  }
 
  T const max_axis = util::max( radius[0], util::max(radius[1], radius[2]) );
  this->_circumRadius = max_axis+0.5*epsilon;
  if constexpr (!PARTICLE) {
    this->_myMin = {
      this->_pos[0] - this->getCircumRadius(),
      this->_pos[1] - this->getCircumRadius(),
      this->_pos[2] - this->getCircumRadius()
    };
    this->_myMax = {
      this->_pos[0] + this->getCircumRadius(),
      this->_pos[1] + this->getCircumRadius(),
      this->_pos[2] + this->getCircumRadius()
    };
    this->init();
  }
}

References olb::util::degreeToRadian(), and olb::util::max().

Here is the call graph for this function:

Member Function Documentation

calcMofiAndMass()

template<typename T , typename S , bool PARTICLE>

Vector< S, 4 > olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::calcMofiAndMass ( const S density )

override

Definition at line 180 of file smoothIndicatorF3D.hh.

{
  Vector<S,3> const radius = _ind.getRadius();
  T const mass = getVolume()*density;
  T const xHalfAxis2 = radius[0]*radius[0];
  T const yHalfAxis2 = radius[1]*radius[1];
  T const zHalfAxis2 = radius[2]*radius[2];
  Vector<S,3> mofi;
  mofi[0] = 0.2*mass*(yHalfAxis2+zHalfAxis2);
  mofi[1] = 0.2*mass*(xHalfAxis2+zHalfAxis2);
  mofi[2] = 0.2*mass*(yHalfAxis2+xHalfAxis2);
  return Vector<S,4>(mofi[0], mofi[1], mofi[2], mass);
}

getIndicator()

template<typename T , typename S , bool PARTICLE>

IndicatorEllipsoid3D< S > & olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::getIndicator

(

)

Definition at line 168 of file smoothIndicatorF3D.hh.

                                                                               {
  return _ind;
}

getVolume()

template<typename T , typename S , bool PARTICLE>

S olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::getVolume ( )

override

Definition at line 173 of file smoothIndicatorF3D.hh.

{
  Vector<S,3> const radius = _ind.getRadius();
  return 4./3.*M_PI*radius[0]*radius[1]*radius[2];
}

References M_PI.

operator()()

template<typename T , typename S , bool PARTICLE>

bool olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::operator() ( T output[], const S input[] )

override

Definition at line 217 of file smoothIndicatorF3D.hh.

{
  Vector<T,3> pos(0.);
  Vector<T,9> rotMatrix = {1, 0, 0, 0, 1, 0, 0, 0, 1};
  if constexpr (!PARTICLE) {
    pos = this->getPos();
    rotMatrix = this->getRotationMatrix();
  }
 
  //1.Calculate distance between point and center of unrotated indicator
  T xDist = input[0] - pos[0];
  T yDist = input[1] - pos[1];
  T zDist = input[2] - pos[2];
 
  //2.Calculate point projected to rotated indicator
  // counter-clockwise rotation by _theta=-theta around center
  T x = pos[0] + rotMatrix[0]*xDist + rotMatrix[3]*yDist + rotMatrix[6]*zDist;
  T y = pos[1] + rotMatrix[1]*xDist + rotMatrix[4]*yDist + rotMatrix[7]*zDist;
  T z = pos[2] + rotMatrix[2]*xDist + rotMatrix[5]*yDist + rotMatrix[8]*zDist;
 
  T a = (x - pos[0]) / (_ind.getRadius()[0] - 0.5*this->getEpsilon() );
  T b = (y - pos[1]) / (_ind.getRadius()[1] - 0.5*this->getEpsilon() );
  T c = (z - pos[2]) / (_ind.getRadius()[2] - 0.5*this->getEpsilon() );
  T aEps = (x - pos[0]) / (_ind.getRadius()[0] + 0.5*this->getEpsilon() );
  T bEps = (y - pos[1]) / (_ind.getRadius()[1] + 0.5*this->getEpsilon() );
  T cEps = (z - pos[2]) / (_ind.getRadius()[2] + 0.5*this->getEpsilon() );
 
  if ( (a*a+b*b+c*c) <= 1. ) {
    output[0] = 1.;
    return true;
  }
  if ( (aEps*aEps+bEps*bEps+cEps*cEps) <= 1. ) {
    // TODO: Here the correct distance to the ellipsoid has to be calculated for smooth transition
    //       For now the epsilon region is taken to be 0.5
    output[0] = .5;
    return true;
  }
  output[0] = 0.;
  return false;
}

signedDistance()

template<typename T , typename S , bool PARTICLE>

const S olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::signedDistance ( const PhysR< S, 3 > input )

override

Definition at line 201 of file smoothIndicatorF3D.hh.

{
  // counter-clockwise rotation by _theta=-theta around the center of mass
  Vector<S,3> p;
  if constexpr(!PARTICLE) {
    // counter-clockwise rotation by _theta=-theta around the current position of the center of mass & translation
    p = util::executeRotation<S,3,true>(input, this->_rotMat, this->getPos());
  }
  else {
    p = input;
  }
  return _ind.signedDistance(p + _ind.getCenter());
}

surfaceNormal()

template<typename T , typename S , bool PARTICLE>

Vector< S, 3 > olb::SmoothIndicatorEllipsoid3D< T, S, PARTICLE >::surfaceNormal ( const Vector< S, 3 > & pos, const S meshSize )

override

Definition at line 195 of file smoothIndicatorF3D.hh.

{
  return _ind.surfaceNormal(pos, meshSize);
}

---

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

• src/functors/analytical/indicator/smoothIndicatorF3D.h
• src/functors/analytical/indicator/smoothIndicatorF3D.hh