<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Dear Anton, 
<div class=""><br class="">
<div class="">I’ve gone through your files. The Maxwell Stress Tensor you are using is not correct. You should use exactly the same expression as the one in magnets.pro, </div>
<div class=""><a href="https://gitlab.onelab.info/doc/tutorials/wikis/Magnetic-forces" class="">https://gitlab.onelab.info/doc/tutorials/wikis/Magnetic-forces</a></div>
<div class="">that is</div>
<div class="">
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span id="LC147" class="line" lang="cpp">// Maxwell stress tensor (valid for both formulations and linear materials</span>
<span id="LC148" class="line" lang="cpp">  <span class="n">TM</span><span class="p">[]</span> <span class="o">=</span> <span class="p">(</span> <span class="n">SquDyadicProduct</span><span class="p">[</span><span class="err">$</span><span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="n">SquNorm</span><span class="p">[</span><span class="err">$</span><span class="mi">1</span><span class="p">]</span> <span class="o">*</span> <span class="n">TensorDiag</span><span class="p">[</span><span class="mf">0.5</span><span class="p">,</span> <span class="mf">0.5</span><span class="p">,</span> <span class="mf">0.5</span><span class="p">]</span> <span class="p">)</span> <span class="o">/</span> <span class="n">mu</span><span class="p">[]</span> <span class="p">;</span></span>
</pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">The 0.5 factor comes from the theory not from the geometry. </pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">See e.g.</pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span style="white-space: pre-wrap;" class="">@ARTICLE{2004JCoAM.168..235H,</span></pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><pre style="word-wrap: break-word; white-space: pre-wrap;" class="">   author = {{Henrotte}, F. and {Hameyer}, K.},
    title = "{Computation of electromagnetic force densities: Maxwell stress tensor vs. virtual work principle}",
  journal = {Journal of Computational and Applied Mathematics},
 keywords = {Electromechanical coupling, Electromagnetic force density, Finite element, Maxwell stress tensor, Virtual work principle},
     year = 2004,
    month = jul,
   volume = 168,
    pages = {235-243},
   adsurl = {<a href="https://ui.adsabs.harvard.edu/abs/2004JCoAM.168..235H" class="">https://ui.adsabs.harvard.edu/abs/2004JCoAM.168..235H</a>},
  adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}</pre><div class="">or any other reference dealing with the MST computation.</div></pre>
<pre class="highlight code" lang="cpp"><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">Besides that, the idea of using a shell for bounding the computational domain has to do with using a transformation in the shell and that you can only apply if the two spheres limiting the domain are concentric. If the center is not (0,0,0), you will have to add the center for applying the transformation. </span></font></pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span class="line" lang="cpp"><span class="p">HTH, </span></span></pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">Ruth</pre>
<pre class="highlight code" lang="cpp" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span class="line" lang="cpp"><span class="p"><br class=""></span></span></pre>
<div class=""><br class="">
</div>
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
—<br class="">
Prof. Ruth V. Sabariego<br class="">
KU Leuven, 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="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<a href="http://www.energyville.be" class="">http://www.energyville.be</a></div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
</div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
</div>
</div>
</div>
</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 26 May 2019, at 22:06, Panteleev, Anton <<a href="mailto:anton.panteleev@haw-hamburg.de" class="">anton.panteleev@haw-hamburg.de</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; direction: ltr; font-family: Tahoma; font-size: 10pt;" class="">
<div class="">dear Ruth,</div>
<div class=""><br class="">
</div>
<div class="">
<pre class="tw-ta tw-data-text tw-text-small" data-placeholder="Übersetzung" id="tw-target-text" data-fulltext="" dir="ltr" style="text-align: left; height: 168px;"><span tabindex="0" lang="en" class="">in my previous letter, namely in the pro. file I made a thousand mistakes and therefore it did not work for me, <br class="">now I have finished this file and it gives me all the forces in all three planes Fx,Fy,Fz). <br class="">I hope that these forces are calculated correctly.<br class=""><br class="">best regards<br class="">Anton Panteleev<br class=""></span></pre>
</div>
</div>
<span id="cid:1E358444-5D17-4407-B98A-D69042A97ADA@esat.kuleuven.be"><P5.geo></span><span id="cid:789A8B84-4D11-46E2-A946-9FE407C9C79D@esat.kuleuven.be"><P5.pro></span><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">getdp
 mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="mailto:getdp@onelab.info" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">getdp@onelab.info</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="http://onelab.info/mailman/listinfo/getdp" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://onelab.info/mailman/listinfo/getdp</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</body>
</html>