Discussion Forum

Posted: 1 decade ago 2011年1月11日 GMT-5 11:05

Hi are you sure about those "square brackets" isnt' it (10^-14)*(exp((ec.normE/10000000)^15)+1). On the other side your function is rather steep, are you sure the binary floating point numeber scheme can represent it correctly ? -- Good luck Ivar

Posted: 1 decade ago 2011年1月12日 GMT-5 10:17

Hi Ivar, Thanks for reply. Yea I have noticed square brackets is not an allowed syntax and hence removed. I have reduced the power of the variable in the function from 15 to 2 to reduce the steepness of the function. Changed the geometry dimension unit (and hence field magnitude reduces accordingly) to take care of the floating point scheme to some extent(if I had understood right). Any better way to take care of floating point scheme?. After the above modifications still not able to execute. The following is the error message. Added a small value to ec.normE in the expression to resolve negative power of zero. Still not working. Any idea?? Attempt to evaluate negative power of zero. - Function: ^ Failed to evaluate temporary symbolic derivative variable. - Variable: mod1.ec.normE@VDN$mod1.Vy - Defined as: ((0.5*((realdot(mod1.ec.Ez,mod1.ec.Ez)+(realdot(mod1.ec.Ey,mod1.ec.Ey)+realdot(mod1.ec.Ex,mod1.ec.Ex)))^-0.5))*((-conj(mod1.ec.Ey))+(-conj(mod1.ec.Ey)))) Failed to evaluate temporary symbolic derivative variable. - Variable: mod1.ec.sigmayy@VDN$mod1.Vy - Defined as: ((((0.5*((realdot(mod1.ec.Ez,mod1.ec.Ez)+(realdot(mod1.ec.Ey,mod1.ec.Ey)+realdot(mod1.ec.Ex,mod1.ec.Ex)))^-0.5))*((-conj(mod1.ec.Ey))+(-conj(mod1.ec.Ey)))))*dmod1.mat1.def.sigma_dec.normE(mod1.ec.normE)) Failed to evaluate temporary symbolic derivative variable. - Variable: mod1.ec.Jiy@VDN$mod1.Vy - Defined as: ((-mod1.ec.sigmayy)+((((((0.5*((realdot(mod1.ec.Ez,mod1.ec.Ez)+(realdot(mod1.ec.Ey,mod1.ec.Ey)+realdot(mod1.ec.Ex,mod1.ec.Ex)))^-0.5))*((-conj(mod1.ec.Ey))+(-conj(mod1.ec.Ey)))))*dmod1.mat1.def.sigma_dec.normE(mod1.ec.normE)))*mod1.ec.Ey)) Failed to evaluate temporary symbolic derivative variable. - Variable: mod1.ec.Jy@VDN$mod1.Vy - Defined as: (((-mod1.ec.sigmayy)+((((((0.5*((realdot(mod1.ec.Ez,mod1.ec.Ez)+(realdot(mod1.ec.Ey,mod1.ec.Ey)+realdot(mod1.ec.Ex,mod1.ec.Ex)))^-0.5))*((-conj(mod1.ec.Ey))+(-conj(mod1.ec.Ey)))))*dmod1.mat1.def.sigma_dec.normE(mod1.ec.normE)))*mod1.ec.Ey))) Failed to evaluate expression. - Expression: d(d((mod1.ec.Jx*test(mod1.Vx)+mod1.ec.Jy*test(mod1.Vy)+mod1.ec.Jz*test(mod1.Vz))*(dvol),{test@1}),mod1.Vy) Regards, Venkat [QUOTE] Hi are you sure about those "square brackets" isnt' it (10^-14)*(exp((ec.normE/10000000)^15)+1). On the other side your function is rather steep, are you sure the binary floating point numeber scheme can represent it correctly ? -- Good luck Ivar [/QUOTE]

Attachments:

• analytical function usage.docx

Posted: 1 decade ago 2011年1月13日 GMT-5 01:47

Hi I believe its because of a scaling issue, and a numerical limitation, try to smoothen even more. Do not forget that ratios are limited to "eps" and rather sqrt(eps) -- Good luck Ivar

Posted: 1 decade ago 2011年1月17日 GMT-5 07:42

Hi Ivar, Thanks for the reply. I have tried various ways and even a simple smooth function for sigma taking care of scaling issue and eps. Still not able to solve. Please check the error message in the model (same as before).Find the model as an attachment(.mph). I have attached simple model for conveying the modelling issues. If we disable analytic fn and give sigma of air as some value say 10E-4 it successfully solves. But whatever may be the sigma fn in terms of E field it does not solve. All the material properties in the attached file are not realistic values but meaningful. Any idea how to resolve this? second issue is about vector addition. For constant sigma, comsol gives the results. Then I had extracted data of ec.Ex,ec.Ey, ec.Ez, ec.normE and ec.Ex+ec.Ey. Here Ex,y and z are vectors.If we add then it should be not simple scalar addition.But for ec.Ex+ec.Ey comsol simply gives output addition of Ex and Ey components like scalars. Where is the confusion?

Attachments:

• uisng analytic function for E-field dependent sigma.mph

Posted: 1 decade ago 2011年1月17日 GMT-5 08:09

Hi I do not have access to COMSOL just now, not in my office, but I'll try to check tomorrow. One thing though you talk about ec.Ex Ec.Ey ... for me in COMSOL notation this is the components of the gradient of E , hence dE(x,y,z)/dx , dE(x,y,z)/dy ? are there a naming interaction between COMSOl internal variables and your defined ones or am I simply missing something ? And Comsol only give you access to scalars, or the separate scalar components of a vector, no ? -- Good luck Ivar

Posted: 1 decade ago 2011年1月24日 GMT-5 02:27

Hi I found 2 way to avoid an error 1) set a better starting potential, if you have 0.5V on the upper electrode and 0=gnd on the base, you can also simply say the Initial conditions voltage across your air slab is 0.5[V]*(1-(z/1[m]) and it will solve. Then why ? the reason 2) is that the conductivity function you have defined is "=0" if there is no voltage field, and you use "0" as initial conditions. The best is to give a little offset to the "air" conductivity "eps" is enough, but "1[S/m]" is probably more realistic By the way, you can leave COMSOL look after the derivative "automatically" in your case -- Good luck Ivar

Posted: 1 decade ago 2011年1月24日 GMT-5 09:00

Hi Ivar, Thanks for the reply. What you had suggested is working. It seems even 1*(z-1) as an initial condition for the potential is also working fine. Last week I have learnt from the support that we need to define initial voltage as some 1e-6*sqrt(x^2+y^2+z^2) to avoid the singularity. Mvh Venkat.

Posted: 1 decade ago 2011年9月23日 GMT-4 09:54

[QUOTE] Hi I found 2 way to avoid an error 1) set a better starting potential, if you have 0.5V on the upper electrode and 0=gnd on the base, you can also simply say the Initial conditions voltage across your air slab is 0.5[V]*(1-(z/1[m]) and it will solve. Then why ? the reason 2) is that the conductivity function you have defined is "=0" if there is no voltage field, and you use "0" as initial conditions. The best is to give a little offset to the "air" conductivity "eps" is enough, but "1[S/m]" is probably more realistic By the way, you can leave COMSOL look after the derivative "automatically" in your case -- Good luck Ivar [/QUOTE] Hi Ivar, thanks for the suggestions, the 0.5V initial voltage also worked well for me. But I don't understand here the usage of that 'z'. What does it mean? Why not just '0.5[V]' but that '0.5[V]*(1-(z/1[m]))' ? thanks again Moyan

Posted: 1 decade ago 2011年9月23日 GMT-4 10:26

Hi "z" is the height parameter a voltage between two electrodes spaced by "1" would be something like V0*(1-z) the voltage is V0 but the voltage at a height of a point in the domain would be V0*(d0-z)/d0 where d0[m] is the electrode separation distance I assumed to be "d0=1"(careful with the offsets where to put your geometry in absolute coordinates though) -- Good luck Ivar

Posted: 8 years ago 2017年4月11日 GMT-4 02:05

Hi Ivar , I am trying fluid structure interaction .. I want to know why the particle is getting deformed as it progresses down the channel .. Shall i attach my file? Can you help?