TR: [Getdp] Re: (no subject)

[Patrick Lefevre]

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>