[Getdp] Magnetostatic 3D

Alain J-A.Martin at sympatico.ca
Wed Apr 4 00:36:51 CEST 2007


Olivier Castany wrote:

>Private message :
>-----------------
>
>Hello,
>
>you posted a message on the GetDP mailing list a few hours ago.
>
>It seems that you had an interesting private discussing with Bernhard 
>Kubicek. 
>
>I am also interested by a working source code with tree-gauge. If you 
>have a working example, could you consider sending it on the mailing 
>list (or putting it in the wiki) ?
>
>Regards,
>
>O.C.
>
>  
>
>>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
>>>      
>>>
>                             
>
>  
>
Hello Olivier,
    My sample code was send to the getdp list on 03/30/2007 and it is 
called ring_magnet..........
and it is sort of working. All communications was done through the list 
as far as I know.
I hope this help.
Alain