olb::initialPsi Struct Reference

#include <phaseFieldCoupling.h>

Collaboration diagram for olb::initialPsi:

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

Definition at line 41 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

using olb::initialPsi::parameters = meta::list<descriptors::INTERFACE_WIDTH>

Definition at line 43 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

template<typename CELL , typename PARAMETERS >

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

inline

Definition at line 50 of file phaseFieldCoupling.h.

  {
    using V = typename CELL::value_t;
    auto w = parameters.template get<descriptors::INTERFACE_WIDTH>();
    V phi = cell.template getFieldComponent<descriptors::STATISTIC>(0);
    V psi = 0;
    if (phi > (0.99)) {
      psi += -4.595120;
    }
    else if (phi < (0.01)) {
      psi += 4.595120;
    }
    else {
      psi += util::log(1/phi-1);
    }
    psi *= -w/4.;
    cell.template setField<descriptors::SCALAR>(psi/util::sqrt(psi*psi+1.));
    cell.template setField<descriptors::PSI>(psi);
    cell.template setField<descriptors::PSI0>(psi);
  }

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

Here is the call graph for this function:

getPriority()

int olb::initialPsi::getPriority ( ) const

inline

Definition at line 45 of file phaseFieldCoupling.h.

                          {
    return 1;
  }

Member Data Documentation

scope

OperatorScope olb::initialPsi::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 42 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h