涡街模拟建议

1、1我计算的是一个二维自维持振荡问题（好多文献都这样说），我采用层流算法也得到了类似的结果，k-e-rng 也可以。而别的模型都不行了，一般都是最后得到一个稳定的解（和文献上说的不同） 。因为雷诺数比较小， 不能确定什么时候转变成紊流，所以想用一个能够计算过渡流动的模型。不知道用 k-e-rng 模型是不是就可以说是准确，因为没有具体的试验数据，是不是可以根据它的计算流场和试验流场相似就确定计算的合理性和准确性呢？多谢多谢hi-key：这种问题的要求比较高， 类似的问题我只算过卡门涡阶的。不过当时是用层流算得。你这个例子里面如果跟湍流模型敏感，我建议你可以尝试以下方法：fzn/5选用其他湍流模型

2、，然后在进出口边界处的湍流相设置时，使用湍流强度和粘性比然后将这两个数值全部给0，再计算时使用绝对压力计算。也许会有变化，也许没用，你可以试试，把结果告诉我。谢谢，另外在所有的湍流模型中k-e-rng 是最适合计算低雷诺数湍流模型的，当然你也可以尝试真正的低雷诺数湍流模型，需要在用户界面中输入命令行激活，至于怎么激活我忘了，哈哈，不好意思，等我查到了给你哈。 f.!z? 流体中文网论坛- 流体力学及相关领域学术问题交流论坛。/.另外判断结果是否正确只能靠实验或者查文献了，流态特征相似只能证明大体上没错，但是精度就不知道了。我计算的是周期性边界条件，和绝对压力有关系吗？我刚才改变了初场的两个湍流

3、变量（不知道是不是你所说的湍流强度和粘性比）计算了一下，发现对结果影响很大，都为零时，没有振荡现象；增大这两个值又会出现不同的流场。绝对压力只是为了使计算更加准确，你也可以用表压计算。re!3湍流的两个变量是入口处的脉动情况，都为0 时跟层流接近但是跟层流不同。你将湍流强度设为5%，粘性比设为 0.01。再试试看，有变化的话，换别的湍流模型再试下。uh 另外周期边界中你设定压降还是流量？流量的指定方式更加容易出现波动。pmy我设定的是流量。刚才我变化了这两个湍流变量，影响比较大。我这段时间一直在试验不同的湍流模型，以找到比较适合的。我试过别的湍流模型，最后都是趋于稳定的。另外，我选用层流和rn

4、g 模型时，最后的振荡基本是一致的（如图，绿色：rng，白色：层流）。我计算的过程是，先用稳态算法算出一个初场，然后改用非稳态算法继续计算 （有人说过这样的方法） ，不知道这个是不是有什么影响？这样计算时，非稳态计算的初场是不同的。换成别的湍流模型时，即使用稳态得到一个类似于r ng 的计算结果，然后改用非稳态算法，最后都是稳定的。上下对称就不会振荡，不对称才会振荡。其他湍流模型的初场都类似于这个不对称流场，但是最终都是稳定了，即上下对称。nuviv按照我的经验，你这样算：p&选用层流模型，通过修改流体的物性来凑雷诺数，然后把流速条的小一些（雷诺数不变） ，然后计算就可以了，所有的松弛

5、因子调到1，连续性方程0.9，用 simplc 算法，离散格式用二阶，应该就没错了。vhspps上面那个图是算淹没射流同志的case ，不知道是不是他想要得结果。有些文献是这样计算的，通过改变粘性来改变雷诺数。但是我现在算的这个东西比较具体了，应该说工质是空气就不会变了。我是改变进口流量来调节雷诺数，松弛因子连续方程取为0.7，其他都一样，采用湍流模型，湍流的松弛因子只要不是1就可以。离散格式换了几个，基本相同。一个问题是，雷诺数大了之后，层流和rng的计算结果就不一样了，层流还是规律性振荡，而rng 的就不太规律了（和文献上类似） 。因为我的计算要涉及到层流流动和过渡流动，但是又不知道什么时

6、候转化，所以想从一开始就使用湍流，以监视流场的变化。d=1 v=1 为了不同的雷诺数，所以在fluent 中改变材料动力粘性系数，但是计算的结果完全相同，甚至于雷诺数到0.5 的时候都会出现涡（层流，定常）而且不同的雷诺数连收敛的步数都是一样的，奇怪哉！模型是不是有问题呢？按照道理应该不一样的，可能的原因1.圆柱饶流在re100多时就有卡门涡街了，应该采用非定常计算方法2.re3900是亚临界雷诺数，尾迹转捩，应该采用湍流模型计算3.网格问题，通常网格在近壁面附近要划分较细呵呵，先就谈这些吧，你给的信息太少了在 fluent 中改变材料动力粘性系数: 检查一下是否再fluid 里面改了各位大侠

7、，看到很多文献中提到涡脱落周期t，而且是保存的1/4t 、1/2t 等各个不同时刻的文件，进行涡脱落分析；我不知道脱落周期如何确定，以及如何保存精确的1/4t 时刻；还有我可否在计算完圆柱绕流后，与文献对比确定计算的准确性，然后将所采用的模型应用到其他形状管子的模拟上；有什么理论指导吗？可以推荐一本好书给我。非常感谢了！周期一般是通过fft 求出的，需要某点的压力或壁面的升力系数。通过计算若干个周期，再用fft 就可以求出涡脱落频率了。知道频率也就知道周期了。然后将1/4t,1/2t 的文件保存下来，可以看看涡，速度等图片的变化情况。不知道 fluent 行不行，具体该怎么做？fft 是什么？

8、傅立叶变换fft 是快速傅立叶变换，fluent 里面应该没有对某点数值的时程图进行傅立叶变换吧（ matlab 里面有这个功能） ，不过我想，要是能把某点的数值的时程图提取出来，不用进行fft，应该也不难判断周期了吧。最近正在算一个圆柱绕流问题，re=200，velocity=1 m/s,desity=1kg/m*3 ot)diameter=1m , viscosity=0.005 kg/m.s,$%n虽然后面有涡街，但是升力系数和阻力系数和实验的或者别人计算的对不上。d0请问各位大虾这是怎么回事？望高手不吝赐教！yi*x（report里的 reference value 我取的是： dep

9、th=1m,velocity=1m/s,zk2area=1 m*2,density=1kg/m*3) 我觉得 reference value中的 area=1m*1,因为你的 diameter=1m，那么结果就应该是2倍或者是 1/2 的关系，可能你的结果没错，只是个放大倍数的问题！我也有类似的问题 ,不同的是 ,我的阻力系数不是一个基本恒定的值,而是随时间震荡变化的值 ,请问这是为什么啊 ? u,现在在用 lbm 做方柱绕流，采用郭照立提出的不可压模型，边界条件也是他的非平衡外推方法。 -%!g计算中粘度系数用以下表达式jkhniu=(2*tau-1)/6*(c*dx) cj方柱绕流的 re

10、数定义为： sypca9re=u*l/niu ，计算中 u1，l1（u 为来流速度， l 为方柱边长） v-y)k计算区域为： 0=x=60l 0=y100,算不出来 votex shedding现象不知道哪位大侠有高招，帮帮忙呀v0我用 flunt 算一个立方体绕流问题，立方体为600mm*600mm*200mm ，利用对称性，取 600mm300mm200mm，计算域取10200mm4500mm1000mm，划分了 180 多万个四面体网格。网格用gambit 生成。 b6l导入 fluent 中，采用层流模型，设置了入口的速度边界条件，如x=12m/s，进行计算。结果发现得出的结构表面的

11、压力系数的数量级为10的二次方， 这显然是不对的。不知可能在那里出现问题。可否有人指点一下。万分感谢。大家好，在看圆柱绕流的文章时，看到一个对称边界条件，不太理解，这个对称边界条件是如何定义的，有什么物理意义呢？&a,j请不吝赐教。谢谢。对称边界条件很常见的.+w顾名思义 , 就是流场中的一个对称面, 面两边的流动型态是完全对称的. 在求解的时候为了节约计算量, 经常只解一半 , 而把对称面设为对称边界条件.tg(t由于对称性 , 对称边界条件有以下特点:,d+1. 穿过对称面的通量 (法向速度 )为 0%/2. 对称面法线方向的导数为0b2a3. 不同于壁面 , 对称边界条件容许有切

12、向速度和切线方向的变化. 3 楼的回答的正确，但是在圆柱绕流中不能用对称边界，因为要是re较大的话，会产生涡街，不符合三个条件。p大家好，圆柱绕流数值模拟的时候，时间步长如何选择合适呢？有个老师告诉我通常取无量纲的时间步长为0.002-0.01，为什么这样选择？时间步长小应该来说计算结果更精确，但是不是由此会导致收敛速度慢，容易发散？l 在圆柱绕流的数值模拟中，除了低雷诺数绕流， 哪些情况会导致无漩涡脱落、无分离流动现象出现？ 望大虾不吝赐教，以助后学者进步，谢谢。时间步长小 , 每个步长收敛速度会变快，也更不容易发散 . 但总计算量可能变大 . 别人算一个周期用10 步, 你用 100 步当

13、然计算时间可能变长了.9=步长太长不行 , 可能导致无漩涡脱落 , 或看不到脱落过程 .edadbp一个周期至少20 个时间步吧 . 18x时间步长太小也没必要. 简单看看漩涡脱落 ,一个周期 100 个时间步足够了 ._x&s上次我算着玩的时候, 不是用了 20 就是 40 个时间步 . 漩涡脱落没问题 .要想算出漩涡脱落一般要用高阶精度的格式.r 时间步长小 , 每个步长收敛速度会变快，也更不容易发散f7c可是我用很小的时间步，每时间步内的结果却不容易收敛，最后还没有涡脱，不知怎么回事。自己编的程序还是用商业的?okj自己编的话就很难说了. 在别人的基础上改编的，呵呵sk 这实际上

14、如同实验的采样频率，为了扑作到某一个绕流特征，每一个周期内必须采集一定数量的样本，如：无量纲时间步长0.01 对应 100 次，不同的绕流特征需要的采样频率有一定的差别，因此有0.010.002 一说。 ，9pwh2if there is vortex shedding does not only depend on the numerical parameters. for the flow over circular cilinder, you need to know how large the reynolds number is. if re0 时速度为 u0+u,z另, 10e+5

15、 级别的雷诺数的绕流完全不能用层流分析么么? 我们在比较高的流速时,得到了脱体现象 , x!4f在圆柱边界有一层反向流动,但后面的涡 ku6不是交替释放的啊 , 也怀疑是否是湍流破坏了涡结构.+边界层太细了好像也不好，加我qq 把，我也作者方面的工作；啊啊； ；279536328#为什么用层流的模型可以得到周期性的涡旋释放呢,d在你给定的雷诺数范围内，流动是非定常，层流。用湍流模型算层流一般不行。ke1、解释型 udf 要求函数中不能包含子函数，有可能使这个原因。zxxp2、动网格求解中的锯齿形残差是正常的。yno3、计算过程中，某几个单元温度越界没有关系，只要它能回来，别发散就行，这个原因较

16、多，一般来说不是网格的问题，比如在求解超音速流动时，使用segregate算法时就容易出现。 k4、 分离求解器可以用来求解可压超音速流，但要求结构网格同时网格间距跳跃不大，马赫数小于 2.5 时，将松弛因子全部设为0.1，是可以计算的。 马赫数高了就完蛋了。还有计算结果跟couple 的一样。 gt这些是我的经验，也不一定十分准确，仅供参考?3我现在算的是弹丸出炮口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每个时间步很难收敛,我就换成分离的 ,而且我亚松弛因子调的很大,压力有 0.9,动量为 0.8,其他都是 1.我用的是三角形的非结构网格,气流速度有 1500多,马赫数达到

17、了 4 左右,可是在每个时间步上它收敛很快,一般 16步左右就收敛了 .想问下你 ,象我这样设那么大的亚松弛因子可以吗?还有就想问你 ,象我这样的非稳态怎么看每一个时间步是否收敛了?我现在就是看残差图7nih)? 流体中文网论坛- 流体力学及相关领域学术问题交流论坛。=!直觉不太对，不过如果你的时间步长很小也许就可以避过每个时间步长内发散的可能性，如果时间步长比较大就不对了哈，另外炮弹出口如果直线运动用layering，比remesh好多了，又准确又快。aps，膛内外压力差多大？mpa 级？u:% 是啊,有 60 多 mpa.我时间步长为0.000001s4 还要请教楼主 ,想我这样的非稳态计

18、算怎么样看每一步是否收敛了?看残差吗 ? 可以采用稳态算法中的办法，监控某点的某个变量值，要在流场变化区域内，如喷口位置。如果平了就好了。4我计算车室内流场，出风口很多，回风口（排风口）只有一个，采用hex core网格、rng k-e 模型，计算过程监视出口（回风口）流量，计算到300 步的时候差不多成了一条直线，而残差图也差不多成了直线，接着迭代300 步，残差图还是下不去，也就是不能降低到低于1e-3的水平，在 report 里面看 net imbalance 在 0.1%左右，请问这样能不能算收敛了？z#i压力我用的 standard ，其他的是二阶迎风，有什么办法能使残差降到1e-3

19、 以下，谢谢?bn有时候残差下不来是正常的，尤其是连续性方程，一般说来x,y,z速度收敛了，只剩连续性方程没收敛，其他监控变量已水平，可以认为计算收敛（包含传热时除外） 。也就是说你这个问题可以算收敛了，要想将连续方程残差降下来要改网格的，或者改变松弛因子造成假收敛，当然此法不可取。5wormxie u你可以看下 fluent6.1 的第 19 个例子，是个欧拉模型的例子。可以参考一下大家好， 圆柱绕流数值模拟的时候，时间步长如何选择合适呢？有个老师告诉我通常取无量纲的时间步长为0.002-0.01，为什么这样选择？时间步长小应该来说计算结果更精确，但是不是由此会导致收敛速度慢，容易发散？在圆

20、柱绕流的数值模拟中，除了低雷诺数绕流，哪些情况会导致无漩涡脱落、无分离流动现象出现？望大虾不吝赐教，以助后学者进步，谢谢。时间步长小 , 每个步长收敛速度会变快，也更不容易发散. 但总计算量可能变大. 别人算一个周期用10 步, 你用 100 步当然计算时间可能变长了. 步长太长不行 , 可能导致无漩涡脱落, 或看不到脱落过程. 一个周期至少20 个时间步吧 . 时间步长太小也没必要. 简单看看漩涡脱落,一个周期100 个时间步足够了. 上次我算着玩的时候, 不是用了20 就是 40 个时间步 . 漩涡脱落没问题. 要想算出漩涡脱落一般要用高阶精度的格式. 自己编的程序还是用商业的? 自己编的话就很难说了. 我上面提到的那位老师说选取无量纲的时间步为0.002-0.01，我看了很多圆柱绕流的数值模拟的文章，也基本如此，为什么会是这样呢？anyone 说， “时间步长太小也没必要. 简单看看漩涡脱落,一个周期 10