CFD中文网

我不是大神，老菜了😂
不好意思，找不到了。
我搞错了，compressibleInterFoam也不是做气泡溃灭的吧。cavitatingFoam这种带空化的才可以吧。

不做😂
之前学习的时候用compressibleInterFoam练手算了一个，各方面设置也不是很严格

@吴小帅 前两天刚好看到这个知识，它叫虚假流动，和界面的算法有关。你可以再找找相关的讨论。我也不是很了解这方面的内容。希望能帮到你。
https://www.cfd-china.com/topic/2990/fluent-vof-耦合level-set方法时的density-correction-抑制虚拟流动
https://zhuanlan.zhihu.com/p/12448562051

再贴一下Pope的Turbulent flows书中的channel flow得到的壁面速度分布，里面的系数是5.2，而且还提到了拟合的差距应该在5%即4.94 ~ 5.46之间。

@coolhhh
一些槽道流DNS的文献里面应该有系数方面的讨论吧，我倒是没深入学习过这方面的内容。

目标$u_\tau$和实际$u_\tau$之间的差距，取决于槽道驱动的方式。因为计算之前基本上是不知道最后能算出来多少的，如果给一个固定压力梯度值，最后的结果不一定一致。但是可以在计算过程中不断地调整驱动的压力梯度值，到目标摩擦雷诺数之后，大概就平衡了，最终就会得到比如1000和550的摩擦雷诺数。
我计算的时候，550的使用了理论公式计算出来的压力梯度，下面这个：

学流体的小明 在 LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？ 中说：

方法四：
在完全发展的槽道流中
$$ \frac{{\partial p}}{{\partial x}} \times h = {\tau _w} $$

得到的结果就是550，零点几的差距。

计算$U^+$和$y^+$我都是用实际的$u_\tau$的。不过看情况喽，有理由的情况下用理论的$u_\tau$也可以。

我认为首先确实是在分辨率不够的情况下，channel395是算不好的。我后续算过摩擦雷诺数1000和550的雷诺数，在网格分辨率足够的情况下，对数区的速度都是符合$\kappa=0.4,\ -C=5.5$对应的对数关系的。
如果按照改变常数的思路去，得到的-C大概为7，应该不在以前的研究者们推荐的那个范围了。

12k网格的结果和15k网格的结果你都发一下。
map之后你画个速度，视觉上看起来就应该完全一样，只是15k网格实际上的cell中心的数据来自于插值嘛。

边界条件应该没啥问题，就是接触角边界条件，我用的时候可能会搞反0°和180°，你自己最好确定一下。
这种网格，大网格和小网格之间的变化是2倍的变化，有些地方生成的可能也不是很好。反正你现在是二维，要不就全加密到最小的尺寸？

我算气液两相流的时候也确实遇到了这样的问题，某些情况下气泡内部的气体速度明显比外部液体的大，但没你差距这么明显。我感觉我那个是比较合理的。
感觉和网格分辨率有关。你这样的情况，最小尺寸的那些固壁圆，网格分辨率有多大？有20个网格/直径？

编程解决的，直接改求解器

直接mapFields就可以吧？我用过，是可以把自适应加密之后的网格，映射到没有加密的规整网格上的。
如果是说机理上面可不可以，那其实是损失了加密之后那些网格的信息。你用另外一个计算模型也是要自适应网格吗？那可以直接把自适应加密之后的网格也拷贝过去。

@黒洛克 提取完之后就是foam格式的，算是ascii编码的数据，文本文件就能读。
我是用matlab处理的。
用tecplot的话，你试试vtk格式的？再不成，就自己拿ascii的写一个tecplot能处理的文件。

我不太懂网格
只是感觉还是得从网格着手吧

我自己的经验：同一块代码，可以看看不同进程的执行耗时，看看是哪个进程比较慢，也许就拖累了进度。
我也就使用reduce、gather、scatter这些进程间通信的函数。执行到这些函数的时候，快的进程会停下来等慢的进程就绪吧。

在controlDict中使用#includeFunc wallShearStress就可以得到wallShearStress的场数据，应该不是东岳老师说的那个原因。
wallShearStress只在边界的patch上面有数据，提取一条线的数据容易选到没有数据的内部网格上？所以没有文件输出吧。
用surfaces提取面上的数据好了

#includeEtc "caseDicts/postProcessing/visualization/surfaces.cfg"

fields       (wallShearStress);

surfaces
(
    bottomWall
    {
        type patch;
        patches ("bottomWall");
        interpolate     false;
        invariant       true;  // Unaffected by mesh motion
        surfaceFormat   foam; //ensight vtk foam
    }
)

迭代到1000次是没收敛，和网格、边界条件、时间步长都有关系。
10mm量级的网格，要画出10um的狭缝，网格感觉就很难画了。得保证狭缝那里有足够的分辨率吧，看看网格？

@风云5091 我们这些都是槽道流里面的实际坐标，除以特征长度，得到的无量纲化坐标y+。
你如果只做验证的话，OF输出的yPlus足够了，应该可以在paraview显示出来。
OF输出的yPlus确实只存在于那个边界patch上。

是的，praview这里会有问题，我也遇到过。应该是paraview没有正确识别noSlip边界条件的原因。

h大小会影响到槽道的无量纲参数，比如摩擦雷诺数$Re_\tau=u_\tau h / \nu$，你要算湍流的话，网格要跟着$Re_\tau$去画。当然不能简单地缩放。

我不理解这句话

你说的这个是绝对速度，看相对速度的话，和转轴没多大关系

旋转中心设置为-2 m，计算域还是0到0.04m，计算域距离旋转中心就很远了啊，旋转产生的速度$v=\omega \times r$也应该很大的吧。