@bestucan 谢谢老师,这个问题我再去认真了解一下并尝试尝试。
Latest posts made by Caijinjin
-
RE: 关于multiphaseEulerFoam运算分散相体积分数问题
@caijinjin 老师,这是我的离散格式和求解器设置
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; limited cellLimited Gauss linear 1; } divSchemes { default none; "div\(phi,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss vanLeer; "div\(alphaRhoPhi.*,U.*\)" Gauss linearUpwind limited; "div\(phi.*,U.*\)" Gauss linearUpwind limited; //"div\(alphaRhoPhi.*,Yi\)" Gauss linearUpwind limited; "div\(alphaRhoPhi.*,(h|e).*\)" Gauss linearUpwind limited; "div\(alphaPhi.*,f.*\)" Gauss linearUpwind limited; "div\(alphaRhoPhi.*,K.*\)" Gauss linearUpwind limited; "div\(alphaRhoPhi.*,\(p\|thermo:rho.*\)\)" Gauss linearUpwind limited; "div\(alphaRhoPhi.*,(k|epsilon|omega).*\)" Gauss linearUpwind limited; "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; } wallDist { method meshWave; nRequired yes; }
求解器设置
solvers { "alpha.*" { nAlphaCorr 1; nAlphaSubCycles 1; } droplets { nCorr 1; tolerance 1e-4; scale true; solveOnFinalIterOnly true; sourceUpdateInterval 20; } p_rgh { solver GAMG; smoother DIC; tolerance 1e-8; relTol 0.01; maxIter 20; minIter 2; } p_rghFinal { $p_rgh; relTol 0; } "U.*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; minIter 1; } "(e|h).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1; relTol 0; minIter 0; maxIter 0; } "(k|epsilon|omega).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; minIter 1; } "f.*" { solver PBiCGStab; preconditioner DILU; tolerance 1e-6; relTol 0; } } PIMPLE { nOuterCorrectors 3; nCorrectors 1; nNonOrthogonalCorrectors 0; pRefPoint (0.0276 0 0); pRefValue 101325; } relaxationFactors { equations { ".*" 0.7; } }
-
RE: 关于multiphaseEulerFoam运算分散相体积分数问题
@李东岳 老师,我这个网格是三维的,图上展示了二维,用的结构网格,of9计算。目前发现我只要把流速设的大一点,就不会有这种问题出现,但是流速变小就不行了。另外,网格壁面不加密且数量少的时候它也不报错,自己也还在找原因中。
-
关于multiphaseEulerFoam运算分散相体积分数问题
目前本人在用multiphaseEulerFoam求解器求解孔板中的油水两相流动,孔板架构如图
模拟过程中水是分散相且体积分数设为了0.005,运算时发现有时候运算会报错,计算过程显示Mules算法求解相分数时,不管连续相还是分散相,最小相分数都为0,偶尔为负。请问各位老师和同学们有没有了解这一问题的?不吝赐教。
-
RE: 离散法求解群体平衡方程中的程序含义
我自己暂时的理解是,程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目,也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的,程序里面为什么还要再乘以i颗粒体积fi.x()?这样子的话量纲就不是每秒了。
-
离散法求解群体平衡方程中的程序含义
在离散法求解群体平衡方程中,程序文件中颗粒聚并生成项(birthByCoalescence)和文献中的给的式子形式不太一样,请问哪位大佬可以帮忙解释以下程序中为什么还要多乘以颗粒的直径?
方程里面η在程序里面应该对应gamma,指的是颗粒的分配体积,c对应程序里的coalescenceRate -
RE: 并行计算速度没有变快
@李东岳 我台式是AMD 3600 6核, 工作站AMD EPYC 7502 32核,平时用都快不少,就这个算例基本没有加速效果,可能和网格有关系吗?
-
并行计算速度没有变快
最近画了个釜式反应器,然后用的是结构化网格,一共50多万个六面体,使用自己台式机分四块并行计算和工作站分24块计算作比较,发现计算速度是一样的,请问这是什么原因?按理说速度最起码要提升一点啊,一点没提升是什么问题?网格是这样比较规整的