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.

E distribution of nonlinear material under DC

Please login with a confirmed email address before reporting spam

I want to know the E distribution of nonlinear material under DC.
Electric resistivity is a function of electric field in the material.
For example: Rou=3*E^12*exp(-0.215*E)

I established a simple rectangle model, with DC supply and nonlinear equation in the subdomain setting.
But error came out as below,

Messages: Undefined value found solvermodel

Stack trace:
at solvermodel.cpp, row 1837, (undefined value found in the equation residual vector, There are 417 degrees of freedom giving NaN/Inf in the vector for the variable V at coordinates: .......................)

At the first time, I thought maybe the boundary setting did some wrong. But to my model, no difficult setting in it.
So what should I do to solve this problem and get the answer?

Thanks a lot.

Janet


4 Replies Last Post 2010年10月4日 GMT-4 05:10

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月1日 GMT-4 11:30
Janet,

did you make sure you used the right variable in your expression? Check variable names according to the application mode you use. The expression you showed in the posting will definitely not work.

Best regards
Edgar
Janet, did you make sure you used the right variable in your expression? Check variable names according to the application mode you use. The expression you showed in the posting will definitely not work. Best regards Edgar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月1日 GMT-4 22:14
Hi, Janet,

I think maybe your definition of electric field E is not correct. If like, you can post your mph file here. We can try to solve this problem.
Hi, Janet, I think maybe your definition of electric field E is not correct. If like, you can post your mph file here. We can try to solve this problem.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月4日 GMT-4 01:31
Thanks both of you.
Actually the electric field E I typed in model is `normE_embc`.

And I attched my model here. (A basic model)
Hope you can get me out of the problem.

Janet
Thanks both of you. Actually the electric field E I typed in model is `normE_embc`. And I attched my model here. (A basic model) Hope you can get me out of the problem. Janet


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月4日 GMT-4 05:10
Janet,

I think you must check your expression. Your boundary settings generate a field of 710 kV/m. The exponential term in your expression exp(0.0152*710000) will result in a numerical error. This is a huge number. Also I assume, that your denominator should be E1*1e16 and not E1^16. 710000^16 is something like 1e93.
I presume you have an estimate what value the conductivity should be. Check if your expression gets you close to this estimate.
Your expression is numerically extremely unstable.

Best regards
Edgar

Janet, I think you must check your expression. Your boundary settings generate a field of 710 kV/m. The exponential term in your expression exp(0.0152*710000) will result in a numerical error. This is a huge number. Also I assume, that your denominator should be E1*1e16 and not E1^16. 710000^16 is something like 1e93. I presume you have an estimate what value the conductivity should be. Check if your expression gets you close to this estimate. Your expression is numerically extremely unstable. Best regards Edgar

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.