<div dir="ltr"><div><div><div><div>Dear Christophe,<br><br></div>thank you for recommending me to look at the magnetodynamic problem. As I see A-V formalism is used there and the Current_2D constraint is applied to the electric scalar potential part i FunctionSpace  Hregion_u_Mag_2D. However I don't have scalar potential in my particular problem. I use the magnetic vector potential only. <br><br></div>Do you think it is possible to apply the same constraint to my problem too?<br><br></div>Best Regards<br><br></div>Peter<br><div><div><div> <br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 13, 2015 at 6:26 PM, Christophe Geuzaine <span dir="ltr"><<a href="mailto:cgeuzaine@ulg.ac.be" target="_blank">cgeuzaine@ulg.ac.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
You should simply impose the total current (including the induced currents) as a global quantity. See<br>
<br>
  <a href="http://www.geuz.org/getdp/doc/texinfo/getdp.html#Magnetodynamic-problem" rel="noreferrer" target="_blank">http://www.geuz.org/getdp/doc/texinfo/getdp.html#Magnetodynamic-problem</a><br>
<br>
for a complete example (cf. the "Current_2D" constraint).<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
> On 13 Aug 2015, at 16:47, Peter Kis <<a href="mailto:md2z34@gmail.com">md2z34@gmail.com</a>> wrote:<br>
><br>
> Dear Frederic,<br>
><br>
> 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>
> Best Regards<br>
><br>
> Peter<br>
><br>
> On Thu, Aug 13, 2015 at 4:40 PM, Frederic Trillaud Pighi <<a href="mailto:ftrillaudp@pumas.iingen.unam.mx">ftrillaudp@pumas.iingen.unam.mx</a>> wrote:<br>
> Dear Peter,<br>
><br>
> You may adjust your TimeFunction to your need in the constraint. It is interpolated during the resolution.<br>
><br>
> Best,<br>
><br>
> Frederic<br>
><br>
>> On Aug 13, 2015, at 9:34 AM, Peter Kis <<a href="mailto:md2z34@gmail.com">md2z34@gmail.com</a>> wrote:<br>
>><br>
>> Dear List Members,<br>
>><br>
>> 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 current_factor,<br>
>><br>
>> Function{<br>
>>     ......<br>
>>     current_factor = 1.;<br>
>> }<br>
>><br>
>> which appears in the formulation as a scaling factor for the current density:<br>
>><br>
>> 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{ [-current_factor*Dof{js},{Aphi}];<br>
>>       In Omega_coil; Jacobian JVol; Integration I1;}<br>
>>     }<br>
>>   }<br>
>> }<br>
>><br>
>> and I'd like to update this parameter during the time iteration:<br>
>><br>
>> 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>
>>             current_factor = 1/1.7;  // This has no effect!!<br>
>>             Generate[Syst]; Solve[Syst]; SaveSolution[Syst];<br>
>>         }<br>
>>     }<br>
>>   }<br>
>> }<br>
>><br>
>> 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>
>> Thanks for your help in advance.<br>
>><br>
>><br>
>> Peter<br>
>><br>
>> _______________________________________________<br>
>> getdp mailing list<br>
>> <a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
>> <a href="http://www.geuz.org/mailman/listinfo/getdp" rel="noreferrer" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br>
><br>
><br>
> _______________________________________________<br>
> getdp mailing list<br>
> <a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
> <a href="http://www.geuz.org/mailman/listinfo/getdp" rel="noreferrer" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Prof. Christophe Geuzaine<br>
University of Liege, Electrical Engineering and Computer Science<br>
<a href="http://www.montefiore.ulg.ac.be/~geuzaine" rel="noreferrer" target="_blank">http://www.montefiore.ulg.ac.be/~geuzaine</a><br>
<br>
<br>
<br>
</font></span></blockquote></div><br></div>