olb::BouzidiAdeDirichletPostProcessor Class Reference

#include <setBouzidiBoundary.h>

Collaboration diagram for olb::BouzidiAdeDirichletPostProcessor:

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

Definition at line 101 of file setBouzidiBoundary.h.

Member Function Documentation

apply()

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

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

inline

Definition at line 110 of file setBouzidiBoundary.h.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    const auto q = x_b.template getFieldPointer<descriptors::BOUZIDI_DISTANCE>();
    const auto phi_d = x_b.template getFieldPointer<descriptors::BOUZIDI_ADE_DIRICHLET>();
    V f = V{1/3};       // D2Q5
    if (descriptors::q<DESCRIPTOR>() == 7) {
      f = V{0.25};      // D3Q7
    }
    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
        auto source_d = phi_d[iPop] * f;                                        // source term set by the dirichlet condition
        auto t_i = descriptors::t<V,DESCRIPTOR>(iPop);
        auto t_iopp = descriptors::t<V,DESCRIPTOR>(iPop_opposite);
 
        x_b[iPop_opposite] = (q[iPop] <= V{0.5}) // cut is closer to the fluid cell
                           * (V{-2} * q[iPop] * (x_s[iPop] + t_i) + (V{2} * q[iPop] - V{1}) * (x_b[iPop] + t_i) + source_d)
                           + (q[iPop] >  V{0.5}) // cut is closer to the solid cell
                           * (V{-1} / (V{2} * q[iPop]) * (x_s[iPop] + t_i) + (V{1} - V{1} / (V{2} * q[iPop])) * (x_f[iPop_opposite] + t_iopp) + (V{1} / (V{2} * q[iPop])) * source_d)
                           - t_iopp;
      }
      // 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 source_d = phi_d[iPop] * f;
        auto t_i = descriptors::t<V,DESCRIPTOR>(iPop);
        auto t_iopp = descriptors::t<V,DESCRIPTOR>(iPop_opposite);
        x_b[iPop_opposite] = -(x_b[iPop] + t_i) + source_d - t_iopp;
      }
    }
  }

getPriority()

int olb::BouzidiAdeDirichletPostProcessor::getPriority ( ) const

inline

Definition at line 105 of file setBouzidiBoundary.h.

                          {
    return -1;
  }

Member Data Documentation

scope

constexpr OperatorScope olb::BouzidiAdeDirichletPostProcessor::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 103 of file setBouzidiBoundary.h.

---

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

• src/boundary/setBouzidiBoundary.h