olb::normGradPsi Struct Reference

#include <phaseFieldCoupling.h>

Collaboration diagram for olb::normGradPsi:

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

Definition at line 73 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

using olb::normGradPsi::parameters = meta::list<descriptors::EPSILON>

Definition at line 75 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

template<typename CELL , typename PARAMETERS >

void olb::normGradPsi::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 82 of file phaseFieldCoupling.h.

  {
    using V = typename CELL::value_t;
    using DESCRIPTOR = typename CELL::descriptor_t;
 
    auto psi = cell.template getField<descriptors::PSI>();
    auto psi0 = cell.template getField<descriptors::PSI0>();
    auto epsilon = parameters.template get<descriptors::EPSILON>();
    Vector<V,DESCRIPTOR::d> gradPsi{};
    V normGradPsi = 1.;
    if (util::fabs(psi) <= epsilon) {
      V a = psi - cell.neighbor({-1,0}).template getField<descriptors::PSI>();
      V b = cell.neighbor({1,0}).template getField<descriptors::PSI>() - psi;
      V c = psi - cell.neighbor({0,-1}).template getField<descriptors::PSI>();
      V d = cell.neighbor({0,1}).template getField<descriptors::PSI>() - psi;
      V 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()

int olb::normGradPsi::getPriority ( ) const

inline

Definition at line 77 of file phaseFieldCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

scope

OperatorScope olb::normGradPsi::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 74 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h