olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal > Struct Template Reference

#include <boundaryPostProcessors3D.h>

Collaboration diagram for olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal >:

Public Member Functions
int	getPriority () const
template<CONCEPT(Cell) CELL>
void	apply (CELL &cell) any_platform

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

Detailed Description

template<typename T, typename DESCRIPTOR, int xNormal, int yNormal, int zNormal>
struct olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal >

Definition at line 111 of file boundaryPostProcessors3D.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , int xNormal, int yNormal, int zNormal>

template<CONCEPT(Cell) CELL>

void olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal >::apply

(

CELL &

cell

)

Definition at line 240 of file boundaryPostProcessors3D.hh.

{
  using namespace olb::util::tensorIndices3D;
 
  auto& dynamics = cell.getDynamics();
 
  T rho100 = cell.neighbor({-1*xNormal, -0*yNormal, -0*zNormal}).computeRho();
  T rho010 = cell.neighbor({-0*xNormal, -1*yNormal, -0*zNormal}).computeRho();
  T rho001 = cell.neighbor({-0*xNormal, -0*yNormal, -1*zNormal}).computeRho();
  T rho200 = cell.neighbor({-2*xNormal, -0*yNormal, -0*zNormal}).computeRho();
  T rho020 = cell.neighbor({-0*xNormal, -2*yNormal, -0*zNormal}).computeRho();
  T rho002 = cell.neighbor({-0*xNormal, -0*yNormal, -2*zNormal}).computeRho();
 
  T rho = (T)4/(T)9 * (rho001 + rho010 + rho100) - (T)1/(T)9 * (rho002 + rho020 + rho200);
 
  T dx_u[DESCRIPTOR::d], dy_u[DESCRIPTOR::d], dz_u[DESCRIPTOR::d];
  fd::DirectedGradients3D<T, DESCRIPTOR, 0, xNormal, 0, true>::interpolateVector(dx_u, cell);
  fd::DirectedGradients3D<T, DESCRIPTOR, 1, yNormal, 0, true>::interpolateVector(dy_u, cell);
  fd::DirectedGradients3D<T, DESCRIPTOR, 2, zNormal, 0, true>::interpolateVector(dz_u, cell);
 
  T dx_ux = dx_u[0];
  T dy_ux = dy_u[0];
  T dz_ux = dz_u[0];
  T dx_uy = dx_u[1];
  T dy_uy = dy_u[1];
  T dz_uy = dz_u[1];
  T dx_uz = dx_u[2];
  T dy_uz = dy_u[2];
  T dz_uz = dz_u[2];
  T omega = dynamics.getOmegaOrFallback(std::numeric_limits<T>::signaling_NaN());
  T sToPi = - rho / descriptors::invCs2<T,DESCRIPTOR>() / omega;
  T pi[util::TensorVal<DESCRIPTOR >::n];
  pi[xx] = (T)2 * dx_ux * sToPi;
  pi[yy] = (T)2 * dy_uy * sToPi;
  pi[zz] = (T)2 * dz_uz * sToPi;
  pi[xy] = (dx_uy + dy_ux) * sToPi;
  pi[xz] = (dx_uz + dz_ux) * sToPi;
  pi[yz] = (dy_uz + dz_uy) * sToPi;
 
  // Computation of the particle distribution functions
  // according to the regularized formula
  T u[DESCRIPTOR::d];
  cell.computeU(u);
 
  for (int iPop = 0; iPop < DESCRIPTOR::q; ++iPop) {
    cell[iPop] = dynamics.computeEquilibrium(iPop,rho,u)
               + equilibrium<DESCRIPTOR>::template fromPiToFneq<T>(iPop, pi);
  }
}

References olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, orthogonal >::interpolateVector().

Here is the call graph for this function:

getPriority()

template<typename T , typename DESCRIPTOR , int xNormal, int yNormal, int zNormal>

int olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal >::getPriority ( ) const

inline

Definition at line 114 of file boundaryPostProcessors3D.h.

                          {
    return 1;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , int xNormal, int yNormal, int zNormal>

constexpr OperatorScope olb::OuterVelocityCornerProcessor3D< T, DESCRIPTOR, xNormal, yNormal, zNormal >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 112 of file boundaryPostProcessors3D.h.

---

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

• src/boundary/boundaryPostProcessors3D.h
• src/boundary/boundaryPostProcessors3D.hh