comsol计算出现错误

试图计算负数的实数平方根。

  - 函数: sqrt

无法计算算子。

  - 算子: nif

  - 几何: geom1

  - 边界: 3, 16

无法计算表达式。

  - 表达式: nif(comp1.spf.d_w_plus_d<(10/(down(comp1.nitf1.Pr)^0.33333333333333331483)),nojac(comp1.spf.d_w_plus_d*down(comp1.nitf1.Pr)),comp1.spf.d_w_plus_d<sqrt((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt)),nojac((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))-(500/(comp1.spf.d_w_plus_d^2))),nojac(((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))+((down(comp1.nitf1.Prt)*log(comp1.spf.d_w_plus_d))/comp1.spf.kappav))-(((0.5*down(comp1.nitf1.Prt))*(1+log((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt))))/comp1.spf.kappav)))

无法计算表达式。

  - 表达式: if(comp1.spf.hasWF_d,nif(comp1.spf.d_w_plus_d<(10/(down(comp1.nitf1.Pr)^0.33333333333333331483)),nojac(comp1.spf.d_w_plus_d*down(comp1.nitf1.Pr)),comp1.spf.d_w_plus_d<sqrt((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt)),nojac((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))-(500/(comp1.spf.d_w_plus_d^2))),nojac(((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))+((down(comp1.nitf1.Prt)*log(comp1.spf.d_w_plus_d))/comp1.spf.kappav))-(((0.5*down(comp1.nitf1.Prt))*(1+log((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt))))/comp1.spf.kappav))),0)

无法计算变量。

  - 变量: comp1.nitf1.Tplus_d

  - 定义为: if(comp1.spf.hasWF_d,nif(comp1.spf.d_w_plus_d<(10/(down(comp1.nitf1.Pr)^0.33333333333333331483)),nojac(comp1.spf.d_w_plus_d*down(comp1.nitf1.Pr)),comp1.spf.d_w_plus_d<sqrt((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt)),nojac((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))-(500/(comp1.spf.d_w_plus_d^2))),nojac(((15*(down(comp1.nitf1.Pr)^0.66666666666666662966))+((down(comp1.nitf1.Prt)*log(comp1.spf.d_w_plus_d))/comp1.spf.kappav))-(((0.5*down(comp1.nitf1.Prt))*(1+log((comp1.spf.kappav*1000)/down(comp1.nitf1.Prt))))/comp1.spf.kappav))),0)

无法计算表达式。

  - 表达式: if(comp1.spf.hasWF_d,(((-comp1.nitf1.qwfc_d)*down(comp1.ht.Cp))*(comp1.ht.Td-comp1.nitf1.TdWF))/sqrt(1.0E-6+(max(comp1.nitf1.Tplus_d,0)^2)),0)

无法计算变量。

  - 变量: comp1.nitf1.qwf_d

  - 定义为: if(comp1.spf.hasWF_d,(((-comp1.nitf1.qwfc_d)*down(comp1.ht.Cp))*(comp1.ht.Td-comp1.nitf1.TdWF))/sqrt(1.0E-6+(max(comp1.nitf1.Tplus_d,0)^2)),0)

无法计算表达式。

  - 表达式: d(((-comp1.nitf1.qwf_d)*test(comp1.nitf1.TdWF-comp1.ht.Td))*comp1.ht.d,{test@1})

无法计算表达式的雅可比矩阵。

  - 表达式: ((-comp1.nitf1.qwf_d)*test(comp1.nitf1.TdWF-comp1.ht.Td))*comp1.ht.d

无法计算表达式的雅可比矩阵。

  - 表达式: (if(comp1.spf.hasWF_d,-comp1.nitf1.qwf_d*test(comp1.nitf1.TdWF-comp1.ht.Td)*comp1.ht.d,0))*(dvol_spatial)