Verlet dynamics for particles aware of their DYNAMIC_STATE. More...

#include <particleDynamicsBase.h>

Inheritance diagram for olb::particles::dynamics::ParticleDetachmentDynamics< T, PARTICLETYPE >:

Collaboration diagram for olb::particles::dynamics::ParticleDetachmentDynamics< T, PARTICLETYPE >:

Public Member Functions
	ParticleDetachmentDynamics (SolidBoundary< T, PARTICLETYPE::d > &solidBoundary, Vector< T, PARTICLETYPE::d > &mainFlowDirection, T tiltThreshold=0.3 *M_PI)
	Constructor.

void	process (Particle< T, PARTICLETYPE > &particle, T timeStepSize) override
	Procesisng step.

Public Member Functions inherited from olb::particles::dynamics::VerletParticleDynamics< T, PARTICLETYPE, conditions::active_particles >
	VerletParticleDynamics ()
	Constructor.

void	process (Particle< T, PARTICLETYPE > &particle, T timeStepSize) override
	Procesisng step.

Public Member Functions inherited from olb::particles::dynamics::ParticleDynamics< T, PARTICLETYPE >
virtual	~ParticleDynamics ()
	Destructor: virtual to enable inheritance.

std::string &	getName ()
	read and write access to name

std::string const &	getName () const
	read only access to name

Detailed Description

template<typename T, typename PARTICLETYPE>
class olb::particles::dynamics::ParticleDetachmentDynamics< T, PARTICLETYPE >

Verlet dynamics for particles aware of their DYNAMIC_STATE.

For now, manly used for wall-flow application
- source: https://gitlab.com/openlb/olb/-/blob/private/nicolas/main4/apps/nicolas/A_WallFlow/wallFlow/README.md
- partially assuming axis aligned surfaces and provision of main flow direction

Definition at line 252 of file particleDynamicsBase.h.

Constructor & Destructor Documentation

◆ ParticleDetachmentDynamics()

template<typename T , typename PARTICLETYPE >

olb::particles::dynamics::ParticleDetachmentDynamics< T, PARTICLETYPE >::ParticleDetachmentDynamics	(	SolidBoundary< T, PARTICLETYPE::d > &	solidBoundary,
		Vector< T, PARTICLETYPE::d > &	mainFlowDirection,
		T	tiltThreshold = 0.3*M_PI )

Constructor.

Definition at line 373 of file particleDynamicsBase.hh.

  : _solidBoundary(solidBoundary), _mainFlowDirection(mainFlowDirection), _tiltThreshold(tiltThreshold)
{
  this->getName() = "ParticleDetachmentDynamics";
  static_assert(PARTICLETYPE::template providesNested<descriptors::SURFACE,descriptors::ANGLE>(),
                "Field SURFACE:ANGLE has to be provided");
  static_assert(PARTICLETYPE::template providesNested<descriptors::DYNBEHAVIOUR,descriptors::DETACHING>(),
                "Field DYNBEHAVIOUR:DETACHING has to be provided");
}

References olb::particles::dynamics::ParticleDynamics< T, PARTICLETYPE >::getName().

Here is the call graph for this function:

Member Function Documentation

◆ process()

template<typename T , typename PARTICLETYPE >

void olb::particles::dynamics::ParticleDetachmentDynamics< T, PARTICLETYPE >::process	(	Particle< T, PARTICLETYPE > &	particle,
		T	timeStepSize )

overridevirtual

Procesisng step.

Implements olb::particles::dynamics::ParticleDynamics< T, PARTICLETYPE >.

Definition at line 386 of file particleDynamicsBase.hh.

{
  using namespace particles::access;
  using namespace descriptors;
 
  //Check for valid particles condition
  doWhenMeetingCondition<T,PARTICLETYPE,conditions::valid_particles>( particle,[&](){
    //Check current dynamic state
    bool detaching = isDetaching( particle );
    //State: Normal motion
    if (!detaching){
      //Execute process of VerletParticleDynamcis (valid, not detaching, active)
      VerletParticleDynamics<T,PARTICLETYPE,conditions::active_particles>
        ::process(particle,timeStepSize);
    }
    //State: Detaching
    else{
      //Check adhesion threshold (set active, if passed once)
      bool isAdhering = interaction::checkAdhesion( _solidBoundary, _mainFlowDirection, particle );
      //Perform detachment dynamics (valid, detaching, adhering)
      if (!isAdhering){
        //Calculate angular acceleration
        auto angularAcceleration = getAngAcceleration( particle );
        //Verlet algorithm
        particles::dynamics::velocityVerletRotation(
          particle, timeStepSize, timeStepSize*timeStepSize, angularAcceleration );
        //Check if rotation matrix provided and update
        if constexpr ( providesRotationMatrix<PARTICLETYPE>() ) {
          updateRotationMatrix( particle );
        }
        //Handle detachment from surface
        interaction::handleDetachment( _solidBoundary, _mainFlowDirection, particle );
        //Reevaluate state: Check, whether state has to be changed (e.g. detachment finished)
        interaction::evaluateDetachmentState( _solidBoundary, particle, _tiltThreshold );
      } //if (!isAdhering)
    }
  }); //doWhenMeetingCondition<T,PARTICLETYPE,conditions::valid_particles>
}

References olb::particles::interaction::checkAdhesion(), olb::particles::interaction::evaluateDetachmentState(), olb::particles::interaction::handleDetachment(), olb::particles::dynamics::updateRotationMatrix(), and olb::particles::dynamics::velocityVerletRotation().