<div dir="ltr"><div class="gmail_default" style="color:rgb(0,0,0)">For the record, I found out my mistake : the coil was set as a conductor. However, if the current is set, it should be in the non-conductors domain.<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">Best regards.<br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>--<br>Gilles VOGT<br></div>Ingénieur généraliste de l'École Centrale de Lille, docteur en génie électrique<br><div>+33 (0)6 25 98 72 02</div></div></div></div>
<br><div class="gmail_quote">2015-03-17 15:19 GMT+01:00 Gilles Vogt <span dir="ltr"><<a href="mailto:vogtgilles@gmail.com" target="_blank">vogtgilles@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="color:rgb(0,0,0)">Hi,<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">I am using getdp to solve a 2D magneto-thermic problem based on the wiki example (<a href="http://onelab.info/wiki/Magnetodynamics_with_cohomology_conditions" target="_blank">http://onelab.info/wiki/Magnetodynamics_with_cohomology_conditions</a>).<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">At the moment, I am using the A-V formulation, which is in my opinion easier than the T-Ω formulation. However, I have some trouble about the heating power inside my coil domain.<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">In post-pro, my evaluation of the ohmic losses is right in every domain :<br><span style="font-family:monospace,monospace">sigma[]*((Re[{ur}]-Im[w*{a}])*(Re[{ur}]-Im[w*{a}]) + (Im[{ur}]+Re[w*{a}])*(Im[{ur}]+Re[w*{a}]))</span><br></div><div class="gmail_default" style="color:rgb(0,0,0)">where <span style="font-family:monospace,monospace">w=2*pi*f</span>, <span style="font-family:monospace,monospace">{ur}</span> the electric potential and <span style="font-family:monospace,monospace">{a}</span> the magnetic vector potential<br>(which is the same as <span style="font-family:monospace,monospace">sigma[]*SquNorm[Dt[{a}]+{ur}]</span>).<br><br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">However, the following formulation does not produce good results when the coil is included in <span style="font-family:monospace,monospace">Omega_c2 :</span><br><span style="font-family:monospace,monospace">      Galerkin { [ -sigma[]*((Re[{ur}]-Im[w*{a}])*(Re[{ur}]-Im[w*{a}]) + (Im[{ur}]+Re[w*{a}])*(Im[{ur}]+Re[w*{a}]))  , {t} ];<br>      In Omega_c2; Integration CurlCurl; Jacobian Vol;  }</span><br></div><div class="gmail_default" style="color:rgb(0,0,0)">For some reason, the temperature increases far too much in the coil even if the ohmics losses are still OK when post-processed.<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">I have been able to work around this behavior with :<br><span style="font-family:monospace,monospace">Galerkin { [ -qVol[]  , {t} ];<br>      In ***COIL DOMAIN***; Integration CurlCurl; Jacobian Vol;  }</span><br><span style="font-family:monospace,monospace">Galerkin { [ -sigma[]*((Re[{ur}]-Im[w*{a}])*(Re[{ur}]-Im[w*{a}]) + (Im[{ur}]+Re[w*{a}])*(Im[{ur}]+Re[w*{a}]))  , {t} ];<br>      In ***OTHER DOMAINS***; Integration CurlCurl; Jacobian Vol;  }</span><br></div><div class="gmail_default" style="color:rgb(0,0,0)">where qVol=(I/section)^2/sigma...<br></div><div class="gmail_default" style="color:rgb(0,0,0)"><br></div><div class="gmail_default" style="color:rgb(0,0,0)">My results seem right, but I am a bit vexed because I cannot understand where is my mistake... Can anyone pinpoint my error ?<br><br></div><div class="gmail_default" style="color:rgb(0,0,0)">PS : obviously <span style="font-family:monospace,monospace"><br>Galerkin { [ -1/sigma[]* (Re[{js}]*Re[{js}]+Im[{js}]*Im[{js}]), {t} ];</span></div><div class="gmail_default" style="color:rgb(0,0,0)"><span style="font-family:monospace,monospace">In ***COIL DOMAIN***; Integration CurlCurl; Jacobian Vol;  }<br></span></div><div class="gmail_default" style="color:rgb(0,0,0)"><span style="font-family:arial,helvetica,sans-serif">doesn't work either...</span><span style="font-family:monospace,monospace"><br><br></span></div><div class="gmail_default" style="color:rgb(0,0,0)">Best regards,<br clear="all"></div><div><div><div dir="ltr"><div>--<br>Gilles VOGT<br></div><br></div></div></div>
</div>
</blockquote></div><br></div>