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

#include <phaseFieldCoupling.h>

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

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

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 xNormal, int yNormal>
struct olb::normGradPsiBoundary2D< T, DESCRIPTOR, xNormal, yNormal >

Definition at line 129 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

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

using olb::normGradPsiBoundary2D< T, DESCRIPTOR, xNormal, yNormal >::parameters = meta::list<descriptors::EPSILON>

Definition at line 131 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

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

template<typename CELL , typename PARAMETERS >

void olb::normGradPsiBoundary2D< T, DESCRIPTOR, xNormal, yNormal >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 138 of file phaseFieldCoupling.h.

  {
    auto psi = cell.template getField<descriptors::PSI>();
    auto psi0 = cell.template getField<descriptors::PSI0>();
    auto epsilon = parameters.template get<descriptors::EPSILON>();
    Vector<T,DESCRIPTOR::d> gradPsi{};
    T normGradPsi = 1.;
    if (util::fabs(psi) <= epsilon) {
      T a = psi - cell.neighbor({-1,0}).template getField<descriptors::PSI>();
      T b = cell.neighbor({1,0}).template getField<descriptors::PSI>() - psi;
      T c = psi - cell.neighbor({0,-1}).template getField<descriptors::PSI>();
      T d = cell.neighbor({0,1}).template getField<descriptors::PSI>() - psi;
      if (xNormal == 1) b = a;
      else if (xNormal == -1) a = b;
      if (yNormal == 1) d = c;
      else if (yNormal == -1) c = d;
      T gradX, gradY;
      if (psi0 > 0) {
        if (a<0) a = 0;
        if (b>0) b = 0;
        if (c<0) c = 0;
        if (d>0) d = 0;
        gradX = a;
        if (b*b>a*a) gradX = b;
        gradY = c;
        if (d*d>c*c) gradY = d;
        normGradPsi = util::sqrt(gradX*gradX+gradY*gradY);
      }
      else if (psi0 < 0) {
        if (a>0) a = 0;
        if (b<0) b = 0;
        if (c>0) c = 0;
        if (d<0) d = 0;
        gradX = a;
        if (b*b>a*a) gradX = b;
        gradY = c;
        if (d*d>c*c) gradY = d;
        normGradPsi = util::sqrt(gradX*gradX+gradY*gradY);
      }
      else {}
      cell.template setField<descriptors::NORMGRADPSI>(normGradPsi);
    } else {
      cell.template setField<descriptors::NORMGRADPSI>(normGradPsi);
    }
  }

References olb::util::fabs(), and olb::util::sqrt().

Here is the call graph for this function:

getPriority()

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

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

inline

Definition at line 133 of file phaseFieldCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

scope

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

OperatorScope olb::normGradPsiBoundary2D< T, DESCRIPTOR, xNormal, yNormal >::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 130 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h