: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中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的

[2] [1] ##00 [38] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:? [38] #1 Foam::sigFpe::sigHandler(int) at at ??:???:? [1] [2] ##11 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int) at ??:? [38] #2 ? at ??:? [2] #2 ? at ??:? [1] #2 ? in "/lib64/libpthread.so.0" [38] #3 ? in "/lib64/libpthread.so.0" [2] #3 ? in "/lib64/libpthread.so.0" [1] #3 ? in "/lib64/libm.so.6" [38] #4 ? in "/lib64/libm.so.6" [2] #4 ? in "/lib64/libm.so.6" [1] #4 ? in "/lib64/libm.so.6" [38] #5 ? in "/lib64/libm.so.6" [2] #5 ? in "/lib64/libm.so.6" [1] #5 ? in "/lib64/libm.so.6" in "/lib64/libm.so.6" in "/lib64/libm.so.6" _pmiu_daemon(SIGCHLD): [NID 00126] [c0-0c1s15n2] [Tue Jul 9 10:43:27 2019] PE RANK 38 exit signal Floating point exception _pmiu_daemon(SIGCHLD): [NID 00082] [c0-0c1s4n2] [Tue Jul 9 10:43:27 2019] PE RANK 1 exit signal Floating point exception [NID 00126] 2019-07-09 10:43:27 Apid 8697814: initiated application termination

I have tried it again, the same problem occurs.
Last time i have used shear gradient boundary, and this time i just use zeroGraident. However, ....

I guess. the problem is the signal floating exception???

Tks

@队长别开枪

谢谢您的建议，感觉这个方案可行。 我尽早试试，给您答复。。。

谢谢

把你的新库通过controlDict挂载一下，在controlDict里面添加

libs ( "libnew1.so" "libnew2.so" );

有没有老哥用过SPARTA这个软件，有没有什么好的入门读物推荐呀？谢谢啦（不知道帖子该放哪就发到这里了）

@夏雨天 你把7.55-57三个方程加起来，就会发现他们用的是C-N格式，三者加和是7.59。在书上179页最底下也说啦 :duang:

@东岳 谢谢东岳老师的关心，后来还是有点进展，不过进度不是很大，现在所走的路就是修改求解器，文章还没有出来。等出来之后我会详细说下这个问题。

看起来只不过是模拟准确不准确的问题哦？可能下图模拟出来就是这样啊