olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ > Struct Template Reference

#include <setSlipBoundary3D.hh>

Collaboration diagram for olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ >:

Public Member Functions
int	getPriority () const
template<typename CELL , typename V = typename CELL::value_t>
void	apply (CELL &x_b) any_platform

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

Detailed Description

template<typename, typename, int NX, int NY, int NZ>
struct olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ >

Definition at line 36 of file setSlipBoundary3D.hh.

Member Function Documentation

apply()

template<typename , typename , int NX, int NY, int NZ>

template<typename CELL , typename V = typename CELL::value_t>

void olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ >::apply ( CELL & x_b )

inline

Definition at line 45 of file setSlipBoundary3D.hh.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    int reflectionPop[DESCRIPTOR::q];
    int mirrorDirection0;
    int mirrorDirection1;
    int mirrorDirection2;
    int mult = 2 / (NX*NX + NY*NY + NZ*NZ);
    reflectionPop[0] =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)*NX + descriptors::c<DESCRIPTOR>(iPop,1)*NY + descriptors::c<DESCRIPTOR>(iPop,2)*NZ;
      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*NX;
          mirrorDirection1 = descriptors::c<DESCRIPTOR>(iPop,1) - mult*scalarProduct*NY;
          mirrorDirection2 = descriptors::c<DESCRIPTOR>(iPop,2) - mult*scalarProduct*NZ;
        }
 
        // 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) {
        //do reflection
        x_b[iPop] = x_b[reflectionPop[iPop]];
      }
    }
  }

getPriority()

template<typename , typename , int NX, int NY, int NZ>

int olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ >::getPriority ( ) const

inline

Definition at line 39 of file setSlipBoundary3D.hh.

                          {
    return -1;
  }

Member Data Documentation

scope

template<typename , typename , int NX, int NY, int NZ>

constexpr OperatorScope olb::FullSlipBoundaryPostProcessor3D< typename, typename, NX, NY, NZ >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 37 of file setSlipBoundary3D.hh.

---

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

• src/boundary/setSlipBoundary3D.hh