olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ > Struct Template Reference

#include <setNewSlipBoundary3D.hh>

Collaboration diagram for olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ >:

Public Member Functions
int	getPriority () const
template<typename CELL >
void	apply (CELL &cell) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCell

Detailed Description

template<typename T, typename DESCRIPTOR, int discreteNormalX, int discreteNormalY, int discreteNormalZ>
struct olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ >

Definition at line 32 of file setNewSlipBoundary3D.hh.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , int discreteNormalX, int discreteNormalY, int discreteNormalZ>

template<typename CELL >

void olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ >::apply ( CELL & cell )

inline

Definition at line 41 of file setNewSlipBoundary3D.hh.

                                     {
    int mirrorDirection0;
    int mirrorDirection1;
    int mirrorDirection2;
    int mult = 2 / (discreteNormalX*discreteNormalX + discreteNormalY*discreteNormalY + discreteNormalZ*discreteNormalZ);
    int reflectionPop[DESCRIPTOR::q] = {0};
    for (int iPop = 1; iPop < DESCRIPTOR::q; iPop++) {
      reflectionPop[iPop] = 0;
      // iPop are the directions which pointing into the fluid, discreteNormal is pointing outwarts
      int scalarProduct = descriptors::c<DESCRIPTOR>(iPop,0)*discreteNormalX + descriptors::c<DESCRIPTOR>(iPop,1)*discreteNormalY + descriptors::c<DESCRIPTOR>(iPop,2)*discreteNormalZ;
      if (scalarProduct < 0) {
        // bounce back for the case discreteNormalX = discreteNormalY = discreteNormalZ = 1, that is mult=0
        if (mult == 0) {
          mirrorDirection0 = -descriptors::c<DESCRIPTOR>(iPop,0);
          mirrorDirection1 = -descriptors::c<DESCRIPTOR>(iPop,1);
          mirrorDirection2 = -descriptors::c<DESCRIPTOR>(iPop,2);
        }
        else {
          mirrorDirection0 = descriptors::c<DESCRIPTOR>(iPop,0) - mult*scalarProduct*discreteNormalX;
          mirrorDirection1 = descriptors::c<DESCRIPTOR>(iPop,1) - mult*scalarProduct*discreteNormalY;
          mirrorDirection2 = descriptors::c<DESCRIPTOR>(iPop,2) - mult*scalarProduct*discreteNormalZ;
        }
        // run through all lattice directions and look for match of direction
        for (int i = 1; i < DESCRIPTOR::q; i++) {
          if (descriptors::c<DESCRIPTOR>(i,0)==mirrorDirection0
              && descriptors::c<DESCRIPTOR>(i,1)==mirrorDirection1
              && descriptors::c<DESCRIPTOR>(i,2)==mirrorDirection2) {
            reflectionPop[iPop] = i;
            break;
          }
        }
      }
    }
    for (int iPop = 1; iPop < DESCRIPTOR::q; iPop++) {
      if (reflectionPop[iPop]!=0) {
        cell[iPop] = cell[reflectionPop[iPop]];
      }
    }
  }

getPriority()

template<typename T , typename DESCRIPTOR , int discreteNormalX, int discreteNormalY, int discreteNormalZ>

int olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ >::getPriority ( ) const

inline

Definition at line 36 of file setNewSlipBoundary3D.hh.

                          {
    return 0;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , int discreteNormalX, int discreteNormalY, int discreteNormalZ>

constexpr OperatorScope olb::SlipBoundaryPostProcessor3D< T, DESCRIPTOR, discreteNormalX, discreteNormalY, discreteNormalZ >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 34 of file setNewSlipBoundary3D.hh.

---

The documentation for this struct was generated from the following file:

• src/boundary/setNewSlipBoundary3D.hh