olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR > Class Template Reference

This class interpolates missing f_i from values near the boundary to get a more stable outflow condition for the density. More...

#include <advectionDiffusionBoundaryPostProcessor3D.h>

Inheritance diagram for olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >:

Collaboration diagram for olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >:

Public Member Functions
	ConvectionBoundaryProcessor3D (int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, int discreteNormalX_, int discreteNormalY_, int discreteNormalZ_)
int	extent () const override
	Extent of application area (0 for purely local operations)
int	extent (int whichDirection) const override
	Extent of application area along a direction (0 or 1)
void	process (BlockLattice< T, DESCRIPTOR > &blockLattice) override
	Execute post-processing step.
void	processSubDomain (BlockLattice< T, DESCRIPTOR > &blockLattice, int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override
	Execute post-processing step on a sublattice.
Public Member Functions inherited from olb::PostProcessor3D< T, DESCRIPTOR >
	PostProcessor3D ()
virtual	~PostProcessor3D ()
std::string &	getName ()
	read and write access to name
std::string const &	getName () const
	read only access to name
int	getPriority () const
	read only access to priority

Additional Inherited Members
Protected Attributes inherited from olb::PostProcessor3D< T, DESCRIPTOR >
int	_priority

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >

This class interpolates missing f_i from values near the boundary to get a more stable outflow condition for the density.

It is assumed that the next two cells in direction of this f_i have viable values.

Definition at line 43 of file advectionDiffusionBoundaryPostProcessor3D.h.

Constructor & Destructor Documentation

ConvectionBoundaryProcessor3D()

template<typename T , typename DESCRIPTOR >

olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >::ConvectionBoundaryProcessor3D	(	int	x0_,
		int	x1_,
		int	y0_,
		int	y1_,
		int	z0_,
		int	z1_,
		int	discreteNormalX_,
		int	discreteNormalY_,
		int	discreteNormalZ_ )

Definition at line 35 of file advectionDiffusionBoundaryPostProcessor3D.hh.

  : x0(x0_), x1(x1_), y0(y0_), y1(y1_), z0(z0_), z1(z1_)
{
  OLB_PRECONDITION(x0==x1 || y0==y1 || z0==z1);
 
  this->getName() = "ConvectionBoundaryProcessor3D";
 
  interpolationPop[0] = 0;
  for (int iPop = 1; iPop < DESCRIPTOR::q; iPop++) {
    interpolationPop[iPop] = 0;
    // find incoming iPop from material 0
    if (descriptors::c<DESCRIPTOR>(iPop,0)*discreteNormalX + descriptors::c<DESCRIPTOR>(iPop,1)*discreteNormalY + descriptors::c<DESCRIPTOR>(iPop,2)*discreteNormalZ > 0) {
      // check for material number of neighbours has to be one level higher
      interpolationPop[iPop] = 1;
    }
  }
}

References olb::PostProcessor3D< T, DESCRIPTOR >::getName(), and OLB_PRECONDITION.

Here is the call graph for this function:

Member Function Documentation

extent() [1/2]

template<typename T , typename DESCRIPTOR >

int olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >::extent ( ) const

inlineoverridevirtual

Extent of application area (0 for purely local operations)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 48 of file advectionDiffusionBoundaryPostProcessor3D.h.

  {
    return 0;
  }

extent() [2/2]

template<typename T , typename DESCRIPTOR >

int olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >::extent ( int direction ) const

inlineoverridevirtual

Extent of application area along a direction (0 or 1)

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 52 of file advectionDiffusionBoundaryPostProcessor3D.h.

  {
    return 0;
  }

process()

template<typename T , typename DESCRIPTOR >

void olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >::process ( BlockLattice< T, DESCRIPTOR > & blockLattice )

overridevirtual

Execute post-processing step.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 92 of file advectionDiffusionBoundaryPostProcessor3D.hh.

{
  processSubDomain(blockLattice, x0, x1, y0, y1, z0, z1);
}

processSubDomain()

template<typename T , typename DESCRIPTOR >

void olb::ConvectionBoundaryProcessor3D< T, DESCRIPTOR >::processSubDomain ( BlockLattice< T, DESCRIPTOR > & blockLattice, int x0_, int x1_, int y0_, int y1_, int z0_, int z1_ )

overridevirtual

Execute post-processing step on a sublattice.

Implements olb::PostProcessor3D< T, DESCRIPTOR >.

Definition at line 56 of file advectionDiffusionBoundaryPostProcessor3D.hh.

{
  int newX0, newX1, newY0, newY1, newZ0, newZ1;
  if ( util::intersect (
         x0, x1, y0, y1, z0, z1,
         x0_, x1_, y0_, y1_, z0_, z1_,
         newX0, newX1, newY0, newY1, newZ0, newZ1 ) ) {
 
    auto& pop = blockLattice.template getField<descriptors::POPULATION>();
 
#ifdef PARALLEL_MODE_OMP
    #pragma omp parallel for
#endif
    for (int iPop = 1; iPop < DESCRIPTOR::q; ++iPop) {
      if (interpolationPop[iPop] != 0) {
        for (int iX=newX0; iX<=newX1; ++iX) {
          for (int iY=newY0; iY<=newY1; ++iY) {
            for (int iZ=newZ0; iZ<=newZ1; ++iZ) {
              const auto c = descriptors::c<DESCRIPTOR>(iPop);
              if (blockLattice.isInside(iX+c[0],iY+c[0],iZ+c[0]) && blockLattice.isInside(iX+2*c[0],iY+2*c[0],iZ+2*c[0])) {
                //do reflection
                T v = 0.5 * (  pop[iPop][blockLattice.getCellId(iX+  c[0],iY+  c[1],iZ+  c[2])]
                             + pop[iPop][blockLattice.getCellId(iX+2*c[0],iY+2*c[1],iZ+2*c[2])]);
                pop[iPop][blockLattice.getCellId(iX,iY,iZ)] = v;
              }
            }
          }
        }
      }
    }
  }
}

References olb::BlockStructureD< D >::getCellId(), olb::util::intersect(), and olb::BlockStructureD< D >::isInside().

Here is the call graph for this function:

---

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

• src/boundary/postprocessor/advectionDiffusionBoundaryPostProcessor3D.h
• src/boundary/postprocessor/advectionDiffusionBoundaryPostProcessor3D.hh