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

#include <robinBoundaryLatticePostProcessor3D.h>

Collaboration diagram for olb::robinBoundaryLatticePostProcessor3Dother< 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::robinBoundaryLatticePostProcessor3Dother< T, DESCRIPTOR, Normal1, Normal2, Normal3 >

Definition at line 113 of file robinBoundaryLatticePostProcessor3D.h.

Member Typedef Documentation

parameters

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

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

Definition at line 116 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::robinBoundaryLatticePostProcessor3Dother< T, DESCRIPTOR, Normal1, Normal2, Normal3 >::apply ( CELL & cell, PARAMETERS & parameters )

inline

define variables, names match the characters in the formula in thesis

necessary procedure to calculate B

final calculation of the new populations

Definition at line 122 of file robinBoundaryLatticePostProcessor3D.h.

                                                             {
 
    T omega = parameters.template get<descriptors::OMEGA>();
    auto v = cell.template getField<descriptors::VELOCITY>();
    auto a = cell.template getField<descriptors::G>(); //get robin boundary coefficients
    T a1 = a[0];
    T a2 = a[1];
    T a3 = a[2];
 
    const int direction = abs(Normal1*0 + Normal2*1 + Normal3*2); // ==0,1,2
    const int orientation = -(Normal1 + Normal2 + Normal3); // ==+-1
    Vector<int,3> n(-Normal1, -Normal2, -Normal3); //normal points into domain
    T NdotV = n[0]*v[0] + n[1]*v[1] + n[2]*v[2];
    T w = descriptors::t<T,DESCRIPTOR>(1); //lattice weight = 1/8
    T gamma = -descriptors::invCs2<T,DESCRIPTOR>() * omega; //a factor used in scheme
 
    Vector<T,1> unknownIndices;
    for(int iPop = 0; iPop<DESCRIPTOR::q; iPop++){
      if(descriptors::c<DESCRIPTOR>(iPop)[direction] == orientation){
        unknownIndices[0] = iPop;
      }
    }
 
    T chi = (a2!=0.) ? 1/ (gamma*a2) : (1/omega)/(1/omega-0.5);
    T k = a1; //reaction rate
    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 contribution
      cell[iPop] = alpha * (cell[descriptors::opposite<DESCRIPTOR>(iPop)] + w) + beta - w; //set new populations
      B += 2*(cell[descriptors::opposite<DESCRIPTOR>(iPop)]+w); //restore B
    }
  }

References olb::abs(), and olb::GenericVector< T, D, IMPL >::size().

Here is the call graph for this function:

getPriority()

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

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

inline

Definition at line 117 of file robinBoundaryLatticePostProcessor3D.h.

                          {
    return 0;
  }

Member Data Documentation

scope

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

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

staticconstexpr

Definition at line 115 of file robinBoundaryLatticePostProcessor3D.h.

---

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

• src/boundary/postprocessor/robinBoundaryLatticePostProcessor3D.h