From 0c5f36f40a92863d17b5c18ad3b999e79121f5cb Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 18 Jun 2018 15:00:47 +0100
Subject: [PATCH] ENH: Move stabilisation to steady-state part of LES motorBike
tutorial. See #856.
---
.../motorBike/motorBike/system/controlDict | 68 ++-----------------
.../LES/motorBike/motorBike/system/fvSchemes | 2 +-
.../motorBike/system/stabilizationSchemes | 0
.../simpleFoam/motorBike/system/controlDict | 1 -
.../simpleFoam/motorBike/system/fvSchemes | 2 +-
5 files changed, 9 insertions(+), 64 deletions(-)
rename tutorials/incompressible/{simpleFoam => pisoFoam/LES/motorBike}/motorBike/system/stabilizationSchemes (100%)
diff --git a/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/controlDict b/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/controlDict
index 452037f221d..e20cbd30fb7 100644
--- a/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/controlDict
+++ b/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/controlDict
@@ -49,68 +49,14 @@ runTimeModifiable true;
functions
{
- streamLines
- {
- type streamLine;
- writeControl writeTime;
- setFormat vtk;
- trackForward true;
- fields (p U);
- lifeTime 10000;
- nSubCycle 5;
- cloud particleTracks;
- seedSampleSet
- {
- type uniform;
- axis x;
- start (-1.001 1e-07 0.0011);
- end (-1.001 1e-07 1.0011);
- nPoints 20;
- }
- }
- cuttingPlane
- {
- type surfaces;
- libs ("libsampling.so");
- writeControl writeTime;
- surfaceFormat vtk;
- fields (p U);
- interpolationScheme cellPoint;
- surfaces
- (
- yNormal
- {
- type cuttingPlane;
- planeType pointAndNormal;
- pointAndNormalDict
- {
- point (0 0 0);
- normal (0 1 0);
- }
- interpolate true;
- }
- );
- }
- forces
- {
- type forceCoeffs;
- libs ("libforces.so");
- writeControl timeStep;
- writeInterval 1;
- patches ("motorBike.*");
- rho rhoInf;
- log true;
- rhoInf 1;
- liftDir (0 0 1);
- dragDir (1 0 0);
- CofR (0.72 0 0);
- pitchAxis (0 1 0);
- magUInf 20;
- lRef 1.42;
- Aref 0.75;
- }
-
+ // Demonstrate runtime postprocessing
+ #include "streamLines"
+ #include "cuttingPlane"
#include "runtimePostProcessing"
+
+ // Do some other functionObjects
+ #include "forceCoeffs"
+ #include "stabilizationSchemes"
}
diff --git a/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/fvSchemes b/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/fvSchemes
index 88088155947..ecb07d042b2 100644
--- a/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/fvSchemes
+++ b/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/fvSchemes
@@ -30,7 +30,7 @@ gradSchemes
divSchemes
{
default none;
- div(phi,U) bounded Gauss linearUpwindV grad(U);
+ div(phi,U) bounded Gauss localBlended linearUpwindV grad(U) linear;
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
diff --git a/tutorials/incompressible/simpleFoam/motorBike/system/stabilizationSchemes b/tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/stabilizationSchemes
similarity index 100%
rename from tutorials/incompressible/simpleFoam/motorBike/system/stabilizationSchemes
rename to tutorials/incompressible/pisoFoam/LES/motorBike/motorBike/system/stabilizationSchemes
diff --git a/tutorials/incompressible/simpleFoam/motorBike/system/controlDict b/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
index 30048b9ae62..92992749fa9 100644
--- a/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
+++ b/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
@@ -51,7 +51,6 @@ functions
#include "cuttingPlane"
#include "forceCoeffs"
#include "ensightWrite"
- #include "stabilizationSchemes"
}
diff --git a/tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes b/tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes
index 88d86cc6309..12fdd994df6 100644
--- a/tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes
+++ b/tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes
@@ -28,7 +28,7 @@ gradSchemes
divSchemes
{
default none;
- div(phi,U) bounded Gauss localBlended linearUpwindV grad(U) linear;
+ div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
--
GitLab