涡街模拟建议.docx

1、涡街模拟建议1我计算的是一个二维自维持振荡问题（好多文献都如此说），我采纳层流算法也取得了类似的结果，k-e-rng也能够。而别的模型都不行了，一样都是最后取得一个稳固的解（和文献上说的不同）。因为雷诺数比较小，不能确信何时转变成紊流，因此想用一个能够计算过渡流动的模型。不明白用k-e-rng模型是不是就能够够说是准确，因为没有具体的实验数据，是不是能够依照它的计算流场和实验流场相似就确信计算的合理性和准确性呢？多谢多谢Hi-key：这种问题的要求比较高，类似的问题我只算过卡门涡阶的。只是那时是用层流算得。你那个例子里面若是跟湍流模型灵敏，我建议你能够尝试以下方式：FzN/5选用其他湍流模型，

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

3、明白是不是你所说的湍流强度和粘性比）计算了一下，发觉对结果阻碍专门大，都为零时，没有振荡现象；增大这两个值又会显现不同的流场。绝对压力只是为了使计算加倍准确，你也能够用表压计算。Re!3湍流的两个变量是入口处的脉动情形，都为0时跟层流接近可是跟层流不同。你将湍流强度设为5%，粘性比设为。再碰运气，有转变的话，换别的湍流模型再试下。UH另外周期边界中你设定压降仍是流量？流量的指定方式加倍容易显现波动。pMy我设定的是流量。适才我转变了这两个湍流变量，阻碍比较大。我这段时刻一直在实验不同的湍流模型，以找到比较适合的。我试过别的湍流模型，最后都是趋于稳固的。另外，我选用层流和rng模型时，最后的振荡

4、大体是一致的（如图，绿色：rng，白色：层流）。我计算的进程是，先用稳态算法算出一个初场，然后改用非稳态算法继续计算（有人说过如此的方式），不明白那个是不是有什么阻碍？如此计算时，非稳态计算的初场是不同的。换成别的湍流模型时，即利用稳态取得一个类似于r ng的计算结果，然后改用非稳态算法，最后都是稳固的。上下对称就可不能振荡，不对称才会振荡。其他湍流模型的初场都类似于那个不对称流场，可是最终都是稳固了，即上下对称。nUvIv依照我的体会，你如此算：P&选用层流模型，通过修改流体的物性来凑雷诺数，然后把流速条的小一些（雷诺数不变），然后计算就能够够了，所有的松弛因子调到1，持续性方程，用simp

5、lc算法，离散格式用二阶，应该就没错了。vHSpPS上面那个图是算淹没射流同志的case，不明白是不是他想要得结果。有些文献是如此计算的，通过改变粘性来改变雷诺数。可是我此刻算的那个东西比较具体了，应该说工质是空气就可不能变了。我是改变入口流量来调剂雷诺数，松弛因子持续方程取为，其他都一样，采纳湍流模型，湍流的松弛因子只要不是1就能够够。离散格式换了几个，大体相同。一个问题是，雷诺数大了以后，层流和rng的计算结果就不一样了，层流仍是规律性振荡，而rng的就不太规律了（和文献上类似）。因为我的计算要涉及到层流流动和过渡流动，可是又不明白何时转化，因此想从一开始就利用湍流，以监视流场的转变。D=

6、1 V=1 为了不同的雷诺数，因此在Fluent中改变材料动力粘性系数， 可是计算的结果完全相同，乃至于雷诺数到的时候都会显现涡（层流，定常） 而且不同的雷诺数连收敛的步数都是一样的，奇怪哉！模型是不是有问题呢？依照道理应该不一样的，可能的缘故1.圆柱饶流在re100多时就有卡门涡街了，应该采纳非定常计算方式3900是亚临界雷诺数，尾迹转捩，应该采纳湍流模型计算3.网格问题，通常网格在近壁面周围要划分较细 呵呵，先就谈这些吧，你给的信息太少了在Fluent中改变材料动力粘性系数: 检查一下是不是再fluid里面改了列位大侠，看到很多文献中提到涡脱落周期T，而且是保留的1/4T、1/2T等各个不

7、同时刻的文件，进行涡脱落分析；我不明白脱落周期如何确信，和如何保留精准的1/4T时刻； 还有我可否在计算完圆柱绕流后，与文献对照确信计算的准确性，然后将所采纳的模型应用到其他形状管子的模拟上；有什么理论指导吗？能够推荐一本好书给我。超级感激了！周期一样是通过FFT求出的，需要某点的压力或壁面的升力系数。通过计算假设干个周期，再用FFT就能够够求出涡脱落频率了。明白频率也就明白周期了。然后将1/4T,1/2T的文件保留下来，能够看看涡，速度等图片的转变情形。不明白FLUENT行不行，具体该怎么做？FFT是什么？傅立叶变换FFT是快速傅立叶变换，fluent里面应该没有对某点数值的时程图进行傅立叶

8、变换吧（matlab里面有那个功能），只是我想，若是能把某点的数值的时程图提掏出来，不用进行FFT，应该也不难判定周期了吧。最近正在算一个圆柱绕流问题，Re=200，velocity=1 m/s,desity=1kg/m*3ot)diameter=1m , viscosity= kg/,$%n尽管后面有涡街，可是升力系数和阻力系数和实验的或他人计算的对不上。D0请问列位大虾这是怎么回事？望高手不吝赐教！yi*x（Report里的reference value 我取的是：depth=1m,velocity=1m/s,ZK2area=1 m*2,density=1kg/m*3)我感觉referen

9、ce value中的area=1m*1,因为你的diameter=1m，那么结果就应该是2倍或是1/2的关系，可能你的结果没错，只是个放大倍数的问题！我也有类似的问题,不同的是,我的阻力系数不是一个大体恒定的值,而是随时刻震荡转变的值,请问这是什么缘故啊?u, 此刻在用lbm做方柱绕流，采纳郭照立提出的不可压模型，边界条件也是他的非平稳外推方式。-%!G 计算中粘度系数用以下表达式jKh niu=(2*tau-1)/6*(c*dx) cj 方柱绕流的Re数概念为：syPcA9 Re=U*L/niu，计算中U1，L1（U为来流速度，L为方柱边长）v-y)K 计算区域为： 0=x=60L 0=y1

10、00,算不出来Votex Shedding现象不明白哪位大侠有高招，帮帮忙呀V0我用FLUNT算一个立方体绕流问题，立方体为600mm*600mm*200mm，利用对称性，取600mm300mm200mm，计算域取10200mm4500mm1000mm，划分了180多万个四面体网格。网格用GAMBIT生成。b6L导入FLUENT中，采纳层流模型，设置了入口的速度边界条件，如X=12m/s，进行计算。结果发感觉出的结构表面的压力系数的数量级为10的二次方，这显然是不对的。不知可能在那里显现问题。可否有人指点一下。万分感激。大伙儿好，在看圆柱绕流的文章时，看到一个对称边界条件，不太明白得，那个对称

11、边界条件是如何概念的，有什么物理意义呢？&A,j请不吝赐教。谢谢。对称边界条件很常见的.+W顾名思义, 确实是流场中的一个对称面, 面两边的流动型态是完全对称的. 在求解的时候为了节约计算量, 常常只解一半, 而把对称面设为对称边界条件.TG(t由于对称性, 对称边界条件有以下特点:,D+1. 穿过对称面的通量(法向速度)为0%/2. 对称面法线方向的导数为0b2a3. 不同于壁面, 对称边界条件允许有切向速度和切线方向的转变.3楼的回答的正确，可是在圆柱绕流中不能用对称边界，因为若是Re较大的话，会产生涡街，不符合三个条件。p大伙儿好，圆柱绕流数值模拟的时候，时刻步长如何选择适合呢？有个教师

12、告知我通常取无量纲的时刻步长为，什么缘故如此选择？时刻步长小应该来讲计算结果更精准，可是不是由此会致使收敛速度慢，容易发散？l在圆柱绕流的数值模拟中，除低雷诺数绕流，哪些情形会致使无漩涡脱落、无分离流动现象显现？望大虾不吝赐教，以助后学者进步，谢谢。时刻步长小, 每一个步长收敛速度会变快，也更不容易发散. 但共计算量可能变大. 他人算一个周期用10步, 你用100步固然计算时刻可能变长了.9=步长太长不行, 可能致使无漩涡脱落, 或看不到脱落进程.eDaDbP一个周期至少20个时刻步吧. 18x时刻步长过小也没必要. 简单看看漩涡脱落,一个周期100个时刻步足够了._x&s上次我算着玩的时候,

13、 不是用了20确实是40个时刻步. 漩涡脱落没问题.要想算出漩涡脱落一样要用高阶精度的格式.r时刻步长小, 每一个步长收敛速度会变快，也更不容易发散F7c可是我用很小的时刻步，每时刻步内的结果却不容易收敛，最后尚未涡脱，不知怎么回事。自己编的程序仍是用商业的?OKJ自己编的话就很难说了.在他人的基础上改编的，呵呵sk这事实上犹如实验的采样频率，为了扑作到某一个绕流特点，每一个周期内必需搜集必然数量的样本，如：无量纲时刻步长对应100次，不同的绕流特点需要的采样频率有必然的不同，因此有一说。，9pWh2If there is vortex shedding does not only depen

14、d 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 级别的雷诺数的绕流 完全不能用层流分析么么?咱们在比较高的流速时,取得了脱表现象, x!4f在圆柱边界有一层反向流动,但后面的涡kU6不是交替释放的啊, 也疑心是不是是湍流破坏了涡结构.+边界层太细了仿佛也不行，加我QQ把，我也作者方面的工作；啊啊；8#什么缘故用层流的模型能够取得周期性的涡旋释放呢,d在你给定

15、的雷诺数范围内，流动是非定常，层流。用湍流模型算层流一样不行。kE一、说明型udf要求函数中不能包括子函数，有可能使那个缘故。ZxXp二、动网格求解中的锯齿形残差是正常的。YNO3、计算进程中，某几个单元温度越界没有关系，只要它能回来，别发散就行，那个缘故较多，一样来讲不是网格的问题，比如在求解超音速流动时，利用segregate算法时就容易显现。K4、分离求解器能够用来求解可压超音速流，但要求结构网格同时网格间距跳跃不大，马赫数小于时，将松弛因子全数设为，是能够计算的。马赫数高了就完蛋了。还有计算结果跟couple的一样。Gt这些是我的体会，也不必然十分准确，仅供参考?3我此刻算的是弹丸出炮

16、口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每一个时刻步很难收敛,我就换成份离的,而且我亚松弛因子调的专门大,压力有,动量为,其他都是1.我用的是三角形的非结构网格,气流速度有1500多,马赫数达到了4左右,可是在每一个时刻步上它收敛专门快,一样16步左右就收敛了.想问下你,象我如此设那么大的亚松弛因子能够吗?还有就想问你,象我如此的非稳态怎么看每一个时刻步是不是收敛了?我此刻确实是看残差图7nih)流体中文网论坛 - 流体力学及相关领域学术问题交流论坛。=!直觉不太对，只是若是你的时刻步长很小或许就能够够避过每一个时刻步长内发散的可能性，若是时刻步长比较大就不对了哈，另外炮弹

17、出口若是直线运动用layering，比remesh好多了，又准确又快。APS，膛内外压力差多大？Mpa级？U:%是啊,有60多mpa.我时刻步长为还要请教楼主,想我如此的非稳态计算怎么样看每一步是不是收敛了?看残差吗?能够采纳稳态算法中的方法，监控某点的某个变量值，要在流场转变区域内，如喷口位置。若是平了就好了。4我计算车室内流场，出风口很多，回风口（排风口）只有一个，采纳hex core网格、RNG k-e 模型，计算进程监视出口（回风口）流量，计算到300步的时候差不多成了一条直线，而残差图也差不多成了直线，接着迭代300步，残差图仍是下不去，也确实是不能降低到低于1e-3的水平，在rep

18、ort里面看net imbalance 在%左右，请问如此能不能算收敛了？Z#i压力我用的standard ，其他的是二阶迎风，有什么方法能使残差降到1e-3以下，谢谢?bN有时候残差下不来是正常的，尤其是持续性方程，一样说来x,y,z速度收敛了，只剩持续性方程没收敛，其他监控变量已水平，能够以为计算收敛（包括传热时除外）。也确实是说你那个问题能够算收敛了，要想将持续方程残差降下来要改网格的，或改变松弛因子造成假收敛，固然此法不可取。5wormxieu你能够看下的第19个例子，是个欧拉模型的例子。能够参考一下大伙儿好，圆柱绕流数值模拟的时候，时刻步长如何选择适合呢？有个教师告知我通常取无量纲的

19、时刻步长为，什么缘故如此选择？时刻步长小应该来讲计算结果更精准，可是不是由此会致使收敛速度慢，容易发散？在圆柱绕流的数值模拟中，除低雷诺数绕流，哪些情形会致使无漩涡脱落、无分离流动现象显现？望大虾不吝赐教，以助后学者进步，谢谢。时刻步长小, 每一个步长收敛速度会变快，也更不容易发散. 但共计算量可能变大. 他人算一个周期用10步, 你用100步固然计算时刻可能变长了.步长太长不行, 可能致使无漩涡脱落, 或看不到脱落进程.一个周期至少20个时刻步吧. 时刻步长过小也没必要. 简单看看漩涡脱落,一个周期100个时刻步足够了.上次我算着玩的时候, 不是用了20确实是40个时刻步. 漩涡脱落没问题.

20、要想算出漩涡脱落一样要用高阶精度的格式.自己编的程序仍是用商业的?自己编的话就很难说了.我上面提到的那位教师说选取无量纲的时刻步为，我看了很多圆柱绕流的数值模拟的文章，也大体如此，什么缘故会是如此呢？anyone说，“时刻步长过小也没必要. 简单看看漩涡脱落,一个周期100个时刻步足够了，上次我算着玩的时候, 不是用了20确实是40个时刻步. 漩涡脱落没问题.”如此做，漩涡脱落是没问题了吧，可是计算绕流阻力会可不能有较大阻碍呢？这事实上犹如实验的采样频率，为了扑作到某一个绕流特点，每一个周期内必需搜集必然数量的样本，如：无量纲时刻步长对应100次，不同的绕流特点需要的采样频率有必然的不同，因此

21、有一说。， 此刻在用lbm做方柱绕流，采纳郭照立提出的不可压模型，边界条件也是他的非平稳外推方式。 计算中粘度系数用以下表达式 niu=(2*tau-1)/6*(c*dx) 方柱绕流的Re数概念为： Re=U*L/niu，计算中U1，L1（U为来流速度，L为方柱边长） 计算区域为： 0=x=60L 0=y100,算不出来Votex Shedding现象 以下内容由 wsl2006 在 2006年12月07日 07:11pm 时添加 不明白哪位大侠有高招，帮帮忙呀想做一个绕圆柱产生卡门涡街的动画，但在求解进程中老是取得对称的流场。明白的人都知道那个确信不对，我想问问高手什么地址需要注意，可能是我什么地址设置错了!1Ks我用的是ideal-gasF;jk-e湍流模型=边界都是远场边界条件，来流速Ma= 表压1atm-O还有什么地址可能错了？n你是不是在流场初始化的时候也是对称的哦？gB0你试着初始化流场时不对称，看看4 可是卡门涡街本身没有要求非对称来流亚？Dk 我说的是流场初始化，不能是对称的，但能够是对称的来流，你能够看一下李万平那本书$e