#include <phaseFieldCoupling.h>

Collaboration diagram for olb::AllenCahnOutletCoupling:

Classes
struct	NORMAL

struct	SIGMA

struct	W

Public Types
using	parameters = meta::list<SIGMA,W,NORMAL>

Public Member Functions
int	getPriority () const

template<typename CELL , typename PARAMETERS >
void	apply (CELL &cells, PARAMETERS &parameters) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCellWithParameters

Detailed Description

Definition at line 382 of file phaseFieldCoupling.h.

Member Typedef Documentation

◆ parameters

using olb::AllenCahnOutletCoupling::parameters = meta::list<SIGMA,W,NORMAL>

Definition at line 387 of file phaseFieldCoupling.h.

Member Function Documentation

◆ apply()

template<typename CELL , typename PARAMETERS >

void olb::AllenCahnOutletCoupling::apply	(	CELL &	cells,
		PARAMETERS &	parameters )

inline

Definition at line 394 of file phaseFieldCoupling.h.

  {
    using V = typename CELL::template value_t<names::NavierStokes>::value_t;
    using DESCRIPTOR = typename CELL::template value_t<names::NavierStokes>::descriptor_t;
 
    auto& cellNS = cells.template get<names::NavierStokes>();
    auto& cellAC = cells.template get<names::Component1>();
 
    auto w = parameters.template get<W>();
    auto sigma = parameters.template get<SIGMA>();
    Vector<int, DESCRIPTOR::d> normal = parameters.template get<NORMAL>();
    auto rho = cellNS.template getField<descriptors::RHO>();
    auto nablaRho = cellNS.template getField<descriptors::NABLARHO>();
 
    V phi = cellAC.template getFieldComponent<descriptors::STATISTIC>(0);
    auto cellAC_n1 = cellAC.neighbor({-normal[0],-normal[1]});
    auto cellNS_n1 = cellNS.neighbor({-normal[0],-normal[1]});
    V phi_n1 = cellAC_n1.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 = cellAC.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>();
    }
    V laplacePhi_n1 = 0;
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      const V phi_i_n1 = cellAC_n1.neighbor(descriptors::c<DESCRIPTOR>(iPop)).template getFieldComponent<descriptors::STATISTIC>(0);
      laplacePhi_n1 += 2*(phi_i_n1 - phi_n1) * descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>();
    }
 
    // Computation of thermodynamic pressure and update populations
    V k = 1.5*sigma*w;
    V beta = 12*sigma/w;
    V mu = 4*beta*phi*(phi-1.)*(phi-0.5)-k*laplacePhi;
    V mu_n1 = 4*beta*phi_n1*(phi_n1-1.)*(phi_n1-0.5)-k*laplacePhi_n1;
 
    V u[DESCRIPTOR::d] {};
    //cellAC_n1.computeU(u);
    V rhoU = u[0]*nablaRho[0]+u[1]*nablaRho[1];
    V p_n1 = cellNS_n1.computeRho();
    V p_n2 = cellNS.neighbor({-2*normal[0],-2*normal[1]}).computeRho();
    V p_old = cellNS.computeRho();
    V p_out = (p_n1 + (mu+mu_n1)/2*(phi-phi_n1));//0.8*(0.25*sigma*curv*(2*phi-1)*((2*phi-1)*(2*phi-1)-3)+0.5*sigma*curv);
    //V p_star = p_out;
    if (phi >= 0.95) {
      p_out = 0;
    }
    //V p_star = p_out-3./10.*rhoU;
    cellNS.defineRho(p_out);
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      cellNS[iPop] += descriptors::t<V,DESCRIPTOR>(iPop) * descriptors::invCs2<V,DESCRIPTOR>() * (p_out-p_old);
    }
  }

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

Here is the call graph for this function:

◆ getPriority()

int olb::AllenCahnOutletCoupling::getPriority ( ) const

inline

Definition at line 389 of file phaseFieldCoupling.h.

                          {
    return 0;
  }

Member Data Documentation

◆ scope

OperatorScope olb::AllenCahnOutletCoupling::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 383 of file phaseFieldCoupling.h.

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

src/dynamics/phaseFieldCoupling.h

Classes

Public Types

Public Member Functions

Static Public Attributes