Rename sMinus to omegaMinus for consistency
In Palabos the relaxation rate is called "omega", in the TRT there is a change of notation omega->s. This leads to a mixed notation that can make confusion. Moreover "guarantees that the wall is exactly half-way between nodes" is true only in the case of walls parallel to the lattice. Proposition:
- rename sMinus and sPlus to omegaMinus and omegaPlus
- add a omegaMinus getter
- add a method to convert omegaMinus to Magic Parameter and vice versa
[palabos] src/complexDynamics/trtDynamics.h (Lines 45-58)
class TRTdynamics : public IsoThermalBulkDynamics<T,Descriptor> {
public:
/* *************** Construction / Destruction ************************ */
// sMinus is the relaxation time associated to the anti-symmetric part.
// If a value of 0 is provided for sMinus, this parameter will default
// to a value depending on omega, which guarantees that the wall is
// exactly half-way between nodes, according to Eq. 11 in Pan et al 2006:
// https://doi.org/10.1016/j.compfluid.2005.03.008
// The original formulation of the model is probably in Ginzburg and
// Adler 1994:
// https://doi.org/10.1051/jp2:1994123.
// It corresponds to a "magic parameter" of 3/16.
TRTdynamics(T omega_, T sMinus_ = T());
TRTdynamics(HierarchicUnserializer& unserializer);