functor to get pointwise momentum exchange on local lattice (block level) More...

#include <latticeMomentumExchangeForce.h>

Inheritance diagram for olb::BlockLatticeMomentumExchangeForceLocal< T, DESCRIPTOR, PARTICLETYPE, useTorque >:

Collaboration diagram for olb::BlockLatticeMomentumExchangeForceLocal< T, DESCRIPTOR, PARTICLETYPE, useTorque >:

Public Member Functions
	BlockLatticeMomentumExchangeForceLocal (BlockLattice< T, DESCRIPTOR > &blockLattice, const BlockGeometry< T, DESCRIPTOR::d > &blockGeometry, particles::ParticleSystem< T, PARTICLETYPE > &particleSystem, const UnitConverter< T, DESCRIPTOR > &converter)

bool	operator() (T output[], const int input[]) override

Detailed Description

template<typename T, typename DESCRIPTOR, typename PARTICLETYPE, bool useTorque = false>
class olb::BlockLatticeMomentumExchangeForceLocal< T, DESCRIPTOR, PARTICLETYPE, useTorque >

functor to get pointwise momentum exchange on local lattice (block level)

Definition at line 147 of file latticeMomentumExchangeForce.h.

Constructor & Destructor Documentation

◆ BlockLatticeMomentumExchangeForceLocal()

template<typename T , typename DESCRIPTOR , typename PARTICLETYPE , bool useTorque>

olb::BlockLatticeMomentumExchangeForceLocal< T, DESCRIPTOR, PARTICLETYPE, useTorque >::BlockLatticeMomentumExchangeForceLocal	(	BlockLattice< T, DESCRIPTOR > &	blockLattice,
		const BlockGeometry< T, DESCRIPTOR::d > &	blockGeometry,
		particles::ParticleSystem< T, PARTICLETYPE > &	particleSystem,
		const UnitConverter< T, DESCRIPTOR > &	converter )

Definition at line 259 of file latticeMomentumExchangeForce.hh.

  : BlockLatticePhysF<T,DESCRIPTOR>(blockLattice,converter,
      useTorque ? utilities::dimensions::convert<DESCRIPTOR::d>::rotation : DESCRIPTOR::d),
    _blockLattice(blockLattice),
    _blockGeometry(blockGeometry),
    _particleSystem(particleSystem)
{
  this->getName() = "localMomentumExchange";
}

References olb::getName().

Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename DESCRIPTOR , typename PARTICLETYPE , bool useTorque>

bool olb::BlockLatticeMomentumExchangeForceLocal< T, DESCRIPTOR, PARTICLETYPE, useTorque >::operator()	(	T	output[],
		const int	input[] )

override

Definition at line 274 of file latticeMomentumExchangeForce.hh.

{
  using namespace descriptors;
  const unsigned D = DESCRIPTOR::d;
  const unsigned Drot = utilities::dimensions::convert<D>::rotation;
 
  for (unsigned iDim=0; iDim<D; ++iDim) {
    output[iDim] = 0.;
  }
 
  // iterate over all particles in _particleSystem
  int _iP0 = 0;
  for (std::size_t iP=_iP0; iP!=_particleSystem.size(); ++iP) {
    auto particle = _particleSystem.get(iP);
    auto circumRadius = particle.template getField<SURFACE,SINDICATOR>()->getCircumRadius();
    auto position = particle.template getField<GENERAL,POSITION>();
 
    LatticeR<DESCRIPTOR::d> start, end;
    if ( particles::getBlockParticleIntersection(this->_blockGeometry, T(1.)/this->_converter.getPhysDeltaX(), start, end,
                                      position, circumRadius) ){
 
      Vector<T,D>tmpForce(0.);
      PhysR<T,D> lever(0.);
      if ( particles::resolved::momentumExchangeAtSurfaceLocation(tmpForce.data(),
            lever, input, this->_blockGeometry,
            this->_blockLattice, this->_converter, particle) ){
 
        //Shift centre of mass (for oblique rotation) if offset provided
        if constexpr ( PARTICLETYPE::template providesNested<SURFACE,COR_OFFSET>() ){
          lever -= Vector<T,D> ( particle.template getField<SURFACE,COR_OFFSET>() );
        }
 
        //Calculate torque
        const auto torque = particles::dynamics::torque_from_force<D,T>::calculate( tmpForce, lever );
 
        //Add force and torque to output
        if constexpr(!useTorque){
          for (unsigned iDim=0; iDim<D; ++iDim) {
            output[iDim] = tmpForce[iDim];
          }
        } else {
          for (unsigned iRot=0; iRot<Drot; ++iRot) {
            output[iRot] = torque[iRot];
          }
        }
      }
    }
  }
  return true;
}

References olb::particles::dynamics::torque_from_force< D, T >::calculate(), olb::Vector< T, Size >::data(), olb::particles::getBlockParticleIntersection(), and olb::particles::resolved::momentumExchangeAtSurfaceLocation().