Source[kth].ref() += (1 - kth)/(kth + 1)* pow(a[n], kth + 2)*g*w[n];

上面要改为

Source[kth].ref() += (1.0 - kth)/(kth + 1.0)* pow(a[n], kth + 2)*g*w[n];

否则 (1 - kth)/(kth + 1)在kth大于1的时候为0

东岳大哥您好，算几十万的网格是没有问题的，我扩展了虚拟机的内存到64g和虚拟内存到10g都不行。 。。也是奇怪

请问在OpenFOAM中，如何定义一个轴对称模型来仿真Taylor Couette flow？如何设定wedge边界？谢谢大神！

@mr-h 在 标量场壁面梯度的处理方式 中说：

在温度方程里的拉普拉斯项fvm::laplacian(ka,T) 中 存在ka的梯度的计算。

此处Ka为volScalarField，那么fvm::laplacian(ka,T)是不是等价于(fvm::grad(ka)&fvm::grad(T)) + ka*fvm::laplacian(T);这两种表达方式的意义一样吗？

@东岳 是类似与MultiRegionFoam的方法吗？存在多套网格分别对应不同的region？这种方法我没试过。我用的求解器是interfoam，现在只是想的是从一套网格上生成一个subMesh，这样绝大部分的网格节点是对应的，求解标量C方程所需要的数据可以通过原始的网格map过来。

图片
我最近在使用reactingParcelFilmFoam计算水滴与结构表面撞击并形成水膜的过程中，我想往模型中添加一个水滴的收集系数，这个系数需要知道水滴在撞击到结构时的位置坐标和初始的位置坐标，请问有人知道该怎么确定这两个坐标吗？

从代码来看，e表示回弹系数，mu表示一种损耗。e越大，回弹越完全，mu越大，速度被消耗的越大。这些值影响你颗粒撞击壁面的回弹行为。
你可以手动计算一下，看看U有什么变化

case PatchInteractionModel<CloudType>::itRebound: { keepParticle = true; active = true; vector nw; vector Up; this->owner().patchData(p, pp, nw, Up); // Calculate motion relative to patch velocity U -= Up; scalar Un = U & nw; vector Ut = U - Un*nw; if (Un > 0) { U -= (1.0 + e_)*Un*nw; } U -= mu_*Ut; // Return velocity to global space U += Up; break; }

AFAIK, you are the only one working in this area. Any update?

@东岳 好的，明白啦，我试一下。谢谢东岳老师。

@东岳 老师，您好。我目前在CFD中文网只发了一封求助帖。但因帖子中涉及到我们课题组提出的在某方面研究的新结构模型，目前还未申请专利和发表论文。为防止我们的结构创意被别人抢先用去，麻烦您把我目前唯一的这一封帖子给删了。发帖时是我有欠考虑，由此给您带来的麻烦，我十分抱歉！！谢谢您，麻烦了！！！

已删除

我在用interFoam时，笛卡尔网格可以得到合理的计算结果。但换成非结构网格（prism之类）后，计算虽没有发散，但界面有问题，请问有什么需要改进的地方吗？@东岳 @tidedrinker

无动网格，parcel量200000，k-epsilon，sprayFoam，spraycloud，solver和边界未进行任何其他编译

@东岳 好滴:xinlei: :xiezuoye:

同问，重心的位置怎么求，，，

@东岳 谢谢老师解答~

@东岳 感谢东岳老师的回答。

我参考了相关文献，low-Mach-number assumption 应该是忽略了压力波动，使得密度、焓等变量与压力波动无关。

那是不是意味着问题还是属于可压缩的呢？论文里面提到：

With these assumptions, the compressible energy equation is greatly simplified so that the resulting low-Machnumber equations are represented in tensorial notations as follows:

然后就是上面的三个控制方程。

那上面有关DNS的问题老师了解吗？

谢谢！

把文件里面的 simulationtype 设置成 laminar，网格足够密，时间步长足够小，是不是就算是DNS呢？

should be quasi-DNS

代码里面turbulence->divDevRhoReff(U)以及turbulence->alphaEff()这些湍流相关的项，会不会对DNS直接求解NS方程的方法造成影响？

如果设置为层流，turbulence->divDevRhoReff(U)直接返回层流项
因此不会有影响

@wwzhao 谢谢，sorry，没有仔细搜索，建议删帖

人找到了

@东岳 在 可压缩状态方程 中说：

可能是液体的压力和压强关系较大

我试过了，计算时这个方程可以考虑液体可压缩性并且结果也可以接受，本来我还以为有什么流体力学的理论给出的这个公式呢，谢谢东岳老师~