（求助）请教流体性质设置的问题

hanzhang654321

大家好，我在做热流耦合的问题，是流体受热在自重力情况下的温度分布。遇到一个问题，已经困惑几天。

我的问题是，如果设置流体的密度不变，即FLDATA7,PROT,DENS,liquid，fldata8,nomi,dens,999.8
则程序可以运行，可以得到结果（FLDATA1,SOLU,TEMP,1  FLDATA1,SOLU,flow,1   同时计算温度场，流体场）

但是如果设置密度可变，即加上一条FLDATA10,COF2,DENS,-0.416，

则程序一直出现“unsymmetric cg method failed catastrophically
coefficient matrix has a negative diagonal. probably indicates divergent solution.”的错误提示。

请问怎么处理才使其能够计算呢？

拜托哪位朋友知道的帮忙解答，不胜感谢。