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

#include <phaseFieldCoupling.h>

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

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

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::IsoPhaseFieldCurvedWallProcessor2D< T, DESCRIPTOR, xNormal, yNormal >

Definition at line 816 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

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

using olb::IsoPhaseFieldCurvedWallProcessor2D< T, DESCRIPTOR, xNormal, yNormal >::parameters = meta::list<descriptors::INTERFACE_WIDTH>

Definition at line 818 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

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

template<typename CELL , typename PARAMETERS >

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

inline

Definition at line 825 of file phaseFieldCoupling.h.

  {
    auto phi_s = cell.template getField<descriptors::STATISTIC>();
    auto phi_b = cell.neighbor({-xNormal,-yNormal}).template getFieldComponent<descriptors::STATISTIC>(0);
    if (phi_b >= 0.995) {
      phi_s[0] = 1.;
    }
    else if (phi_b <= 0.005) {
      phi_s[0] = 0.;
    }
    else {
      auto theta = cell.template getField<descriptors::THETA>();
      auto w = parameters.template get<descriptors::INTERFACE_WIDTH>();
      T z_z0 = w/2.*util::atanh(2.*phi_b-1.);
      T cx = util::cos(-theta)*xNormal-util::sin(-theta)*yNormal;
      T cy = util::sin(-theta)*xNormal+util::cos(-theta)*yNormal;
      Vector<T,DESCRIPTOR::d> c{cx,cy};
      Vector<T,DESCRIPTOR::d> n{T(xNormal),T(yNormal)};
      c = util::normalize(c);
      T d = c[0]*n[0]+c[1]*n[1];
      phi_s[0] = 0.5*(1+util::tanh(2.*(z_z0+d)/w));
    }
    cell.template setField<descriptors::STATISTIC>(phi_s);
  }

References olb::util::atanh(), olb::util::cos(), olb::util::normalize(), olb::util::sin(), and olb::util::tanh().

Here is the call graph for this function:

getPriority()

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

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

inline

Definition at line 820 of file phaseFieldCoupling.h.

                          {
    return 1;
  }

Member Data Documentation

scope

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

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

staticconstexpr

Definition at line 817 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h