olb::MorsiDragModel< T, Lattice, Particle > Class Template Reference

Class to compute the standard drag coefficient as in Morsi and Alexander (1972). More...

#include <dragModels3D.h>

Inheritance diagram for olb::MorsiDragModel< T, Lattice, Particle >:

Collaboration diagram for olb::MorsiDragModel< T, Lattice, Particle >:

Public Member Functions
	MorsiDragModel (UnitConverter< T, Lattice > &converter)
	Constructor.
virtual T	operator() (Particle< T > *p, T latticeVelF[], T latticeVelP[], int globicFull[], T continuousPhaseFraction=1.) override
	Returns the scalar drag coefficient. globicFull = { globic, latticeRoundedP[0, ..., 2] }.
std::string	name () override
Public Member Functions inherited from olb::DragModelBase< T, Lattice, Particle >
	DragModelBase (UnitConverter< T, Lattice > &converter)
	Constructor.

Additional Inherited Members
Protected Attributes inherited from olb::DragModelBase< T, Lattice, Particle >
UnitConverter< T, Lattice > &	_converter
	Returns the scalar drag coefficient to overload.
Protected Attributes inherited from olb::DragModel< T, Particle >
std::shared_ptr< ParticleReynoldsNumber< T, Particle > >	_reP
	Functional to compute particle Reynolds number.

Detailed Description

template<typename T, typename Lattice, template< typename V > class Particle>
class olb::MorsiDragModel< T, Lattice, Particle >

Class to compute the standard drag coefficient as in Morsi and Alexander (1972).

Definition at line 81 of file dragModels3D.h.

Constructor & Destructor Documentation

MorsiDragModel()

template<typename T , typename Lattice , template< typename V > class Particle>

olb::MorsiDragModel< T, Lattice, Particle >::MorsiDragModel

(

UnitConverter< T, Lattice > &

converter

)

Constructor.

Definition at line 60 of file dragModels3D.hh.

  : DragModelBase<T,Lattice,Particle>(converter)
{
  this->_reP = std::make_shared<NewtonianParticleReynoldsNumber<T,Lattice,Particle> > (this->_converter);
}

References olb::DragModelBase< T, Lattice, Particle >::_converter, and olb::DragModel< T, Particle >::_reP.

Member Function Documentation

name()

template<typename T , typename Lattice , template< typename V > class Particle>

std::string olb::MorsiDragModel< T, Lattice, Particle >::name ( )

inlineoverridevirtual

Implements olb::DragModel< T, Particle >.

Reimplemented in olb::PowerLawMorsiDragModel< T, Lattice, Particle >.

Definition at line 87 of file dragModels3D.h.

{return "MorsiDragModel";}

operator()()

template<typename T , typename Lattice , template< typename V > class Particle>

T olb::MorsiDragModel< T, Lattice, Particle >::operator() ( Particle< T > * p, T latticeVelF[], T latticeVelP[], int globicFull[], T continuousPhaseFraction = 1. )

overridevirtual

Returns the scalar drag coefficient. globicFull = { globic, latticeRoundedP[0, ..., 2] }.

Implements olb::DragModel< T, Particle >.

Definition at line 67 of file dragModels3D.hh.

{
  T physVelRelative = this->_converter.getPhysVelocity (
                        util::sqrt( util::pow(latticeVelF[0] - latticeVelP[0],2) +
                                    util::pow(latticeVelF[1] - latticeVelP[1],2) +
                                    util::pow(latticeVelF[2] - latticeVelP[2],2) ) );
 
  T ReP = continuousPhaseFraction * this->_reP->operator() (p, physVelRelative, globicFull);
 
  T a[3] = {T(), T(), T()};
  if (ReP < 0.1) {
    a[0] = 0.0;
    a[1] = 24.0;
    a[2] = 0.0;
  }
  else if (ReP < 1.0) {
    a[0] = 3.69;
    a[1] = 22.73;
    a[2] = 0.0903;
  }
  else if (ReP < 10.0) {
    a[0] = 1.222;
    a[1] = 29.16667;
    a[2] =-3.8889;
  }
  else if (ReP < 100.0) {
    a[0] = 0.6167;
    a[1] = 46.5;
    a[2] =-116.67;
  }
  else if (ReP < 1000.0) {
    a[0] = 0.3644;
    a[1] = 98.33;
    a[2] =-2778;
  }
  else if (ReP < 5000.0) {
    a[0] = 0.357;
    a[1] = 148.62;
    a[2] =-4.75e4;
  }
  else if (ReP < 10000.0) {
    a[0] = 0.46;
    a[1] =-490.546;
    a[2] = 57.87e4;
  }
  else {
    a[0] = 0.5191;
    a[1] =-1662.5;
    a[2] = 5.4167e6;
  }
 
  return ( a[0] + a[1]/ReP + a[2]/(ReP*ReP) );
}

References 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:

• src/particles/subgrid3DLegacyFramework/twoWayCouplings/dragModels3D.h
• src/particles/subgrid3DLegacyFramework/twoWayCouplings/dragModels3D.hh