<div dir="ltr"><div><div><div>Hi getDP users,<br><br></div>I am about to implement a single turn coil with nonlinear sigma using AV formalism for the coil region and A formalism for the surrounding air domain.<br><br></div>I use a global term to prescribe the total current in the eddy current domain, which works like a charm for the linear case, while the time iteration is done with the DtDof{} operator. When I replace that term with a manual implementation of the backward Euler and put that into the equation, it becomes divergent even for constant sigma.<br><br></div>I followed Michael's instructions to implement Picard iteration:   <br><div><div><div><div><br><a href="http://www.geuz.org/pipermail/getdp/2012/001533.html">http://www.geuz.org/pipermail/getdp/2012/001533.html</a><br><br></div><div>But using Picard iteration for the nonlinear problem is my final goal. Now I just want to make the linear case work with manual backward Euler.<br></div><div><br></div><div>The symptoms are the following despite using constant sigma:<br><ol><li>the vector potential I get is slightly different from the DtDof{} case, but not too much<br></li><li>the gradV becomes much larger by one order of magnitude<br></li><li>the total current, which is imposed by global terms reaches 1700 A and it should be 100 A</li></ol><p>These are formulations:</p><p><font size="1"><span style="font-family:monospace,monospace">Formulation{ <br>  {Name AV_lin; Type FemEquation;<br>   Quantity{<br>     {Name Aphi; Type Local; NameOfSpace Afs;}<br>     {Name gradV; Type Local; NameOfSpace gradVfs;}<br>     {Name I; Type Global; NameOfSpace gradVfs [I]; }<br>     {Name U; Type Global; NameOfSpace gradVfs [U]; }<br>   }<br>   Equation{<br>    Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];<br>      In Omega; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ DtDof[sigma_lin[]*Dof{Aphi},{Aphi}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [sigma_lin[]*Dof{gradV},{Aphi}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ DtDof[sigma_lin[]*Dof{Aphi},{gradV}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [sigma_lin[]*Dof{gradV},{gradV}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    GlobalTerm { [ Dof{I} , {U} ]; In Omega_coil; }<br>   }<br>  }<br><br>  {Name AV_Picard; Type FemEquation;<br>   Quantity{<br>     {Name Aphi; Type Local; NameOfSpace Afs;}<br>     {Name gradV; Type Local; NameOfSpace gradVfs;}<br>     {Name Aphipast; Type Local; NameOfSpace Apastfs;}<br>     {Name gradVpast; Type Local; NameOfSpace gradVpastfs;}<br>     {Name I; Type Global; NameOfSpace gradVfs [I]; }<br>     {Name U; Type Global; NameOfSpace gradVfs [U]; }<br>   }<br>   Equation{<br>    Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];<br>      In Omega; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [(sigma_lin[]*Dof{Aphi} - sigma_lin[]*{Aphipast})/$DTime,{Aphi}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [sigma_lin[]*Dof{gradV},{Aphi}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [(sigma_lin[]*Dof{Aphi} - sigma_lin[]*{Aphipast})/$DTime,{gradV}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [sigma_lin[]*Dof{gradV},{gradV}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br><br>    GlobalTerm { [ Dof{I} , {U} ]; In Omega_coil; }<br>   }<br>  }}</span></font></p><p><br></p><p>Could you take a look at my code? Maybe you can spot something.</p><p>Thank for your help in advance.</p><p>Best Regards</p><p>Peter<br></p></div></div></div></div></div>