计算出现负相率和负湍动能
-
如图，是一个搅拌釜内的液固两相流，用的twoPhaseEulerFoam,计算过程中逐渐发散，出现了很多负值，请问这个可能是什么原因导致的？对网格进行了质量检查，网格为非结构网格，这个非正交性的值是偏大还是正常范围？
fvSolution设置为：（为了不求解温度，我将h/e的迭代数设为0，是否对计算有影响？）solvers { "alpha.*" { nAlphaCorr 1; nAlphaSubCycles 2; smoothLimiter 0.1; implicitPhasePressure yes; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-9; relTol 0; minIter 1; } p_rgh { solver GAMG; smoother DIC; tolerance 1e-8; relTol 0; } p_rghFinal { $p_rgh; relTol 0; } "U.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-5; relTol 0; minIter 1; } "(h|e).*" { solver smoothSolver; smoother symGaussSeidel; // tolerance 1e-6; relTol 0; minIter 0; maxIter 0; } "Theta.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-6; relTol 0; minIter 1; } "(k|epsilon).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-5; relTol 0; minIter 1; } } PIMPLE { nOuterCorrectors 3; nCorrectors 1; nNonOrthogonalCorrectors 0; } relaxationFactors { equations { ".*" 1; } }希望各位大佬帮忙解答，感谢！
-
另外附上我的离散格式，出现负值是否因为离散格式设置错误？未包含有界限制？
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default none; "div\(phi,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss vanLeer; "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaRhoPhi.*,(h|e).*\)" Gauss limitedLinear 1; "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,p\)" Gauss limitedLinear 1; div(alphaRhoPhi.particles,Theta.particles) Gauss limitedLinear 1; "div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1; div((((alpha.water*thermo:rho.water)*nuEff.water)*dev2(T(grad(U.water))))) Gauss linear; div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I))) Gauss linear; } laplacianSchemes { default Gauss linear uncorrected; bounded Gauss linear uncorrected; } interpolationSchemes { default linear; } snGradSchemes { default uncorrected; bounded uncorrected; } -
