Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

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.

How to solve the problem"Attempt to evaluate non-integral power of negative number. Function: ^"

Please login with a confirmed email address before reporting spam

Hello,

I am solving a Drift diffusion model for plasma using general form PDE. However, I am facing the following error,

Failed to find consistent initial values.
Attempt to evaluate non-integral power of negative number.
Function: ^
Failed to evaluate variable.
Variable: mod1.uth, Defined as: (((mod1.Te*(k B const*8))/(pi*mod1.me))^0.5)
Failed to evaluate expression.
Expression: d((0.5*mod1.uth*mod1.u*test(mod1.u))*(dvol),{test0})
Last time step is not converged.


Here, Te= (2/3)*(n/ne)
n, ne is the dependent variables for the problem and their initial ratio is 3. So it is positive. And also k_B_const is Boltzman constant is a constant, me is the mass of electron also a positive quantity.

Can you tell me why it happens?

Thanks is advance.

Rajib

2 Replies Last Post 2017年6月15日 GMT-4 18:09
COMSOL Moderator

Hello Rajib MD

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.


Please login with a confirmed email address before reporting spam

Posted: 9 years ago 2015年12月28日 GMT-5 14:26
Hi,

I'm also wondering about that same behavoiur by using:
T^2*8.227e-07-0.0004936*T+0.153 for defining heat source

Thanks in advance
Niklas

Hi, I'm also wondering about that same behavoiur by using: T^2*8.227e-07-0.0004936*T+0.153 for defining heat source Thanks in advance Niklas

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 2017年6月15日 GMT-4 18:09
Hello Guys,

I'm having a problem to model an adsorption column using Freundlich isotherm in coefficient form PDE module. My initial value of u1 is zero. Could you give me some help?

Attempt to evaluate negative power of zero.
- Function: ^
Failed to evaluate expression.
Expression: d((-(comp1.u2t+S*comp1.u2-9.468679245283018E-5 (0.158*comp1.u1)^0.09)*test(comp1.u2))*(dvol),{test2})
Failed to evaluate Jacobian of expression.
Expression: (-(comp1.u2t+S*comp1.u2-9.468679245283018E-5*(0.158*comp1.u1)^0.09)*test(comp1.u2))*(dvol)
Last time step is not converged.
- Feature: Time-Dependent Solver 1
Regards,

João
Hello Guys, I'm having a problem to model an adsorption column using Freundlich isotherm in coefficient form PDE module. My initial value of u1 is zero. Could you give me some help? Attempt to evaluate negative power of zero. - Function: ^ Failed to evaluate expression. Expression: d((-(comp1.u2t+S*comp1.u2-9.468679245283018E-5 (0.158*comp1.u1)^0.09)*test(comp1.u2))*(dvol),{test2}) Failed to evaluate Jacobian of expression. Expression: (-(comp1.u2t+S*comp1.u2-9.468679245283018E-5*(0.158*comp1.u1)^0.09)*test(comp1.u2))*(dvol) Last time step is not converged. - Feature: Time-Dependent Solver 1 Regards, João

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.