[Getdp] Magnetostatic 3D

Alain J-A.Martin at sympatico.ca
Tue Apr 3 04:59:16 CEST 2007


Kubicek Bernhard wrote:

>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.
>
>  
>
Thank you for this clarification of the vector-potential problem which 
was puzzling. If I understand well, the best way to visualize the 
magnetic field lines is by using the streamlines plugin?

>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
>  
>
Thank you also for the very interesting paper.
Alain