From 91ea5ddd77bb1479de65baf8145157a81d4690cd Mon Sep 17 00:00:00 2001
From: niikonst <niikonst@blogin7.usr.hlrn.de>
Date: Fri, 10 Dec 2021 13:52:13 +0100
Subject: [PATCH] add acoustic scaling
---
apps/cpu/ViskomatXL/viskomat.cpp | 13 ++++++++++---
1 file changed, 10 insertions(+), 3 deletions(-)
diff --git a/apps/cpu/ViskomatXL/viskomat.cpp b/apps/cpu/ViskomatXL/viskomat.cpp
index 8bbcf6e19..c88840742 100644
--- a/apps/cpu/ViskomatXL/viskomat.cpp
+++ b/apps/cpu/ViskomatXL/viskomat.cpp
@@ -69,9 +69,16 @@ void bflow(string configname)
//double nuLB = OmegaLB * R * 1e3 * R * 1e3 / Re;
- double dx = deltax * 1e-3;
+ // double dx = deltax * 1e-3;
+ // double nuLB = OmegaLB * (R / dx)*(R / dx) / Re;
+
+ double dx = 1.0 * 1e-3;
double nuLB = OmegaLB * (R / dx)*(R / dx) / Re;
+ //acustic scaling
+ OmegaLB /= 2.0;
+ nuLB *= 2.0;
+
SPtr<LBMUnitConverter> conv = SPtr<LBMUnitConverter>(new LBMUnitConverter());
//SPtr<LBMUnitConverter> conv = SPtr<LBMUnitConverter>(new LBMUnitConverter(1, 1461, 970, 1e3));
//UBLOG(logINFO, conv->toString());
@@ -405,9 +412,9 @@ void bflow(string configname)
SPtr<UbScheduler> stepGhostLayer(new UbScheduler(1));
SPtr<Calculator> calculator(new BasicCalculator(grid, stepGhostLayer, endTime));
calculator->addCoProcessor(npr);
- calculator->addCoProcessor(fp);
+ //calculator->addCoProcessor(fp);
calculator->addCoProcessor(fp2);
- calculator->addCoProcessor(writeMQCoProcessor);
+ //calculator->addCoProcessor(writeMQCoProcessor);
//calculator->addCoProcessor(writeThixotropicMQCoProcessor);
calculator->addCoProcessor(restartCoProcessor);
--
GitLab