OpenLB 1.7
Loading...
Searching...
No Matches
Public Member Functions | List of all members
olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR > Class Template Reference

#include <schillerNaumannDragForce3D.h>

+ Inheritance diagram for olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >:
+ Collaboration diagram for olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >:

Public Member Functions

 SchillerNaumannDragForce3D (SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > &getVel, T dynVisc, T physDensity)
 
 ~SchillerNaumannDragForce3D () override
 
void applyForce (typename std::deque< PARTICLETYPE< T > >::iterator p, int pInt, ParticleSystem3D< T, PARTICLETYPE > &psSys) override
 
- 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::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >

Definition at line 39 of file schillerNaumannDragForce3D.h.

Constructor & Destructor Documentation

◆ SchillerNaumannDragForce3D()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >
olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::SchillerNaumannDragForce3D ( SuperLatticeInterpPhysVelocity3D< T, DESCRIPTOR > & getVel,
T dynVisc,
T physDensity )

Definition at line 38 of file schillerNaumannDragForce3D.hh.

38 :
39 Force3D<T, PARTICLETYPE>(), _getVel(getVel), _dynVisc(dynVisc), _physDensity(physDensity)
40{
41// this->_name = "SchillerNaumannDragForce3D";
42}

◆ ~SchillerNaumannDragForce3D()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >
olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::~SchillerNaumannDragForce3D ( )
inlineoverride

Definition at line 45 of file schillerNaumannDragForce3D.h.

45{};

Member Function Documentation

◆ applyForce()

template<typename T , template< typename U > class PARTICLETYPE, typename DESCRIPTOR >
void olb::SchillerNaumannDragForce3D< T, PARTICLETYPE, DESCRIPTOR >::applyForce ( typename std::deque< PARTICLETYPE< T > >::iterator p,
int pInt,
ParticleSystem3D< T, PARTICLETYPE > & psSys )
overridevirtual

Implements olb::Force3D< T, PARTICLETYPE >.

Definition at line 45 of file schillerNaumannDragForce3D.hh.

48{
49
50 Vector<T, 3> force = {T(0), T(0), T(0)} ;
51
52 T fluidVel[3] = {0., 0., 0.};
53 _getVel(fluidVel, &p->getPos()[0], p->getCuboid());
54
55 T fluidVelAbs = util::sqrt(std::pow(fluidVel[0], 2.) + util::pow(fluidVel[1], 2.) + util::pow(fluidVel[2], 2.));
56 T partVelAbs = util::sqrt(std::pow(p->getVel()[0], 2.) + util::pow(p->getVel()[1], 2.) + util::pow(p->getVel()[2], 2.));
57 T particleRe = 2. * p->getRad() * abs(fluidVelAbs - partVelAbs) * _physDensity / _dynVisc;
58 T coeffSN = 1. + 0.15 * util::pow(particleRe, 0.687);
59
60 for (int i = 0; i < 3; i++) {
61
62 force[i] = -1. * 6 * M_PI * p->getRad() * (p->getVel()[i] - fluidVel[i]) * _dynVisc * coeffSN;
63 p->getForce()[i] += force[i] ;
64 }
65
66}
#define M_PI
cpu::simd::Pack< T > sqrt(cpu::simd::Pack< T > value)
Definition pack.h:100
cpu::simd::Pack< T > pow(cpu::simd::Pack< T > base, cpu::simd::Pack< T > exp)
Definition pack.h:112
std::enable_if_t< std::is_arithmetic< T >::type::value, T > abs(T x)
Definition util.h:396

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

+ Here is the call graph for this function:

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