CFD中文网

能否根据一些散点的坐标xyz，即这些散点的风压值，生成paraview可以读取的格式让其能插值显示曲面风压？

现在有一个类似冷却塔的空间曲面的模型以及它的表面散点风压的实验数据，请问下有什么办法能利用paraview处理这种数据得到三维曲面风压云图吗，因为python上面似乎找不到能直接根据散点坐标生成这种任意空心空间曲面云图的

@lizhisongsjtu 又测试了一下，似乎对于blockMesh生成的网格都不会出问题，SHM生成的网格有一定程度会出问题，然后pointwise生成的网格就都会出问题报这种错误

@lizhisongsjtu 感谢分享，请问下在of7下同样编译成功，但运行Rotex -latestTime会出现这种报错是什么原因呢

@xpqiu 谢谢您！我大致明白了，去试一试

您好，可以请教下binData中设置什么含义吗，如果我想沿圆柱展向z方向取各个截面三分力系数

用SHM对一个旋转过的立方体网格生成时，stl文件小立方体表面的网格还是按照原有坐标系下加密，立方体的角点网格贴体效果不是太好

请问下有哪个参数可以控制网格在一定区域内加密按旋转后的方向加密吗像下面这样
或是有没有其他优化贴体网格的建议？

顶一下，现在有方法解决了吗？最好各个时刻的数据能保存到一个文件中便于python读取分析

@coolhhh 感谢！

@coolhhh 感谢回复！另外您引的帖子中没有找到大库朗数算例的网盘诶？

@coolhhh 感谢回复！这两天做了两个算例

1. 我把坐标轴改成了您建议的，同时加上了最近两天的使用壁面函数（WMLES，y+在20+左右）的结果如下图。

因为算例设置的似乎对比性不大，只能看出用不用壁面函数确实像您说的对风压模拟影响不大。

2. 模拟出口的压力出口指的是？入口用的是师兄提出的改进的CDRFG，是经过通量修正的。

3. 对于LESy+1.2时，为保证最大库朗数=1，时间步长取得很小在0.000025左右，所以计算的很慢十天左右才一个结果

4. 几个结果确实都显示侧面漩涡脱落模拟不足，导致侧面脉动风压小，平均负压绝对值偏小。

5. 总的模拟时长在10s，取的是后5-7s统计的数据

然后我还有几个问题想请教下各位师兄

1. 这种侧面漩涡脱落模拟不足的主要原因是什么呢，网格、亚格子模型或是继续尝试别的离散格式

2. 采用pimpleFoam能否加大时间步长库朗数提高计算结果还是否准确？因为我看是有人在大库朗数下进行一些模拟，不知道对于风压模拟有没有采用过大库朗数计算的，但是我似乎库朗数一旦大于5就立刻发散停止？

3. 我采用的是nutUSalding壁面函数，壁面函数模化近壁面的速度分布（湍流粘性分布？）得到近壁面切应力，那壁面切应力和风压的关系是什么呢？也即这与风压模拟的准确与否有何关联？

4. 对于入口的通量修正像Kim, Y., Castro, I.P., Xie, Z.T., 2013. Divergence-free turbulence inflow conditions for large-eddy simulations with incompressible flow solvers. Comput. Fluids 84, 56–68.这篇文章提到的说基本不会改变入流的统计特性，但是我尝试了下无论是NSRFG还是CDRFG经过通流修正后都会引起入口统计参数改变尤其是顺风积分尺度和功率谱如下图，未通量修正前与目标值对比良好。不知道您有没有这种情况。
   
   猜测原因应该是生成的风速入口本身的通量时变太大导致达不到修正系数的这个精度只达到了1±0.1的精度！ Gungor, A.G., Sillero, J.A., Jimenez, J., 2012. Pressure statistics from direct simulation of turbulent boundary layer. Hawaii Seventh International Conference on Computational Fluid Dynamics (ICCFD7). (Accessed July 2012), 9-13.
   所以RFG系列方法本身生成时是不是没有考虑尽可能保证入口总通量守恒，以及零散度条件与通量守恒应该没关系吧？

5. 还有就是像这篇文章说的[Simulation of approaching boundary layer flow and wind loads on high-rise buildings by wall-modeled LES]相对于在入口二维平面生成零散度入口来说，保持入口通量守恒是不是重要一些？

不好意思问题可能有点多，盼复！

@coolhhh
感谢各位师兄回复，空计算域风场已经调整过应该还算对得上，但是最终算的结果尤其是脉动风压对不太上（以及侧面的平均风压），如上上图整体偏小，算例近壁面y+平均值在1.2左右，库朗数最大值始终在1左右。现在考虑可能的原因是对流项离散格式和亚格子模型的选取问题？或者说可能是因为TPU除了Iu剖面的其他统计量与我自己的假设不同导致的结果差异？另外一个问题想请教各位师兄的是因为网格密和时间不长小所以每次调整参数比较慢，我看有的文章会采用较大的时间步长比如0.001，但不知道当取大时间步长以及壁面函数比如前面师兄的nutUSaldingWallFunction结果尤其是脉动风压结果还能否计算准确呢？

@李东岳 李老师，我湍流模型用的是Smagorinsky，Cs参考一个师兄论文取得0.1

@HITSC30 您好，我想请问下您的这里的图示是否使用了Spalding壁面函数呢？我做了一个算例均值表现似乎还行但是脉动结果差的有点多，如下图。我没有使用壁面函数，但保证了近壁面y+约为1.2左右，亚格子模型选用的标准的SGS，对流项离散格式选取的linearupwind，据我的尝试似乎对流项格式会对结果影响较大？想问下您的对流项格式最终选取的是什么呢？以及想请教下哪些设置可能会对风压结果尤其是脉动风压结果影响较大呢？

@HITSC30 在 钝体建筑扰流的大涡模拟 中说：

@HITSC30

@一朵苍凉 您好，请问如何将压力云图离散化并使等值线作为其颜色分界线呢？

@coolhhh 好的好的，感谢回复！我去学习一下

@xzdbjx 您好！我看您使用过NSRFG方法导入过OF中计算，请问您在瞬态计算中有遇到过入口附近有非常大的虚拟压力波动的问题吗？我想验证下出现这种情况是否是算法的问题。我看您对湍流强度平均风剖面以及功率谱进行了验证，我这里湍流强度、平均风剖面以及功率谱这些与风速相关的统计量都还正常，就是涉及到脉动压强就会出问题。

@sunss 您好，感谢分享您的入口合成方法，我在使用您的程序生成的风速导入到空流场时似乎遇到了一些文章中提到的入口不满足divergence-free所导致的很大的虚拟压力波动（见后压力方差图），这会对放入建筑时建筑的表面脉动风载产生很大的影响，我看很多采用NSRFG的文章都没有提到这个事情，不知道您有没有遇到这种情况以及是如何解决的呢？

@HITSC30 好的，感谢！

@HITSC30 感谢回复！没有看到上文有网盘连接嘞，只看到了您最初的设置，请问您最终是如何解决这一问题了呢？设置了更小的时间步长吗，还是更换了网格，以及有没有出现十分大的压力波动现象呢？