olb::BouzidiVelocityPostProcessor Class Reference

Post processor for the velocity Bouzidi boundary. More...

#include <setBouzidiBoundary.h>

Collaboration diagram for olb::BouzidiVelocityPostProcessor:

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

Post processor for the velocity Bouzidi boundary.

Definition at line 148 of file setBouzidiBoundary.h.

Member Function Documentation

apply()

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

void olb::BouzidiVelocityPostProcessor::apply ( CELL & x_b )

inline

Definition at line 157 of file setBouzidiBoundary.h.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    const auto q = x_b.template getFieldPointer<descriptors::BOUZIDI_DISTANCE>();
    const auto veloCoeff = x_b.template getFieldPointer<descriptors::BOUZIDI_VELOCITY>();
    for (int iPop=1; iPop < descriptors::q<DESCRIPTOR>(); ++iPop) {
      // update missing population if valid bouzidi distance
      if (q[iPop] > 0) {
        const auto c = descriptors::c<DESCRIPTOR>(iPop);
        const int iPop_opposite = descriptors::opposite<DESCRIPTOR>(iPop);
        auto x_s = x_b.neighbor(c);                                             // solid side neighbor
        auto x_f = x_b.neighbor(descriptors::c<DESCRIPTOR>(iPop_opposite));     // fluid side neighbor opposite to the missing population
        auto veloTerm = veloCoeff[iPop] * (descriptors::t<V,DESCRIPTOR>(iPop)) * (descriptors::invCs2<V,DESCRIPTOR>());
 
        x_b[iPop_opposite] = (q[iPop] <= V{0.5}) // cut is closer to the fluid cell
                           * (V{2} * q[iPop] * x_s[iPop] + (V{1} - V{2} * q[iPop]) * x_b[iPop] - V{2} * veloTerm)
                           + (q[iPop] >  V{0.5}) // cut is closer to the solid cell
                           * (V{0.5} / q[iPop] * x_s[iPop] + V{0.5} * (V{2} * q[iPop] - V{1}) / q[iPop] * x_f[iPop_opposite] - V{1}/q[iPop] * veloTerm);
      }
      // if intersection point is on the cell then fall back to full-way bounce back
      else if (q[iPop] == V{0}) {
        const int iPop_opposite = descriptors::opposite<DESCRIPTOR>(iPop);
        auto veloTerm = veloCoeff[iPop] * (descriptors::t<V,DESCRIPTOR>(iPop)) * (descriptors::invCs2<V,DESCRIPTOR>());
        x_b[iPop_opposite] = x_b[iPop] - V{2} * veloTerm;
      }
    }
  }

getPriority()

int olb::BouzidiVelocityPostProcessor::getPriority ( ) const

inline

Definition at line 152 of file setBouzidiBoundary.h.

                          {
    return -1;
  }

Member Data Documentation

scope

constexpr OperatorScope olb::BouzidiVelocityPostProcessor::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 150 of file setBouzidiBoundary.h.

---

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

• src/boundary/setBouzidiBoundary.h