OpenFOAM テトラメッシュ（四面体）　流体解析

以前OpenFOAMでfluent用のメッシュをインポートする方法を

以下で紹介しました。

openfoamquest.hatenablog.jp

作成した条件でテトラメッシュを解析するとANSYS fluentと

計算結果が異なったり、発散して計算できなかったりしました。

そこで、OpenFOAMでテトラメッシュを使用しても

fluentと計算結果がかなり近くなった方法を紹介します。

OpenFOAM-v2006、simpleFoam　乱流モデルｋ-εモデルを想定

（1）systemフォルダにあるfvSchemesを編集

//////////////////////////////////////////////////////////////////////////////////

FoamFile{
　　version 2.0;
　　format ascii;
　　class dictionary;
　　object fvSchemes;
}
ddtSchemes{
   default steadyState;
}
gradSchemes{
   default faceMDLimited Gauss linear 0.333;
}
divSchemes{
   default none;
   div(phi,U) bounded Gauss linearUpwindV grad(U);
   div((nuEff*dev2(T(grad(U))))) Gauss linear;
   div(phi,k) bounded Gauss upwind;
   div(phi,epsilon) bounded Gauss upwind;
}
laplacianSchemes{
   default Gauss linear limited 0.333;
}
interpolationSchemes{
   default linear;
}
snGradSchemes{
   default limited 0.333;
}
fluxRequired{
   default no;
   pcorr ;
   p;
}

//////////////////////////////////////////////////////////////////////////////////

（1）systemフォルダにあるfvSolutionを編集

FoamFile{
   version 2.0;
   format ascii;
   class dictionary;
   object fvSolution;
}
solvers{
   p{
     solver GAMG;
     tolerance 1e-06;
     relTol 0.01;
     smoother FDIC;
     nPreSweeps 1;
     nPostSweeps 2;
     nFinestSweeps 2;
     cacheAgglomeration on;
     agglomerator faceAreaPair;
     nCellsInCoarsestLevel 1000;
     mergeLevels 1;
   }
   U{
   type coupled;
   solver PBiCCCG;
   preconditioner DILU;
     tolerance (1e-6 1e-6 1e-6);
     relTol (0.01 0.01 0.01);
   }
   k{
     solver PBiCG;
     preconditioner DILU;
     tolerance 1e-06;
     relTol 0.01;
   }
   epsilon{
     solver PBiCG;
     preconditioner DILU;
     tolerance 1e-06;
     relTol 0.01;
   }
}
SIMPLE{
   nNonOrthogonalCorrectors 3;