[Getdp] Magnetostatic 3D

[Kubicek Bernhard]

Dear all,
especially Olivier..
It seems my last mail to the list was rejected, probably because of a 300kb paper I half-legaly attached (and which is now excluded).
And because of Alain's reply, I did not believe it was necessary to repost the content. However, if this mail comes through, you will recieve all the files that were transferred.

The mag. vec. potential it giving me literary nightmares at the moment, as I am working on it 9+ hours a day for the last weeks, and still cannot get it to work on my real geometries, although it simply _has_ to for us to being able to present at Compumag in Achen. On hexahedral meshes, I have not been able to create anything useful yet. But I have constructed a test case which I will present in a mail that is being composed within the hour...

very nice greetings,
 Bernhard

PS: I just found that I accidentially had also attached the 1MB preprocessing file. Stupid me.

-----Ursprüngliche Nachricht-----
Von: Kubicek Bernhard 
Gesendet: Montag, 2. April 2007 18:23
An: 'Alain'; getdp at geuz.org
Betreff: Re: [Getdp] Magnetostatic 3D


Hello Alain,

after fooling around a bit with your files (and finding your dual-constraint in the formulation quite useful for my stuff aswell !) I managed to calculate a (probably) useful magnet field. However, you cannot expect to get a physically meaningful vector-potential due to the limitations of the tree-gauge, only the derived magnetic field is meaningful.

See e.g. the following attached paper:
Ticar, I.; Biro, O.; Preis, K.: Vector potential expanded by edge basis functions associated with loops on finite element facets. IEEE transactions on magnetics 38 (2002) , S. 437-440

What I did, basically, to make your files work somehow is to make the mesh finer, change nbfill in solverpar to 20 (I tend to use something between the max/average bandwith of the inital matrix), and uncomment the output to the mag. field.

Personally, I would expand the calculation region further along the axis of the core.

Now, if I only could make my own calculations work (electric current based magnetic vec. potential in 3d on a _hexahedral_ grid)....

Very nice greetings,
 Bernhard





-----Ursprüngliche Nachricht-----
Von: getdp-bounces at geuz.org [mailto:getdp-bounces at geuz.org] Im Auftrag von Alain
Gesendet: Mittwoch, 28. März 2007 22:55
An: getdp at geuz.org
Betreff: [Getdp] Magnetostatic 3D


Hi all,
I have tried everything I could find but cannot solve to my satisfaction 
any 3D magnetostatic cases either from the list archive or the wiki 
examples. The only case where I get a good 3D solution is the 
magnetodynamic case CoilAboveCrack but only for high frequencies. Is their something specific I have to follow to get good results?
e.g. Basis Function   using EdgesOf or NodesOf  or  ?
             Form0, Form1, Form1P  or ?
       Constraints  :  Gauge or scalar ?
Is the formulation good?
Did anyone succeeded in getting a good solution for example the simple 
cube magnet describe somewhere in this list?
Thanks,
Alain
_______________________________________________
getdp mailing list
getdp at geuz.org
http://www.geuz.org/mailman/listinfo/getdp
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Integration_Lib.pro
Type: application/octet-stream
Size: 1237 bytes
Desc: Integration_Lib.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070404/4f712ad1/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MagSta_a_3Dp.pro
Type: application/octet-stream
Size: 3189 bytes
Desc: MagSta_a_3Dp.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070404/4f712ad1/attachment-0001.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ring_magnet.geo
Type: application/octet-stream
Size: 4229 bytes
Desc: ring_magnet.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070404/4f712ad1/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ring_magnet.pro
Type: application/octet-stream
Size: 2382 bytes
Desc: ring_magnet.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070404/4f712ad1/attachment-0002.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solver.par
Type: application/octet-stream
Size: 5059 bytes
Desc: solver.par
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070404/4f712ad1/attachment.par>