求教采用同位网格simple方法计算翼形分离流的问题

runheart

我参考陶老师的《计算传热学的近代进展》p197－非结构化网格上的simple算法。编写了计算翼型的二维程序（不可压缩流)。不过对原方法做了一点修改：1.结构网格计算。2.单独推导压力的计算方程。再用压力修正方程对速度压力修正。计算发现有些问题想求教一下这里的同仁。
1，对接边界M－N应当怎样处理？
   TDMA方法交替迭代求解，即从A-B联立求解，再D-C；G-F；D-E求解。这样尾迹区的求解被分成了两部分，对网格对接边界M-N处采用了插值边界条件，初步地计算发现这种尾迹区分离式求解会导致对接边界两侧压力不连续，并且在高Re数下导致振荡发散。后来，我又参考了RHIE的计算程序，RHIE对尾迹区整体求解，即TDMA方法对C-E联立求解取代C-D,D-E的分离求解,并且RHIE的注释强调了此问题，不过我的计算程序采用这样的跌代求解立刻就发散。上图计算的压力云图是对接边界采用插值边界条件计算的，在边界连接处压力还是不很连续。我该怎样处理对接边界？
2，如果不计湍流模型，只用层流计算（用quick格式），对naca0012翼型来说，Re数和攻角能达到多大？我的计算程序在Re=10000左右对15度攻角都能算。而Re再高就发散了。同仁有没有相关的计算经验，使Re=10e6对15度攻角采用层流方程也能计算？

supersonic

参考DUNS的多块网格情况下的ADI系数矩阵构造