CFD中文网

@yxg2018 谢谢，我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域，计算依旧会报错，提示Divergence detected in AMG solver: x-momentum，我将迭代次数增加并改小了动量的松弛因子，计算还是发散。

@yxg2018 我现在动网格是采用弹簧更新和网格重构，将圆柱边界设置为刚体运动

各位老师，我在做亚临界雷诺数下二维圆柱的涡激振动（y向），目前模拟的工况是Re=5000，圆柱直径0.0381m，速度0.132m/s，约化速度为8，质量比2.6，阻尼因子0.0036，振动方程通过编写四阶龙格库塔的UDF导入到fluent中，振动方程进行了无量纲化。

模型如下，上游10D，下游20D，宽度上下各15D，圆柱附近采用四边形网格，其余为三角形网格。第一层网格高度0.0023D，时间步长0.0014s，湍流模型k-wSST。

在低雷诺数时基于编好的UDF进行过计算（此时的模型不一样，直径0.0016m），与实验数据接近，但在亚临界时，计算开始一段时间后圆柱位移一直增大，直到出现了负体积，此时的位移已接近10D。自己更改过时间步长，但是问题也是一样，想问下有老师做过类似的模拟吗？这个问题排查了很久还是没能解决，谢谢大家！
编写的UDF如下。

@random_ran 好的谢谢老师

@random_ran 老师可以帮忙看看吗？我后面也换了k-omega SST模型，结果也相差不多，看了一下也并没有漩涡脱落产生。

谢谢东岳老师

大家好，我最近在模拟Re=10000时，二维圆柱的升力系数和阻力系数。采用的湍流模型为k-omega，网格划分采用非结构化网格，边界层采用四边形网格，圆柱壁面y+值在1附近。换过不同的网格，时间步长等，阻力系数为1.0-1.2左右，与文献可以对应上，但是升力系数模拟出来在1.25左右，文献数据为0.4-0.6附近。大家能帮着看看问题出在哪里嘛？fluent里再选择k-omega模型后还可以选择低雷诺数修正或者自由剪切修正，我在计算的时候并没有勾选，大家知道这个一般怎么用吗？

各位老师好，我在画圆柱扰流网格的时候想在圆柱周围加边界层，再在一定范围内加密，该怎么实现呢？我现在试了两种方法。

1.分成了两个区域，但是由于后面圆柱需要运动，当运动范围超过紫色的部分时，就会报错。两个区域中间使用了interior，静止时没影响，但网格运动应该是全局网格，如果有interior其实是不正确的。

2.采用了密度盒，但是该方法要求网格划分采用Patch Independent，这样就不会遵守我设置的全局网格参数，而且设置的圆柱边界层网格也生成不了，直接全部都是三角形网格。

想请问各位老师，应该如何实现呢？

或者有老师知道如何将ICEM的模型输出为.stp格式吗？

在ICEM中建立模型如下图

通过两个圆面和圆柱侧面形成一个名为soild的part。ICEM生成.stl文件后导入到Design Modeler中。

但是在DM中，圆面和圆柱侧面均变成了属性为soild的part。
请问在ICEM中如何设置part可以使得圆柱的三个面在DM中以一个soild形式生成？

@东岳 嗯嗯，确实是尺度问题，谢谢东岳老师！

@ibelief 确实是scale的问题，谢谢老师！

各位老师，我在用fluent模拟圆柱二维扰流时，采用瞬态求解器，雷诺数200，但是无法看到漩涡脱落现象，加大雷诺数后也没有变化，只能在圆柱周围看到少量的涡量线。请问老师们问题可能出现在什么地方呢？网格和相关设置都是跟着书里算例走的，所以不清楚问题在哪里，谢谢!

大家好，我用OpenFOAM计算三维风力机的转矩和转矩系数时，结果和文献相差大约10倍，风速10m/s，转速31.4rad/s，风机叶片高度3m，旋转轴高4m，风轮半径1m。计算转矩时对整个风机求解。大家可以帮忙看看问题出在哪里吗？谢谢！ （用gambit划分网格，圆柱区域是滑移区。） （fvschemes文件） （controldict文件）

@dzw05 湍流模型中有一个运动黏度nu，正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值，由于之后的计算我都用的是这个错误的湍流模型文件，导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。

@dzw05 老师我排查出来是我的湍流模型文件里面的参数设置出了问题，所以导致阻力系数一直存在问题，现在已经解决了。十分感谢老师这几天的帮助！

@dzw05 老师，我修改了边界条件，对网格也进行了一些修改，但是重新计算后阻力系数还是基本没变。会不会是采用的SA模型本身对于这种情况模拟就不是很好？

@东岳 老师您说的分辨率指的是gambit中划分网格时这个方向的successive ratio吗？

@dzw05 老师，我这两天重新划分网格的时候发现之前计算过程中郎库数非常大，现在我根据郎库数为1重新设定了时间步长，但阻力系数依旧不正确。我在gambit网格质量检查的时候，远场质量为0,翼型附近质量为0.1-0.2。不知道这样的网格质量可以吗，或者还有什么其他的可能问题吗（比如说我的nut，nutlida值会不会设定有问题），因为我现在修改了很多次网格，但是阻力系数总是维持在10以上。谢谢老师！

@dzw05 老师您好，我的第一层网格高度大约为0.18mm，昨天跟您说错了。刚才重新划分网格后，checkMesh显示正常，但是进行计算后升力系数和阻力系数都与实验数据有偏差。老师我把我的网格文件，controlDict文件和初始条件传上来，您能帮我看看问题出在哪里吗？谢谢老师！4_1523608055301_U.pdf 3_1523608055301_p.pdf 2_1523608055301_nuTilda.pdf 1_1523608055301_nut.pdf 0_1523608055300_controlDict.pdf [网格](链接: https://pan.baidu.com/s/129VWe7pswJfTMWmG4oHxrg 密码: x7g4)