olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 > Struct Template Reference

#include <robinBoundaryLatticePostProcessor3D.h>

Collaboration diagram for olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >:

Public Types
using	parameters = meta::list<descriptors::OMEGA>

Public Member Functions
int	getPriority () const
template<typename CELL , typename PARAMETERS >
void	apply (CELL &cell, PARAMETERS &parameters) any_platform

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

Detailed Description

template<typename T, typename DESCRIPTOR, int Normal1, int Normal2, int Normal3>
struct olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >

Definition at line 245 of file robinBoundaryLatticePostProcessor3D.h.

Member Typedef Documentation

parameters

template<typename T , typename DESCRIPTOR , int Normal1, int Normal2, int Normal3>

using olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >::parameters = meta::list<descriptors::OMEGA>

Definition at line 248 of file robinBoundaryLatticePostProcessor3D.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , int Normal1, int Normal2, int Normal3>

template<typename CELL , typename PARAMETERS >

void olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >::apply ( CELL & cell, PARAMETERS & parameters )

inline

define variables, names match the characters in the formula

calculation of the two sums used in the MengCurvedCorr scheme

define variables, names match the characters in the formula

necessary procedure to calculate B

final calculation of the new populations

Definition at line 254 of file robinBoundaryLatticePostProcessor3D.h.

                                                             {
    bool schemeSwitch = false; //switch between two schemes; true=MengCurvedCorr, false=Ju2020
 
    // get external values
    T omega = parameters.template get<descriptors::OMEGA>();
    auto v = cell.template getField<descriptors::VELOCITY>(); //velocity field
    auto a = cell.template getField<descriptors::G>(); //field containing coefficients for the robin boundary condition
    T a1 = a[0]; T a2 = a[1]; T a3 = a[2];
    std::vector<int> n(3);
    n[0] = -Normal1; n[1] = -Normal2; n[2] = -Normal3; //change normal to point into domain
 
    constexpr auto unknownIndices = util::subIndexOutgoing3DonCorners<DESCRIPTOR,Normal1,Normal2,Normal3>(); //function returns the indices of missing populations
    assert(unknownIndices.size() == 3);
 
    if(schemeSwitch){
      //--- MengCurvedCorr ---//
      T w = descriptors::t<T,DESCRIPTOR>(1); //lattice weight = 1/8
      T gamma = -descriptors::invCs2<T,DESCRIPTOR>() * omega;
      T b1 = a1-a2*gamma*(n[0]*v[0]+n[1]*v[1]+n[2]*v[2]);
      T b2 = a2*gamma;
 
      T sum1 = cell.computeRho();
      T sum2 = 0;
      for (unsigned iPop : unknownIndices) {
        sum1 -= (cell[iPop]+w + cell[descriptors::opposite<DESCRIPTOR>(iPop)]+w);
        sum2 += cell[descriptors::opposite<DESCRIPTOR>(iPop)] + w;
      }
      T beta = (a3-b1*sum1+2*b2*sum2)/(unknownIndices.size()*w*(b1+b2));
 
      for (unsigned iPop : unknownIndices) {
        cell[iPop] = w * beta - (cell[descriptors::opposite<DESCRIPTOR>(iPop)] + w) - w; //set new populations
      }
    } else{
      //--- Ju2020 ---//
      T w = descriptors::t<T,DESCRIPTOR>(1); //lattice weight = 1/8
      T chi = (a2!=0.) ? -1/ (descriptors::invCs2<T,DESCRIPTOR>()*a2*omega) : (1/omega)/(1/omega-0.5);
      T k = a1; //reaction rate
      T NdotV = n[0]*v[0] + n[1]*v[1] + n[2]*v[2];
      T A = unknownIndices.size()-1;
      T alpha = (-chi*k+NdotV+2*w*(1-A)) / (chi*k-NdotV+2*w*(1+A)); //bounced-back part
 
      T B = 0;
      for (unsigned iPop : unknownIndices) {
        B += 2*(cell[descriptors::opposite<DESCRIPTOR>(iPop)]+w);
      }
 
      for (unsigned iPop : unknownIndices) {
        B -= 2*(cell[descriptors::opposite<DESCRIPTOR>(iPop)]+w); //final value of B
        T beta = (2*w*chi*a3+2*w*B) / (chi*k-NdotV+2*w*(1+A)); //reverse reaction part
        cell[iPop] = alpha * (cell[descriptors::opposite<DESCRIPTOR>(iPop)] + w) + beta - w; //set new populations
        B += 2*(cell[descriptors::opposite<DESCRIPTOR>(iPop)]+w); //restore B for next iteration
      }
    }
  }

getPriority()

template<typename T , typename DESCRIPTOR , int Normal1, int Normal2, int Normal3>

int olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >::getPriority ( ) const

inline

Definition at line 249 of file robinBoundaryLatticePostProcessor3D.h.

                          {
    return 1;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , int Normal1, int Normal2, int Normal3>

constexpr OperatorScope olb::robinBoundaryExtendedPostProcessor3DCorners< T, DESCRIPTOR, Normal1, Normal2, Normal3 >::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 247 of file robinBoundaryLatticePostProcessor3D.h.

---

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

• src/boundary/postprocessor/robinBoundaryLatticePostProcessor3D.h