求帮分析全位势程序（full potential equation求解

shiyongjie

因研究需要，程序中部分使用全位势方法求解。根据相关资料编制求解二维翼型的全位势方程。全位势在遇到激波问题时，需要增加人工粘性修正（对密度离散采用逆风格式）。编好的程序在计算亚音速流时收敛速度和计算结果都比较不错。但是在出现强激波的时候总是出错。程序检查了很多变，并未发现任何问题。现将一个典型计算算例在不同迭代步骤上的计算结果给出，请了解的人帮忙分析指点下，谢谢了。
马赫数0.8 ，迎角2度
以下四个图分别是90步，100步，110步，120步的计算结果。刚开始时，计算捕捉到了激波，但是随着计算的进行，激波前和激波后开始出现发散现象。请问这是什么原因呢？

[ 本帖最后由 shiyongjie 于 2009-11-4 11:25 编辑 ]

shiyongjie

没有人知道吗？

xinboy

你的方法没用过！是不是边界条件的问题啊？
我现在也在算这个翼型，我发现远场边界离翼型太近的话，需要进行涡修正。
不知道你的cl是多少？
另外，我现在在搜寻一本外文文献，“Test Cases for Inviscid Flow Field Methods”
不知道你有没有啊？
如果没有和我一起找吧，这是一个标准算例，可以拿自己的结果与之对比！