涡街模拟建议.docx
- 文档编号:16495480
- 上传时间:2023-07-14
- 格式:DOCX
- 页数:14
- 大小:26.76KB
涡街模拟建议.docx
《涡街模拟建议.docx》由会员分享,可在线阅读,更多相关《涡街模拟建议.docx(14页珍藏版)》请在冰点文库上搜索。
涡街模拟建议
涡街模拟建议
流体模拟
1.我计算的是一个二维自维持振荡问题(好多文献都这样说),我采用层流算法也得到了类似的结果,k-e-rng也可以。
而别的模型都不行了,一般都是最后得到一个稳定的解(和文献上说的不同)。
因为雷诺数比较小,不能确定什么时候转变成紊流,所以想用一个能够计算过渡流动的模型。
不知道用k-e-rng模型是不是就可以说是准确,因为没有具体的试验数据,是不是可以根据它的计算流场和试验流场相似就确定计算的合理性和准确性呢?
多谢多谢
Hi-key:
这种问题的要求比较高,类似的问题我只算过卡门涡阶的。
不过当时是用层流算得。
你这个例子里面如果跟湍流模型敏感,我建议你可以尝试以下方法:
FzN/5[
选用其他湍流模型,然后在进出口边界处的湍流相设置时,使用湍流强度和粘性比然后将这两个数值全部给0,再计算时使用绝对压力计算。
也许会有变化,也许没用,你可以试试,把结果告诉我。
谢谢,另外在所有的湍流模型中k-e-rng是最适合计算低雷诺数湍流模型的,当然你也可以尝试真正的低雷诺数湍流模型,需要在用户界面中输入命令行激活,至于怎么激活我忘了,哈哈,不好意思,等我查到了给你哈。
f.!
Z流体中文网论坛--流体力学及相关领域学术问题交流论坛。
/.另外判断结果是否正确只能靠实验或者查文献了,流态特征相似只能证明大体上没错,但是精度就不知道了。
我计算的是周期性边界条件,和绝对压力有关系吗?
我刚才改变了初场的两个湍流变量(不知道是不是你所说的湍流强度和粘性比)计算了一下,发现对结果影响很大,都为零时,没有振荡现象;增大这两个值又会出现不同的流场。
绝对压力只是为了使计算更加准确,你也可以用表压计算。
Re\!
3
湍流的两个变量是入口处的脉动情况,都为0时跟层流接近但是跟层流不同。
你将湍流强度设为5%,粘性比设为0.01。
再试试看,有变化的话,换别的湍流模型再试下。
UH另外周期边界中你设定压降还是流量?
流量的指定方式更加容易出现波动。
pMy
我设定的是流量。
刚才我变化了这两个湍流变量,影响比较大。
我这段时间一直在试验不同的湍流模型,以找到比较适合的。
我试过别的湍流模型,最后都是趋于稳定的。
另外,我选用层流和rng模型时,最后的振荡基本是一致的(如图,绿色:
rng,白色:
层流)。
我计算的过程是,先用稳态算法算出一个初场,然后改用非稳态算法继续计算(有人说过这样的方法),不知道这个是不是有什么影响?
这样计算时,非稳态计算的初场是不同的。
换成别的湍流模型时,即使用稳态得到一个类似于rng的计算结果,然后改用非稳态算法,最后都是稳定的。
上下对称就不会振荡,不对称才会振荡。
其他湍流模型的初场都类似于这个不对称流场,但是最终都是稳定了,即上下对称。
nUvIv
按照我的经验,你这样算:
^P
选用层流模型,通过修改流体的物性来凑雷诺数,然后把流速条的小一些(雷诺数不变),然后计算就可以了,所有的松弛因子调到1,连续性方程0.9,用simplc算法,离散格式用二阶,应该就没错了。
vHSp
PS上面那个图是算淹没射流同志的case,不知道是不是他想要得结果。
流体模拟
有些文献是这样计算的,通过改变粘性来改变雷诺数。
但是我现在算的这个东西比较具体了,应该说工质是空气就不会变了。
我是改变进口流量来调节雷诺数,松弛因子连续方程取为0.7,其他都一样,采用湍流模型,湍流的松弛因子只要不是1就可以。
离散格式换了几个,基本相同。
一个问题是,雷诺数大了之后,层流和rng的计算结果就不一样了,层流还是规律性振荡,而rng的就不太规律了(和文献上类似)。
因为我的计算要涉及到层流流动和过渡流动,但是又不知道什么时候转化,所以想从一开始就使用湍流,以监视流场的变化。
D=1V=1
为了不同的雷诺数,所以在Fluent中改变材料动力粘性系数,
但是计算的结果完全相同,甚至于雷诺数到0.5的时候都会出现涡(层流,定常)而且不同的雷诺数连收敛的步数都是一样的,奇怪哉!
!
模型是不是有问题呢?
按照道理应该不一样的,可能的原因1.圆柱饶流在re100多时就有卡门涡街了,应该采用非定常计算方法2.re3900是亚临界雷诺数,尾迹转捩,应该采用湍流模型计算3.网格问题,通常网格在近壁面附近要划分较细
呵呵,先就谈这些吧,你给的信息太少了
在Fluent中改变材料动力粘性系数:
检查一下是否再fluid里面改了
各位大侠,看到很多文献中提到涡脱落周期T,而且是保存的1/4T、1/2T等各个不同时刻的文件,进行涡脱落分析;我不知道脱落周期如何确定,以及如何保存精确的1/4T时刻;
还有我可否在计算完圆柱绕流后,与文献对比确定计算的准确性,然后将所采用的模型应用到其他形状管子的模拟上;有什么理论指导吗?
可以推荐一本好书给我。
非常感谢了!
流体模拟
周期一般是通过FFT求出的,需要某点的压力或壁面的升力系数。
通过计算若干个周期,再用FFT就可以求出涡脱落频率了。
知道频率也就知道周期了。
然后将1/4T,1/2T的文件保存下来,可以看看涡,速度等图片的变化情况。
不知道FLUENT行不行,具体该怎么做?
FFT是什么?
?
?
?
傅立叶变换
FFT是快速傅立叶变换,fluent里面应该没有对某点数值的时程图进行傅立叶变换吧(matlab里面有这个功能),不过我想,要是能把某点的数值的时程图提取出来,不用进行FFT,应该也不难判断周期了吧。
最近正在算一个圆柱绕流问题,Re=200,velocity=1m/s,desity=1kg/m**3ot)diameter=1m,viscosity=0.005kg/m.s,$%n
虽然后面有涡街,但是升力系数和阻力系数和实验的或者别人计算的对不上。
D0请问各位大虾这是怎么回事?
望高手不吝赐教!
yi*x
(Report里的referencevalue我取的是:
depth=1m,velocity=1m/s,Z[K2
area=1m**2,density=1kg/m**3)
我觉得referencevalue中的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,计算中U=1,L=1(U为来流速度,L为方柱边长)v-y)K计算区域为:
0=x=60L0=y=8LY
在计算中我固定Re(100)、tau(0.65)和dx(1/30),然后确定c,开始我用tau=1.6,然后发现Re最大只能到30,后来我设tau=0.65,我发现可以计算出Re=100,但是不知为何,始终得不到涡脱落的结果,仍然只是定常的对涡现象,不知道问题出在哪里。
我想请问一下,无论设置tau和c为多少,只要保证Re数一定,是否计算的流场也应该一定?
[+b"
谢谢!
流体模拟
我算的是圆柱绕流,同样存在这个问题.@Ur_
即使Re100,算不出来VotexShedding现象
不知道哪位大侠有高招,帮帮忙呀V0
我用FLUNT算一个立方体绕流问题,立方体为600mm*600mm*200mm,利用对称性,取600mm×300mm×200mm,计算域取*****mm×4500mm×1000mm,划分了180多万个四面体网格。
网格用GAMBIT生成。
b6`L
导入FLUENT中,采用层流模型,设置了入口的速度边界条件,如X=12m/s,进行计算。
结果发现得出的结构表面的压力系数的数量级为10的二次方,这显然是不对的。
不知可能在那里出现问题。
可否有人指点一下。
万分感谢。
大家好,在看圆柱绕流的文章时,看到一个对称边界条件,不太理解,这个对称边界条件是如何定义的,有什么物理意义呢?
A,j
请不吝赐教。
谢谢。
对称边界条件很常见的.+W
顾名思义,就是流场中的一个对称面,面两边的流动型态是完全对称的.在求解的时候为了节约计算量,经常只解一半,而把对称面设为对称边界条件.TG(t由于对称性,对称边界条件有以下特点:
D+'
1.穿过对称面的通量(法向速度)为0%/
2.对称面法线方向的导数为0b2a
3.不同于壁面,对称边界条件容许有切向速度和切线方向的变化.
3楼的回答的正确,但是在圆柱绕流中不能用对称边界,因为要是Re较大的话,会产生涡街,不符合三个条件。
p
流体模拟
大家好,圆柱绕流数值模拟的时候,时间步长如何选择合适呢?
有个老师告诉我通常取无量纲的时间步长为0.002-0.01,为什么这样选择?
\
时间步长小应该来说计算结果更精确,但是不是由此会导致收敛速度慢,容易发散?
l在圆柱绕流的数值模拟中,除了低雷诺数绕流,哪些情况会导致无漩涡脱落、无分离流动现象出现?
[
望大虾不吝赐教,以助后学者进步,谢谢。
时间步长小,每个步长收敛速度会变快,也更不容易发散.但总计算量可能变大.别人算一个周期用10步,你用100步当然计算时间可能变长了.@~9=
步长太长不行,可能导致无漩涡脱落,或看不到脱落过程.eDaDbP
一个周期至少20个时间步吧.'18x
时间步长太小也没必要.简单看看漩涡脱落,一个周期100个时间步足够了._xs上次我算着玩的时候,不是用了20就是40个时间步.漩涡脱落没问题.
要想算出漩涡脱落一般要用高阶精度的格式.
r
"时间步长小,每个步长收敛速度会变快,也更不容易发散"F7c
可是我用很小的时间步,每时间步内的结果却不容易收敛,最后还没有涡脱,不知怎么回事。
自己编的程序还是用商业的?
OKJ
自己编的话就很难说了.
在别人的基础上改编的,呵呵sk
这实际上如同实验的采样频率,为了扑作到某一个绕流特征,每一个周期内必须采集一定数量的样本,如:
无量纲时间步长0.01对应100次,不同的绕流特征需要的采样频率有一定的差别,因此有0.01~0.002一说。
,9pWh2
Ifthereisvortexsheddingdoesnotonlydependonthenumericalparameters.Fortheflowovercircularcilinder,youneedtoknowhowlargetheReynoldsnumberis.IfRe1,nowaytoobservethevortexshedding.
CheckBatchlor'stexttoseeathowlargeRe#theflowcanbecomeunsteady.
流体模拟
我的毕业论文有关于剪切流绕过圆柱的计算,坐标系固连于圆柱上,远场来流是一个剪切流,z0时速度为U0+u,z0时速度为U0-u。
用Fortran编程,采用有限差分法计算。
请问边界条件该如何提,网格又该如何划分?
有人用fluent做过圆柱绕流的模拟吗?
交流一下。
)
godbless@http:
//.cnks
我正在做钝体绕流的问题,请你把问题说的具体一些。
我做过一次圆柱绕流的问题,你有什么问题吗?
(}.
俺做过几次,不过不是用软件算的。
湍流模型采用的k-emodel,近壁区采用两层壁面函数,但高雷诺数结果不是太令人满意。
我的同学就是做圆柱绕流的,较果不怎么好!
Re=100,water,网格数*****结构网格i;
还有,reportrefferrence里的参数需要设置吗?
F{BkM
谢谢!
*-
当然要设置I
作的是海洋立管的绕流,用2D,转化为平面问题,直径0.3米,试验了1m/sQBi
0.5m/s等流速,雷诺数当然比较大了.#
网格也作了边界层,但怎么也得不到边界层脱体,也看不见周期性的涡旋释放.升力系数,曳力系数振荡后为零
最终的目的是获得周期性的升力,并与立管固有频率对比.
非常感谢,我是大连理工造船系的,谢谢!
!
又试验了较小的速度,0.05m/s,结果在相同网格条件下,用层流就分析出了脱体涡,并获得了周期性的升力和阻力系数;.N
但若用湍流模型则没有脱体现象
我这两天也在算这个,把你的cas,date贴上来研究一下了!
`
好的,我下载下来看看!
/
是不是边界层设定有些问题啊,lz检查一下
流体模拟
我用的是结构化网格的计算了好久,我用K-E和RSM模型计算了,但也没有出现边界层脱体现象!
难道是网格还不够密!
为什么用层流的模型可以得到周期性的涡旋释放呢,在论坛中看见以前也有人问过这个问题,l另,10E+5级别的雷诺数的绕流完全不能用层流分析么么?
我们在比较高的流速时,得到了脱体现象,x!
4f}
在圆柱边界有一层反向流动,但后面的涡kU6@
不是交替释放的啊,也怀疑是否是湍流破坏了涡结构.+
边界层太细了好像也不好,加我QQ把,我也作者方面的工作;
啊啊;;***-*****8#
为什么用层流的模型可以得到周期性的涡旋释放呢,d
~~~在你给定的雷诺数范围内,流动是非定常,层流。
用湍流模型算层流一般不行。
k另,10E+5级别的雷诺数的绕流完全不能用层流分析么么?
~~~~~一般情况是湍流了,若用层流来算,应该是DNS,这个级别的DNS的资源需要太大了。
D
版主你好,非常感谢你的回答,但是我还有问题不明白,关于海洋立管得绕流`有如下介绍
Re在40到150内是周期性交替泄放的层流涡旋@Z[t
Re在300到3乘10的五次方是周期性交替泄放的紊流涡旋,也是工程上最关心的流速范围ya
这样0.05m/s,也是紊流啊,另外所说的层流涡旋需要用层流模型分析,紊流涡旋需要用紊流模型分析,对么r"+et
另外,我用的是6.1.22版本壁面模型怎么没有双层区模型啊,却是加强的壁面处理
另外两个是标准壁函数和非平衡壁面函数,什么原因呢s
是不是我问得有不合理的地方呢,谢谢版主{
我来说两句,好像有很多人关心这个问题,也包括香港地区的朋友,我把我的经验告诉大家:
N3a4
1)网格需要足够密~K\b
2)一般的湍流模型也能算出来,最好用大涡模拟或分离涡模拟(*****.2版本比较成熟)pU
3)时间步长要合理,不能太大,如果太大就算不出来了,而且不同的雷诺数也需要不一样。
众所周知,St数与Re数相关,St中有包含涡脱落的周期(或频率),如果时间步长大于此周期将观察不到,但是如果你的时间步长太小又会出现很长计算时间也观察不到周期现象。
HWO/
4)如果上述方法还不能奏效,你可以在计算时激发它一下,比如初场不均匀,然后经过计算一般就会出现稳定的周期现象。
这一点很重要的喔,这可是压箱底的经验喔
流体模拟
1)我们的网格画不好没,尤其是圆柱边界层附近的网格需要加密,我们画不出来结构化的四边网格,不知道怎么把流场区域分成可以画的块.比如加上径向的虚边,再布设不均匀边网格点,总难以成功,简直心灰意冷啊,比有限元网格难画多了!
|'ELuv
2)我们只有6.1.22版2DFLUENT好象没有LES啊,3D的有啊,分离涡模拟是什么方法啊,中科大教材中算例,是用RNG模型加上双模型壁面得到的大雷诺数结果,我们的IPqIK
软件只有3个壁面处理标准,非平衡,加强的?
?
也是困惑啊gE%:
_
希望有人加以指点啊,谢谢!
!
f
我用的也是fluent6.1.22版,我个人认为3个壁面处理标准,非平衡,加强的中的加强的就是双模型壁面了,但是我计算后还是没有得到想要的结果,看来还是步长和网格的问题了!
1)网格问题我们解决了,但还是没有得到想要的大雷诺数绕流的周期性释放的湍流涡旋,升力系数和阻力系数曲线在震荡后都变成直线了;u
2)我们仅试验了采用如下方法扰动网格,先把上部速度改变一点,算一段时间后再改回来继续算,不知道能不能相当于初始化.f
3)关于"为加速旋涡的分离,流场初始化时给流场加一个扰动,上半部分速度为0.005m/s,下半部分速度为0,(需要用到用户定义场函数功能),然后在初始化后如果速度后patch中调用该函数,对整个流场进行初始化"k
我们还不会用啊,前面的问题也不太清楚,希望有人指点一二,谢谢!
jf
这样的高手应该提携一下众兄弟啊,可以认识一下吗,我的QQ***-*****3
2.一共四个问题:
nE
1、解释型udf要求函数中不能包含子函数,有可能使这个原因。
ZxXp
2、动网格求解中的锯齿形残差是正常的。
YNO
3、计算过程中,某几个单元温度越界没有关系,只要它能回来,别发散就行,这个原因较多,一般来说不是网格的问题,比如在求解超音速流动时,使用segregate算法时就容易出现。
K`
4、分离求解器可以用来求解可压超音速流,但要求结构网格同时网格间距跳跃不大,马赫数小于2.5时,将松弛因子全部设为0.1,是可以计算的。
马赫数高了就完蛋了。
还有计算结果跟couple的一样。
Gt
这些是我的经验,也不一定十分准确,仅供参考?
3.我现在算的是弹丸出炮口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每个时间步很难收敛,我就换成分离的,而且我亚松弛因子调的很大,压力有0.9,动量为0.8,其他都是1.我用的是三角形的非结构网格,气流速度有1500多,马赫数达到了4左右,可是在每个时间步上它收敛很快,一般16步左右就收敛了.想问下你,象我这样设那么大的亚松弛因子可以吗?
还有就想问你,象我这样的非稳态怎么看每
流体模拟
一个时间步是否收敛了?
我现在就是看残差图7nih)
流体中文网论坛--流体力学及相关领域学术问题交流论坛。
={!
直觉不太对,不过如果你的时间步长很小也许就可以避过每个时间步长内发散的可能性,如果时间步长比较大就不对了哈,另外炮弹出口如果直线运动用layering,比remesh好多了,又准确又快。
A~
PS,膛内外压力差多大?
Mpa级?
U:
%
是啊,有60多mpa.我时间步长为0.000001s4
还要请教楼主,想我这样的非稳态计算怎么样看每一步是否收敛了?
看残差吗?
可以采用稳态算法中的办法,监控某点的某个变量值,要在流场变化区域内,如喷口位置。
如果平了就好了。
4.我计算车室内流场,出风口很多,回风口(排风口)只有一个,采用hexcore网格、RNGk-e模型,计算过程监视出口(回风口)流量,计算到300步的时候差不多成了一条直线,而残差图也差不多成了直线,接着迭代300步,残差图还是下不去,也就是不能降低到低于1e-3的水平,在report里面看netimbalance在0.1%左右,请问这样能不能算收敛了?
Z#i
压力我用的standard,其他的是二阶迎风,有什么办法能使残差降到1e-3以下,谢谢?
[bN
有时候残差下不来是正常的,尤其是连续性方程,一般说来x,y,z速度收敛了,只剩连续性方程没收敛,其他监控变量已水平,可以认为计算收敛(包含传热时除外)。
也就是说你这个问题可以算收敛了,要想将连续方程残差降下来要改网格的,或者改变松弛因子造成假收敛,当然此法不可取。
5.wormxieu
你可以看下fluent6.1的第19个例子,是个欧拉模型的例子。
可以参考一下
大家好,圆柱绕流数值模拟的时候,时间步长如何选择合适呢?
有个老师告诉我通常取无量纲的时间步长为0.002-0.01,为什么这样选择?
时间步长小应该来说计算结果更精确,但是不是由此会导致收敛速度慢,容易发散?
在圆柱绕流的数值模拟中,除了低雷诺数绕流,哪些情况会导致无漩涡脱落、无分离流动现象出现?
望大虾不吝赐教,以助后学者进步,谢谢。
时间步长小,每个步长收敛速度会变快,也更不容易发散.但总计算量可能变大.别人算一个周期用10步,你用100步当然计算时间可能变长了.
步长太长不行,可能导致无漩涡脱落,或看不到脱落过程.
一个周期至少20个时间步吧.
时间步长太小也没必要.简单看看漩涡脱落,一个周期100个时间步足够了.
上次我算着玩的时候,不是用了20就是40个时间步.漩涡脱落没问题.
流体模拟
要想算出漩涡脱落一般要用高阶精度的格式.
自己编的程序还是用商业的?
自己编的话就很难说了.
我上面提到的那位老师说选取无量纲的时间步为0.002-0.01,我看了很多圆柱绕流的数值模拟的文章,也基本如此,为什么会是这样呢?
anyone说,“时间步长太小也没必要.简单看看漩涡脱落,一个周期100个时间步足够了,上次我算着玩的时候,不是用了20就是40个时间步.漩涡脱落没问题.”
这样做,漩涡脱落是没问题了吧,但是计算绕流阻力会不会有较大影响呢?
这实际上如同实验的采样频率,为了扑作到某一个绕流特征,每一个周期内必须采集一定数量的样本,如:
无量纲时间步长0.01对应100次,不同的绕流特征需要的采样频率有一定的差别,因此有0.01~0.002一说。
,
现在在用lbm做方柱绕流,采用郭照立提出的不可压模型,边界条件也是他的非平衡外推方法。
计算中粘度系数用以下表达式
niu=(2*tau-1)/6*(c*dx)
方柱绕流的Re数定义为:
Re=U*L/niu,计算中U=1,L=1(U为来流速度,L为方柱边长)
计算区域为:
0=x=60L0=y=8L
在计算中我固定Re(100)、tau(0.65)和dx(1/30),然后确定c,开始我用tau=1.6,然后发现Re最大只能到30,后来我设tau=0.65,我发现可以计算出Re=100,但是不知为何,始终得不到涡脱落的结果,仍然只是定常的对涡现象,不知道问题出在哪里。
我想请问一下,无论设置tau和c为多少,只要保证Re数一定,是否计算的流场也应该一定?
我算的是圆柱绕流,同样存在这个问题
即使Re100,算不出来VotexShedding现象
[以下内容由wsl2006在2021年12月07日07:
11pm时添加]
不知道哪位大侠有高招,帮帮忙呀
想做一个绕圆柱产生卡门涡街的动画,但在求解过程中总是得到对称的流场。
知道的人都晓得这个肯定不对,我想问问高手什么地方需要注意,可能是我什么地方设置错了!
`1Ks我用的是ideal-gasF;j
k-e湍流模型=
边界都是远场边界条件,来流速Ma=0.05表压1atm-@O
还有什么地方可能错了?
?
n
流体模拟
你是不是在流场初始化的时候也是对称的哦?
gB0
你试着初始化流场时不对称,看看4
但是卡门涡街本身没有要求非对称来流亚?
Dk
我说的是流场初始化,不能是对称的,但可以是对称的来流,^你可以看一下李万平那本书$e
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 模拟 建议