olb::LaddMomentumExchange< T, Lattice > Class Template Reference

Using Ladd mechanism. More...

#include <twoWayHelperFunctionals.h>

Inheritance diagram for olb::LaddMomentumExchange< T, Lattice >:

Collaboration diagram for olb::LaddMomentumExchange< T, Lattice >:

Public Member Functions
	LaddMomentumExchange (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, std::shared_ptr< SuperLatticeInterpDensity3Degree3D< T, Lattice > > interpLatticeDensity, std::shared_ptr< SuperLatticeInterpPhysVelocity3D< T, Lattice > > interpLatticeVelocity)
	Constructor.
virtual bool	operator() (T gF[], T latticeVelF[], T latticeVelP[], T physPosP[], int latticeRoundedP[], int globic) override
	Computes the momentum transfer from fluid to particle.
Public Member Functions inherited from olb::TwoWayHelperFunctional< T, Lattice >
virtual	~TwoWayHelperFunctional ()

Additional Inherited Members
Protected Member Functions inherited from olb::TwoWayHelperFunctional< T, Lattice >
	TwoWayHelperFunctional (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice)
	Constructor.
Protected Attributes inherited from olb::TwoWayHelperFunctional< T, Lattice >
UnitConverter< T, Lattice > &	_converter
SuperLattice< T, Lattice > &	_sLattice
std::shared_ptr< SuperLatticeInterpDensity3Degree3D< T, Lattice > >	_interpLatticeDensity
std::shared_ptr< SuperLatticeInterpPhysVelocity3D< T, Lattice > >	_interpLatticeVelocity

Detailed Description

template<typename T, typename Lattice>
class olb::LaddMomentumExchange< T, Lattice >

Using Ladd mechanism.

Definition at line 76 of file twoWayHelperFunctionals.h.

Constructor & Destructor Documentation

LaddMomentumExchange()

template<typename T , typename Lattice >

olb::LaddMomentumExchange< T, Lattice >::LaddMomentumExchange	(	UnitConverter< T, Lattice > &	converter,
		SuperLattice< T, Lattice > &	sLattice,
		std::shared_ptr< SuperLatticeInterpDensity3Degree3D< T, Lattice > >	interpLatticeDensity,
		std::shared_ptr< SuperLatticeInterpPhysVelocity3D< T, Lattice > >	interpLatticeVelocity )

Constructor.

Definition at line 89 of file twoWayHelperFunctionals.hh.

  : TwoWayHelperFunctional<T, Lattice>(converter, sLattice)
{
  this->_interpLatticeDensity = interpLatticeDensity;
  this->_interpLatticeVelocity = interpLatticeVelocity;
}

References olb::TwoWayHelperFunctional< T, Lattice >::_interpLatticeDensity, and olb::TwoWayHelperFunctional< T, Lattice >::_interpLatticeVelocity.

Member Function Documentation

operator()()

template<typename T , typename Lattice >

bool olb::LaddMomentumExchange< T, Lattice >::operator() ( T gF[], T latticeVelF[], T latticeVelP[], T physPosP[], int latticeRoundedP[], int globic )

overridevirtual

Computes the momentum transfer from fluid to particle.

Implements olb::TwoWayHelperFunctional< T, Lattice >.

Definition at line 101 of file twoWayHelperFunctionals.hh.

{
  T physLatticeL = this->_converter.getConversionFactorLength();
 
  // force density gF
  gF[0] = T();
  gF[1] = T();
  gF[2] = T();
 
  T fiPop = T();
  T sp = T(); // dot product for particle velocity
  T faPos[3] = {T(), T(), T()}; // fAlphaPosition = particle position
  T fbPos[3] = {T(), T(), T()}; // fBetaPosition = neighbor position to particle position in direction iPop
 
  T fa[Lattice::q] = { T() }; // fAlpha = interpolated density to fAlphaPosition
  T fb[Lattice::q] = { T() }; // fBeta = interpolated density to fBetaPosition
  T lFU[3] = {T(), T(), T()};
 
  // runs through all q discrete velocity directions
  for (unsigned iPop = 0; iPop < Lattice::q; ++iPop) {
    // physical position on neighbor to get pre-streaming collision part
    faPos[0] = physPosP[0] + physLatticeL * descriptors::c<Lattice>(iPop,0);
    faPos[1] = physPosP[1] + physLatticeL * descriptors::c<Lattice>(iPop,1);
    faPos[2] = physPosP[2] + physLatticeL * descriptors::c<Lattice>(iPop,2);
    // Lagrange interpolated polynomial to get density on particle position
    this->_interpLatticeDensity->operator() (fa, faPos, globic);
 
    // physical position on neighbor to get pre-streaming collision part
    fbPos[0] = physPosP[0] - physLatticeL * descriptors::c<Lattice>(iPop,0);
    fbPos[1] = physPosP[1] - physLatticeL * descriptors::c<Lattice>(iPop,1);
    fbPos[2] = physPosP[2] - physLatticeL * descriptors::c<Lattice>(iPop,2);
    // Lagrange interpolated polynomial to get density on particle position
    this->_interpLatticeDensity->operator() (fb, fbPos, globic);
 
    // fiPop = density on fBetaPosition in direction iPop
    fiPop = fb[descriptors::opposite<Lattice >(iPop)];
    // Get f_l of the boundary cell
    // add density fAlphaL of opposite direction to iPop
    fiPop -= fa[iPop];
 
    // physical velocity
    lFU[0] = -descriptors::c<Lattice>(iPop,0) * fiPop;
    lFU[1] = -descriptors::c<Lattice>(iPop,1) * fiPop;
    lFU[2] = -descriptors::c<Lattice>(iPop,2) * fiPop;
 
    // point product
    sp = descriptors::c<Lattice>(iPop,0) * latticeVelP[0] + descriptors::c<Lattice>(iPop,1) * latticeVelP[1]
         + descriptors::c<Lattice>(iPop,2) * latticeVelP[2];
    sp *= 2. * descriptors::invCs2<T,Lattice>()
          * descriptors::t<T,Lattice>(iPop);
 
    // external force density that acts on particles
    gF[0] += (lFU[0] - descriptors::c<Lattice>(iPop,0) * (sp));
    gF[1] += (lFU[1] - descriptors::c<Lattice>(iPop,1) * (sp));
    gF[2] += (lFU[2] - descriptors::c<Lattice>(iPop,2) * (sp));
  }
  gF[0] = util::fabs(gF[0]);
  gF[1] = util::fabs(gF[1]);
  gF[2] = util::fabs(gF[2]);
 
  return true;
}

References olb::util::fabs().

Here is the call graph for this function:

---

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

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