olb::FullSlipBoundaryPostProcessor2D< typename, typename, NX, NY > Struct Template Reference

#include <setSlipBoundary2D.hh>

Collaboration diagram for olb::FullSlipBoundaryPostProcessor2D< typename, typename, NX, NY >:

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>
struct olb::FullSlipBoundaryPostProcessor2D< typename, typename, NX, NY >

Definition at line 35 of file setSlipBoundary2D.hh.

Member Function Documentation

apply()

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

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

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

inline

Definition at line 44 of file setSlipBoundary2D.hh.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    int reflectionPop[DESCRIPTOR::q];
    int mirrorDirection0;
    int mirrorDirection1;
    int mult = 2 / (NX*NX + NY*NY);
    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;
      if ( scalarProduct < 0) {
        // bounce back for the case discreteNormalX = discreteNormalY = 1, that is mult=1
        if (mult == 1) {
          mirrorDirection0 = -descriptors::c<DESCRIPTOR>(iPop,0);
          mirrorDirection1 = -descriptors::c<DESCRIPTOR>(iPop,1);
        }
        else {
          mirrorDirection0 = descriptors::c<DESCRIPTOR>(iPop,0) - mult*scalarProduct*NX;
          mirrorDirection1 = descriptors::c<DESCRIPTOR>(iPop,1) - mult*scalarProduct*NY;
        }
 
        // 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) {
            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 olb::FullSlipBoundaryPostProcessor2D< typename, typename, NX, NY >::getPriority ( ) const

inline

Definition at line 38 of file setSlipBoundary2D.hh.

                          {
    return -1;
  }

Member Data Documentation

scope

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

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

staticconstexpr

Definition at line 36 of file setSlipBoundary2D.hh.

---

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

• src/boundary/setSlipBoundary2D.hh