请教二维翼形定常计算的问题

runheart

那位同仁做过二维翼形的流场数值模拟，请教几个问题：）
     我做的算例是naca0012，在零度攻角下的不可压流场模拟。自己做的程序，采用同位
网格上的simple方法，用动量插值法解决的压力速度耦合问题，采用的quick对流离散格式
。用的C型网格。
     出现的问题是，1）在翼形尾迹里速度竟然加快了，而且很大！后来对尾迹网格稀疏
后，在程序收敛后速度加快现象不存在了。我想请教一下这个问题的原因是什么？能否在
一开始计算时就能保证尾迹中速度是亏损的？
                   2）大家都知道，尾迹处也是C型网格上下部分的衔接处，这里要给出
边界条件的。我是给了速度的周期边界调，即：u下边界＝u上边界，v下边界＝v上边界，
压力的边界由速度来保证（用连续方程来得到压力方程时在边界上直接引入速度），而压
力梯度的计算却要给这部分的边界条件，我是直接外推了，这样做是否正确？

paradize

我做过可压流NACA0012翼型跨声速流场的模拟，用的是O型网格，不存在尾迹边界处理。
我觉得如果采用C型网格，边界条件不应给“u下边界＝u上边界，v下边界＝v上边界”，而应该给定压力“P下边界＝P上边界”，速度上下平均。
不知可否？

大地震

我的一位朋友用拟压缩性方法计算不可压流时，也出过类似情况！压力分布尾缘处出现明显上敲。究其原因，未讨论出结果。只是后来听他说尾迹区的网格蔬密对结果有影响，但消除不掉加速现象！

runheart

多谢两位了：
    我把压力的边条这样改了，程序很快就发散了，不知和simple方法有没有什么关系。
    在尾迹中用粗网格是能够消除加速现象的，不过要等很收敛了才行，这样会影响到湍流粘性的计算的。