哪位做过激波抖振计算

junior

我用超临界翼型进行计算，翼面发生激波抖振，但是随计算时间增长，抖振幅度越来越小，谁知道怎么判断这个抖振究竟是物理的还是数值的？

通流

是激波本身在震动吗？
这个问题不太容易。数值与物理都可能造成激波震动。我还想问问搞CFD可靠性的人怎么看这个问题的。

junior

是激波在动，导致升力系数、阻力系数振荡，类似于正弦波，但振幅逐渐衰减，所以不清楚是物理的还是数值的。

通流

是非定常解吧？频率是多少？把时间步长设成大大小于周期（如小于0.01周期），看看结果。如果是物理现象，振动频率可能与（1）尾缘的涡脱落频率，或（2）激波后的边界层分离。第一种情况比较好弄，你可以估算一下尾缘涡的频率，看看与你计算的频率是不是很接近。第二种情况我只是听人作报告时说过。到底物理上怎么解释，我也不清楚。
如果是数值问题的话，你可以试试在激波处加密网格，看看是不是对频率有很大影响，如果是的话，那么这种震动是由于数值造成的可能性就很大了。

junior

谢谢版主！应该是数值振荡，不是激波抖振。

通流

希望这是个科学判断，也就是有证据支持的，不是别人，如你的老板，告诉你的。

junior

我觉得这是个科学判断，因为在流动条件没有变化的时候，激波抖振不会随时间衰减。

通流

给你举个例子吧。一个带阻尼的单摆，随着时间推移，它的增幅在衰减。这是个真实的，而且情形与你描写的差不多例子。
换句话说，在你设的流动条件不变时，流体震动的增幅可以增大，可以衰减，也可能不变更复杂的情形是有极限环，甚至有不止一个的极限环。

梅羡仙

流体力学就是个集散微粒学，你从集散微粒的角度出发，一切皆迎刃而解。

junior

用一阶格式算的时候后缘没有旋涡脱落，用二阶格式算的时候有旋涡脱落，不过应该没有极线环那么复杂。

通流

极限环是很常见的。其实任何有限的震动，在振幅稳定后的情况都可以看成是极限环。

按照你的判断，到底有没有涡的脱落？

回到你的激波震动的问题上来。计算最后成了什么样了？振幅减小，最后消失了吗？测了衰减率了吗？

[ 本帖最后由 通流 于 2010-2-20 21:04 编辑 ]