请问有哪位大侠用chemkin模拟过超音速喷管中的燃烧吗？

cyclone0907

燃料在燃烧室中燃烧后（此时压力5MPa，燃烧后温度3431.46k）进入拉法尔喷管，小弟想要用chemkin4.0中的柱塞流反应器模型模拟燃气在喷管中流动时各燃烧产物（包括自由基）浓度的变化情况。但是模拟过程中遇到一些问题：对于一个给定形状的喷管（喷管直径随轴向距离而变），柱塞流反应器的进口速度大于一定程度后便会出现错误  TASK INTERRUPTED: TROUBLE IN THE TRANSIENT SOLVER DDASPK；THE NONLINEAR SYSTEM SOLVER CANNOT CONVERGE。如果进口速度比该值小的时候虽然可以计算，但是燃气在喷管中并没有达到超声速，根据速度、压力和温度分布可以看出扩张段中是一个扩压的过程，说明此时喉部没有达到声速。难道chemkin的PFR模型不能计算跨声速吗？
我试着将收敛段与扩张段分开计算。同样采用PFR模型，在计算收敛段时，进口速度大于一定程度后便会出现错误。计算扩张段时，进口速度在一个特定区间内的时候便会出现错误。小于该区间的进口速度可以计算，但是燃气在喷管中并没有达到超声速，扩张段中是一个扩压的过程。大于该区间的进口速度也可以计算，而且燃气在喷管中达到了超声速，扩张段中是膨胀的过程。
现在的问题是为什么在声速附近的时候程序无法继续计算？这个跟喷管形状有关系吗？计算中的喷管面积比很大，达到100.
有没有哪位大神算过此类问题，小弟在此谢过了！