summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-01-17 12:45:25 +0100
committerAdrian Kummerlaender2019-06-24 15:17:09 +0200
commit5b8bad24546c0c01b11113fcd9b212bc633b6757 (patch)
treef3fcbc316783cfd17730c1057fd6b6c4707c8e17
parent69f40a510a39ba561d290ad07ccbd3f0c3786a7c (diff)
downloadgrid_refinement_openlb-5b8bad24546c0c01b11113fcd9b212bc633b6757.tar
grid_refinement_openlb-5b8bad24546c0c01b11113fcd9b212bc633b6757.tar.gz
grid_refinement_openlb-5b8bad24546c0c01b11113fcd9b212bc633b6757.tar.bz2
grid_refinement_openlb-5b8bad24546c0c01b11113fcd9b212bc633b6757.tar.xz
grid_refinement_openlb-5b8bad24546c0c01b11113fcd9b212bc633b6757.zip
Fix non-equilibrium scaling factor
-rw-r--r--src/refinement/coupler2D.hh3
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>