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

Abstact class for all the local forward-coupling models, viz., momentum coupling from fluid to particle. More...

#include <forwardCouplingModels3D.h>

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

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

Public Member Functions
virtual bool	operator() (Particle< T > *p, int globic) override
	Class operator to apply the coupling, for overload.

Protected Member Functions
	LocalBaseForwardCouplingModel (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, SuperGeometry< T, 3 > &sGeometry, std::shared_ptr< DragModel< T, Particle > > dragModel)
	Constructor.
Protected Member Functions inherited from olb::BaseForwardCouplingModel< T, Lattice, Particle >
	BaseForwardCouplingModel (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, SuperGeometry< T, 3 > &sGeometry, std::shared_ptr< DragModel< T, Particle > > dragModel)
	Constructor.
Protected Member Functions inherited from olb::ForwardCouplingModel< T, Particle >
	ForwardCouplingModel (SuperGeometry< T, 3 > &sGeometry, std::shared_ptr< DragModel< T, Particle > > dragModel)
	Constructor.

Additional Inherited Members
Protected Attributes inherited from olb::BaseForwardCouplingModel< T, Lattice, Particle >
UnitConverter< T, Lattice > &	_converter
SuperLattice< T, Lattice > &	_sLattice
std::shared_ptr< SuperLatticeInterpDensity3Degree3D< T, Lattice > >	_interpLatticeDensity
std::shared_ptr< SuperLatticeInterpPhysVelocity3D< T, Lattice > >	_interpLatticeVelocity
std::shared_ptr< TwoWayHelperFunctional< T, Lattice > >	_momentumExchange
Protected Attributes inherited from olb::ForwardCouplingModel< T, Particle >
SuperGeometry< T, 3 > &	_sGeometry
std::shared_ptr< DragModel< T, Particle > >	_dragModel

Detailed Description

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

Abstact class for all the local forward-coupling models, viz., momentum coupling from fluid to particle.

Input parameters in attice units.

Definition at line 78 of file forwardCouplingModels3D.h.

Constructor & Destructor Documentation

LocalBaseForwardCouplingModel()

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

olb::LocalBaseForwardCouplingModel< T, Lattice, Particle >::LocalBaseForwardCouplingModel ( UnitConverter< T, Lattice > & converter, SuperLattice< T, Lattice > & sLattice, SuperGeometry< T, 3 > & sGeometry, std::shared_ptr< DragModel< T, Particle > > dragModel )

protected

Constructor.

Definition at line 67 of file forwardCouplingModels3D.hh.

  : BaseForwardCouplingModel<T,Lattice,Particle>(converter, sLattice, sGeometry, dragModel)
{}

Member Function Documentation

operator()()

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

bool olb::LocalBaseForwardCouplingModel< T, Lattice, Particle >::operator() ( Particle< T > * p, int globic )

overridevirtual

Class operator to apply the coupling, for overload.

Getting the particle and its containing cell's position

Computing drag force in dimensionless units

Computing physical drag force

Updating the particle

Implements olb::ForwardCouplingModel< T, Particle >.

Definition at line 76 of file forwardCouplingModels3D.hh.

{
  T physPosP[3] = { p->getPos()[0],
                    p->getPos()[1],
                    p->getPos()[2]
                  }; // particle's physical position
  // particle's dimensionless position, rounded at neighbouring voxel
  int latticeRoundedPosP[3] = {0, 0, 0};
  this->_sLattice.getCuboidGeometry().get(globic).getLatticeR (
    latticeRoundedPosP, physPosP );
  // { globic, latticeRoundedP[0, ..., 2] }
  int globicFull[4] = { globic,
                        latticeRoundedPosP[0],
                        latticeRoundedPosP[1],
                        latticeRoundedPosP[2]
                      };
 
  // Particle's velocity
  T physVelP[3] = { p->getVel()[0],
                    p->getVel()[1],
                    p->getVel()[2]
                  }; // Physical
  // Lattice
  T latticeVelP[3] = { this->_converter.getLatticeVelocity(physVelP[0]),
                       this->_converter.getLatticeVelocity(physVelP[1]),
                       this->_converter.getLatticeVelocity(physVelP[2])
                     }; // particle's dimensionless velocity
 
  // Lattice's velocity at particle's location
  T physVelF[3] = {T(), T(), T()}; // Physical
  this->_interpLatticeVelocity->operator() (physVelF, physPosP, globic);
  // Lattice
  T latticeVelF[3] = { this->_converter.getLatticeVelocity(physVelF[0]),
                       this->_converter.getLatticeVelocity(physVelF[1]),
                       this->_converter.getLatticeVelocity(physVelF[2])
                     }; // Lattice's dimensionless velocity at particle's location
 
  // Computing fluid-particle momentum transfer
  T gF[3] = {T(), T(), T()}; // force density gF
  this->_momentumExchange->operator() (gF, latticeVelF, latticeVelP, physPosP, latticeRoundedPosP, globic);
 
  // Computing drag coefficient
  T Cd = this->_dragModel->operator()(p, latticeVelF, latticeVelP, globicFull);
#ifdef VERBOSE
  std:: cout << Cd
             << " physPosP=(" << physPosP[0] << ", " << physPosP[1] << ", " << physPosP[2] << ") "
             << " physVelP=(" << physVelP[0] << ", " << physVelP[1] << ", " << physVelP[2] << ") "
             << " physVelF=(" << physVelF[0] << ", " << physVelF[1] << ", " << physVelF[2] << ") "
             << std::endl;
#endif
  /*
  if (Cd > 100.)
    throw std::range_error ( "LocalBaseForwardCouplingModel::operator(). Cd="+std::to_string(Cd)
                      + "\nphysVelP=(" + std::to_string(physVelP[0]) + ", " + std::to_string(physVelP[1]) + ", " + std::to_string(physVelP[2]) + ")"
                      + "\nphysVelF=(" + std::to_string(physVelF[0]) + ", " + std::to_string(physVelF[1]) + ", " + std::to_string(physVelF[2]) + ")"
                      + "\n");
  */
 
  T latticePRad = p->getRad() / this->_converter.getConversionFactorLength();
  T latticeForceP[3] = { .5 * Cd * M_PI*util::pow(latticePRad,2) * gF[0] * (latticeVelF[0] - latticeVelP[0]),
                         .5 * Cd * M_PI*util::pow(latticePRad,2) * gF[1] * (latticeVelF[1] - latticeVelP[1]),
                         .5 * Cd * M_PI*util::pow(latticePRad,2) * gF[2] * (latticeVelF[2] - latticeVelP[2])
                       }; // dimensionless force acting on the particle
 
  std::vector<T> physForceP(3, T()); // physical force acting on the particle
  physForceP[0] = latticeForceP[0] * this->_converter.getConversionFactorForce();
  physForceP[1] = latticeForceP[1] * this->_converter.getConversionFactorForce();
  physForceP[2] = latticeForceP[2] * this->_converter.getConversionFactorForce();
 
  p->setStoreForce(physForceP);
  return true;
}

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

Here is the call graph for this function:

---

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

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