找回密码
 注册
查看: 237|回复: 2

关于前台阶绕流问题的壁面边界条件的设置

[复制链接]
发表于 2019-11-28 21:11:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?注册

x
本帖最后由 Sin_Luo 于 2019-11-29 08:31 编辑

最近写论文要添加前台阶问题的算例,被壁面边界的设置难住了。我尝试了两种处理方式:

第一种:根据对称性设置虚网格,然后正常计算壁面上的数据,计算完之后使垂直方向的速度为0,并修正切向速度使速度大小不变。简单地说,若1是壁面,我设置虚网格0,-1,-2与内侧的2,3,4对称,然后照常使用数值格式计算和更新1上的数据,只不过每一步都做修正,使垂直方向的速度为0,并修正切向速度使速度大小不变。
第二种:同样照常使用数值格式计算和更新1上的数据,只不过在设置虚网格时,使0和1对称,-1和2对称,-2和3对称。相当于真正的壁面在0和1之间,是隐性的。


我用两种方式处理,第一种和文献的图像更接近,并且激波反射的位置不随网格变化,但不稳定,使格式参数的选取受限。第二种更稳定,能用非常低耗散的格式计算,但激波反射的位置跟文献有些许差异,激波反射的位置随网格质量变化,网格变密后位置趋于一致。


诸位是怎么设置无粘流壁面边界条件的?我这两种处理对不对?求教高手!

QQ图片20191128203121.png
这是第一种处理的结果
QQ图片20191128205441.png
这是第二种处理的结果,反射位置靠后
QQ图片20191128210200.png
这是第二种处理,网格加密后的结果

发表于 4 天前 | 显示全部楼层

回帖奖励 +2 流金金钱

两种格式应该都是合理的,第一种是有限差分法的思路,第二种是有限体积法的思路,但是在第一种方法处理壁面问题时不知你是否是直接令法向速度等于零,如果是的话则可能导入比较大的数值误差引起不稳定,安德森的书上推荐一种局部膨胀/压缩的办法进行修正,可能效果会更好一些。
回复 支持 1 反对 0

使用道具 举报

 楼主| 发表于 3 天前 | 显示全部楼层
周华 发表于 2019-12-4 11:34
两种格式应该都是合理的,第一种是有限差分法的思路,第二种是有限体积法的思路,但是在第一种方法处理壁面 ...

谢谢老师!问题应该就在这里,我直接使法向速度为0了。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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