大佬们帮忙看看这个是什么问题？

caigen

快疯了。调试了快2周了，完全不知到怎么搞。

找不到一致的初始值。
试图计算负数的实数平方根。
Function: sqrt
未能计算变量。
Variable: comp2.chem.eta C, Defined as: (((2.669E-6*sqrt((unit mol cf*(comp2.chem.M C*(comp2.chem.T*1000)))/(unit K cf*unit kg cf)))*unit Pa cf)*unit s cf)/(((((1.16145/((comp2.chem.T/comp2.chem.epsilonkb C)^0.14874))+(0.52487*exp((comp2.chem.T*(-0.7732))/comp2.chem.epsilonkb C)))+(2.16178*exp((comp2.chem.T*(-2.4378700000000002))/comp2.chem.epsilonkb C)))+((((4.998E-40*((comp2.chem.mu C/(unit C cf*unit m cf))^4))*(unit kg cf^2))*(unit m cf^10))/((comp2.chem.T*(comp2.chem.epsilonkb C*(((1.0E10*comp2.chem.sigma C)^6)*(k B const^2))))*(unit s cf^4))))*(((1.0E10*comp2.chem.sigma C)/unit m cf)^2))
未能计算变量。
Variable: comp2.chem.eta, Defined as: (((comp2.chem.eta C/(1+((((((comp2.chem.c CH4*((1+(((comp2.chem.eta C/comp2.chem.eta CH4)^0.5)*((comp2.chem.M CH4/comp2.chem.M C)^0.25)))^2))/((1+(comp2.chem.M C/comp2.chem.M CH4))^0.5))+((comp2.chem.c H2*((1+(((comp2.chem.eta C/comp2.chem.eta H2)^0.5)*((comp2.chem.M H2/comp2.chem.M C)^0.25)))^2))/((1+(comp2.chem.M C/comp2.chem.M H2))^0.5)))+((comp2.chem.c a*((1+(((comp2.chem.eta C/comp2.chem.eta a)^0.5)*((comp2.chem.M a/comp2.chem.M C)^0.25)))^2))/((1+(comp2.chem.M C/comp2.chem.M a))^0.5)))*0.35355339059327379)/(comp2.chem.c*max(eps,comp2.chem.c C/comp2.chem.c)))))+(comp2.chem.eta CH4/(1+((((((comp2.chem.c C*((1+(((comp2.chem.eta CH4/comp2.chem.eta C)^0.5)*((comp2.chem.M C/comp2.chem.M CH4)^0.25)))^2))/((1+(comp2.chem.M CH4/comp2.chem.M C))^0.5))+((comp2.chem.c H2*((1+(((comp2.chem.eta CH4/comp2.chem.eta H2)^0.5)*((comp2.chem.M H2/comp2.chem.M CH4)^0.25)))^2))/((1+(comp2.chem.M CH4/comp2.chem.M H2))^0.5)))+((comp2.chem.c a*((1+(((comp2.chem.eta CH4/comp2.chem.eta a)^0.5)*((comp2.chem.M a/comp2.chem.M CH4)^0.25)))^2))/((1+(comp2.chem.M CH4/comp2.chem.M a))^0.5)))*0.35355339059327379)/(comp2.chem.c*max(eps,comp2.chem.c CH4/comp2.chem.c))))))+(comp2.chem.eta H2/(1+((((((comp2.chem.c C*((1+(((comp2.chem.eta H2/comp2.chem.eta C)^0.5)*((comp2.chem.M C/comp2.chem.M H2)^0.25)))^2))/((1+(comp2.chem.M H2/comp2.chem.M C))^0.5))+((comp2.chem.c CH4*((1+(((comp2.chem.eta H2/comp2.chem.eta CH4)^0.5)*((comp2.chem.M CH4/comp2.chem.M H2)^0.25)))^2))/((1+(comp2.chem.M H2/comp2.chem.M CH4))^0.5)))+((comp2.chem.c a*((1+(((comp2.chem.eta H2/comp2.chem.eta a)^0.5)*((comp2.chem.M a/comp2.chem.M H2)^0.25)))^2))/((1+(comp2.chem.M H2/comp2.chem.M a))^0.5)))*0.35355339059327379)/(comp2.chem.c*max(eps,comp2.chem.c H2/comp2.chem.c))))))+(comp2.chem.eta a/(1+((((((comp2.chem.c C*((1+(((comp2.chem.eta a/comp2.chem.eta C)^0.5)*((comp2.chem.M C/comp2.chem.M a)^0.25)))^2))/((1+(comp2.chem.M a/comp2.chem.M C))^0.5))+((comp2.chem.c CH4*((1+(((comp2.chem.eta a/comp2.chem.eta CH4)^0.5)*((comp2.chem.M CH4/comp2.chem.M a)^0.25)))^2))/((1+(comp2.chem.M a/comp2.chem.M CH4))^0.5)))+((comp2.chem.c H2*((1+(((comp2.chem.eta a/comp2.chem.eta H2)^0.5)*((comp2.chem.M H2/comp2.chem.M a)^0.25)))^2))/((1+(comp2.chem.M a/comp2.chem.M H2))^0.5)))*0.35355339059327379)/(comp2.chem.c*max(eps,comp2.chem.c a/comp2.chem.c)))))
未能计算变量。
Variable: comp2.spf.mu, Defined as: comp2.chem.eta
未能计算变量。
Variable: comp2.spf.tauNSViscTerm
未能计算表达式。
Expression: nojac(1/sqrt((((comp2.spf.rho^2)*max(comp2.spf.dt2Inv,0))+(4*comp2.spf.tauNSViscTerm))+((4*(comp2.spf.rho^2))*comp2.spf.emetric u)))
未能计算变量。
Variable: comp2.spf.tau m
未能计算变量。
Variable: comp2.spf.u primer
未能计算变量。
Variable: comp2.ht.u primer
未能计算变量。
Variable: comp2.ht.res prime
未能计算表达式。
Expression: nojac(max(((0.5*sqrt(max(comp2.ht.res prime^2,0)*max((comp2.ht.Cp^2)*((((comp2.Tr*comp2.ht.Jijrr)+(comp2.Tz*comp2.ht.Jijrz))^2)+(((comp2.Tr*comp2.ht.Jijzr)+(comp2.Tz*comp2.ht.Jijzz))^2)),0)))-(max(comp2.ht.res prime^2,0)*comp2.ht.tau prime))/max((comp2.ht.Cp^2)*((((comp2.Tr*comp2.ht.Jijrr)+(comp2.Tz*comp2.ht.Jijrz))^2)+(((comp2.Tr*comp2.ht.Jijzr)+(comp2.Tz*comp2.ht.Jijzz))^2)),eps),0))
未能计算临时符号衍生变量。
Variable: comp2.ht.crosswindVDN${test12}
未能计算表达式。
Expression: d((2*comp2.ht.crosswind*pi*r)*(dvol spatial),{test12})
未能计算表达式的 Jacobian 矩阵。
Expression: (2*comp2.ht.crosswind*pi*r)*(dvol spatial)
最后一个时间步不收敛。

caigen

某个地方的值给的太大了。我是算化学反应的，化学反应速率太快了。