diff options
Fix non-equilibrium scaling factor
Diffstat (limited to 'src/refinement')
| -rw-r--r-- | src/refinement/coupler2D.hh | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/refinement/coupler2D.hh b/src/refinement/coupler2D.hh index 2218c28..91a7db8 100644 --- a/src/refinement/coupler2D.hh +++ b/src/refinement/coupler2D.hh @@ -33,7 +33,8 @@ namespace olb { template <typename T, template<typename> class DESCRIPTOR> T Coupler2D<T,DESCRIPTOR>::getScalingFactor() const { - return 4. - _coarse.getConverter().getLatticeRelaxationFrequency(); + const T coarseTau = _coarse.getConverter().getLatticeRelaxationTime(); + return (coarseTau - 0.25) / coarseTau; } template <typename T, template<typename> class DESCRIPTOR> |