用twoPhaseEulerFoam求解水砂输送的问题
-
@东岳 非常感谢东岳老师这么晚还在回复问题!!!希望您也能休息好!!!目前有三个问题:
(1)我在residuals字典里面写明了U.particles(U.particle监控不了),才出现的N/A这个问题;
(2)模拟结果表明砂堤没有沉积(注入的固相百分数就为0.15,而且停止注入固相后,固相整体被往前推动),不知道是不是颗粒动理论模型设置的问题
(3)液相、固相注入速度都是(0.5 0 0),但计算过程中底部达到4m/s+,这是否合理(液相的internalFiled也给的0.5,感觉液-固混合液比纯液跑得还快了):
(边界给的固定值(0 0 0),插值看不到结果 )希望您有空的时候能帮忙解答,不胜感谢!
-
@东岳 老师,您看:
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; div((((alpha.particles*thermo:rho.particles)*nuEff.particles)*dev2(T(grad(U.particles))))) Gauss linear;//改层流后多出 } laplacianSchemes { default Gauss linear uncorrected; bounded Gauss linear uncorrected; } interpolationSchemes { default linear; }
-
@李东岳 东岳老师,在kineticTheoryModel.C中有这样一段代码:
Foam::tmp<Foam::volScalarField> Foam::RASModels::kineticTheoryModel::k() const { NotImplemented; return nut_; } Foam::tmp<Foam::volScalarField> Foam::RASModels::kineticTheoryModel::epsilon() const { NotImplemented; return nut_; }
但是感觉固相湍流开关打开的话,其湍流黏度完全是按照如下公式计算的:
(其中a为固相,b为液相)
也就是固相被看作连续相后不存在真实的k和epsilon哇?