找回密码
 注册
楼主: lixl-imech

发布OpenCFD-EC2D; 更新 Ver 1.5.4 ; 三维软件使用协议

  [复制链接]
发表于 2011-1-30 13:49:09 | 显示全部楼层
是不是我算的问题由于会出现涡街,不能用对称算?

非定常,我觉得有些问题也是可以用对称来算的。

谢谢您的回答
发表于 2011-1-30 14:16:16 | 显示全部楼层
我想说的就是因为非定常,中线上不对称。
发表于 2011-1-30 16:00:00 | 显示全部楼层
发表于 2011-2-1 09:29:01 | 显示全部楼层
李老师,那个压力远场边界能不能给个公式说明,
有些地方没看懂,
看了这本书后Computational Fluid Dynamics: principle and application ,感觉亚音速入口和出口的远场条件很难控制。

[ 本帖最后由 ch06 于 2011-2-13 01:44 编辑 ]
 楼主| 发表于 2011-2-25 15:42:53 | 显示全部楼层
原帖由 ch06 于 2011-2-1 09:29 发表
李老师,那个压力远场边界能不能给个公式说明,
有些地方没看懂,
看了这本书后Computational Fluid Dynamics: principle and application ,感觉亚音速入口和出口的远场条件很难控制。


亚音速情况下,扰动波双向传播。因此,无论入口条件还是出口条件都比超音速情况复杂。 以亚声速出口边界条件为例进行说明,亚声速入口条件类似。
假设b点为远场出口的边界点,d点为紧邻它的上游内点。
如果流动是超声速的,边界条件就简单了,直接令b点处的物理量(密度、速度、压力)等于d点的值即可。由于是超声速流,扰动不会向上游传播,因此b点处的边界条件误差不会影响计算域内部。
如果是b点处的流动是亚声速的,边界处的扰动会传播到计算域内部。物理模型需要给定一个边界条件。 通常情况下,给定无穷远压力(即背压)pb。 这样,b点处的压力就确定了。 b点处其他量(密度d, 速度u,v,w)由上游值及pb来确定。一个简单的模型就是J. Blazek书283页的公式 (8.23):  b点的密度为: rb=rd+(pb-pd)/(c0*c0)
这个公式的含义是:
如果计算域内部(d点)的压力刚好等于背压pb, 则流动既不膨胀也不压缩,出口(b)点的密度为d点的密度rd;
如果内点的压力高于背压(pd>pb), 则流动在出口要膨胀,出口的密度会低于内部;
如果内点的压力小于背压(pd<pb), 则流动在出口要被压缩,出口的密度会高于内部;
《空气动力学》或《流体力学》喷口流动规律里面会有类似介绍。
应当指出,(8.23)式只是一个近似公式。但当出口远离计算域中心时,其误差对计算域内部影响不大。

[ 本帖最后由 lixl-imech 于 2011-2-25 15:44 编辑 ]
发表于 2011-2-25 18:37:48 | 显示全部楼层
感谢李老师的耐心解答,我再试一下,

[ 本帖最后由 ch06 于 2011-2-25 12:09 编辑 ]
发表于 2011-2-25 18:39:00 | 显示全部楼层
李老师请教如果出现下面的情况,应该如何处理

激波在边界附近时,可能刚开始是亚声速出口,满足下面的情况
else if (Ma_n .gt. 0.d0) then  !  亚声速出口
       pb=p_out
            db=d1+(p_out-p1)/(c1*c1)
            ub=u1+(p1-p_out)/(d1*c1)*n1
            vb=v1+(p1-p_out)/(d1*c1)*n2
            p2=2.d0*pb-p1 ; d2=2.d0*db-d1 ; u2=2.d0*ub-u1 ; v2=2.d0*vb-v1

Ma_n .gt. 0,但如果激波离边界非常近,可能存在p2=2.d0*pb-p1 <0的情况(即激波在p1和p2之间,或在p1点上,此时p1很大)

[ 本帖最后由 ch06 于 2011-2-25 12:15 编辑 ]
发表于 2011-2-27 10:52:34 | 显示全部楼层
我在模拟这两幅图时,就遇到了上面的情况
QQ截图未命名.jpg
发表于 2011-2-27 10:54:27 | 显示全部楼层
采用自由出流,总会有一些波在边界形成干扰
采用远场,当亚音速入口和出口时,经常会出现负压
 楼主| 发表于 2011-2-28 13:56:26 | 显示全部楼层
原帖由 ch06 于 2011-2-27 10:54 发表
采用自由出流,总会有一些波在边界形成干扰
采用远场,当亚音速入口和出口时,经常会出现负压


边界上波的反射不太容易消除。 比较简便的方法是把计算域取得大一些,出口边界远离感兴趣的区域,这样即使边界上有反射,对感兴趣的区域影响也不大。 另外,计算域取大以后,远离中心的区域网格比较稀疏,大的网格间距会产生较大的数值耗散,可以吸收掉反射波。
发表于 2011-2-28 16:17:55 | 显示全部楼层
好的,我扩大一下计算域试试。
发表于 2011-3-20 20:09:31 | 显示全部楼层
请问输出的残差文件中4个最大,4个rms,这4个值都是什么的呀?u,v,p?
发表于 2011-3-25 12:35:58 | 显示全部楼层
多谢李老师!
 楼主| 发表于 2011-3-26 10:36:43 | 显示全部楼层
原帖由 shirazbj 于 2011-3-20 20:09 发表
请问输出的残差文件中4个最大,4个rms,这4个值都是什么的呀?u,v,p?

输出的4个值是守恒型方程的残差。 最大残差即所有控制体单元里残差的最大值;均方根残差是所有控制体单元残差的均方根值。
“残差”是积分型N-S方程的右端项,其物理含义为净流量。 即该时间步内流入控制体单元的净质量、净动量(包括x 以y方向的动量)及净能量。
发表于 2011-3-30 20:55:47 | 显示全部楼层
李老师,这个程序可以算非结构网格的吗
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表