From 5b6017f8979a0b78150a8ea4465467ab3f415aac Mon Sep 17 00:00:00 2001
From: Christoph Rettinger <christoph.rettinger@fau.de>
Date: Fri, 12 Apr 2019 15:58:34 +0200
Subject: [PATCH] Fixed and updated logging in force on sphere benchmark
---
.../ForcesOnSphereNearPlaneInShearFlow.cpp | 17 ++++++++++++-----
1 file changed, 12 insertions(+), 5 deletions(-)
diff --git a/apps/benchmarks/ForcesOnSphereNearPlaneInShearFlow/ForcesOnSphereNearPlaneInShearFlow.cpp b/apps/benchmarks/ForcesOnSphereNearPlaneInShearFlow/ForcesOnSphereNearPlaneInShearFlow.cpp
index 52e9bf90a..ea62c617f 100644
--- a/apps/benchmarks/ForcesOnSphereNearPlaneInShearFlow/ForcesOnSphereNearPlaneInShearFlow.cpp
+++ b/apps/benchmarks/ForcesOnSphereNearPlaneInShearFlow/ForcesOnSphereNearPlaneInShearFlow.cpp
@@ -385,7 +385,10 @@ void initializeCouetteProfile( const shared_ptr< StructuredBlockStorage > & bloc
}
void logFinalResult( const std::string & fileName, real_t Re, real_t wallDistance, real_t diameter, uint_t numLevels,
- real_t dragCoeff, real_t liftCoeff, real_t timestep )
+ uint_t domainLength, uint_t domainWidth, uint_t domainHeight,
+ real_t dragCoeff, real_t dragCoeffRef,
+ real_t liftCoeff, real_t liftCoeffRef,
+ uint_t timestep, real_t nonDimTimestep )
{
WALBERLA_ROOT_SECTION()
{
@@ -393,12 +396,16 @@ void logFinalResult( const std::string & fileName, real_t Re, real_t wallDistanc
file.open( fileName.c_str(), std::ofstream::app );
file << Re << "\t " << wallDistance << "\t " << diameter << "\t " << numLevels << "\t "
- << dragCoeff << "\t " << liftCoeff << "\t " << timestep
+ << domainLength << "\t " << domainWidth << "\t " << domainHeight << "\t "
+ << dragCoeff << "\t " << dragCoeffRef << "\t "<< std::abs(dragCoeff-dragCoeffRef) / dragCoeffRef << "\t "
+ << liftCoeff << "\t " << liftCoeffRef << "\t "<< std::abs(liftCoeff-liftCoeffRef) / liftCoeffRef << "\t "
+ << timestep << "\t " << nonDimTimestep
<< "\n";
file.close();
}
}
+
//////////
// MAIN //
//////////
@@ -420,7 +427,7 @@ void logFinalResult( const std::string & fileName, real_t Re, real_t wallDistanc
* - Domain size [x, y, z] = [48 x 16 x 8 ] * diameter = [L(ength), W(idth), H(eight)]
* - horizontally periodic, bounded by two planes in z-direction
* - top plane is moving with constant wall velocity -> shear flow
- * - sphere is placed in the vicinity of the bottom plane at [ L/2 + xOffset, W72 + yOffset, dist * diameter]
+ * - sphere is placed in the vicinity of the bottom plane at [ L/2 + xOffset, W/2 + yOffset, dist * diameter]
* - distance of sphere center to the bottom plane is crucial parameter
* - viscosity is adjusted to match specified Reynolds number = shearRate * diameter * wallDistance / viscosity
* - dimensionless drag and lift forces are evaluated and written to logging file
@@ -799,8 +806,8 @@ int main( int argc, char **argv )
if ( !zeroShearTest )
{
std::string resultFileName( baseFolderLogging + "/ResultForcesNearPlane.txt");
- logFinalResult(resultFileName, ReynoldsNumberShear, normalizedWallDistance, diameter, numberOfLevels,
- logger->getLiftCoefficient(), logger->getLiftCoefficient(), real_c(timestep * lbmTimeStepsPerTimeLoopIteration) / physicalTimeScale );
+ logFinalResult(resultFileName, ReynoldsNumberShear, normalizedWallDistance, diameter, numberOfLevels, domainSize[0], domainSize[1], domainSize[2],
+ logger->getDragCoefficient(), dragCorrelationZeng, logger->getLiftCoefficient(), liftCorrelationZeng, timestep, real_c(timestep) / physicalTimeScale );
}
return EXIT_SUCCESS;
--
GitLab