COMSOL Forums: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/
Most recent forum messagesFri, 18 Apr 2014 12:40:15 +0000COMSOL Forums: Thermoelectric Effect - Adding the physicshttp://www.comsol.com/shared/images/logos/comsol_logo.gif
http://www.comsol.com/community/forums/general/thread/32340/
Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p117463
Are you using the physics builder guide method?<br />
<br />
To meassure at boundrary, create a boundrary probe and use the insert or replace button in the expression section and find your identifier.Jx Jy or Jz from the drop down tree, choose your used physics then current density. To get current, change the type in the probe setting section to Integral. I write the expressions in variables with aproriate names then export the results to tables through global variable probes to get tables with nice headers.Fri, 18 Apr 2014 12:40:15 +000032340.1397824815.117463Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p117437
Dear Jonatan, <br />
How did you get the current density? Now what I have is voltage but I need current density to calculate the power generation. <br />
Thanks! <br />
JunceThu, 17 Apr 2014 18:03:17 +000032340.1397757797.117437Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p117435
Dear Sergei, <br />
I have a question for the peltier test model: So you will get T and V for the result. How can I get the current (I) in the p and n legs to calculate the power? <br />
Thanks a lot! <br />
JunceThu, 17 Apr 2014 17:55:19 +000032340.1397757319.117435Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109481
Yes, thank you<br />
<br />
Can I ask you what sw you use to write down equation in a fast way?Mon, 28 Oct 2013 17:08:23 +000032340.1382980103.109481Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109479
Dario,<br />
<br />
Yes, this is the right equation. (Please correct my typo for the very last factor d(Vtest)/dy->d(Vtest/dz).<br />
<br />
Regards,<br />
SergeiMon, 28 Oct 2013 16:44:21 +000032340.1382978661.109479Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109477
ok it works perfectly now<br />
<br />
forgive me for my ingenuityMon, 28 Oct 2013 16:40:17 +000032340.1382978417.109477Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109475
Thank you,<br />
<br />
Yes I'm making the assumption that it is isotropic. So the right equation should simply be the one in the attached file?<br />
<br />
Mon, 28 Oct 2013 16:27:58 +000032340.1382977678.109475Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109469
Dario,<br />
<br />
Seebeck weak contribution is placed for current balance, "emdc" interface, in your case.<br />
This term is shown in the attached image.<br />
<br />
There are couple issues with your weak contribution expression:<br />
1. According to Comsol sign convention, Seebeck weak term should be with “minus” sign.<br />
2. Seems you are modeling material with isotropic electric conductivity: sigma_xx=sigma_yy=sigma_zz=sigma. Then, all off-diagonal terms for electric conductivity should be zero. Yet, those terms (containing sigma_xy, sigma_xz, etc.) are non-zero in your expression. <br />
<br />
Regards,<br />
Sergei<br />
Mon, 28 Oct 2013 15:14:24 +000032340.1382973264.109469Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109463
Dear Sergei I'm afraid that I need again your help.<br />
I'm trying to introduce the the Seebeck effect in the 3D model, my equation is the following one:<br />
dVol_emdc*(S*sigma)*((TSx_htgh+TSy_htgh+TSz_htgh)*test(Ex_emdc)+ (TSx_htgh+TSy_htgh+TSz_htgh)*test(Ey_emdc)+(TSx_htgh+TSy_htgh+TSz_htgh)*test(Ez_emdc))<br />
<br />
this equation is accepted by Comsol (not in red)<br />
<br />
I put at 273 K one end of the leg and at 473 the other end, one end is grounded, the other is at a floating potential V.<br />
<br />
When solving there is no voltage generation and all the leg is at 0V <br />
<br />
for the Peltier effect the model works if I put the respective equation in the field "dweack"<br />
<br />
where It is wrong?<br />
where is better to put the equation for the Seebeck?<br />
<br />
Thank you<br />
<br />
DarioMon, 28 Oct 2013 13:45:11 +000032340.1382967911.109463Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109339
Hi Sergei,<br />
<br />
Thanks. In my application, I want to run heat transfer model in steady state and TEC model in transient mode to see change in temperature field after TEC is activated.<br />
<br />
Regards,<br />
SusantThu, 24 Oct 2013 14:08:05 +000032340.1382623685.109339Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109334
Thank you sergei.<br />
<br />
Now i'll try to introduce also the generation effect in the 3D modelThu, 24 Oct 2013 12:21:36 +000032340.1382617296.109334Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109111
Susant,<br />
<br />
It should work without any modifications in the weak contribution node. Just replace Stationary Step by Time Dependent Step. Then click Heat Transfer (ht) root node and make sure that Equation form is selected as “Time dependent”. Then click Electric Currents (ht) root node and make sure that Equation form is selected as “Stationary” (assuming that you are solving AC case).<br />
<br />
Regards,<br />
Sergei<br />
Fri, 18 Oct 2013 16:47:31 +000032340.1382114851.109111Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p109108
Hi Sergei,<br />
<br />
This is with reference to your model, Is there any modification to be done to the model if I want to do time dependent analysis instead of steady state?<br />
<br />
Regards,<br />
SusantFri, 18 Oct 2013 16:12:02 +000032340.1382112722.109108Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108936
<div class="quote"><br />
Thank you again Sergei, can you tell me if this equation is good for a 2d axial symmetry domain?<br />
<br />
dVol_htgh*(S*T)*(Jr_emdc*test(TSr_htgh)+Jz_emdc*test(TSz_htgh)) <br />
<br />
<br />
dario<br />
</div><br />
<br />
Dario,<br />
<br />
Equation seems OK for 2D axial symmetry.<br />
<br />
Regards,<br />
Sergei<br />
Mon, 14 Oct 2013 14:58:23 +000032340.1381762703.108936Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108918
Thank you again Sergei, can you tell me if this equation is good for a 2d axial symmetry domain?<br />
<br />
dVol_htgh*(S*T)*(Jr_emdc*test(TSr_htgh)+Jz_emdc*test(TSz_htgh)) <br />
<br />
<br />
darioMon, 14 Oct 2013 08:02:27 +000032340.1381737747.108918Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108876
<div class="quote"><br />
ok thank you so much!<br />
Therefore also in 2D axial symmetry I don't need to do so?<br />
<br />
You believe that this equation is good as weak term for a 3D model?<br />
dVol_htgh*(S*T)*(-Jx_emdc*test(TSx_htgh)-Jy_emdc*test(TSy_htgh)-Jz_emdc*test(TSz_htgh)) <br />
0<br />
</div><br />
<br />
Dario,<br />
<br />
1. No, 2D plane and 2D axial symmetry are different. Multiplier dVol is Jacobian of coordinate transformation. In 2D plane geometry it reduces to 1 and in 2D axial symmetry case it reduces to 2*pi*r.<br />
<br />
I think, confusion here is due to different notations in 3.5 and 4.3. As far as I remember, in 3.5 multiplier dVol is always present in weak form notation (it was a while ago I use 3.5, though).<br />
<br />
<br />
2. I believe, your weak term is OK for 3D model in version 3.5a.<br />
<br />
Regards,<br />
Sergei<br />
Fri, 11 Oct 2013 10:13:33 +000032340.1381486413.108876Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108865
ok thank you so much!<br />
Therefore also in 2D axial symmetry I don't need to do so?<br />
<br />
You believe that this equation is good as weak term for a 3D model?<br />
dVol_htgh*(S*T)*(-Jx_emdc*test(TSx_htgh)-Jy_emdc*test(TSy_htgh)-Jz_emdc*test(TSz_htgh)) <br />
0<br />
<br />
Fri, 11 Oct 2013 06:51:18 +000032340.1381474278.108865Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108856
<div class="quote"><br />
Dear Sergei,<br />
<br />
I'd like to ask you some questions:<br />
<br />
1)why in the 2D model you gave me you write:<br />
(-k_ht*Tx+P*Jx_dc)*test(Tx)+(-k_ht*Ty+P*Jy_dc)*test(Ty) without volume integration?<br />
<br />
2)where you have declared "P" ?<br />
<br />
Attached you find my 2D axial symmetric model, there I can find only a heat release at one junction but no heat absorption at the other. Can you explain me where I am wrong and why?<br />
<br />
Note that my file is more than 10 MB so I was forced to split it in 3 parts with winrar. Also i was forced to rename *.rar files as *.mph files so you have to rename the files in *.rar (otherwise the forum does not allow me to attach)<br />
<br />
Hope this is not too boring for you...<br />
<br />
Thank you soo much!<br />
</div><br />
<br />
Dario,<br />
<br />
1. I didn’t include volume integration because out-of-plane thickness was assumed to be one in that particular 2D example.<br />
<br />
2. Peltier coefficient “P” is defined under Model Variables on the domains level.<br />
<br />
3. I was not able to open your *.rar files.<br />
<br />
Regards,<br />
Sergei<br />
Thu, 10 Oct 2013 17:59:26 +000032340.1381427966.108856Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108834
<div class="quote"><br />
<div class="quote"><br />
Thank you so much Sergei, yesterday I managed to do it by myself in 2D axyal symmetric and 3D thanks to the pdf file you have attached. With the 3D model I make a comparison with my previous model in PDE form obtaining exactly the same result. Now I can also check if i did it correctly the 2D model.<br />
</div><br />
<br />
Dear Jonathan,<br />
<br />
Could you tell me how did you add the weak equation to your 3D model? or How did you make it work in a 3D model?<br />
Thank you!<br />
<br />
Ying<br />
</div><br />
<br />
I'm Dario not Jonathan :D<br />
<br />
This is the equation I put in my 3D model:<br />
<br />
dVol_htgh*(S*T)*(-Jx_emdc*test(TSx_htgh)-Jy_emdc*test(TSy_htgh)-Jz_emdc*test(TSz_htgh)) <br />
<br />
please wait for Sergei validation<br />
<br />
<br />
Thu, 10 Oct 2013 08:36:32 +000032340.1381394192.108834Re: Thermoelectric Effect - Adding the physics
http://www.comsol.com/community/forums/general/thread/32340/#p108832
Dear Sergei,<br />
<br />
I'd like to ask you some questions:<br />
<br />
1)why in the 2D model you gave me you write:<br />
(-k_ht*Tx+P*Jx_dc)*test(Tx)+(-k_ht*Ty+P*Jy_dc)*test(Ty) without volume integration?<br />
<br />
2)where you have declared "P" ?<br />
<br />
Attached you find my 2D axial symmetric model, there I can find only a heat release at one junction but no heat absorption at the other. Can you explain me where I am wrong and why?<br />
<br />
Note that my file is more than 10 MB so I was forced to split it in 3 parts with winrar. Also i was forced to rename *.rar files as *.mph files so you have to rename the files in *.rar (otherwise the forum does not allow me to attach)<br />
<br />
Hope this is not too boring for you...<br />
<br />
Thank you soo much!<br />
Thu, 10 Oct 2013 08:33:48 +000032340.1381394028.108832