CFD中文网

@chon01 这也只是我的猜测，因为我按了一下计算器，lg和ln函数算出来的结果差距确实比较明显，希望能帮到你吧。

@chon01 我在想有没有可能是他们写错了呢，如果用ln函数算，我认为你的结果是没问题的。但是如果将ln换算成底数为10的lg，好像就可以得到你最开始放的那张图里面的数值了

@chon01 我刚刚用计算器初步算了一下，用ln函数算出来当z=0.5时，u=0.88。用lg函数算出来当z=0.5时，u=0.38，跟你给出来的图是能对上的。这就比较奇怪了，是文献中把公式写错了还是咱们用错了。

你好，我请教一下，别人文章中用的ln函数，你写的是log函数，c++中ln函数是这么写的嘛？

@xpqiu 好的，谢谢您！

我最近在使用force函数监测流体作用力的过程中，出现了困惑。

流体作用力的计算应该是对物体表面做压力和粘性力的积分，openFoam在做二维计算时，需要在z方向设置一定的厚度，程序在计算流体作用力时是如何选取宽度的呢？是按照我们设置的厚度去计算，还是针对二维计算会默认取1m？

还有force函数输出的文件中，力的单位是什么？有没有大佬能给我解答一下。谢谢。

@Mina-Lee 好的，谢谢啦。我刚刚试了很久终于试出来了，确实是少给了参数。

@Mina-Lee 您好，可以再请教一下嘛？我扩散方案想试一下exponential，但是他会报错超出EOF，这是什么原因引起的呢？

@Mina-Lee 好的，我再排查一下。

@Mina-Lee 好的，谢谢您。我觉得应该是我扩散方案选的有问题，我再检查检查。

我网格变形是这样的，按照我的想法，整个变形域的网格应该一起变化才对，为什么现在他只压缩第一层网格呢？

@Mina-Lee 您好，我可以再请教您一个问题嘛？我使用了您提供的方法，现在上层网格固定住了，但是我计算过程中直接发散了，我检查了一下，是因为变形网格和静止网格交界处的网格被压缩的很小，而且他只压缩第一层网格，这是为啥呢？还是因为我扩散率选的不对嘛？

@Mina-Lee 谢谢您！

我在使用openFoam动网格的时候，希望能根据文献实现上面的三角形网格不变，下面的正方形网格变形，但是在使用过程中，我发现一旦下面变形区域内的网格变化，便会带动上面的网格一起变形，是否是因为我扩散模型选的有问题，我该如何实现呢？下面是我扩散模型的配置。

@yms00yms 不好意思哈，这我真不好说，那个cyclic边界条件我没用过。

@EricLiu 你用openFoam自带的那个函数去监测y+，别用你自己算的那个，不准。我个人觉得哈

@EricLiu 你的这些设置没有问题，我觉得是网格的问题，你从层流过渡到湍流，网格要更细一些，你算y+了吗？

@EricLiu 我觉得层流和湍流差异很大很正常吧，而且你要模拟液膜被切碎我感觉已经出来了，有没有文献对照一下？看看流速啥的，还有你的湍流边界条件具体咋给的？

@EricLiu 什么问题？

@李东岳 好的，谢谢老师。