From 2f0f99bda79a7b99ddc9ddfa35992fe32d6dfa7d Mon Sep 17 00:00:00 2001
From: Samuel Kemmler <samuel.kemmler@fau.de>
Date: Tue, 20 Sep 2022 18:16:28 +0200
Subject: [PATCH] Fix velocityAveragingFactor
---
apps/showcases/Antidunes/Antidunes.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/apps/showcases/Antidunes/Antidunes.cpp b/apps/showcases/Antidunes/Antidunes.cpp
index ce99f3d5e..54a19f556 100644
--- a/apps/showcases/Antidunes/Antidunes.cpp
+++ b/apps/showcases/Antidunes/Antidunes.cpp
@@ -1123,7 +1123,7 @@ int main(int argc, char** argv)
// add sweep for evaluating the fluid's mean velocity
const std::shared_ptr< Vector3< real_t > > meanVelocity = std::make_shared< Vector3< real_t > >(real_c(0));
- const real_t velocityAveragingFactor = 1_r / (liquidHeight * real_c(domainSize[0]) * real_c(domainSize[1]));
+ const real_t velocityAveragingFactor = 1_r / (h * real_c(domainSize[0]) * real_c(domainSize[1]));
MeanVelocityComputer< AntidunesBoundaryHandling_T, PdfField_T, FlagField_T > meanVelocityComputer(
blockForest, antidunesBoundaryHandling, pdfFieldID, meanVelocity, velocityAveragingFactor);
//timeloop.addFuncAfterTimeStep(meanVelocityComputer, "Evaluator: mean velocity");
--
GitLab