olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR > Class Template Reference

#include <stokesDragForce3D.h>

Inheritance diagram for olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >:

Collaboration diagram for olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >:

Public Member Functions
	StokesDragForce3D (SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &getVel, T dT, T mu)
	Constructor, FluidVelocity, physicalTimeStep, physicalDynamicViscosity.
	StokesDragForce3D (SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &getVel, UnitConverter< T, DESCRIPTOR > const &converter)
	Constructor using values from converter.
	StokesDragForce3D (SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &getVel, UnitConverter< T, DESCRIPTOR > const &converter, T scaleFactor)
	Constructor using values from converter.
	~StokesDragForce3D () override
	Destructor.
void	applyForce (typename std::deque< PARTICLETYPE< T > >::iterator p, int pInt, ParticleSystem3D< T, PARTICLETYPE > &psSys) override
	6 Pi r mu (u_f-u_p)
void	computeForce (int pInt, ParticleSystem3D< T, PARTICLETYPE > *psSys, T force[3])
	Compute Force for subgrid scale particles.
Public Member Functions inherited from olb::Force3D< T, PARTICLETYPE >
	Force3D ()
	Force3D (Force3D< T, PARTICLETYPE > &)
	Force3D (const Force3D< T, PARTICLETYPE > &)
virtual	~Force3D ()

Additional Inherited Members
Protected Attributes inherited from olb::Force3D< T, PARTICLETYPE >
OstreamManager	clout

Detailed Description

template<typename T, template< typename U > class PARTICLETYPE, typename DESCRIPTOR>
class olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >

Definition at line 44 of file stokesDragForce3D.h.

Constructor & Destructor Documentation

StokesDragForce3D() [1/3]

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::StokesDragForce3D	(	SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &	getVel,
		T	dT,
		T	mu )

Constructor, FluidVelocity, physicalTimeStep, physicalDynamicViscosity.

Definition at line 37 of file stokesDragForce3D.hh.

  : Force3D<T, PARTICLETYPE>(),
    _getVel(getVel),
    _mu(mu)
{
  _C1 = 6. * M_PI * _mu * dT ;
  _dTinv = 1. / dT;
  _scaleFactor = 1.;
}

References M_PI.

StokesDragForce3D() [2/3]

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::StokesDragForce3D	(	SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &	getVel,
		UnitConverter< T, DESCRIPTOR > const &	converter )

Constructor using values from converter.

Definition at line 48 of file stokesDragForce3D.hh.

  : Force3D<T, PARTICLETYPE>(),
    _getVel(getVel)
{
  //implicit formulation
  _C1 = 6. * M_PI * converter.getPhysViscosity() * converter.getPhysDensity() * converter.getConversionFactorTime();
  _mu = converter.getPhysViscosity() * converter.getPhysDensity();
  // explicit formulation
  //  _C1 = 6. * M_PI * converter.getDynamicViscosity();
  _dTinv = 1. / converter.getConversionFactorTime();
  _scaleFactor = 1. ;
}

References olb::UnitConverter< T, DESCRIPTOR >::getConversionFactorTime(), olb::UnitConverter< T, DESCRIPTOR >::getPhysDensity(), olb::UnitConverter< T, DESCRIPTOR >::getPhysViscosity(), and M_PI.

Here is the call graph for this function:

StokesDragForce3D() [3/3]

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::StokesDragForce3D	(	SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &	getVel,
		UnitConverter< T, DESCRIPTOR > const &	converter,
		T	scaleFactor )

Constructor using values from converter.

Definition at line 62 of file stokesDragForce3D.hh.

  : Force3D<T, PARTICLETYPE>(),
    _getVel(getVel),
    _scaleFactor(scaleFactor)
{
  //implicit formulation
  _C1 = 6. * M_PI * converter.getPhysViscosity() * converter.getPhysDensity() * converter.getConversionFactorTime();
  _mu = converter.getPhysViscosity() * converter.getPhysDensity();
  // explicit formulation
  //  _C1 = 6. * M_PI * converter.getDynamicViscosity();
  _dTinv = 1. / converter.getConversionFactorTime();
}

References olb::UnitConverter< T, DESCRIPTOR >::getConversionFactorTime(), olb::UnitConverter< T, DESCRIPTOR >::getPhysDensity(), olb::UnitConverter< T, DESCRIPTOR >::getPhysViscosity(), and M_PI.

Here is the call graph for this function:

~StokesDragForce3D()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::~StokesDragForce3D ( )

inlineoverride

Destructor.

Definition at line 54 of file stokesDragForce3D.h.

{}

Member Function Documentation

applyForce()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

void olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::applyForce ( typename std::deque< PARTICLETYPE< T > >::iterator p, int pInt, ParticleSystem3D< T, PARTICLETYPE > & psSys )

overridevirtual

6 Pi r mu (u_f-u_p)

Implements olb::Force3D< T, PARTICLETYPE >.

Definition at line 77 of file stokesDragForce3D.hh.

{
  T fluidVel[3] = {0., 0., 0.};
 
  //implicit formulation
  _getVel(fluidVel, &p->getPos()[0], p->getCuboid());
  fluidVel[0] *= _scaleFactor;
  fluidVel[1] *= _scaleFactor;
  fluidVel[2] *= _scaleFactor;
 
  T c = _C1 * p->getRad() * p->getInvMass();
  T C2 = 1. / (1. + c);
 
  // p->getVel() is particle velocity of the last time step
  // formulation of new force with particle velocity of the last time step with
  // implicit Euler
  p->getForce()[0] += p->getMass() * _dTinv
                      * ((c * fluidVel[0] + p->getVel()[0]) * C2 - p->getVel()[0]);
  p->getForce()[1] += p->getMass() * _dTinv
                      * ((c * fluidVel[1] + p->getVel()[1]) * C2 - p->getVel()[1]);
  p->getForce()[2] += p->getMass() * _dTinv
                      * ((c * fluidVel[2] + p->getVel()[2]) * C2 - p->getVel()[2]);
}

computeForce()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >

void olb::StokesDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::computeForce	(	int	pInt,
		ParticleSystem3D< T, PARTICLETYPE > *	psSys,
		T	force[3] )

Compute Force for subgrid scale particles.

Definition at line 104 of file stokesDragForce3D.hh.

{
  T fluidVel[3] = {0., 0., 0.};
 
  _getVel(fluidVel, &psSys->operator[](pInt).getPos()[0], psSys->operator[](pInt).getCuboid());
 
  T c = _C1 * psSys->operator[](pInt).getRad() * psSys->operator[](pInt).getInvMass();
  T C2 = 1. / (1. + c);
  T mass = psSys->operator[](pInt).getMass();
  std::vector<T> vel = psSys->operator[](pInt).getVel();
  force[0] = mass * _dTinv * ((c * fluidVel[0] + vel[0]) * C2 - vel[0]);
  force[1] = mass * _dTinv * ((c * fluidVel[1] + vel[1]) * C2 - vel[1]);
  force[2] = mass * _dTinv * ((c * fluidVel[2] + vel[2]) * C2 - vel[2]);
}

---

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

• src/particles/subgrid3DLegacyFramework/forces/stokesDragForce3D.h
• src/particles/subgrid3DLegacyFramework/forces/stokesDragForce3D.hh