<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>Dear Prof. Sabariego and All,</div><div><br></div><div>I looked into the new indheat.pro benchmark and change several lines to use AV formulation, MagDynAV formulation.</div><div><br></div><div>The diff file is this.</div><div><br></div><div>325c325,326<br>< { Name h; Type Local; NameOfSpace HSpace; }<br>---<br>> { Name a; Type Local; NameOfSpace ASpace; }<br>> { Name e; Type Local; NameOfSpace ESpace; }<br>333c334<br>< Galerkin { [ -0.5/sigma[]*<h>[Re[{d h}]*Re[{d h}] + Im[{d h}]*Im[{d h}]], {t} ];<br>---<br>> Galerkin { [ -0.5/sigma[]*<a>[Re[Dt[-{a}]+{e}]*Re[Dt[-{a}]+{e}] + Im[Dt[-{a}]+{e}]*Im[Dt[-{a}]+{e}]], {t} ];<br>378c379<br>< { Name A; NameOfFormulation MagDynTO;<br>---<br>> { Name A; NameOfFormulation MagDynAV;<br>478c479<br>< { Name p; Value{ Local{ [ 1./sigma[]*( Re[{d h}]*Re[{d h}] + Im[{d h}]*Im[{d h}] ) ] ;<br>---<br>> { Name p; Value{ Local{ [ -0.5/sigma[]*<a>[Re[Dt[-{a}]+{e}]*Re[Dt[-{a}]+{e}] + Im[Dt[-{a}]+{e}]*Im[Dt[-{a}]+{e}]] ] ;<br><br></div><div>The simulation results are very different from the original TO formulation.</div><div>My concerning points are</div><div> * In AV formulation, “a" and “e" are both complex numbers so <a>[] should be something like <a,e>[].</div><div> * The operator “Dt” would work as expected in the TheDyn formulation, it is transient formulation.</div><div><br></div><div>I tried <a,e>[], but this causes a syntax error.</div><div>I have a feeling I am closing to my goal. I appreciate your kind helps.</div><div>Thank you very much in advance.</div><div><br></div><div>Best Regards,</div><div><br></div><br><div><div>2017/10/18 0:04、Ruth Vazquez Sabariego <<a href="mailto:ruth.sabariego@kuleuven.be">ruth.sabariego@kuleuven.be</a>> のメール：</div><br class="Apple-interchange-newline"><blockquote type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Dear ABE Hiroshi, </div>
<div class=""><br class="">
</div>
<div class="">Using the T-O or the A-V formulation is just a choice, that most of the time depends on the data we have. </div>
<div class="">Refining the mesh, you should observe convergence of the results.</div>
<div class=""><br class="">
</div>
<div class="">As you’ve done, the coupling between the EM and the thermal problem is done via the Joule losses. </div>
<div class="">These losses are calculated from the frequency domain solution (steady state) of the AV-formulation, and therefore they are an average value.</div>
<div class="">The thermal problem is then solve in the time domain. This is possible thanks to the difference in time constants of both problems. </div>
<div class=""><br class="">
</div>
<div class="">The coupling term can be written as:</div>
<div class=""><br class="">
</div>
<div class="">Galerkin { [ -0.5*sigma[] *<a>[ SquNorm[Dt[{a}]+{d v}] ], {t} ];</div>
<div class=""> In DomainC; Integration II; Jacobian Vol; }</div>
<div class=""><br class="">
</div>
<div class="">where <a> indicates that the operation between square brackets is to be done with complex numbers even if the thermal formulation is real.</div>
<div class=""><br class="">
</div>
<div class="">This term is exactly the same as yours (with a factor 0.5, that I think is missing, to check!) if you also indicate there that the quantities are complex, i.e.</div>
<div class="">
<blockquote type="cite" class="">Galerkin { [ -1./sigma[]*( <a>[ <br class="">
<span class="Apple-tab-span" style="white-space: pre;"></span> Re[-(Dt[{a}]+{d v})]*<br class="">
<span class="Apple-tab-span" style="white-space: pre;"></span> Re[-(Dt[{a}]+{d v})]+<br class="">
<span class="Apple-tab-span" style="white-space: pre;"></span> Im[-(Dt[{a}]+{d v})]*<br class="">
<span class="Apple-tab-span" style="white-space: pre;"></span> Im[-(Dt[{a}]+{d v})]] ), {t} ];<br class="">
<span class="Apple-tab-span" style="white-space: pre;"></span>In DomainC; Integration II; Jacobian Vol; }</blockquote>
</div>
<div class=""><br class="">
</div>
<div class="">If you do not indicate that the quantities are complex, the imaginary part is disregarded in the time domain thermal formulation. </div>
<div class=""><br class="">
</div>
<div class="">Regards, </div>
<div class="">Ruth</div>
<div class=""><br class="">
</div>
<div class="">PS: I am correcting the formulation in the benchmarks. </div>
<div class=""><br class="">
</div>
<br class="">
<div class="">
<div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
—<br class="">
Prof. Ruth V. Sabariego<br class="">
KU Leuven <br class="">
Dept. Electrical Engineering ESAT/Electa, EnergyVille<br class="">
<a href="http://www.esat.kuleuven.be/electa" class="">http://www.esat.kuleuven.be/electa</a></div>
<div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<a href="http://www.energyville.be/" class="">http://www.energyville.be</a></div>
<div style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<br class="">
</div>
<div style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Free software: <a href="http://gmsh.info/" class="">http://gmsh.info</a> | <a href="http://getdp.info/" class="">http://getdp.info</a> | <a href="http://onelab.info/" class="">http://onelab.info</a><br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
</div>
</div>
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">On 17 Oct 2017, at 11:04, ABE Hiroshi <<a href="mailto:habe36@gmail.com" class="">habe36@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Dear All,<br class="">
<br class="">
I am working on a coupling analysis of magnetodynamics and thermal dynamics. Referring to “indheat” sample, I build a model.<br class="">
It uses A-V formulation regarding Magnetodynamics, and I would like to couple the electric current in the thermal formulation.<br class="">
<br class="">
They are:<br class="">
<br class="">
Formulation {<br class="">
<br class="">
{ Name MagStaDyn_av_js0_3D ; Type FemEquation ;<br class="">
Quantity {<br class="">
{ Name a ; Type Local ; NameOfSpace HSpace ; }<br class="">
{ Name v ; Type Local ; NameOfSpace USpace ; }<br class="">
}<br class="">
<br class="">
Equation {<br class="">
Galerkin { [ nu[] * Dof{d a} , {d a} ] ;<br class="">
In Domain ; Jacobian Vol ; Integration II ; }<br class="">
Galerkin { DtDof[ sigma[] * Dof{a} , {a} ] ;<br class="">
In DomainC ; Jacobian Vol ; Integration II ; }<br class="">
Galerkin { [ sigma[] * Dof{d v} , {a} ] ;<br class="">
In DomainC ; Jacobian Vol ; Integration II ; }<br class="">
<br class="">
Galerkin { [ -js0[], {a} ] ;<br class="">
In DomainS ; Jacobian Vol ; Integration II ; }<br class="">
<br class="">
<br class="">
Galerkin{ DtDof[ sigma[] * Dof{a}, {d v} ] ;<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In DomainC ; Jacobian Vol ; Integration II ; }<br class="">
Galerkin{ [ sigma[] * Dof{d v} , {d v} ] ;<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In DomainC ; Jacobian Vol ; Integration II ; }<br class="">
<br class="">
}<br class="">
}<br class="">
<br class="">
{ Name Thermal; Type FemEquation;<br class="">
Quantity {<br class="">
{ Name t; Type Local; NameOfSpace TSpace; }<br class="">
{ Name a; Type Local; NameOfSpace HSpace; }<br class="">
{ Name v; Type Local; NameOfSpace USpace; }<br class="">
}<br class="">
Equation {<br class="">
Galerkin { [ K[] * Dof{d t}, {d t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In DomainC; Integration II; Jacobian Vol; }<br class="">
Galerkin { DtDof [ rho[]*Cp[] * Dof{t}, {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In DomainC; Integration II; Jacobian Vol; }<br class="">
Galerkin { [ -1./sigma[]*( <br class="">
<span class="Apple-tab-span" style="white-space:pre"></span> Re[-(Dt[{a}]+{d v})]*<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span> Re[-(Dt[{a}]+{d v})]+<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span> Im[-(Dt[{a}]+{d v})]*<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span> Im[-(Dt[{a}]+{d v})]), {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In DomainC; Integration II; Jacobian Vol; }<br class="">
<br class="">
Galerkin { [ Ht[]*Dof{t}, {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In Skin_ECore; Jacobian Sur ; Integration II ; }<br class="">
Galerkin { [-Ht[]*Tamb[], {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In Skin_ECore; Jacobian Sur ; Integration II ; }<br class="">
Galerkin { [ sigma_sb*Ep[]*(Dof{t})^4, {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In Skin_ECore; Jacobian Sur ; Integration II ; }<br class="">
Galerkin { [ -sigma_sb*Ep[]*(Tamb[])^4, {t} ];<br class="">
<span class="Apple-tab-span" style="white-space:pre"></span>In Skin_ECore; Jacobian Sur ; Integration II ; }<br class="">
<br class="">
}<br class="">
}<br class="">
}<br class="">
<br class="">
The MagStaDyn_av_js0_3D formulation gives a resonable solution but Thermal gives weird results.<br class="">
<br class="">
I know the “indheat” example uses T-O formulation for coupling analysis. Are there any reasons to take T-O formulation instead of A-V formulation?<br class="">
Any ways for A-V formulation?<br class="">
<br class="">
Thank you so much.<br class="">
<br class="">
Best,<br class="">
<br class="">
ABE Hiroshi<br class="">
from Tokorozawa, JAPAN<br class="">
<br class="">
<br class="">
_______________________________________________<br class="">
getdp mailing list<br class="">
<a href="mailto:getdp@onelab.info" class="">getdp@onelab.info</a><br class="">
<a href="http://onelab.info/mailman/listinfo/getdp">http://onelab.info/mailman/listinfo/getdp</a><br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; font-size: 12px; border-spacing: 0px;"><div>ABE Hiroshi</div><div> from Tokorozawa, JAPAN</div></span>
</div>
<br></body></html>