pimplefoam在瞬态模拟过程中库朗数突然增大然后报错停止
-
最近在使用openFOAM对下肢血管进行血流动力学模拟,为了更准确所以在入口处使用了随时间变化的边界条件,但计算一段时间后,库朗数会突然增大,然后程序报错停止,有没有大神可以帮忙看看应该如何修改设置。
版本:openfoam 10
求解器:pimpleFoam
U文件,入口速度方向是z轴的负方向boundaryField { WALL { type noSlip; } INLET { type uniformFixedValue; uniformValue { type tableFile; file "v.txt"; } value (0 0 -0.032180); } OUTLET1 { type zeroGradient; } OUTLET2 { type zeroGradient; }
controlDict文件
application pimpleFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 0.7; deltaT 1e-5; writeControl adjustableRunTime; writeInterval 0.01; purgeWrite 0;
fvSolution文件
solvers { p { solver GAMG; tolerance 1e-7; relTol 0.1; preconditioner DIC; smoother symGaussSeidel; nCellsInCoarsestLevel 10; mergeLevels 1; } pFinal { $p; relTol 0.1; } U { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.1; } UFinal { $U; relTol 0.1; } }
fvSchemes文件
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; }