olb::LiangSinglePostProcessor Struct Reference

#include <phaseFieldCoupling.h>

Collaboration diagram for olb::LiangSinglePostProcessor:

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 451 of file phaseFieldCoupling.h.

Member Typedef Documentation

parameters

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

Definition at line 453 of file phaseFieldCoupling.h.

Member Function Documentation

apply()

template<typename CELL , typename PARAMETERS >

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

inline

Definition at line 460 of file phaseFieldCoupling.h.

  {
    using V = typename CELL::value_t;
    using DESCRIPTOR = typename CELL::descriptor_t;
 
    V phi = cell.template getFieldComponent<descriptors::STATISTIC>(0);
    Vector<V,DESCRIPTOR::d> gradPhi{};
    V laplacePhi = 0;
 
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      const V phi_i = cell.neighbor(descriptors::c<DESCRIPTOR>(iPop)).template getFieldComponent<descriptors::STATISTIC>(0);
      gradPhi += phi_i * descriptors::c<DESCRIPTOR>(iPop) * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>();
      laplacePhi += 2*(phi_i - phi) * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>();
    }
    auto w = parameters.template get<descriptors::INTERFACE_WIDTH>();
 
    // Computation and storage of forces
    Vector<V,DESCRIPTOR::d> forceAC{};
    Vector<V,DESCRIPTOR::d> old_phiU{};
    Vector<V,DESCRIPTOR::d> n{};
    Vector<V,DESCRIPTOR::d> phiU{};
    V u[DESCRIPTOR::d] {};
    cell.computeU(u);
    V gradPhiSqr = 0;
    for (int iD=0; iD < DESCRIPTOR::d; ++iD) {
      gradPhiSqr += gradPhi[iD]*gradPhi[iD];
      phiU[iD] = phi*u[iD];
    }
    if (gradPhiSqr >= 1e-28) {
      n += gradPhi/util::sqrt(gradPhiSqr);
    }
 
    old_phiU = cell.template getField<descriptors::OLD_PHIU>();
    cell.template setField<descriptors::OLD_PHIU>(phiU);
    V lambda = 4*phi*(1.-phi)/w;
    forceAC += (phiU - old_phiU) + lambda*n/descriptors::invCs2<V,DESCRIPTOR>();
    cell.template setField<descriptors::FORCE>(forceAC);
  }

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

Here is the call graph for this function:

getPriority()

int olb::LiangSinglePostProcessor::getPriority ( ) const

inline

Definition at line 455 of file phaseFieldCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

scope

OperatorScope olb::LiangSinglePostProcessor::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 452 of file phaseFieldCoupling.h.

---

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

• src/dynamics/phaseFieldCoupling.h