From b4c303907f6fbf8b5fccb3620e65c588e93b17ad Mon Sep 17 00:00:00 2001
From: Konstantin Kutscher <kutscher@irmb.tu-bs.de>
Date: Tue, 23 Jan 2018 17:27:54 +0100
Subject: [PATCH] implemented variable bulk viscosity
---
source/Applications/DLR-F16-Solid/f16.cpp | 4 +-
...ibleCumulant4thOrderViscosityLBMKernel.cpp | 37 +++++++++++++------
...ssibleCumulant4thOrderViscosityLBMKernel.h | 8 ++--
3 files changed, 34 insertions(+), 15 deletions(-)
diff --git a/source/Applications/DLR-F16-Solid/f16.cpp b/source/Applications/DLR-F16-Solid/f16.cpp
index 1670b12df..479658d94 100644
--- a/source/Applications/DLR-F16-Solid/f16.cpp
+++ b/source/Applications/DLR-F16-Solid/f16.cpp
@@ -154,6 +154,7 @@ void run(string configname)
//dynamicPointerCast<CompressibleCumulantLBMKernel>(kernel)->setRelaxationParameter(CompressibleCumulantLBMKernel::NORMAL);
SPtr<LBMKernel> kernel = SPtr<LBMKernel>(new CompressibleCumulant4thOrderViscosityLBMKernel());
+ dynamicPointerCast<CompressibleCumulant4thOrderViscosityLBMKernel>(kernel)->setBulkViscosity(nuLB*1.0e4);
kernel->setNX(std::array<int,3>{{blockNx[0], blockNx[1], blockNx[2]}});
SPtr<BCProcessor> bcProc;
@@ -251,7 +252,8 @@ void run(string configname)
RefineCrossAndInsideGbObjectBlockVisitor refVisitorCylinderL2(refCylinderL2, level);
grid->accept(refVisitorCylinderL2);
- SPtr<GbObject3D> refBoxL2(new GbCuboid3D(0.015, 0.0, -0.04, 1.100, 0.1, 0.04));
+ //SPtr<GbObject3D> refBoxL2(new GbCuboid3D(0.015, 0.0, -0.04, 1.100, 0.1, 0.04));
+ SPtr<GbObject3D> refBoxL2(new GbCuboid3D(0.015, 0.0, -0.04, 0.5, 0.1, 0.04));
if (myid==0) GbSystem3D::writeGeoObject(refBoxL2.get(), pathOut+"/geo/refBoxL2", WbWriterVtkXmlASCII::getInstance());
RefineCrossAndInsideGbObjectBlockVisitor refVisitorBoxL2(refBoxL2, level);
grid->accept(refVisitorBoxL2);
diff --git a/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.cpp b/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.cpp
index 63ab6e24b..9284224cd 100644
--- a/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.cpp
+++ b/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.cpp
@@ -12,6 +12,7 @@
CompressibleCumulant4thOrderViscosityLBMKernel::CompressibleCumulant4thOrderViscosityLBMKernel()
{
this->compressible = true;
+ bulkViscosity = 0;
}
//////////////////////////////////////////////////////////////////////////
CompressibleCumulant4thOrderViscosityLBMKernel::~CompressibleCumulant4thOrderViscosityLBMKernel(void)
@@ -39,14 +40,15 @@ SPtr<LBMKernel> CompressibleCumulant4thOrderViscosityLBMKernel::clone()
kernel->setIndex(ix1, ix2, ix3);
kernel->setDeltaT(deltaT);
- if (bulkOmegaToOmega)
+ if (bulkViscosity != 0)
{
- dynamicPointerCast<CompressibleCumulant4thOrderViscosityLBMKernel>(kernel)->OxxPyyPzz = collFactor;
- }
+ OxxPyyPzz = 1.0/(3.0*bulkViscosity/deltaT+0.5);
+ }
else
{
- dynamicPointerCast<CompressibleCumulant4thOrderViscosityLBMKernel>(kernel)->OxxPyyPzz = one;
+ OxxPyyPzz = one;
}
+ dynamicPointerCast<CompressibleCumulant4thOrderViscosityLBMKernel>(kernel)->OxxPyyPzz = this->OxxPyyPzz;
return kernel;
}
@@ -101,10 +103,18 @@ void CompressibleCumulant4thOrderViscosityLBMKernel::calculate()
int maxX3 = bcArrayMaxX3-ghostLayerWidth;
LBMReal omega = collFactor;
- LBMReal OxyyPxzz = eight*(-two+omega)*(one+two*omega)/(-eight-fourteen*omega+seven*omega*omega);//one;
- LBMReal OxyyMxzz = eight*(-two+omega)*(-seven+four*omega)/(fiftysix-fifty*omega+nine*omega*omega);//one;
- LBMReal Oxyz = twentyfour*(-two+omega)*(-two-seven*omega+three*omega*omega)/(fourtyeight+c152*omega-c130*omega*omega+twentynine*omega*omega*omega);
+ //LBMReal OxyyPxzz = eight*(-two+omega)*(one+two*omega)/(-eight-fourteen*omega+seven*omega*omega);//one;
+ //LBMReal OxyyMxzz = eight*(-two+omega)*(-seven+four*omega)/(fiftysix-fifty*omega+nine*omega*omega);//one;
+ //LBMReal Oxyz = twentyfour*(-two+omega)*(-two-seven*omega+three*omega*omega)/(fourtyeight+c152*omega-c130*omega*omega+twentynine*omega*omega*omega);
+ LBMReal OxyyPxzz = 8.0*(omega-2.0)*(OxxPyyPzz*(3.0*omega-1.0)-5.0*omega)/(8.0*(5.0-2.0*omega)*omega+OxxPyyPzz*(8.0+omega*(9.0*omega-26.0)));
+ LBMReal OxyyMxzz = 8.0*(omega-2.0)*(omega+OxxPyyPzz*(3.0*omega-7.0))/(OxxPyyPzz*(56.0-42.0*omega+9.0*omega*omega)-8.0*omega);
+ LBMReal Oxyz = 24.0*(omega-2.0)*(4.0*omega*omega+omega*OxxPyyPzz*(18.0-13.0*omega)+OxxPyyPzz*OxxPyyPzz*(2.0+omega*(6.0*omega-11.0)))/(16.0*omega*omega*(omega-6.0)-2.0*omega*OxxPyyPzz*(216.0+5.0*omega*(9.0*omega-46.0))+OxxPyyPzz*OxxPyyPzz*(omega*(3.0*omega-10.0)*(15.0*omega-28.0)-48.0));
+
+ //LBMReal A = (four + two*omega - three*omega*omega) / (two - seven*omega + five*omega*omega);
+ //LBMReal B = (four + twentyeight*omega - fourteen*omega*omega) / (six - twentyone*omega + fiveteen*omega*omega);
+ LBMReal A = (4.0*omega*omega+2.0*omega*OxxPyyPzz*(omega-6.0)+OxxPyyPzz*OxxPyyPzz*(omega*(10.0-3.0*omega)-4.0))/((omega-OxxPyyPzz)*(OxxPyyPzz*(2.0+3.0*omega)-8.0*omega));
+ LBMReal B = (4.0*omega*OxxPyyPzz*(9.0*omega-16.0)-4.0*omega*omega-2.0*OxxPyyPzz*OxxPyyPzz*(2.0+9.0*omega*(omega-2.0)))/(3.0*(omega-OxxPyyPzz)*(OxxPyyPzz*(2.0+3.0*omega)-8.0*omega));
for (int x3 = minX3; x3 < maxX3; x3++)
{
@@ -695,8 +705,7 @@ void CompressibleCumulant4thOrderViscosityLBMKernel::calculate()
//CUMcbb += wadjust * (-CUMcbb);
//////////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////
- LBMReal A = (four + two*omega - three*omega*omega) / (two - seven*omega + five*omega*omega);
- LBMReal B = (four + twentyeight*omega - fourteen*omega*omega) / (six - twentyone*omega + fiveteen*omega*omega);
+
//////////////////////////////////////////////////////////////////////////
@@ -1042,11 +1051,17 @@ double CompressibleCumulant4thOrderViscosityLBMKernel::getCalculationTime()
return timer.getTotalTime();
}
//////////////////////////////////////////////////////////////////////////
-void CompressibleCumulant4thOrderViscosityLBMKernel::setBulkOmegaToOmega(bool value)
+void CompressibleCumulant4thOrderViscosityLBMKernel::setBulkViscosity(LBMReal value)
{
- bulkOmegaToOmega = value;
+ bulkViscosity = value;
}
+//////////////////////////////////////////////////////////////////////////
+//void CompressibleCumulant4thOrderViscosityLBMKernel::setBulkOmegaToOmega(bool value)
+//{
+// bulkOmegaToOmega = value;
+//}
+
//void CompressibleCumulant4thOrderViscosityLBMKernel::setViscosityFlag(bool vf)
//{
// viscosityFlag = vf;
diff --git a/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.h b/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.h
index 96bd3288e..0a22436fc 100644
--- a/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.h
+++ b/source/VirtualFluidsCore/LBM/CompressibleCumulant4thOrderViscosityLBMKernel.h
@@ -22,7 +22,8 @@ public:
virtual void calculate();
virtual SPtr<LBMKernel> clone();
double getCalculationTime() override;
- void setBulkOmegaToOmega(bool value);
+ void setBulkViscosity(LBMReal value);
+ //void setBulkOmegaToOmega(bool value);
//void setViscosityFlag(bool vf);
protected:
virtual void initDataSet();
@@ -42,8 +43,9 @@ protected:
LBMReal forcingX3;
// bulk viscosity
- bool bulkOmegaToOmega;
- LBMReal OxxPyyPzz;
+ //bool bulkOmegaToOmega;
+ LBMReal OxxPyyPzz; //omega2 (bulk viscosity)
+ LBMReal bulkViscosity;
//bool viscosityFlag;
//LBMReal OxyyPxzz;
--
GitLab