FSI + kOmegaSST turbulence model

perpendicular-flap/fluid-openfoam/0/k

@@ -0,0 +1,51 @@
+FoamFile
+{
+    version 2.0;
+    format ascii;
+    class volScalarField;
+    object k;
+}
+
+dimensions [ 0 2 -2 0 0 0 0 ];
+
+internalField uniform 5;
+
+boundaryField
+{
+
+    inlet
+    {
+        type fixedValue;
+        value $internalField;
+    }
+
+    outlet
+    {
+        type inletOutlet;
+        inletValue $internalField;
+        value $internalField;
+    }
+
+    flap
+    {
+        type kqRWallFunction;
+        value $internalField;
+    }
+
+    frontAndBack
+    {
+        type empty;
+    }
+
+    upperWall
+    {
+        type zeroGradient;
+    }
+
+    lowerWall
+    {
+        type zeroGradient;
+    }
+
+}
+

perpendicular-flap/fluid-openfoam/0/nut

@@ -0,0 +1,52 @@
+FoamFile
+{
+    version 2.0;
+    format ascii;
+    class volScalarField;
+    object k;
+}
+
+dimensions [ 0 2 -1 0 0 0 0 ];
+
+internalField uniform 5;
+
+boundaryField
+{
+
+    inlet
+    {
+        type calculated;
+        value uniform 0;
+    }
+
+    outlet
+    {
+        type calculated;
+        value uniform 0;
+    }
+
+    flap
+    {
+        type nutkWallFunction;
+        value uniform 0;
+    }
+
+    frontAndBack
+    {
+        type empty;
+    }
+
+    upperWall
+    {
+        type calculated;
+        value uniform 0;
+    }
+
+    lowerWall
+    {
+        type calculated;
+        value uniform 0;
+    }
+
+}
+

perpendicular-flap/fluid-openfoam/0/omega

@@ -0,0 +1,51 @@
+FoamFile
+{
+    version 2.0;
+    format ascii;
+    class volScalarField;
+    object k;
+}
+
+dimensions [ 0 0 -1 0 0 0 0 ];
+
+internalField uniform 10;
+
+boundaryField
+{
+
+    inlet
+    {
+        type fixedValue;
+        value $internalField;
+    }
+
+    outlet
+    {
+        type inletOutlet;
+        inletValue $internalField;
+        value $internalField;
+    }
+
+    flap
+    {
+        type omegaWallFunction;
+        value $internalField;
+    }
+
+    frontAndBack
+    {
+        type empty;
+    }
+
+    upperWall
+    {
+        type zeroGradient;
+    }
+
+    lowerWall
+    {
+        type zeroGradient;
+    }
+
+}
+

perpendicular-flap/fluid-openfoam/constant/transportProperties

@@ -8,4 +8,4 @@ FoamFile
 
 transportModel  Newtonian;
 
-nu              nu [ 0 2 -1 0 0 0 0 ] 1;
+nu              nu [ 0 2 -1 0 0 0 0 ] 1e-5;

perpendicular-flap/fluid-openfoam/constant/turbulenceProperties

@@ -6,4 +6,11 @@ FoamFile
     object      turbulenceProperties;
 }
 
-simulationType laminar;
+simulationType RAS;
+
+RAS
+{
+    RASModel kOmegaSST;
+    turbulence on;
+    printCoeffs on;
+}

perpendicular-flap/fluid-openfoam/system/fvSchemes

@@ -20,6 +20,9 @@ divSchemes
 {
     default         none;
     div(phi,U)      bounded Gauss upwind;
+    div(phi,k)      bounded Gauss limitedLinear 1;
+    div(phi,omega)  bounded Gauss upwind;
+
     div((nuEff*dev2(T(grad(U))))) Gauss linear;
 }
 
@@ -37,3 +40,8 @@ snGradSchemes
 {
     default         corrected;
 }
+
+wallDist
+{
+    method          Poisson;
+}

perpendicular-flap/fluid-openfoam/system/fvSolution

@@ -54,6 +54,36 @@ solvers
         $U;
         relTol          0;
     }
+
+    "(U|epsilon|k|omega)"
+    {
+        solver smoothSolver;
+        smoother symGaussSeidel;
+        tolerance 1e-6;
+        relTol 1e-3;
+        // minIter 100;
+    }
+    "(U|k|epsilon|omega)Final"
+    {
+        $U;
+        tolerance 1e-6;
+        relTol 0;
+        minIter 2;
+    }
+
+    yPsi
+    {
+        solver PBiCGStab;
+        preconditioner DILU;
+        tolerance 1e-5;
+        relTol 0;
+    }
+
+    yPsiFinal
+    {
+        $yPsi
+        relTol 0;
+    }
 }
 
 PIMPLE