超声速湍流平板计算问题

[云飞]

各位前辈和大佬，我在Fortran编程时将SA湍流模型加入原有的层流求解器时，对M=2的超声速平板进行算例验证。
计算求得的边界层涡粘系数量级与参考解一致，但是分布比参考解薄一些，边界层厚度也约为参考解的2/3。
检查了湍流模型的源项没有发现问题，用到的最小壁面距离也是无量纲的，不知道问题可能在哪？希望各位不吝赐教，感谢！

[line]

网格有没有再加密，对流通量精度是多少，一阶精度会比二阶低很多

[云飞]

网格加密到Y+<1，计算结果变化不大，对流通量用二阶MUSCL, AUSM格式