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

#include <smoothIndicatorF2D.h>

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

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

Public Member Functions
	SmoothIndicatorCustom2D (T latticeSpacing, std::shared_ptr< IndicatorF2D< T > > indPtr, Vector< T, 2 > pos, T epsilon, T theta=0.)
Vector< T, 2 >	getLocalCenter ()
S	getArea () override
Vector< T, 2 >	calcMofiAndMass (T rhoP) override
	calculates and returns mofi and mass of the particle (mofi is at index 0 and mass at index 1)
Vector< S, 2 >	surfaceNormal (const Vector< S, 2 > &pos, const S meshSize) override
const S	signedDistance (const Vector< S, 2 > input) override
bool	regardCell (int input[2])
bool	operator() (T output[], const S input[]) override

Detailed Description

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

Definition at line 123 of file smoothIndicatorF2D.h.

Constructor & Destructor Documentation

SmoothIndicatorCustom2D()

template<typename T , typename S , bool PARTICLE>

olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::SmoothIndicatorCustom2D	(	T	latticeSpacing,
		std::shared_ptr< IndicatorF2D< T > >	indPtr,
		Vector< T, 2 >	pos,
		T	epsilon,
		T	theta = 0. )

Definition at line 283 of file smoothIndicatorF2D.hh.

  :_indPtr(indPtr),
   _latticeSpacing(latticeSpacing)
{
  OstreamManager clout(std::cout,"createIndicatorCustom2D");
  this->_name = "custom2D";
  this->_epsilon = epsilon;
  if constexpr(!PARTICLE) {
    this->_pos = pos;         // global position of the local center
    this->_theta = util::degreeToRadian(theta);
    this->init();
  }
 
  initData(*_indPtr);
}

References olb::util::degreeToRadian().

Here is the call graph for this function:

Member Function Documentation

calcMofiAndMass()

template<typename T , typename S , bool PARTICLE>

Vector< T, 2 > olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::calcMofiAndMass ( T rhoP )

override

calculates and returns mofi and mass of the particle (mofi is at index 0 and mass at index 1)

Definition at line 377 of file smoothIndicatorF2D.hh.

{
  // TODO - calculation
  T mofi = 0.;
  T mass;
  unsigned nCells = 0;
  int input[2];
  for (input[0] = 0; input[0] < this->_blockData->getNx(); ++input[0]) {
    const T dx = util::abs(input[0]*_latticeSpacing - this->_center[0]);
    for (input[1] = 0; input[1] < this->_blockData->getNy(); ++input[1]) {
      if (regardCell(input)) {
        const T dy = util::abs(input[1]*_latticeSpacing - this->_center[1]);
        mofi += (dx*dx+dy*dy);
        ++nCells;
      }
    }
  }
  _area = nCells * util::pow(_latticeSpacing, 2);
  mass = rhoP * _area;
  mofi += util::pow(_latticeSpacing, 4)/ 6.0;
  mofi *= mass/nCells;
 
  return Vector<T,2>(mofi,mass);
}

References olb::util::abs(), and olb::util::pow().

Here is the call graph for this function:

getArea()

template<typename T , typename S , bool PARTICLE>

S olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::getArea ( )

override

Definition at line 370 of file smoothIndicatorF2D.hh.

{
  return _area;
}

getLocalCenter()

template<typename T , typename S , bool PARTICLE>

Vector< T, 2 > olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::getLocalCenter

(

)

Definition at line 438 of file smoothIndicatorF2D.hh.

{
  return this->_center;
}

operator()()

template<typename T , typename S , bool PARTICLE>

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

override

Definition at line 485 of file smoothIndicatorF2D.hh.

{
  PhysR<T,2> pos(input[0], input[1]);
  if constexpr(!PARTICLE) {
    pos = util::executeRotation<S,2,true>(pos, this->_rotMat, this->getPos());
  }
  if(norm(pos) < this->_circumRadius) {
    return SmoothIndicatorF2D<T, S, PARTICLE>::operator()(output, input);
  }
  output[0] = T{0};
  return false;
}

References olb::norm().

Here is the call graph for this function:

regardCell()

template<typename T , typename S , bool PARTICLE>

bool olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::regardCell

(

int

input[2]

)

Definition at line 479 of file smoothIndicatorF2D.hh.

{
  return this->_blockData->get(input) < std::numeric_limits<T>::epsilon();
}

signedDistance()

template<typename T , typename S , bool PARTICLE>

const S olb::SmoothIndicatorCustom2D< T, S, PARTICLE >::signedDistance ( const Vector< S, 2 > input )

override

Definition at line 450 of file smoothIndicatorF2D.hh.

{
  PhysR<T,2> position = input;
  if constexpr (!PARTICLE) {
    // translation, counter-clockwise rotation by _theta=-theta around (0/0) and movement from rotation center to local center
    position = util::executeRotation<T,2,true>(input, this->_rotMat, this->getPos());
  }
  // The block data originates in (0,0,0) therefore we translate the input position which is relative to center of mass
  const PhysR<T,2> positionInCache = this->_center + position;
 
  T signedDistance(0.);
  if(_interpolateCache->operator()(&signedDistance, positionInCache.data())) {
    return signedDistance;
  }
 
  // If all points were outside return an estimation instead
  LatticeR<2> latticePosition;
  PhysR<T,2> extraDistance;
  for(unsigned iDim=0; iDim<2; ++iDim) {
    latticePosition[iDim] = util::round( positionInCache[iDim] / this->_latticeSpacing );
    latticePosition[iDim] = util::max(0, latticePosition[iDim]);
    latticePosition[iDim] = util::min(this->_blockData->getExtent()[iDim] - 1, latticePosition[iDim]);
    // The extra distance is always positive because it must be outside the geometry
    extraDistance[iDim] = util::abs(_latticeSpacing * latticePosition[iDim] - positionInCache[iDim]);
  }
  return this->_blockData->get(latticePosition.data()) + norm(extraDistance);
}

References olb::util::abs(), olb::Vector< T, D >::data(), olb::util::max(), olb::util::min(), olb::norm(), and olb::util::round().

Here is the call graph for this function:

surfaceNormal()

template<typename T , typename S , bool PARTICLE>

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

override

Definition at line 444 of file smoothIndicatorF2D.hh.

{
  return _indPtr->surfaceNormal(pos, meshSize);
}

---

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

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