后处理得到momentum thickness

[fangyuanaza]

请教各位老师，学生最近在模拟平板流动，想请教有没有在openfoam内部编程得到mmomentum thickness 而不是提取整个场的数据再在其他地方后处理的资料？

[李东岳]

@fangyuanaza 这个是连续形式，有离散形式的方程么

[fangyuanaza]

@李东岳 谢谢老师回复~ 离散形式应该是

本质上我是想请教如果对某个方向的量进行操作，当对某个量的网格点进行计算是是取parameter[celli], 若x方向不进行操作，有没有类似于cell(i,j)这种表达，可以仅对j的量进行操作？

[luofq-sysu]

@fangyuanaza 请问前辈实现了吗

[fangyuanaza]

@luofq-sysu 实现了，我意识到thickness 只能完成模拟后计算得到，因此，我是按照公式用python计算的。基本想法是在我感兴趣的区域监测了150个垂直于壁面的线，监测密度与速度。根据公式 得到每个位置的thickness,然后连接。具体做法是
1: 用python写一个sampleDict（不想手动写150个位置的监测信息）
2: 用python读取每个位置的信息，根据公式计算thickness.

[luofq-sysu]

@fangyuanaza 谢谢分享。还想请教一下：

1. 这个u，是流向速度的时均值吗？
   2.150个线，取的是不同流向位置吗？计算后流向大约有多少个网格呢？以及展向是统一取中间位置吗

[fangyuanaza]

@luofq-sysu 第一个问题：是的，我取的是流向速度时均值。
第二个问题，是的，是不同流向位置。我当时做的是一个平板流动，平板上流向的网格113个。150个不是一个必须的值，看你算例的情况跟需求。我的是二维流动，因此没有展向的问题。

[Amiu_Cheng]

想请教一下怎么进行无穷的积分？做的圆柱绕流在文献中也看到了 momentum thickness，但是不知道具体实现方法