diff --git a/apps/showcases/FreeSurface/BubblyPoiseuille.prm b/apps/showcases/FreeSurface/BubblyPoiseuille.prm
index c6e7d413c5ef6b864f59f6c521542e6ff825a0e4..cfd43e2791bce37964e859f31a5848471e219c32 100644
--- a/apps/showcases/FreeSurface/BubblyPoiseuille.prm
+++ b/apps/showcases/FreeSurface/BubblyPoiseuille.prm
@@ -85,20 +85,22 @@ VTK
normal;
obstacle_normal;
mapped_flag;
- }
+ }
- inclusion_filters
- {
- // only include liquid and interface cells in VTK output
- //liquidInterfaceFilter;
- }
+ inclusion_filters
+ {
+ // only include liquid and interface cells in VTK output
+ //liquidInterfaceFilter;
+ }
- before_functions
- {
- //ghost_layer_synchronization; // only needed if writing the ghost layer
- }
+ before_functions
+ {
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
+ }
}
+
domain_decomposition
{
writeFrequency 0;
diff --git a/apps/showcases/FreeSurface/CapillaryWave.prm b/apps/showcases/FreeSurface/CapillaryWave.prm
index 8d172cadd60e0b0ed616fb4c724c814a2ca23dde..b73208639ee28a8d65244a181d3f34fc01232084 100644
--- a/apps/showcases/FreeSurface/CapillaryWave.prm
+++ b/apps/showcases/FreeSurface/CapillaryWave.prm
@@ -95,9 +95,11 @@ VTK
before_functions
{
- //ghost_layer_synchronization; // only needed if writing the ghost layer
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
}
}
+
domain_decomposition
{
writeFrequency 100;
diff --git a/apps/showcases/FreeSurface/DamBreakCylindrical.prm b/apps/showcases/FreeSurface/DamBreakCylindrical.prm
index de81d688e114b52e226f3c0b42b3f599b31bf33b..aee8429824efd58e5dff7434b0c689bca8fe75b6 100644
--- a/apps/showcases/FreeSurface/DamBreakCylindrical.prm
+++ b/apps/showcases/FreeSurface/DamBreakCylindrical.prm
@@ -98,10 +98,12 @@ VTK
before_functions
{
- //ghost_layer_synchronization; // only needed if writing the ghost layer
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
}
}
+
domain_decomposition
{
writeFrequency 248;
diff --git a/apps/showcases/FreeSurface/DamBreakRectangular.prm b/apps/showcases/FreeSurface/DamBreakRectangular.prm
index 34fdca4b018ce9af4e4f80cf285ee545845246d2..8c4bc96e76bf906072d09786afcac5e34ab3c09b 100644
--- a/apps/showcases/FreeSurface/DamBreakRectangular.prm
+++ b/apps/showcases/FreeSurface/DamBreakRectangular.prm
@@ -88,20 +88,21 @@ VTK
normal;
obstacle_normal;
mapped_flag;
- }
+ }
- inclusion_filters
- {
- // only include liquid and interface cells in VTK output
- //liquidInterfaceFilter;
- }
-
- before_functions
- {
- //ghost_layer_synchronization; // only needed if writing the ghost layer
- }
+ inclusion_filters
+ {
+ // only include liquid and interface cells in VTK output
+ //liquidInterfaceFilter;
+ }
+ before_functions
+ {
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
+ }
}
+
domain_decomposition
{
writeFrequency 991;
diff --git a/apps/showcases/FreeSurface/DropImpact.prm b/apps/showcases/FreeSurface/DropImpact.prm
index 07afacce7480bf9a0307180e4f5da17c1deebbbe..76580d98f79c39163d38b89a83f9045b5d63c095 100644
--- a/apps/showcases/FreeSurface/DropImpact.prm
+++ b/apps/showcases/FreeSurface/DropImpact.prm
@@ -87,20 +87,21 @@ VTK
normal;
obstacle_normal;
mapped_flag;
- }
+ }
- inclusion_filters
- {
- // only include liquid and interface cells in VTK output
- //liquidInterfaceFilter;
- }
-
- before_functions
- {
- //ghost_layer_synchronization; // only needed if writing the ghost layer
- }
+ inclusion_filters
+ {
+ // only include liquid and interface cells in VTK output
+ //liquidInterfaceFilter;
+ }
+ before_functions
+ {
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
+ }
}
+
domain_decomposition
{
writeFrequency 372;
diff --git a/apps/showcases/FreeSurface/DropWetting.prm b/apps/showcases/FreeSurface/DropWetting.prm
index bdfbe4e1ee632724634ad98b8ae2dfd5a8b1cd09..e10ae92738229e615ab22f8a5bcaf97a03787101 100644
--- a/apps/showcases/FreeSurface/DropWetting.prm
+++ b/apps/showcases/FreeSurface/DropWetting.prm
@@ -85,20 +85,21 @@ VTK
normal;
obstacle_normal;
mapped_flag;
- }
+ }
- inclusion_filters
- {
- // only include liquid and interface cells in VTK output
- //liquidInterfaceFilter;
- }
-
- before_functions
- {
- //ghost_layer_synchronization; // only needed if writing the ghost layer
- }
+ inclusion_filters
+ {
+ // only include liquid and interface cells in VTK output
+ //liquidInterfaceFilter;
+ }
+ before_functions
+ {
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
+ }
}
+
domain_decomposition
{
writeFrequency 100;
diff --git a/apps/showcases/FreeSurface/GravityWave.prm b/apps/showcases/FreeSurface/GravityWave.prm
index 363a59c84abbc665dde520ebe4c1c1c89136874e..0400d62da919924847d344523f17fd32f3323d92 100644
--- a/apps/showcases/FreeSurface/GravityWave.prm
+++ b/apps/showcases/FreeSurface/GravityWave.prm
@@ -95,9 +95,11 @@ VTK
before_functions
{
- //ghost_layer_synchronization; // only needed if writing the ghost layer
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
}
}
+
domain_decomposition
{
writeFrequency 0;
diff --git a/apps/showcases/FreeSurface/MovingDrop.prm b/apps/showcases/FreeSurface/MovingDrop.prm
index 184963db1f00e71bb9a3557987d5415a076de5b4..dff95569f3a87410a6a22a46467ad920d67c169e 100644
--- a/apps/showcases/FreeSurface/MovingDrop.prm
+++ b/apps/showcases/FreeSurface/MovingDrop.prm
@@ -85,20 +85,21 @@ VTK
normal;
obstacle_normal;
mapped_flag;
- }
+ }
- inclusion_filters
- {
- // only include liquid and interface cells in VTK output
- //liquidInterfaceFilter;
- }
-
- before_functions
- {
- //ghost_layer_synchronization; // only needed if writing the ghost layer
- }
+ inclusion_filters
+ {
+ // only include liquid and interface cells in VTK output
+ //liquidInterfaceFilter;
+ }
+ before_functions
+ {
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
+ }
}
+
domain_decomposition
{
writeFrequency 0;
diff --git a/apps/showcases/FreeSurface/RisingBubble.prm b/apps/showcases/FreeSurface/RisingBubble.prm
index 5d483426beb659bfe93db4ee94fe595e29c67255..62f1053d208f5b444173b3146932f1425ee6d07d 100644
--- a/apps/showcases/FreeSurface/RisingBubble.prm
+++ b/apps/showcases/FreeSurface/RisingBubble.prm
@@ -96,9 +96,11 @@ VTK
before_functions
{
- //ghost_layer_synchronization; // only needed if writing the ghost layer
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
}
}
+
domain_decomposition
{
writeFrequency 445;
diff --git a/apps/showcases/FreeSurface/TaylorBubble.prm b/apps/showcases/FreeSurface/TaylorBubble.prm
index ee3d8c9f56179798028ebe4e074988a36515549d..e413afc7f856b1cfb44d6c15395dbab171281b97 100644
--- a/apps/showcases/FreeSurface/TaylorBubble.prm
+++ b/apps/showcases/FreeSurface/TaylorBubble.prm
@@ -98,7 +98,8 @@ VTK
before_functions
{
- //ghost_layer_synchronization; // only needed if writing the ghost layer
+ //ghost_layer_synchronization; // only needed if writing the ghost layer
+ gas_cell_zero_setter; // sets velocity=0 and density=1 all gas cells before writing VTK output
}
}
domain_decomposition
diff --git a/src/lbm/free_surface/VtkWriter.h b/src/lbm/free_surface/VtkWriter.h
index f1ff2b93c6e3fdde748c1cf2b579837f37542cc1..f98549a8ff1b6110d3e128af4424eeaf13e0d0ab 100644
--- a/src/lbm/free_surface/VtkWriter.h
+++ b/src/lbm/free_surface/VtkWriter.h
@@ -120,6 +120,60 @@ void addVTKOutput(const std::weak_ptr< StructuredBlockForest >& blockForestPtr,
std::make_shared< field::communication::PackInfo< VectorField_T > >(obstacleNormalFieldID));
beforeFuncs["ghost_layer_synchronization"] = preVTKComm;
+
+ // set velocity and density to zero in obstacle and gas cells (only for visualization purposes); the PDF values in
+ // these cells are not important and thus not set during the simulation;
+ // only enable this functionality if the non-liquid and non-interface cells are not excluded anyway
+ const auto vtkConfigBlock = config->getOneBlock("VTK");
+ const auto fluidFieldConfigBlock = vtkConfigBlock.getBlock("fluid_field");
+ if (fluidFieldConfigBlock)
+ {
+ auto inclusionFiltersConfigBlock = fluidFieldConfigBlock.getBlock("inclusion_filters");
+
+ // liquidInterfaceFilter limits VTK-output to only liquid and interface cells
+ if (!inclusionFiltersConfigBlock.isDefined("liquidInterfaceFilter"))
+ {
+ class ZeroSetter
+ {
+ public:
+ ZeroSetter(const weak_ptr< StructuredBlockForest >& blockForest, const BlockDataID& pdfFieldID,
+ const ConstBlockDataID& flagFieldID,
+ const typename FreeSurfaceBoundaryHandling_T::FlagInfo_T& flagInfo)
+ : blockForest_(blockForest), pdfFieldID_(pdfFieldID), flagFieldID_(flagFieldID), flagInfo_(flagInfo)
+ {}
+
+ void operator()()
+ {
+ auto blockForest = blockForest_.lock();
+ WALBERLA_CHECK_NOT_NULLPTR(blockForest);
+
+ for (auto blockIt = blockForest->begin(); blockIt != blockForest->end(); ++blockIt)
+ {
+ PdfField_T* const pdfField = blockIt->template getData< PdfField_T >(pdfFieldID_);
+ const FlagField_T* const flagField = blockIt->template getData< const FlagField_T >(flagFieldID_);
+ WALBERLA_FOR_ALL_CELLS_INCLUDING_GHOST_LAYER_XYZ(pdfField, uint_c(1), {
+ const typename PdfField_T::Ptr pdfFieldPtr(*pdfField, x, y, z);
+ const typename FlagField_T::ConstPtr flagFieldPtr(*flagField, x, y, z);
+
+ if (flagInfo_.isGas(*flagFieldPtr) || flagInfo_.isObstacle(*flagFieldPtr))
+ {
+ pdfField->setDensityAndVelocity(pdfFieldPtr.cell(), Vector3< real_t >(real_c(0)),
+ real_c(1.0));
+ }
+ }) // WALBERLA_FOR_ALL_CELLS_INCLUDING_GHOST_LAYER_XYZ
+ }
+ }
+
+ private:
+ weak_ptr< StructuredBlockForest > blockForest_;
+ BlockDataID pdfFieldID_;
+ ConstBlockDataID flagFieldID_;
+ typename FreeSurfaceBoundaryHandling_T::FlagInfo_T flagInfo_;
+ };
+
+ beforeFuncs["gas_cell_zero_setter"] = ZeroSetter(blockForest, pdfFieldID, flagFieldID, flagInfo);
+ }
+ }
};
// add VTK output to timeloop
@@ -130,42 +184,6 @@ void addVTKOutput(const std::weak_ptr< StructuredBlockForest >& blockForestPtr,
timeloop.addFuncBeforeTimeStep(output->second.outputFunction, std::string("VTK: ") + output->first,
output->second.requiredGlobalStates, output->second.incompatibleGlobalStates);
}
-
- // only enable the zerosetter (see below) if the non-liquid and non-interface cells are not excluded anyway
- bool enableZeroSetter = true;
- const auto vtkConfigBlock = config->getOneBlock("VTK");
- const auto fluidFieldConfigBlock = vtkConfigBlock.getBlock("fluid_field");
- if (fluidFieldConfigBlock)
- {
- auto inclusionFiltersConfigBlock = fluidFieldConfigBlock.getBlock("inclusion_filters");
-
- if (inclusionFiltersConfigBlock.isDefined("liquidInterfaceFilter"))
- {
- // liquidInterfaceFilter is defined which limits VTK-output to only liquid and interface cells
- enableZeroSetter = false;
- }
- }
-
- // set velocity and density to zero in obstacle and gas cells (only for visualization purposes); the PDF values in
- // these cells are not important and thus not set during the simulation
- if (enableZeroSetter)
- {
- const auto function = [&](IBlock* block) {
- using namespace free_surface;
- PdfField_T* const pdfField = block->getData< PdfField_T >(pdfFieldID);
- const FlagField_T* const flagField = block->getData< const FlagField_T >(flagFieldID);
- WALBERLA_FOR_ALL_CELLS_INCLUDING_GHOST_LAYER_XYZ(pdfField, uint_c(1), {
- const typename PdfField_T::Ptr pdfFieldPtr(*pdfField, x, y, z);
- const typename FlagField_T::ConstPtr flagFieldPtr(*flagField, x, y, z);
-
- if (flagInfo.isGas(*flagFieldPtr) || flagInfo.isObstacle(*flagFieldPtr))
- {
- pdfField->setDensityAndVelocity(pdfFieldPtr.cell(), Vector3< real_t >(0), real_c(1.0));
- }
- }) // WALBERLA_FOR_ALL_CELLS_INCLUDING_GHOST_LAYER_XYZ
- };
- timeloop.add() << Sweep(function, "VTK: zero-setting");
- }
}
} // namespace free_surface