找回密码
 注册
查看: 2267|回复: 3

分析下压缩拐角算例出现这种结果的原因

[复制链接]
发表于 2018-4-26 13:53:15 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 2012ZGZYY 于 2018-4-26 17:23 编辑

模拟来流马赫3的压缩拐角激波,激波左右值的精确结果应该是:左边(rho,rho*u,rho*v,rho*e)=(1, 3.5496, 0,8.8);右边(rho,rho*u,rho*v,rho*e)=(1.6180, 5.2933, 0.8858, 13.8692)。采用守恒形式欧拉方程的DG格式近似空间项,加入人工粘性。时间推进采用rk3或者隐式牛顿迭代。
人工粘性只加在cell_term上,采用的是经典做法,即粘性大小取决于当地残值,引入两个参数C和β,C是系数控制粘性大小,β控制当地网格尺寸的影响。很不解的是为什么会出现图中那种条带形状的结果,这是激波附近的数值振荡吗?(从数值解结果来看是在激波前后确实存在过冲情况) 说明需要加大粘性? 我按照文献给定的粘性C=0.02,β=0.1就是图示的结果。试着增大C值或β值,发现情况略有改善,即在x方向的动量云图上该现象可能消失,但密度始终都存在那种条带状的现象,不知道为什么。如果粘性加得过大可能导致流场出负值,或者收敛到-2次方就收敛不下去了(粘性加的小的时候能收敛到-9次方量级)。
因为缺少这方面经验,不知道怎么样的结果算是好的结果,不知道别人的计算中有没有遇到过类似情况,怎么解决的,不知道在激波前后过冲多少可以接受。特来此请教。

密度云图

密度云图

动量云图

动量云图

能量云图

能量云图

马赫数云图

马赫数云图

压力云图

压力云图
 楼主| 发表于 2018-4-26 14:54:39 | 显示全部楼层
求交流ing
 楼主| 发表于 2018-4-27 12:55:49 | 显示全部楼层
有没有人有相关经验的??
发表于 2018-7-23 10:02:54 | 显示全部楼层
已经算得很好了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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