olb::ChemPotentialPhaseFieldProcessor Struct Reference

#include <phaseFieldCoupling.h>

Collaboration diagram for olb::ChemPotentialPhaseFieldProcessor:

Public Types
using	parameters = meta::list<descriptors::SCALAR,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 666 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

using olb::ChemPotentialPhaseFieldProcessor::parameters = meta::list<descriptors::SCALAR,descriptors::INTERFACE_WIDTH>

Definition at line 668 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

template<typename CELL , typename PARAMETERS >

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

inline

Definition at line 675 of file phaseFieldCoupling.h.

  {
    using V = typename CELL::value_t;
    using DESCRIPTOR = typename CELL::descriptor_t;
 
    V phi = cell.template getFieldComponent<descriptors::STATISTIC>(0);
    V laplacePhi = 0;
    //TODO: use lattice gradient schemes here from Cahn-Hilliard implementation
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      const V phi_i = cell.neighbor(descriptors::c<DESCRIPTOR>(iPop)).template getField<descriptors::PHIWETTING>();
      laplacePhi += 2*(phi_i - phi) * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>();
    }
 
    auto sigma = parameters.template get<descriptors::SCALAR>();
    auto w = parameters.template get<descriptors::INTERFACE_WIDTH>();
 
    // Computation and storage of chemical potential
    V k = 1.5*sigma*w;
    V beta = 12*sigma/w;
    V mu = 4*beta*phi*(phi-1.)*(phi-0.5)+0.25*(phi<0)*phi-k*laplacePhi;
    cell.template setField<descriptors::CHEM_POTENTIAL>(mu);
  }

References olb::descriptors::c(), olb::descriptors::invCs2(), and olb::descriptors::t().

Here is the call graph for this function:

getPriority()

int olb::ChemPotentialPhaseFieldProcessor::getPriority ( ) const

inline

Definition at line 670 of file phaseFieldCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

scope

OperatorScope olb::ChemPotentialPhaseFieldProcessor::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 667 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h