请教大家一个问题， 我正在用twofluidmodel计算solid-water flow。但是我们有自己的模型，需要在原有的基础上改写代码，将solid的effictive visicosity等于water的。其中有一个phasepressure model，可以在这里面改写phase2.nueff_, 并且包含一些header文件，但是我这样写了之后，提示没有declaration。请问怎么解决这个问题呢？

解决了吗？最近在计算周期性的平板层流，感觉温度场和速度场传不进去

Hessian matrix

@莫晓柒 我现在用openfoam6.0中的rhoPimpleFoam，开启滑移网格计算也非常容易发散，后来无奈转为rhoPimpleFoam+MRF才可以计算。你可以把你的设置文件贴出来看下，这样很难看出问题。

不好人为干预，否则不符合物理

好的，谢谢你。我自己在看看相关资料。

@cccrrryyy 好的，谢谢！

https://turbmodels.larc.nasa.gov/spalart.html#sa
可参考这个

我的理解，对于可压缩流来说，求解的变量都是Favre平均量（密度除外，密度是时均密度）。所以，UPrime2Mean得到的也是基于Favre平均的脉动量，因为对于求解器来说它并不知道时均或者Favre平均，你让它算什么它就算什么。我看你贴出来的图的纵轴用的是两撇，这个一般就是Favre平均之后的脉动，因为一般时均后的脉动用一撇表示。

从这个角度讲，如果实验提供了Favre平均的相关统计量，我认为直接用模拟的结果去比对是没有问题的。如果实验只提供了时均的相关统计量，直接用模拟的结果去比对严格来讲不太正确。这里我也还没完全搞明白，看看有没有人来讨论了。

@东岳 谢谢东岳老师

https://www.jianguoyun.com/p/DdvrhJMQ9s3ZBhjw7p4B

可以看一下最后的方程

QQ截图20200113110636.png
不好意思，第一次发帖求助，排版有些不规整，重新上传图一（速度场）如上图所示。

好像是自己定义的壁面函数与foam自带的壁面函数重名了，改了名字后就没问题了。

nutWallFuntion是一个抽象类，它有什么作用？我看其他nutwall函数都是从它继承来的，这个抽象类是必须的吗？它具体是如何发挥作用的？有没有人能解释一下抽象类或推荐一些参考资料，谢谢。

@我是河滩 输出边界上的值列表可看出
1.JPG
2.JPG

@cccrrryyy 好的！谢谢~:xinxin:

请问最后你是怎么设置的初始场啊？我也是用LES直接上来步长实在是太小了，算起来一个case非常久。不知道可否用RANS先算然后转换到LES

@东岳 东岳老师，我可能没有表达清楚，您的意思是我的0文件夹中有一个O2的文件，然后这个C就是个符号而已，随意写，我在使用的时候就用C这个符号就好。
同理只要我0文件夹中有一个叫做C（碳）的文件，随意定义一个符号a，把O2换成C就可以用a作为碳的浓度去计算？

@我是河滩 是的。

@Zephyr QQ截图20200108095030.png
:xinxin: :xinxin2: :xinxin3:
果然可以。。以前还没注意过这个小齿轮按钮:xiexie: