[Getdp] Resolution with a parameter

Peter Kis md2z34 at gmail.com
Fri Aug 14 12:12:17 CEST 2015


Dear Christophe,

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.

Do you think it is possible to apply the same constraint to my problem too?

Best Regards

Peter


On Thu, Aug 13, 2015 at 6:26 PM, Christophe Geuzaine <cgeuzaine at ulg.ac.be>
wrote:

>
> You should simply impose the total current (including the induced
> currents) as a global quantity. See
>
>   http://www.geuz.org/getdp/doc/texinfo/getdp.html#Magnetodynamic-problem
>
> for a complete example (cf. the "Current_2D" constraint).
>
>
> > On 13 Aug 2015, at 16:47, Peter Kis <md2z34 at gmail.com> wrote:
> >
> > Dear Frederic,
> >
> > 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.
> >
> > Regarding your suggestion, is it possible to define such constraint
> which depends on the solution (Aphi)?
> >
> > Best Regards
> >
> > Peter
> >
> > On Thu, Aug 13, 2015 at 4:40 PM, Frederic Trillaud Pighi <
> ftrillaudp at pumas.iingen.unam.mx> wrote:
> > Dear Peter,
> >
> > You may adjust your TimeFunction to your need in the constraint. It is
> interpolated during the resolution.
> >
> > Best,
> >
> > Frederic
> >
> >> On Aug 13, 2015, at 9:34 AM, Peter Kis <md2z34 at gmail.com> wrote:
> >>
> >> Dear List Members,
> >>
> >> I'm just wondering if it's possible to scale the current density from
> the time iteration loop.
> >>
> >> Here I introduced a parameter current_factor,
> >>
> >> Function{
> >>     ......
> >>     current_factor = 1.;
> >> }
> >>
> >> which appears in the formulation as a scaling factor for the current
> density:
> >>
> >> Formulation{ // OK
> >>   {Name VectorPotentialFormulation; Type FemEquation;
> >>    Quantity{
> >>      {Name Aphi; Type Local; NameOfSpace Afs;}
> >>      {Name js; Type Local; NameOfSpace jfs;}
> >>    }
> >>    Equation{
> >>     Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];
> >>       In Omega; Jacobian JVol; Integration I1;}
> >>
> >>     Galerkin{ DtDof[sigma[]*Dof{Aphi},{Aphi}];
> >>       In Omega; Jacobian JVol; Integration I1;}
> >>
> >>     Galerkin{ [-current_factor*Dof{js},{Aphi}];
> >>       In Omega_coil; Jacobian JVol; Integration I1;}
> >>     }
> >>   }
> >> }
> >>
> >> and I'd like to update this parameter during the time iteration:
> >>
> >> Resolution{
> >>   {Name Solution;
> >>     System{
> >>       {Name Syst; NameOfFormulation VectorPotentialFormulation;}
> >>     }
> >>     Operation{
> >>         InitSolution[Syst];SaveSolution[Syst];
> >>         TimeLoopTheta[t0,tmax,dt,theta] {
> >>             current_factor = 1/1.7;  // This has no effect!!
> >>             Generate[Syst]; Solve[Syst]; SaveSolution[Syst];
> >>         }
> >>     }
> >>   }
> >> }
> >>
> >> but this has no impact on the solution. Apparently the value of
> current_vector is still considered one as it's defined in Functions.
> >>
> >> Is there any way to update this parameter from the time loop?
> >>
> >> Thanks for your help in advance.
> >>
> >>
> >> Peter
> >>
> >> _______________________________________________
> >> getdp mailing list
> >> getdp at geuz.org
> >> http://www.geuz.org/mailman/listinfo/getdp
> >
> >
> > _______________________________________________
> > getdp mailing list
> > getdp at geuz.org
> > http://www.geuz.org/mailman/listinfo/getdp
>
> --
> Prof. Christophe Geuzaine
> University of Liege, Electrical Engineering and Computer Science
> http://www.montefiore.ulg.ac.be/~geuzaine
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150814/8b6db15a/attachment.html>