olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 > Struct Template Reference

#include <robinBoundaryLatticePostProcessor3D.h>

Collaboration diagram for olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 >:

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 Plane, int Normal1, int Normal2>
struct olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 >

Definition at line 169 of file robinBoundaryLatticePostProcessor3D.h.

Member Typedef Documentation

parameters

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

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

Definition at line 172 of file robinBoundaryLatticePostProcessor3D.h.

Member Function Documentation

apply()

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

template<typename CELL , typename PARAMETERS >

void olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 >::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 178 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>();
    auto a = cell.template getField<descriptors::G>();
    T a1 = a[0];
    T a2 = a[1];
    T a3 = a[2];
 
    constexpr auto unknownIndices = util::subIndexOutgoing3DonEdges<DESCRIPTOR,Plane,Normal1,Normal2>();
    assert(unknownIndices.size() == 2);
    auto c0=descriptors::c<DESCRIPTOR>(unknownIndices[0]);
    auto c1=descriptors::c<DESCRIPTOR>(unknownIndices[1]);
    Vector<int,3> n(c0[0]+c1[0], c0[1]+c1[1], c0[2]+c1[2]); //normal points into domain
 
 
    if(schemeSwitch){
      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{
      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 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::GenericVector< T, D, IMPL >::size().

Here is the call graph for this function:

getPriority()

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

int olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 >::getPriority ( ) const

inline

Definition at line 173 of file robinBoundaryLatticePostProcessor3D.h.

                          {
    return 2;
  }

Member Data Documentation

scope

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

constexpr OperatorScope olb::robinBoundaryExtendedPostProcessor3DEdges< T, DESCRIPTOR, Plane, Normal1, Normal2 >::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 171 of file robinBoundaryLatticePostProcessor3D.h.

---

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

• src/boundary/postprocessor/robinBoundaryLatticePostProcessor3D.h