<div dir="ltr"><div><div><div>Dear Frederic, <br><br></div>my final goal is to calculate the current (I = int_coil(js - sigma[]*Dt[{Aphi})) at every time instance and update the current density in order to achieve a certain current value. <br><br>Regarding your suggestion, is it possible to define such constraint which depends on the solution (Aphi)?<br><br></div>Best Regards<br><br></div>Peter <br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 13, 2015 at 4:40 PM, Frederic Trillaud Pighi <span dir="ltr"><<a href="mailto:ftrillaudp@pumas.iingen.unam.mx" target="_blank">ftrillaudp@pumas.iingen.unam.mx</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Peter,<div><br></div><div>You may adjust your TimeFunction to your need in the constraint. It is interpolated during the resolution.</div><div><br></div><div>Best,</div><div><br></div><div>Frederic </div><div><br><div><blockquote type="cite"><div><div class="h5"><div>On Aug 13, 2015, at 9:34 AM, Peter Kis <<a href="mailto:md2z34@gmail.com" target="_blank">md2z34@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr"><div><div><div><div><div><div>Dear List Members,<br><br></div>I'm just wondering if it's possible to scale the current density from the time iteration loop. <br><br>Here I introduced a parameter <b>current_factor</b>, <br><br><span style="font-family:monospace,monospace">Function{ <br>    ......<br>    <b>current_factor</b> = 1.;
<br>}
</span><br><br></div>which appears in the formulation as a scaling factor for the current density:<br><br><span style="font-family:monospace,monospace">Formulation{ // OK
<br>  {Name VectorPotentialFormulation; Type FemEquation;
<br>   Quantity{
<br>     {Name Aphi; Type Local; NameOfSpace Afs;}
<br>     {Name js; Type Local; NameOfSpace jfs;}
<br>   }
<br>   Equation{
<br>    Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];
<br>      In Omega; Jacobian JVol; Integration I1;}
<br> <br>    Galerkin{ DtDof[sigma[]*Dof{Aphi},{Aphi}];
<br>      In Omega; Jacobian JVol; Integration I1;}
<br> <br>    Galerkin{ [-<b>current_factor</b>*Dof{js},{Aphi}];
<br>      In Omega_coil; Jacobian JVol; Integration I1;}
<br>    }
<br>  }
<br>}
</span><br><br></div>and I'd like to update this parameter during the time iteration:<br><br><span style="font-family:monospace,monospace">Resolution{
<br>  {Name Solution;
<br>    System{
<br>      {Name Syst; NameOfFormulation VectorPotentialFormulation;}
<br>    }
<br>    Operation{
<br>        InitSolution[Syst];SaveSolution[Syst];
<br>        TimeLoopTheta[t0,tmax,dt,theta] {<br>            <b>current_factor</b> = 1/1.7; 
// This has no effect!!<br>            Generate[Syst]; Solve[Syst]; SaveSolution[Syst];
<br>        }
<br>    }
<br>  }
<br>}</span><br><br></div>but this has no impact on the solution. Apparently the value of current_vector is still considered one as it's defined in Functions. <br><br>Is there any way to update this parameter from the time loop? <br><br></div>Thanks for your help in advance.<br><br><br></div>Peter<br><div><div><br></div></div></div></div></div>
_______________________________________________<br>getdp mailing list<br><a href="mailto:getdp@geuz.org" target="_blank">getdp@geuz.org</a><br><a href="http://www.geuz.org/mailman/listinfo/getdp" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br></div></blockquote></div><br></div></div></blockquote></div><br></div>