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

#include <setBouzidiBoundary.h>

Collaboration diagram for olb::BouzidiPostProcessor:

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 zero-velocity Bouzidi boundary.

Definition at line 67 of file setBouzidiBoundary.h.

Member Function Documentation

◆ apply()

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

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

inline

Definition at line 76 of file setBouzidiBoundary.h.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    const auto q = x_b.template getFieldPointer<descriptors::BOUZIDI_DISTANCE>();
    for (int iPop = 1; iPop < descriptors::q<DESCRIPTOR>(); ++iPop) {
      // update missing population if valid bouzidi distance
      if (q[iPop] > V{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
 
        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])
                           + (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]);
      }
      // 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);
        x_b[iPop_opposite] = x_b[iPop];
      }
    }
  }

◆ getPriority()

int olb::BouzidiPostProcessor::getPriority ( ) const

inline

Definition at line 71 of file setBouzidiBoundary.h.

                          {
    return -1;
  }

Member Data Documentation

◆ scope

constexpr OperatorScope olb::BouzidiPostProcessor::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 69 of file setBouzidiBoundary.h.

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

src/boundary/setBouzidiBoundary.h