| 查看: 2608 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者syf1180将赠送您 50 个金币 | ||
[求助]
求助 出现无法更新分离步骤 1 的残差。 试图计算负数的实数平方根。这种情况如何处理
|
||
|
- 特征: 稳态求解器 2 (sol1/s2) 无法更新分离步骤 1 的残差。 试图计算负数的实数平方根。 - 函数: sqrt 无法计算算子。 - 算子: nif - 几何: geom1 - 边界: 1-2, 4, 6 无法计算表达式。 - 表达式: nif(comp1.hmnf.d_w_plus_d<(10/(down(comp1.hmnf.Pr)^0.33333333333333331483)),nojac(comp1.hmnf.d_w_plus_d*down(comp1.hmnf.Pr)),comp1.hmnf.d_w_plus_d<sqrt((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt)),nojac((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))-(500/(comp1.hmnf.d_w_plus_d^2))),nojac(((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))+((down(comp1.hmnf.Prt)*log(comp1.hmnf.d_w_plus_d))/comp1.hmnf.kappav))-(((0.5*down(comp1.hmnf.Prt))*(1+log((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt))))/comp1.hmnf.kappav))) 无法计算表达式。 - 表达式: if(comp1.hmnf.hasWF_d,nif(comp1.hmnf.d_w_plus_d<(10/(down(comp1.hmnf.Pr)^0.33333333333333331483)),nojac(comp1.hmnf.d_w_plus_d*down(comp1.hmnf.Pr)),comp1.hmnf.d_w_plus_d<sqrt((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt)),nojac((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))-(500/(comp1.hmnf.d_w_plus_d^2))),nojac(((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))+((down(comp1.hmnf.Prt)*log(comp1.hmnf.d_w_plus_d))/comp1.hmnf.kappav))-(((0.5*down(comp1.hmnf.Prt))*(1+log((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt))))/comp1.hmnf.kappav))),0) 无法计算变量。 - 变量: comp1.hmnf.Tplus_d - 定义为: if(comp1.hmnf.hasWF_d,nif(comp1.hmnf.d_w_plus_d<(10/(down(comp1.hmnf.Pr)^0.33333333333333331483)),nojac(comp1.hmnf.d_w_plus_d*down(comp1.hmnf.Pr)),comp1.hmnf.d_w_plus_d<sqrt((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt)),nojac((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))-(500/(comp1.hmnf.d_w_plus_d^2))),nojac(((15*(down(comp1.hmnf.Pr)^0.66666666666666662966))+((down(comp1.hmnf.Prt)*log(comp1.hmnf.d_w_plus_d))/comp1.hmnf.kappav))-(((0.5*down(comp1.hmnf.Prt))*(1+log((1000*comp1.hmnf.kappav)/down(comp1.hmnf.Prt))))/comp1.hmnf.kappav))),0) 无法计算表达式。 - 表达式: if(comp1.hmnf.hasWF_d,(((-comp1.hmnf.fluid1.qwfc_d)*down(comp1.hmnf.Cp))*(comp1.hmnf.Td-comp1.hmnf.TdWF))/sqrt(1.0E-6+(max(comp1.hmnf.Tplus_d,0)^2)),0) 无法计算变量。 - 变量: comp1.hmnf.fluid1.qwf_d - 定义为: if(comp1.hmnf.hasWF_d,(((-comp1.hmnf.fluid1.qwfc_d)*down(comp1.hmnf.Cp))*(comp1.hmnf.Td-comp1.hmnf.TdWF))/sqrt(1.0E-6+(max(comp1.hmnf.Tplus_d,0)^2)),0) 无法计算表达式。 - 表达式: d(((-comp1.hmnf.fluid1.qwf_d)*test(comp1.hmnf.TdWF-comp1.hmnf.Td))*comp1.hmnf.d,{test@1}) 无法计算表达式的雅可比矩阵。 - 表达式: ((-comp1.hmnf.fluid1.qwf_d)*test(comp1.hmnf.TdWF-comp1.hmnf.Td))*comp1.hmnf.d 无法计算表达式的雅可比矩阵。 - 表达式: (if(comp1.hmnf.hasWF_d,-comp1.hmnf.fluid1.qwf_d*test(comp1.hmnf.TdWF-comp1.hmnf.Td)*comp1.hmnf.d,0))*(dvol_spatial) |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
