olb::ChemicalPotentialCoupling3D Struct Reference

#include <freeEnergyCoupling3D.h>

Collaboration diagram for olb::ChemicalPotentialCoupling3D:

Classes
struct	ALPHA
struct	KAPPA1
struct	KAPPA2
struct	KAPPA3

Public Types
using	parameters = meta::list<ALPHA, KAPPA1, KAPPA2, KAPPA3>

Public Member Functions
template<typename CELLS , typename PARAMETERS >
void	apply (CELLS &cells, PARAMETERS &parameters) any_platform

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

Detailed Description

Definition at line 29 of file freeEnergyCoupling3D.h.

Member Typedef Documentation

parameters

using olb::ChemicalPotentialCoupling3D::parameters = meta::list<ALPHA, KAPPA1, KAPPA2, KAPPA3>

Definition at line 39 of file freeEnergyCoupling3D.h.

Member Function Documentation

apply()

template<typename CELLS , typename PARAMETERS >

void olb::ChemicalPotentialCoupling3D::apply ( CELLS & cells, PARAMETERS & parameters )

inline

Definition at line 42 of file freeEnergyCoupling3D.h.

                                                                {
 
    using V = typename CELLS::template value_t<names::A>::value_t;
    using DESCRIPTOR = typename CELLS::template value_t<names::A>::descriptor_t;
 
    // Get the cell of the first lattice
    auto& cellA = cells.template get<names::A>();
 
    // Get the cell of the second lattice
    auto& cellB = cells.template get<names::B>();
 
    V rhoA = cellA.computeRho();
    V rhoB = cellB.computeRho();
 
    V densitySum = rhoA + rhoB;
    V densityDifference = rhoA - rhoB;
 
    V kappa1 = parameters.template get<KAPPA1>();
    V kappa2 = parameters.template get<KAPPA2>();
 
    V term1 = 0.125 * kappa1 * (densitySum)
                  * (densitySum-1.) * (densitySum-2.);
 
    V term2 = 0.125 * kappa2 * (densityDifference)
                  * (densityDifference-1.) * (densityDifference-2.);
 
 
    V term3 = 0.;
 
 
    V laplaceRho1 = -24.0 * rhoA;
    V laplaceRho2 = -24.0 * rhoB;
    V laplaceRho3 = 0;
 
    for(int iPop=1; iPop<DESCRIPTOR::q; iPop++) {
 
      auto nextCellA = cellA.neighbor(descriptors::c<DESCRIPTOR>(iPop));
      auto nextCellB = cellB.neighbor(descriptors::c<DESCRIPTOR>(iPop));
 
 
      V nextRhoA = nextCellA.computeRho();
      V nextRhoB = nextCellB.computeRho();
 
      if(iPop == 1 || iPop == 2 || iPop == 3 || iPop == 10 || iPop == 11 || iPop == 12) {
        laplaceRho1 += (2 * nextRhoA);
        laplaceRho2 += (2 * nextRhoB);
      } else {
        laplaceRho1 += nextRhoA;
        laplaceRho2 += nextRhoB;
      }
    }
 
    laplaceRho1 *= 1.0 / 6.0;
    laplaceRho2 *= 1.0 / 6.0;
 
    auto alpha = parameters.template get<ALPHA>();
 
    cellA.template setField<descriptors::CHEM_POTENTIAL>(term1 + term2
            + 0.25*alpha*alpha*( (kappa2 - kappa1) * laplaceRho2
                                 +(kappa2 + kappa1) * (laplaceRho3 - laplaceRho1) ));
 
    cellB.template setField<descriptors::CHEM_POTENTIAL>(term1 - term2
            + 0.25*alpha*alpha*( (kappa2 - kappa1) * (laplaceRho1 - laplaceRho3)
                                 -(kappa2 + kappa1) * laplaceRho2 ));
 
  }

Member Data Documentation

scope

constexpr OperatorScope olb::ChemicalPotentialCoupling3D::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 32 of file freeEnergyCoupling3D.h.

---

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

• src/dynamics/freeEnergyCoupling3D.h