分析下压缩拐角算例出现这种结果的原因
本帖最后由 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次方量级）。
因为缺少这方面经验，不知道怎么样的结果算是好的结果，不知道别人的计算中有没有遇到过类似情况，怎么解决的，不知道在激波前后过冲多少可以接受。特来此请教。 求交流ing 有没有人有相关经验的？？ 已经算得很好了 :lol
页:
[1]