@qjh888 您好，不知道您如何更新的solver？方便讲一下吗？最近我也在做相关的事情。

各位CFD朋友，目前我想学习一下firefoam，但是baidu搜到的资料实在有限，如果方便的话能否提供一些算例资料给我，谢谢大家。

在求解可压缩流动时，出现了如下情况：导致最终结果发散，可能是哪个参数的问题呢？

rhoEqn max/min : 1.21419 0

@东岳 请问东岳老师，可以用速度基求解器比如rhoCentralFoam求解这个算例么？还得加个温度去解压力？
没算过DNS，望指点！谢谢！

@lllwonderliquid 已更新simpleFoam解析 http://dyfluid.com/simplefoam.html

已经解决了。谢谢您。可能是自己安装的经验不够吧。

滤镜-time annodate
即可

:tishizi:

谢谢，谢谢。

@程迪 程老师您好，请问您知道如何把Prt定义成某个表达式吗？

东岳老师好，大家好：
我是正在学习OpenFOAM的新手
最近在学习模拟湍流时运行pisoFoam下自带算例porousBlockage时出现错误：

Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present Creating finite volume options from "constant/fvOptions" Selecting finite volume options model type explicitPorositySource Source: porosity1 - selecting cells using cellZone porousBlockage --> FOAM FATAL ERROR: Cannot find cellZone porousBlockage Valid cellZones are 0() From function void Foam::fv::cellSetOption::setCellSet() in file cellSetOption/cellSetOption.C at line 142. FOAM exiting

请问是什么原因呢？

此外，还想问一下，如果使用标准k-epsilon模型（RAS算例下）模拟湍流，进出口使用什么条件比较合适呢？

如果想要设置两块旋转域（对转），dynamicMeshDict应该怎么编写呢？我像下面那样编写，结果只有sector_part2旋转。

dynamicFvMesh dynamicMotionSolverFvMesh; motionSolverLibs ("libfvMotionSolvers.so"); motionSolver solidBody; cellZone sector_part1; solidBodyMotionFunction rotatingMotion; origin (0 0 0); axis (0 0 1); omega 112.573737; // rad/s cellZone sector_part2; solidBodyMotionFunction rotatingMotion; origin (0 0 0); axis (0 0 1); omega -112.573737; // rad/s

正确的编写形式应该是怎样的呢？

@星星星星晴

:wolaile: 没太细看哈不过

我需要return dClass，PUr和PUz。

你可以把这三个打包成一个struct或者class，而不是三个哦，然后返回这个struct或者class，差不多这种

class test { public: scalar A_; scalar B_; } test returnFunction(test& f) { f.A_ = 1; f.B_ = 2; return f; }

想问一下，针对复杂几何，大家都用什么软件画然后导入啊？

还没关注过这个问题 :xinlei:

对于使用KomegaSST模型作低雷诺数的湍流管流模拟，当y+值小于1时，是不是可以认为网格足够精细，无需再细化？
另外k omega nut三个量大概在哪个范围内合理？我的结果大概是以下：

k 10的负3次方

omega 最大50w+，最小500左右
nut 10的负10到负5次方

这里整理了一些 找对象 的诀窍:xinxin:

应该是Sb.boundaryField() 返回const reference
Sb.boundaryField().ref() 返回 reference
根据我的经验是这样,如果错误的话把编译错误贴上来
认真看看代码很简单的,编译器的错误已经告诉你所有信息了!

《A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM》一文2.6节描述了nutWallFunction，nutAtmRoughWallFunction中nutw公式与之相同，其他量有改变，我再看看。

最近计算一个方腔流动遇到一些问题，这个流动沿流向（x方向）是周期性边界，yz方向是壁面，流动固定x方向的平均速度，自动求解压力梯度，在采用rans计算时是正常的，可以计算出一个合理的结果

pressure gradient = 0.179097

但是当采用PANS（自己编译的湍流模型）计算的时候，压力梯度会变得非常小，导致计算结果很差，

pressure gradient = 0.0903267

自己编译的模型可以在channel395中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的