olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z > Class Template Reference

PostProcessors for the chemical potential boundary condition in the free energy model. More...

#include <boundaryPostProcessors3D.h>

Collaboration diagram for olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >:

Public Member Functions
int	getPriority () const
template<CONCEPT(Cell) CELL>
void	apply (CELL &cell) any_platform

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

Detailed Description

template<typename T, typename DESCRIPTOR, int NORMAL_X, int NORMAL_Y, int NORMAL_Z>
class olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >

PostProcessors for the chemical potential boundary condition in the free energy model.

The chemical potentials on the boundary are set equal to the chemical potential on the fluid cell normal to the boundary. This is necessary because the coupling between the lattices requires the calculation of the gradient of the chemical potential.

It would be preferable if these were implemented as a lattice coupling that ran between the chemical potential and force lattice couplings. However there is no access to the discrete normals in lattice couplings.

Definition at line 208 of file boundaryPostProcessors3D.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

template<CONCEPT(Cell) CELL>

void olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::apply

(

CELL &

cell

)

Definition at line 529 of file boundaryPostProcessors3D.hh.

                                                                                                     {
 
  cell.template setField<descriptors::CHEM_POTENTIAL>(
    cell.neighbor({-NORMAL_X,-NORMAL_Y,-NORMAL_Z}).template getField<descriptors::CHEM_POTENTIAL>());
 
  T rho0 = cell.computeRho();
  T rho1 = cell.neighbor({-NORMAL_X,-NORMAL_Y,-NORMAL_Z}).computeRho();
 
  cell.template setField<descriptors::CHEM_POTENTIAL>(
    cell.template getField<descriptors::CHEM_POTENTIAL>() + (rho1 / rho0 - 1) / descriptors::invCs2<T,DESCRIPTOR>());
}

getPriority()

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

int olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::getPriority ( ) const

inline

Definition at line 212 of file boundaryPostProcessors3D.h.

                          {
    return 0;
  }

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , int NORMAL_X, int NORMAL_Y, int NORMAL_Z>

constexpr OperatorScope olb::FreeEnergyChemPotBoundaryProcessor3DA< T, DESCRIPTOR, NORMAL_X, NORMAL_Y, NORMAL_Z >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 210 of file boundaryPostProcessors3D.h.

---

The documentation for this class was generated from the following files:

• src/boundary/boundaryPostProcessors3D.h
• src/boundary/boundaryPostProcessors3D.hh