[Getdp] Problem with description of magnetisation

Peter Fouquet pf231 at cam.ac.uk
Tue Oct 2 18:42:11 CEST 2001


Dear Christophe,

thank you very much for the help.
Following this, I immediately have three further questions:

1) How can I define hc[] as a vector quantity? I found only scalar
   constants in the getdp documentation. Did I understand anything
   wrong?

2) Is a loop of surfaces already a 'forbidden' geometry from a 
   topological point of view?

3) Is it possible to print a function or constant in the 
   postprocessing (for checking purposes)?
   When I tried this, the program complained, that it was no 
   discrete quantity...

Hope to hear from you soon. :-)

Peter

Patrick Dular wrote:
> 
> Hi Peter,
> 
> For this kind of problem (static fields, no source currents), I think
> the most efficient solution strategy is to use a scalar magnetic
> potential. You can then compute the magnetization as a post-processing
> quantity. In the linear case, this formulation is indeed the one given
> in the demo file of the getdp distribution, i.e.
> 
> Formulation {
>   { Name MagSta_phi ; Type FemEquation ;
>     Quantity {
>       { Name phi ; Type Local ; NameOfSpace Hgrad_phi ; }
>     }
>     Equation {
>       Galerkin { [ - mu[] * Dof{d phi} , {d phi} ] ;
>                  In Domain ; Jacobian JVol ; Integration I1 ; }
>       Galerkin { [ - mu[] * hc[] , {d phi} ] ;
>                  In Domain_M ; Jacobian JVol ; Integration I1 ; }
>     }
>   }
> }
> 
> where 'hc[]' is the volume source term (the coercitive magnetic field in
> your magnets). This formulation is only valid if 'Domain' is
> topologically trivial (no loops): for its construction, see
> http://www.geuz.org/getdp/doc/slides/getdp-18.html,
> http://www.geuz.org/getdp/doc/slides/getdp-25.html (and following), with
> 'j=0'. If the materials are weakly nonlinear, just write 'mu[{d phi}]'
> instead of 'mu[]' in the equations, and use an 'IterativeLoop' in the
> resolution. If the materials are strongly nonlinear, you may want to use
> a Newton method, which requires the definition of an additional term in
> the formulation. Tell me if this is the case: I'll show you how to
> construct it.
> 
> Christophe
> 
> pf231 at cam.ac.uk wrote:
> >
> > Dear all,
> >
> > Thank you, first of all, to Christophe and Patrick for the brilliant
> > getdp program.
> >
> > I am presently working on a 3D FEM simulation of the transition region
> > from a permanent magnetic hexapole field to a dipole field. Being no
> > expert in FEM simulations at all I find myself in a problem, where I
> > need some help, which I hope you can provide me with:
> >
> > I want to describe the system in a vector potential approach with the
> > edge finite elements with the tree gauging, that you sketched in the
> > manual. The problem comes in, when I try to describe the magnetisation.
> >
> > Which function space can I use for it? A vector field didn?t work,
> > because I couldn?t define a curl on it (curl (magnetisation)=0), but
> > otherwise I don?t know how to get it into the Galerkin equation (if I
> > use it as a constant) or to define it properly (if I use a 1form space
> > for the magnetisation).
> >
> > If you could give me a direction to solve the problem, I would be
> > extremely grateful to you.
> >
> > Thank you very much in advance, best regards
> >
> > Peter Fouquet
> >
> > _______________________________________________
> > GetDP mailing list
> > GetDP at geuz.org
> > http://www.geuz.org/mailman/listinfo/getdp
> 
> --
> Christophe Geuzaine
> 
> Tel: 32 (0) 4 366 37 10    http://geuz.org
> Fax: 32 (0) 4 366 29 10    mailto:Christophe.Geuzaine at ulg.ac.be
> _______________________________________________
> GetDP mailing list
> GetDP at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp