CFD中文网

我是用外部导入的网格，计算过Rans，想直接把它加密用于LES计算，大概步骤是怎样的？

mesh refinement好像只是加密局部的网格，如何使用它能加密全部的？它的具体用法怎样的

@wwzhao 好的，非常感谢您！我及时用SST的。 这里y+是指网格的吗？我的是接近1的。

由于我用的是extend版本的of，所以没法升级到3.0以上，要能用omegaWallFunction要做什么改动呢？kqRWallFunction是能用在DES上的，同理的话，如何把omegaWallFunction也链接上？

实在不行，我的omega壁面边界只能用fixedvalue了，那样会影响很大吧？

有没有做DES的，我一点都不懂。DES模拟壁面k,omega之类的是不是也要用壁面函数？ 我使用基于komega的DES，omega使用壁面函数就会运行不了，使用定值或者零梯度就能运行很好。使用定值是不是也行，会不会对结果影响很大？

没人吗？ DES是不是一定要用壁面函数。用fixedvalue行不行？

请问在LES里面没有k,omega等的壁面函数，只有nuSgs啥的，那k和omega之类的wall怎么给边界条件才合适呢？

有没有人声明过LES的k,omega的壁面函数？求分享。

@wwzhao谢谢，已经解决，在extend里面有GlobalPatchesZones，设置这个就能运行了

@wwzhao 我之前打错了，我就是已经使用了preservePatches，最终得到两个周期性边界的nFaces数不一样，从而在decomposepar的时候，出现了uncovered faces founded on master：1625 on slave：0的警告。然后去运行算例两步就发散了。
对于存在cyclic边界的问题，怎么能保证每个processor都一样的nfaces？感觉办不到呀，您的做过的这方面的case能办到吗？

我的case里面存在cyclicggi或者cyclicAMI的边界，decomposepar该怎么办呢？常见的几种我都试过了，出现了两个面网格数目不一样的警告，运算跑两步就发散了。
已经使用了conversePatth了，还是不行，网上说要把周期性的两个面放在一个块里面，怎么办到呢？（不并行运算是没问题的，能跑的）

例如在RASProperties里定义了：

kOmegaSSTCoeffs
{
Prt 0.9;
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
c1 10;
}

怎么在代码里读取这里Prt？

@cfd-china 就是用groovyBC才忽略的为什么。。。fixedvalue是没问题的，每一步都是125000，不会像一直这样的规律增

@cfd-china 我就按照你给的这个运行了几步，还是一样的，这是log文件：0_1461745568682_PyFoamRunner.transonicMRFDyMFoam.logfile

我用的是of-ext-1.6+swak4Foam

@cfd-china 不管我怎么改valueExpression，比如说我直接把valueExpression定义为0，value还是125000. 然后我监测出口压力，压力第一步是125000然后就一直在增大了。 而且我试过不同的式子，就算valueExpression里的式子是错的，出口压力还是一样的增大，计算也不报错。 这样肯定是没有起作用了，如果起作用了，式子错了好歹会报错算不下去了

我模拟叶轮机械出口边界用到径向平衡方程：
为此我在在出口压力边界设置了groovyBC：

OUTLET
{
    type            groovyBC;
variables	"r=sqrt(sqr(coordsX)+sqr(coordsY);Cu=coordsX/r*U.component(vector::Y)-coordsY/r*U.component(vector::X);C@HUB=rho*sqr(Cu)*log(r);";
    valueExpression	"rho*sqr(Cu)*log(r)-C+125000";
value		uniform 125000;
}

但是不知道为什么不起作用，不管式子对错都能算，故意把表达式改为常数0之类的都不行，出口压力一样按照value给的数值按照一样的规律往上升。坐标用pos().x我也试过。问题在哪呢？谢谢！