olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z > Class Template Reference

PostProcessor for the density / velocity outflow boundaries in the free energy model. More...

#include <boundaryPostProcessors3D.h>

Collaboration diagram for olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >:

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 NORMAL_X, int NORMAL_Y, int NORMAL_Z>
class olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >

PostProcessor for the density / velocity outflow boundaries in the free energy model.

The density / order parameters are prescribed to the outflow nodes such that they obey the local-velocity convective boundary condition given in Lou, Gou, Shi (2013).

Definition at line 261 of file boundaryPostProcessors3D.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

template<CONCEPT(Cell) CELL>

void olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::apply

(

CELL &

cell

)

Definition at line 553 of file boundaryPostProcessors3D.hh.

                                                                                               {
 
  T rho, rho0, rho1, u[3];
 
  rho0 = cell.computeRho();
 
  cell.neighbor({-NORMAL_X,-NORMAL_Y,-NORMAL_Z}).computeRhoU(rho1, u);
 
  T uPerp = 0;
 
  Vector<T, 3> normalVec({NORMAL_X,NORMAL_Y,NORMAL_Z});
 
  if (normalVec[2] == 0) {
    if (normalVec[1] == 0) {
      if (normalVec[0] < 0) {
        uPerp = -u[0];
      } else {
        uPerp = u[0];
      }
    } else if (normalVec[0] == 0) {
      if (normalVec[1] < 0) {
        uPerp = -u[1];
      } else {
        uPerp = u[1];
      }
    } else {
      uPerp = util::sqrt(u[0] * u[0] + u[1] * u[1]);
    }
  } else if (normalVec[1] == 0) {
    if (normalVec[0] == 0) {
      if (normalVec[2] < 0) {
        uPerp = -u[2];
      } else {
        uPerp = u[2];
      }
    } else {
      uPerp = util::sqrt(u[0] * u[0] + u[2] * u[2]);
    }
  } else if (normalVec[0] == 0) {
    uPerp = util::sqrt(u[1] * u[1] + u[2] * u[2]);
  } else {
    uPerp = util::sqrt(u[0] * u[0] + u[1] * u[1] + u[2] * u[2]);
  }
 
  rho = (rho0 + uPerp * rho1) / (1. + uPerp);
  cell.defineRho(rho);
}

References olb::util::sqrt().

Here is the call graph for this function:

getPriority()

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

int olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::getPriority ( ) const

inline

Definition at line 265 of file boundaryPostProcessors3D.h.

                          {
    return 0;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

constexpr OperatorScope olb::FreeEnergyConvectiveProcessor3D< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 263 of file boundaryPostProcessors3D.h.

---

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

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