Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
Division by zero error
Posted 2019年5月6日 GMT-4 13:14 Plasma Physics Version 5.4 0 Replies
Please login with a confirmed email address before reporting spam
Hello. I try to model He-Ne plasma. And during the calculation I received the following error:
Division by zero. - Function: / Failed to evaluate temporary symbolic derivative variable. - Variable: comp1.plas.R_wNe_1p@VDN${real@7} - Defined as: (((comp1.plas.c_wecomp1.plas.kf_4)((1/comp1.plas.M_wNe)((comp1.plas.wNe((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p))))))))+(comp1.plas.rho(-exp(real(comp1.WHe_1p)))))))+(((comp1.plas.c_wNecomp1.plas.kf_14)((1/comp1.plas.M_wHe_1p)((comp1.plas.rhoexp(real(comp1.WHe_1p)))+(comp1.plas.wHe_1p((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p)))))))))))+(comp1.plas.c_wHe_1p(comp1.plas.kf_14((1/comp1.plas.M_wNe)((comp1.plas.wNe((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p))))))))+(comp1.plas.rho(-exp(real(comp1.WHe_1p))))))))))+(exp((-comp1.plas.iota)log(comp1.plas.n_wNe_1p))((-comp1.plas.iota)((1/comp1.plas.n_wNe_1p)((1/comp1.plas.M_wNe_1p)((comp1.plas.wNe_1pN_A_const)((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p)))))))))))) Failed to evaluate expression. - Expression: dvol_spatial(r(pi((-2)((d(comp1.WNe_1p,t)(comp1.plas.wNe_1p((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p)))))))))-(comp1.plas.M_wNe_1p((((comp1.plas.c_wecomp1.plas.kf_4)((1/comp1.plas.M_wNe)((comp1.plas.wNe((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p))))))))+(comp1.plas.rho(-exp(real(comp1.WHe_1p)))))))+(((comp1.plas.c_wNecomp1.plas.kf_14)((1/comp1.plas.M_wHe_1p)((comp1.plas.rhoexp(real(comp1.WHe_1p)))+(comp1.plas.wHe_1p((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p)))))))))))+(comp1.plas.c_wHe_1p(comp1.plas.kf_14((1/comp1.plas.M_wNe)((comp1.plas.wNe((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p))))))))+(comp1.plas.rho(-exp(real(comp1.WHe_1p))))))))))+(exp((-comp1.plas.iota)log(comp1.plas.n_wNe_1p))((-comp1.plas.iota)((1/comp1.plas.n_wNe_1p)((1/comp1.plas.M_wNe_1p)((comp1.plas.wNe_1pN_A_const)((1/(R_constcomp1.plas.pes1.minput_temperature))(comp1.plas.pes1.minput_pressure(((-1)/(((((((comp1.plas.wNe/comp1.plas.M_wNe)+(comp1.plas.wNes/comp1.plas.M_wNes))+(comp1.plas.wNe_1p/comp1.plas.M_wNe_1p))+(comp1.plas.wHe/comp1.plas.M_wHe))+(comp1.plas.wHes/comp1.plas.M_wHes))+(comp1.plas.wHe_1p/comp1.plas.M_wHe_1p))^2))(((1/comp1.plas.M_wNe)(-exp(real(comp1.WHe_1p))))+((1/comp1.plas.M_wHe_1p)exp(real(comp1.WHe_1p)))))))))))))))))) Failed to evaluate Jacobian of expression. - Expression: dvol_spatial(r(pi((-2)(((comp1.plas.rhocomp1.plas.wNe_1p)d(comp1.WNe_1p,t))-(comp1.plas.M_wNe_1pcomp1.plas.R_wNe_1p))))) Failed to evaluate Jacobian of expression. - Expression: (-2(comp1.plas.rhocomp1.plas.wNe_1pd(comp1.WNe_1p,t)test(comp1.WNe_1p)-comp1.plas.jr_wNe_1ptest(d(comp1.WNe_1p,r))-comp1.plas.jz_wNe_1ptest(d(comp1.WNe_1p,z))-comp1.plas.R_wNe_1pcomp1.plas.M_wNe_1ptest(comp1.WNe_1p))pir)*(dvol_spatial)
I tried to solve it, but i don't understand, which Ne+ parameter is equal zero. Can you help me? How can understand this expression: comp1.plas.R_wNe_1p@VDN${real@7} I understand, that comp1 - component 1, plas - plasma module. But what is R_wNe_1p@VDN${real@7}?
Hello Даниил Соколов
Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.
If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.