在计算中我固定Re(100)、tau()和dx(1/30),然后确信c,开始我用tau=,然后发觉Re最大只能到30,后来我设tau=,我发觉能够计算出Re=100,可是不知为何,始终得不到涡脱落的结果,仍然只是定常的对涡现象,不明白问题出在哪里。
我想请问一下,不管设置tau和c为多少,只要保证Re数必然,是不是计算的流场也应该必然?
[+b"
谢谢!
我算的是圆柱绕流,一样存在那个问题.@Ur_
即便Re>100,算不出来VotexShedding现象
不明白哪位大侠有高招,帮帮忙呀V0
我用FLUNT算一个立方体绕流问题,立方体为600mm*600mm*200mm,利用对称性,取600mm×300mm×200mm,计算域取10200mm×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
大伙儿好,圆柱绕流数值模拟的时候,时刻步长如何选择适合呢?
有个教师告知我通常取无量纲的时刻步长为,什么缘故如此选择?
\
时刻步长小应该来讲计算结果更精准,可是不是由此会致使收敛速度慢,容易发散?
l在圆柱绕流的数值模拟中,除低雷诺数绕流,哪些情形会致使无漩涡脱落、无分离流动现象显现?
[
望大虾不吝赐教,以助后学者进步,谢谢。
时刻步长小,每一个步长收敛速度会变快,也更不容易发散.但共计算量可能变大.他人算一个周期用10步,你用100步固然计算时刻可能变长了.@~9>=
步长太长不行,可能致使无漩涡脱落,或看不到脱落进程.eDaDbP
一个周期至少20个时刻步吧.'18x
时刻步长过小也没必要.简单看看漩涡脱落,一个周期100个时刻步足够了._x&s
上次我算着玩的时候,不是用了20确实是40个时刻步.漩涡脱落没问题.
要想算出漩涡脱落一样要用高阶精度的格式.r
"时刻步长小,每一个步长收敛速度会变快,也更不容易发散"F7c
可是我用很小的时刻步,每时刻步内的结果却不容易收敛,最后尚未涡脱,不知怎么回事。
自己编的程序仍是用商业的?
OKJ>
自己编的话就很难说了.
在他人的基础上改编的,呵呵sk
这事实上犹如实验的采样频率,为了扑作到某一个绕流特点,每一个周期内必需搜集必然数量的样本,如:
无量纲时刻步长对应100次,不同的绕流特点需要的采样频率有必然的不同,因此有~一说。
,9pWh2
Ifthereisvortexsheddingdoesnotonlydependonthenumericalparameters.Fortheflowovercircularcilinder,youneedtoknowhowlargetheReynoldsnumberis.IfRe<1,nowaytoobservethevortexshedding.
CheckBatchlor'stexttoseeathowlargeRe#theflowcanbecomeunsteady.
我的毕业论文有关于剪切流绕过圆柱的计算,坐标系固连于圆柱上,远场来流是一个剪切流,z>0时速度为U0+u,z<0时速度为U0-u。
用Fortran编程,采纳有限差分法计算。
请问边界条件该如何提,网格又该如何划分?
有人用fluent做过圆柱绕流的模拟吗?
交流一下。
)
我正在做钝体绕流的问题,请你把问题说的具体一些。
我做过一次圆柱绕流的问题,你有什么问题吗?
(}.
俺做过几回,只是不是用软件算的。
湍流模型采纳的k-emodel,近壁区采纳两层壁面函数,但高雷诺数结果不是太令人中意。
我的同窗确实是做圆柱绕流的,较果不怎么好!
Re=100,water,网格数33500结构网格i;
还有,reportrefferrence里的参数需要设置吗?
F{BkM
谢谢!
*-
固然要设置I
作的是海洋立管的绕流, 用2D, 转化为平面问题,直径0.3米,实验了1m/sQBi
0.5m/s等流速,雷诺数固然比较大了.#
网格也作了边界层,但怎么也得不到边界层脱体,也看不见周期性的涡旋释放.升力系数,曳力系数振荡后为零
最终的目的是取得周期性的升力,并与立管固有频率对照.
超级感激,我是大连理工造船系的,谢谢!
!
又实验了较小的速度,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
如此s,也是紊流啊,另外所说的层流涡旋需要用层流模型分析,紊流涡旋需要用紊流模型分析,对么r"+et
另外,我用的是6.1.22版本壁面模型怎么没有双层区模型啊,却是增强的壁面处置
另外两个是标准壁函数和非平稳壁面函数,什么缘故呢s
是不是我问得有不合理的地址呢,谢谢版主{
我来讲两句,仿佛有很多人关切那个问题,也包括香港地域的朋友,我把我的体会告知大伙儿:
N3a4
1)网格需要足够密~K\b
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
如此的高手应该提携一下众兄弟啊,能够熟悉一下吗,我的
2.一共四个问题:
n>E
一、说明型udf要求函数中不能包括子函数,有可能使那个缘故。
ZxXp
二、动网格求解中的锯齿形残差是正常的。
YNO
3、计算进程中,某几个单元温度越界没有关系,只要它能回来,别发散就行,那个缘故较多,一样来讲不是网格的问题,比如在求解超音速流动时,利用segregate算法时就容易显现。
K`
4、分离求解器能够用来求解可压超音速流,但要求结构网格同时网格间距跳跃不大,马赫数小于时,将松弛因子全数设为,是能够计算的。
马赫数高了就完蛋了。
还有计算结果跟couple的一样。
G>t
这些是我的体会,也不必然十分准确,仅供参考?
3.我此刻算的是弹丸出炮口后膛内高温高压气体往外喷的模型.我开始用的是耦合求解器可是每一个时刻步很难收敛,我就换成份离的,而且我亚松弛因子调的专门大,压力有,动量为,其他都是1.我用的是三角形的非结构网格,气流速度有1500多,马赫数达到了4左右,可是在每一个时刻步上它收敛专门快,一样16步左右就收敛了.想问下你,象我如此设那么大的亚松弛因子能够吗?
还有就想问你,象我如此的非稳态怎么看每一个时刻步是不是收敛了?
我此刻确实是看残差图7nih)
©流体中文网论坛--流体力学及相关领域学术问题交流论坛。
={!
直觉不太对,只是若是你的时刻步长很小或许就能够够避过每一个时刻步长内发散的可能性,若是时刻步长比较大就不对了哈,另外炮弹出口若是直线运动用layering,比remesh好多了,又准确又快。
A>~
PS,膛内外压力差多大?
Mpa级?
U:
%
是啊,有60多mpa.我时刻步长为
还要请教楼主,想我如此的非稳态计算怎么样看每一步是不是收敛了?
看残差吗?
能够采纳稳态算法中的方法,监控某点的某个变量值,要在流场转变区域内,如喷口位置。
若是平了就好了。
4.我计算车室内流场,出风口很多,回风口(排风口)只有一个,采纳hexcore网格、RNGk-e模型,计算进程监视出口(回风口)流量,计算到300步的时候差不多成了一条直线,而残差图也差不多成了直线,接着迭代300步,残差图仍是下不去,也确实是不能降低到低于1e-3的水平,在report里面看netimbalance在%左右,请问如此能不能算收敛了?
Z#i
压力我用的standard,其他的是二阶迎风,有什么方法能使残差降到1e-3以下,谢谢?
[bN
有时候残差下不来是正常的,尤其是持续性方程,一样说来x,y,z速度收敛了,只剩持续性方程没收敛,其他监控变量已水平,能够以为计算收敛(包括传热时除外)。
也确实是说你那个问题能够算收敛了,要想将持续方程残差降下来要改网格的,或改变松弛因子造成假收敛,固然此法不可取。
5.wormxieu
你能够看下的第19个例子,是个欧拉模型的例子。
能够参考一下
大伙儿好,圆柱绕流数值模拟的时候,时刻步长如何选择适合呢?
有个教师告知我通常取无量纲的时刻步长为,什么缘故如此选择?
时刻步长小应该来讲计算结果更精准,可是不是由此会致使收敛速度慢,容易发散?
在圆柱绕流的数值模拟中,除低雷诺数绕流,哪些情形会致使无漩涡脱落、无分离流动现象显现?
望大虾不吝赐教,以助后学者进步,谢谢。
时刻步长小,每一个步长收敛速度会变快,也更不容易发散.但共计算量可能变大.他人算一个周期用10步,你用100步固然计算时刻可能变长了.
步长太长不行,可能致使无漩涡脱落,或看不到脱落进程.
一个周期至少20个时刻步吧.
时刻步长过小也没必要.简单看看漩涡脱落,一个周期100个时刻步足够了.
上次我算着玩的时候,不是用了20确实是40个时刻步.漩涡脱落没问题.
要想算出漩涡脱落一样要用高阶精度的格式.
自己编的程序仍是用商业的?
自己编的话就很难说了.
我上面提到的那位教师说选取无量纲的时刻步为,我看了很多圆柱绕流的数值模拟的文章,也大体如此,什么缘故会是如此呢?
anyone说,“时刻步长过小也没必要.简单看看漩涡脱落,一个周期100个时刻步足够了,上次我算着玩的时候,不是用了20确实是40个时刻步.漩涡脱落没问题.”
如此做,漩涡脱落是没问题了吧,可是计算绕流阻力会可不能有较大阻碍呢?
这事实上犹如实验的采样频率,为了扑作到某一个绕流特点,每一个周期内必需搜集必然数量的样本,如:
无量纲时刻步长对应100次,不同的绕流特点需要的采样频率有必然的不同,因此有~一说。
,
此刻在用lbm做方柱绕流,采纳郭照立提出的不可压模型,边界条件也是他的非平稳外推方式。
计算中粘度系数用以下表达式
niu=(2*tau-1)/6*(c*dx)
方柱绕流的Re数概念为:
Re=U*L/niu,计算中U=1,L=1(U为来流速度,L为方柱边长)
计算区域为:
0<=x<=60L0=在计算中我固定Re(100)、tau()和dx(1/30),然后确信c,开始我用tau=,然后发觉Re最大只能到30,后来我设tau=,我发觉能够计算出Re=100,可是不知为何,始终得不到涡脱落的结果,仍然只是定常的对涡现象,不明白问题出在哪里。
我想请问一下,不管设置tau和c为多少,只要保证Re数必然,是不是计算的流场也应该必然?
我算的是圆柱绕流,一样存在那个问题
即便Re>100,算不出来VotexShedding现象
[以下内容由wsl2006在2006年12月07日07:
11pm时添加]
不明白哪位大侠有高招,帮帮忙呀
想做一个绕圆柱产生卡门涡街的动画,但在求解进程中老是取得对称的流场。
明白的人都知道那个确信不对,我想问问高手什么地址需要注意,可能是我什么地址设置错了!
`1Ks
我用的是ideal-gasF;j
k-e湍流模型=
边界都是远场边界条件,来流速Ma=表压1atm-@O
还有什么地址可能错了?
?
n
你是不是在流场初始化的时候也是对称的哦?
gB0
你试着初始化流场时不对称,看看4
可是卡门涡街本身没有要求非对称来流亚?
Dk
我说的是流场初始化,不能是对称的,但能够是对称的来流,^
你能够看一下李万平那本书$e