TR: [Getdp] Re: (no subject)
Patrick Lefevre
Patrick.Lefevre at umh.ac.be
Tue Nov 19 08:58:28 CET 2002
Hello Janvrot, hello Johan,
I'm also trying to solve non linear magnetostatic problems. Before asking
questions, i have one remark. For the phi formulation, I think you don't
have to impose a constraint on the symmetry line. In case.pro, you write :
Dirichlet_phi = Region[ {LineA_Code, LineC_Code} ]; // symetry and
boundary
I think you have to try with
Dirichlet_phi = Region[ {LineC_Code} ]; // boundary
Now the questions about solving a non linear problem with the phi
formulation:
In my problem, magnet's are non linear (Alnico magnets). To give an example
of the method i use for solving see attached file "hlk.geo, hlk.pro,
solver.par". I'm not sure this is the good way to solve the problem, but the
results are better when i compare to magnetic field measure.
Thanks for help
-----Message d'origine-----
De : janvrot at infolink.com.br [mailto:janvrot at infolink.com.br]
Envoyé : mardi 19 novembre 2002 01:21
À : Johan Gyselinck
Cc : GetDP
Objet : Re: [Getdp] Re: (no subject)
Johan,
I did not get the right Galerkin. I suppose the PostProcessing is not
correct either.
Attached is a 2D problem with scalar mag potential that should run with
3D, I guess.
Please
gmsh case -2
getdp case -msh case.msh -solve MagSta_phi_nl -pos phi
on it.
Note that there is no current in the problem.
Janvrot
Johan Gyselinck wrote:
>
> janvrot at infolink.com.br wrote:
>
> > Salut Johan
> >
> > That's it! Otherwise, do you know hot to introduce this non linear
> > term
> > in the scalar magnetic potential formulation? I was running a case
> > with
> > it, and I supose this way I could solve my problem.
>
> You just have to replace dhdb[] and h[] by dbdh[] and b[]
> respectively. This was already partly prepared in the file with the
> tables I sent you a while ago.
>
> Greetings,
>
> Johan
>
>
> --
> Dr. ir. Johan Gyselinck
> University of Liege
> Dept. of Electrical Engineering (ELAP), Institut Montefiore
> Campus du Sart Tilman, Parking P32, Bâtiment B28, B-4000 LIEGE, Belgique
> phone: +32-4-3663732 -- fax: +32-4-3662910 --
http://elap.montefiore.ulg.ac.be/
>
>
--
Janvrot IVM
janvrot at infolink.com.br
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hlk.geo
Type: application/octet-stream
Size: 2310 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20021119/f8eafd95/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solver.par
Type: application/octet-stream
Size: 4891 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20021119/f8eafd95/attachment.par>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hlk.pro
Type: application/octet-stream
Size: 5547 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20021119/f8eafd95/attachment.pro>