CFD中文网

650a2e12-bcea-4ddf-aed1-35dea47389d1-image.png

fluent燃烧问题，ED模型，fluent自带的二步甲烷反应与将二步甲烷反应方程及其置前因子、活化能按同样数值输入后，模拟温度云图差异较大是什么原因

https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C

去这里加个源项不就好了？

// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(alpha, rho, k_) + fvm::div(alphaRhoPhi, k_) - fvm::laplacian(alpha*rho*DkEff(), k_) == alpha*rho*G - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_) - fvm::Sp(alpha*rho*epsilon_/k_, k_) + kSource() // ? + K*k_ );

谢谢东岳老师回复！
试过很多方法，确实没法通过网格划分进行处理，非结构网格也不行；
目前使用的办法是在不改变圆柱和球体公称通径情况下，把角度很小那一部分切除掉，人工做出间隙，这样模型几乎没有变化，而且可以计算，跟您说的办法原理相同。

@ibelief 谢谢，我再看一看。

@田畔的风 这个可能是个bug，但是不确定是不是bug，我自己的双路7702试了一下，关闭一半ccd变成双路64核可以用msmpi识别到两个cpu的物理核心，但是全部128核的情况下只能识别到最多64物理核心，后来我发现如果在后台开一个空的msmpi并行的fluent，再打开一个并行超过64核的就没有问题，只是读取文件的时间有点慢，我试过numa节点设置为0，1，2，4，都是一个情况，只要后台打开一个以msmpi并行的fluent，再打开一个跨越双路u核心的并行就没有问题

今天在CFD-online上看到一个帖子，也在问of2.1.1中的cyclic边界条件，看了一下他遇到的问题跟我一样（包括使用的求解器也一样），所以现在认为应该不是我的Case中设置的问题，应该是openfoam2.1.1中cyclic边界条件本身存在bug。打算从一个高版本的of中copy出cyclic边界条件，再自行编译为用户自定义边界条件，并没有成功。目前还没想到更好的解决办法。

最近在看dynamicRefineFvMesh类，发现他和普通的八叉树网格细化有一些区别，有没有人知道这个算法最早是谁写的？有没有相关的参考文献？

早期的图形界面一般在后面，tty6还是7.
前面的都是终端界面。
后来改了，tty1是加载信息的界面，就是开机的时候一堆字码。最终是登录界面的样子。tty2是登录后的图形界面，往后才是终端界面。
而界面切换的工作量还是有的，如果当前系统有什么占用磁盘读写的程序在卡着，切换会相当慢。甚至切到终端界面，打字母都有一两秒延迟。如果看见没反应持续操作，就会给卡死的系统卡的更死。所以，和平解决方案，就是--切一下，然后等。进top杀死高占用，或者直接重启

大佬们 求教 我做制冷剂流动沸腾 想得到沿着管程的干度
第一种想法：我创建了十个截面，得到各截面的截面含气率，fluent能否得到滑移比，然后通过公式去计算得到干度？

第二种：通过热平衡公式去计算得到？ 但是我通过热平衡公式计算发现算出的干度非常小。

@shifc 不知道你现在解决了没有，我最近做格栅湍流也碰到了同样的问题，我的解决方法是谱空间的复速度对三个坐标轴必须都是共轭对称的。

你好👋

缅怀。。作为CFD萌新，已经无数次地看过这位老教授的名字

@upupsquirrel 你好 请问同样的算例我用大涡fluent时均后得出的resolved reynolds stress和直接用雷诺应力模型得出的结果为啥差了超级多呀?感觉用大涡得出来的结果跟不稳定一样，是我的网格还不够细吗？

@xjwwjl61 有偿有偿啊啊啊

已更新，多谢！

@bestucan 非常感谢老师的提醒，让网格变稀疏之后发现了两处可疑的质量较差的地方。尝试之后发现，有一些长方体网格划的太瘦了，ICEM的网格质量检查没有问题，修改之后OpenFOAM中就也没有问题了。:huahua:

今天青年节，下午放假 :chigua:

这个是理论大坑啊。

如果压力基的求解器，求解的是速度和压力场，那么得用连续性方程来求压力。
如果是密度基的，连续性可以求密度，然后用状态方程求压力。

0_1525756293789_捕获.JPG

什么原因？

ok，感谢。