CFD中文网

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次setFields，给加密后的较密的网格赋值，再去分块。
分块的话也得想想加密的部分会不会在你流场里面移动，你的不同核之间的负载要基本平衡，可能还需要你手动去划分区域。

感谢大佬回复@coolhhh

这两个图分别对应 方法1 和 方法2 的计算结果吗，还是不同方向的结果？

不同方向的结果。$\xi$是$x$方向的间距，$\zeta$是$z$方向的间距。上面我画的都是相关系数，也就是相关函数进行了归一化，所以在间距$\xi$或$\zeta$为0时，相关系数都为1。

看起来对于槽道流这种展向和顺流向是均匀+时间平稳的流动是可以用时间平均之后的。问AI，AI也是这么说的。

我的数据在时间上和空间上也比较多，$\phi(x,z,y=y_0,t)$每个计算网格的值都输出了，时间上也以很小的时间间隔输出了一段时间。无论是时间上的相关系数还是空间上的，计算起来都没什么难度。

@CFD聚沙 我不太理解你的意思。哪个方法呀？
具有周期边界的槽道流初始化的方法，和实际计算时候如果你想得到一个湍流的来流的方法，肯定不一样的。
实际情况下只设置入口边界是不太行的吧，想得到湍流的来流还得发展一段距离。
我印象里论坛里也有好多讨论入口湍流的帖子。

@先先贝
https://cfd-china.com/post/37664
可能理解错你的意思了，我和楼主之前是在另一个地方讨论自适应网格加密之后初始化的，然后他这里又出了这个问题，我之前的操作也遇到过。

@liuxin 借这个楼再问一些相关的问题
相关性的定义，我看到两种。第一种是某个时刻的相关，第二种是在第一种的基础上，在时间上也进行了平均。
第一种定义见Pope, S. B. 2000 Turbulent Flows.

第二种定义我在这个报告中找到 https://www.researchgate.net/publication/304943166_Large-Eddy_Simulation_of_Turbulent_Channel_Flow 。可以看到时间 t 在平均之后消失，而且作者在前面也定义了此方括号包含了时间上的平均。

应该用哪一种互相关函数来证明计算域足够大呢？
Pope书中没有具体的DNS计算设置。
LES报告中的公式以及结果，看起来也都是第二种定义。
这两种定义的结果，画出来确实不太一样。我用$y^+=10$这个平面上的速度脉动来做互相关分析。

上图：某个时刻$t_0$的两点速度脉动的相关函数。可以看到它的起伏还是比较明显的，特别是在展向$z$方向，由于湍流的条带确实存在，随着间距$\zeta$的增大，相关函数是会在0的上下波动。

上图：某一段时间平均之后的两点速度脉动的相关函数。起伏消失，整个相关性曲线很平滑。

目前看起来，Kim, J., Moin, P. & Moser, R. 1987 Turbulence statistics in fully developed channel flow at low Reynolds number. J. Fluid Mech. 177, 133-166.的结果是时间平均之后的两点相关，他们的标识里面也没带时间。
另外想问一下第一张图是那一片文章，有提到过是时间平均的还是某个时刻的吗？

请大佬们指点🙏

@先先贝
串行就算几个时间步，顶多三个，三次自适应加密，网格量已经会增加很多了。算三步消耗的时间也不会太久的。
可能我当时把时间步长也调的很低，1e-8s算三步，算挺快的。

在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？

这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。

楼上说的对，看你的边界层能不能达到半高也就是0.1m的厚度。
看你入口出口边界条件，既然你加压力梯度的话，感觉入口出口也是周期性边界条件。流场最后一定会发展到0.1m的厚度的。所以算的时候最好也把初始化搞好。

我不是大神，老菜了😂
不好意思，找不到了。
我搞错了，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能处理的文件。