我觉得最简单的方法是适当的延长一段出口增加空气域，很多SCI都这样做防止出口的影响

@东岳 谢谢东岳老师，我正在修改弹性模量，让变形小一点。

这样的话，管流的速度不知道为何会持续增涨，然后在入口处变成一个非常大的值后计算崩溃。我之后改成了入口为totalPressure，才能实现稳定的速度，计算也能保持稳定。请问这两个边界条件分别在什么情况下使用呀？

谢谢。

@范准 非常好，感谢分享 :xiexie:

@东岳 @宝丁 哪里哪里，举手之劳而已

参考 http://www.cfd-china.com/topic/933/swak4foam-和-openfoam4-1/2

@郑澄远 你的0文件夹下的U文件有问题。应该是初始场的个数与网格单元数不一致。

@dzw05
k：

dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.081; boundaryField { inlet { type zeroGradient; } outlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type kqRWallFunction; value $internalField; } }

nut

dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value $internalField; } outlet { type calculated; value uniform 0; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type nutkWallFunction; value uniform 0; } }

omega

dimensions [0 0 -1 0 0 0 0]; internalField uniform 0.095; boundaryField { outlet { type zeroGradient; } blades { type omegaWallFunction; value $internalField; } inlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } }

p

dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type zeroGradient; } }

u

dimensions [0 1 -1 0 0 0 0]; internalField uniform (14 0 0); boundaryField { inlet { type fixedValue; value uniform (14 0 0); } outlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type rotatingWallVelocity; origin (0 0 0); axis (0 0 1); omega 0; } }

boundary

inlet { type patch; nFaces 400; startFace 6377766; } outlet { type patch; nFaces 400; startFace 6378166; } frontAndBack { type symmetry; inGroups 1(symmetry); nFaces 4800; startFace 6378566; } blades { type wall; inGroups 1(wall); nFaces 80334; startFace 6383366; } inter1 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 13994; startFace 6463700; matchTolerance 0.0001; transform noOrdering; neighbourPatch inter2; } inter2 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 11572; startFace 6477694; matchTolerance 0.0001; transform noOrdering; neighbourPatch inter1; }

fvSchemes

ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.5; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; }

fvSolution

solvers { "pcorr.*" { solver GAMG; tolerance 0.02; relTol 0; smoother GaussSeidel; } p { $pcorr; tolerance 1e-7; relTol 0.01; } pFinal { $p; tolerance 1e-7; relTol 0; } "(U|k|omega)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-06; relTol 0.1; } "(U|k|omega)Final" { $U; tolerance 1e-06; relTol 0; } cellDisplacement { solver GAMG; tolerance 1e-5; relTol 0; smoother GaussSeidel; } } PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 5; pRefCell 0; pRefPoint 0; pRefValue 0; } relaxationFactors { fields { p 0.2; } equations { "(U|k|omega)" 0.7; "(U|k|omega)Final" 0.7; } } cache { grad(U); }

站在2019年，这个问题就会很有意思了。

首先，从软件的层面，虽然我还咩有机会仔细比对fluent和openfoam的区别，但是从业界的反应来看，开源和商用上纯理论层面没太大区别。只是界面和用户习惯上的问题，这个会是放在企业抉择面前的最大问题，也是目前商业帝国的基础。

企业主要就看性价比。就目前的情势而言，企业在不考虑正版化的情况下，会优先选择fluent。比较好用又不要钱。企业的选择：

使用环境：windows，linux，好像前者简单 人才成本：薪酬、培训费，目前fluent人比较好招 版权费：不存在的。
所以综合下来fluent不赢都难。

但是随着最近正版化越收越紧，微软win10的价格其实非常公道。所以很多企业都会考虑尽可能上正版win10. 而对于of来说，win10是他们的朋友，因为ubuntu内嵌，最近连linux格式都支持了。那安装oe的成本其实很低。所以就平台而言，oe和fluent已经站到了同一起跑线上。

下面的问题就是工程师喜欢哪个了，如果oe的生态链（应用、可视化、培训）整体也上来了，企业的选择就会变成：

使用环境：都是win10，没差 人才成本：薪酬、培训费，fluent大概略低，但差不了多少 版权费：fluent肯定贵

所以，oe的机会就在于怎么把应用做得更好用，类似“傻瓜都能收敛” + 增加用户基数。

which，也就是东岳目前正在做的。

@LiuZhuang http://xiaopingqiu.github.io我基本上看的都是这个人写的总结

各位大佬，请教2个问题

是否可通过paraview获得流场体积流率（或者质量流率）？流场是一个回路，靠泵驱动进行流动，我想知道流场的体积流率 由于OpenFOAM计算不可压缩流动是没有密度rho的，这就导致最后得到的湍流粘度场为运动粘度nut(单位是m2/s)，但是我想获得粘度场（单位为Pa·s），那么nut的值都需要乘以密度的值，paraview只能显示nut的云图，可否用paraview对某个场乘以一个固定值？获得粘度场的云图呢？

望不吝赐教，感谢！

QQ截图20190318163023.jpg
压力云图是这样的。

@东岳 好吧

@东岳 在 OpenFOAM多组分返回热物性 中说：

实际就是rho_

这样！明白了明白了，谢谢。

@东岳 好，谢谢，得啃求解器了，

给大家分享一个编译安装的例子，在Ubuntu18.04上安装OpenFoam2.3.1
https://openfoamwiki.net/index.php/Installation/Linux/OpenFOAM-2.3.1/Ubuntu#Ubuntu_18.04

@cccrrryyy 非常感谢

还有 轴缩放

@LiuZhuang 在 twoPhaseEulerFoam下的bubbleColumn算例的问题 中说：

在outlet边界上气体自由逸出，液体则为壁面（zeroGradient）

好长没太细看:zoule: 不过我觉得你可以设置液体壁面速度为slip试试