Problem with continuity simpleFoam kkl omega airfoils
I am trying to compute lift and drag from 3 airfoils, set up in a triangle in 2D using simpleFoam (OpenFOAM 4.1) and the kkl-omega turbulence model. Upon running the simulation, time step continuity errors occur and my omega, kl and kt blow up. I do not know if I am using the most appropriate fvSolution solvers. I would appreciate if someone could have a look. I have also included my boundary conditions. My polymesh passed the checkMesh fyi, I can submit it if required.
fvSolution: FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-08; relTol 1.e-3; smoother GaussSeidel; nPreSweeps 2; nPostSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } "(U|kt|kl|omega)" { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-10; relTol 1.e-3; } } SIMPLE { nNonOrthogonalCorrectors 2; pRefCell 0; pRefValue 0; } relaxationFactors { p 0.1; U 0.1; kl 0.1; kt 0.1; omega 0.1; } // ************************************************** *********************** // Boundary conditions: U: FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (10.0 0 0); boundaryField { exit { type freestream; freestreamValue uniform (10.0 0 0); } inlet { type freestream; freestreamValue uniform (10.0 0 0); } top { type slip; } bottom { type slip; } airfoila { type fixedValue; value uniform (0 0 0); } airfoilb { type fixedValue; value uniform (0 0 0); } airfoilc { type fixedValue; value uniform (0 0 0); } front { type empty; } back { type empty; } } // ************************************************** *********************** // p: FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } exit { type fixedValue; value uniform 0; } top { type slip; } bottom { type slip; } airfoila { type zeroGradient; } airfoilb { type zeroGradient; } airfoilc { type zeroGradient; } front { type empty; } back { type empty; } } // ************************************************** *********************** // omega: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object omega; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 -1 0 0 0 0]; internalField uniform 1.74964; boundaryField { inlet { type freestream; freestreamValue $internalField; } exit { type freestream; freestreamValue $internalField; } top { type slip; } bottom { type slip; } airfoila { type fixedValue; value uniform 1e-11; } airfoilb { type fixedValue; value uniform 1e-11; } airfoilc { type fixedValue; value uniform 1e-11; } front { type empty; } back { type empty; } } // ************************************************** *********************** // kl: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object kl; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type freestream; freestreamValue $internalField; } exit { type freestream; freestreamValue $internalField; } top { type slip; } bottom { type slip; } airfoila { type fixedValue; value uniform 0; } airfoilb { type fixedValue; value uniform 0; } airfoilc { type fixedValue; value uniform 0; } back { type empty; } front { type empty; } } // ************************************************** *********************** // kt: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object kt; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.015; boundaryField { inlet { type freestream; freestreamValue $internalField; } exit { type freestream; freestreamValue $internalField; } top { type slip; } bottom { type slip; } airfoila { type fixedValue; value uniform 1e-11; } airfoilb { type fixedValue; value uniform 1e-11; } airfoilc { type fixedValue; value uniform 1e-11; } front { type empty; } back { type empty; } } // ************************************************** *********************** // nut: FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } exit { type calculated; value uniform 0; } top { type calculated; value uniform 0; } bottom { type calculated; value uniform 0; } airfoila { type nutLowReWallFunction; value uniform 0; } airfoilb { type nutLowReWallFunction; value uniform 0; } airfoilc { type nutLowReWallFunction; value uniform 0; } front { type empty; } back { type empty; } } // ************************************************** *********************** // |
All times are GMT -4. The time now is 17:40. |