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

Back-coupling is performed on the cell containing the particle and its neighbours within a cube of one lattice spacing, as per in Maier et al. More...

#include <backCouplingModels.h>

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

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

Public Member Functions
	CubicDeltaBackCouplingModel (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, SuperGeometry< T, 3 > &sGeometry, int overlap)
	Constructor.
virtual bool	operator() (Particle< T > *p, int globic, int material, int subCycles=1) override
	Class operator to apply the coupling.
Public Member Functions inherited from olb::BaseNonLocalBackCouplingModel< T, Lattice, Particle >
virtual void	communicate () override
	Communicates POPULATION and FORCE fields if the model is non-local.
	BaseNonLocalBackCouplingModel (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, SuperGeometry< T, 3 > &sGeometry, int overlap)
	Constructor.
Public Member Functions inherited from olb::BaseBackCouplingModel< T, Lattice, Particle >
virtual void	resetExternalField (int material) override
	Resets external field.

Protected Attributes
int	_range = 1
T	_delta [4][4][4] = { T() }
std::shared_ptr< SuperLatticeSmoothDiracDelta3D< T, Lattice > >	_cubicDeltaFunctional
Protected Attributes inherited from olb::BaseBackCouplingModel< T, Lattice, Particle >
UnitConverter< T, Lattice > &	_converter
SuperGeometry< T, 3 > &	_sGeometry
SuperLattice< T, Lattice > &	_sLattice
SuperCommunicator< T, SuperLattice< T, Lattice > >	_commPopulation

Additional Inherited Members
Protected Member Functions inherited from olb::BaseBackCouplingModel< T, Lattice, Particle >
	BaseBackCouplingModel (UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice, SuperGeometry< T, 3 > &sGeometry, int overlap)
	Constructor.

Detailed Description

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

Back-coupling is performed on the cell containing the particle and its neighbours within a cube of one lattice spacing, as per in Maier et al.

(2017). Input parameters in attice units.

Definition at line 127 of file backCouplingModels.h.

Constructor & Destructor Documentation

CubicDeltaBackCouplingModel()

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

olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::CubicDeltaBackCouplingModel	(	UnitConverter< T, Lattice > &	converter,
		SuperLattice< T, Lattice > &	sLattice,
		SuperGeometry< T, 3 > &	sGeometry,
		int	overlap )

Constructor.

Definition at line 162 of file backCouplingModels.hh.

  : BaseNonLocalBackCouplingModel<T,Lattice,Particle>(converter, sLattice, sGeometry, overlap)
{
  _cubicDeltaFunctional = std::make_shared<SuperLatticeSmoothDiracDelta3D<T, Lattice> > (
                            this->_sLattice, this->_converter, this->_sGeometry );
}

References olb::BaseBackCouplingModel< T, Lattice, Particle >::_converter, olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::_cubicDeltaFunctional, olb::BaseBackCouplingModel< T, Lattice, Particle >::_sGeometry, and olb::BaseBackCouplingModel< T, Lattice, Particle >::_sLattice.

Member Function Documentation

operator()()

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

bool olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::operator() ( Particle< T > * p, int globic, int material, int subCycles = 1 )

overridevirtual

Class operator to apply the coupling.

Implements olb::BackCouplingModel< T, Particle >.

Definition at line 174 of file backCouplingModels.hh.

{
  int locIC = this->_sLattice.getLoadBalancer().loc(globic);
 
  // reading the force from the value stored inside the particle
  std::vector<T> physForceP = p->getStoreForce(); // physical force acting on the particle
 
  T latticeForceP[3] = { physForceP[0] / this->_converter.getConversionFactorForce(),
                         physForceP[1] / this->_converter.getConversionFactorForce(),
                         physForceP[2] / this->_converter.getConversionFactorForce()
                       }; // dimensionless force acting on the particle
 
  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 );
 
  // smooth Dirac delta
  this->_cubicDeltaFunctional->operator() (_delta, physPosP, globic);
 
  T tempDelta = T();
  FieldD<T,Lattice,descriptors::FORCE> F; // dimensionless smoothed force
 
  for (int i = -_range; i <= _range; ++i) {
    for (int j = -_range; j <= _range; ++j) {
      for (int k = -_range; k <= _range; ++k) {
        if (this->_sGeometry.getBlockGeometry(locIC).getMaterial(
              latticeRoundedPosP[0] + i, latticeRoundedPosP[1] + j,
              latticeRoundedPosP[2] + k) == material) {
 
          tempDelta = _delta[i + _range][j + _range][k + _range];
 
          F[0] = -latticeForceP[0] * tempDelta / (T)(subCycles);
          F[1] = -latticeForceP[1] * tempDelta / (T)(subCycles);
          F[2] = -latticeForceP[2] * tempDelta / (T)(subCycles);
 
          this->_sLattice.getBlock(locIC).get(
            latticeRoundedPosP[0] + i,
            latticeRoundedPosP[1] + j,
            latticeRoundedPosP[2] + k
          ).template getFieldPointer<descriptors::FORCE>() += F;
        }
      }
    }
  }
  return true;
}

References olb::Vector< T, D >::getBlock().

Here is the call graph for this function:

Member Data Documentation

_cubicDeltaFunctional

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

std::shared_ptr<SuperLatticeSmoothDiracDelta3D<T, Lattice> > olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::_cubicDeltaFunctional

protected

Definition at line 139 of file backCouplingModels.h.

_delta

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

T olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::_delta[4][4][4] = { T() }

protected

Definition at line 138 of file backCouplingModels.h.

{ T() };

_range

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

int olb::CubicDeltaBackCouplingModel< T, Lattice, Particle >::_range = 1

protected

Definition at line 137 of file backCouplingModels.h.

---

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

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