nice getdp results and some questions

[Samuel Kvasnica]

This is a multi-part message in MIME format.
--------------FBA1D314A2CC8D9EDD2009FA
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit

Hi Christophe,

this week a had again some time to play with getdp. All the problems with cyl. axis seem to be
really away after specifying that "Dimension 2" option. Now I'm going to simulate a magnetron
we've constructed few years ago. I made some measurements using an F.W.Bell gaussmeter
to compare the getdp output with reality (well, also to find out Br of used magnets). I
measured
Br and Bz at ~70mm lines in r direction at distance z=2,5,10,20mm from magnetron. After I used
nonlinear BH curve, I started to get quite good results, look at complete.eps. I used just
generic BH curve for steel. I could fit that curve by H=a*(1-exp(-b*B))+c*B in 0-0.4T region.
Unfortunatelly I could not find out how to put it into getdp (how to get unity vector for B,
any hint ? ) so I finally aproximated it by 5th degree polynome and extracted B to get nu[B].
Btw, could you give me an advice for some good source of material property data - e.g.
magnetization curves for standard materials ?
Well, regarding the second attached picture singular.eps, there still seems to be some problem
with interpolation. I generated data for B at line going exactly through iron-air boundary
(magenta color, at z=-6mm  in umag1.geo) and there seems to be some funny noise in Br, Bz is
ok. It seems to correspond to mesh points someway.
The same with line at 3.5mm (cyan), which touches the middle pole, but there's just slight
noise visible.
Do you think it's a bug ?

best regards,

Sam

--
--------------------------------------------------------------
Mag.Samuel Kvasnica
IEMW TU-WIEN
Gusshausstrase 27-29 / 366   email: kvasnica at iaee.tuwien.ac.at
1040 Wien                    tel:   +43-1-58801-36639
Austria                      fax:   +43-1-58801-36695
--------------------------------------------------------------


Content-Type: text/plain; charset=us-ascii;
 name="umag1.geo"
Content-Transfer-Encoding: 7bit
Content-Disposition: inline;
 filename="umag1.geo"

yInt = 350e-3;
yExt = 450e-3;
pInt=7e-3;
pExt = 10e-3;
pDef = 1e-3;

Point(0) = {0,0,0,pDef};

ref2 = 35e-3;
lr = 56e-3;
lh = 2.5e-3;
pr = 10.5e-3;
ph = 29e-3;
Point(7) = {0,ref2,0,pDef};
Point(8) = {0,ref2-ph-lh,0,pDef};
Point(9) = {pr,ref2-ph-lh,0,pDef};
Point(10) = {pr,ref2-lh,0,pDef};
Point(11) = {lr,ref2-lh,0,pDef};
Point(12) = {lr-2e-3,ref2,0,pDef};

ref3 =  ref2 - lh;
mr1 = 36e-3;
mr2 = 54e-3;
mh = 24.5e-3;
Point(13) = {mr1,ref3,0,pDef};
Point(14) = {mr2,ref3,0,pDef};
Point(15) = {mr2,ref3-mh,0,pDef};
Point(16) = {mr1,ref3-mh,0,pDef};

ref4 = ref3-mh;
sr1 = 36e-3;
sr2 = 56e-3;
sh = 6e-3;
/*Point(17) = {sr1,ref4,0,pDef};*/
Point(18) = {sr1,ref4-sh+4e-3,0,pDef};
Point(19) = {sr2-11e-3,ref4-sh+4e-3,0,pDef};
Point(20) = {sr2-11e-3,ref4-sh,0,pDef};
Point(21) = {sr2,ref4-sh,0,pDef};
Point(22) = {sr2,ref4,0,pDef};

ref1 = ref4-sh;
clamp_r1 = 39e-3;
clamp_r2 = 56e-3;
clamp_h = 8e-3;
Point(1) = {clamp_r2-12e-3,ref1-clamp_h,0,pDef};
Point(2) = {clamp_r2,ref1-clamp_h,0,pDef};
/*Point(3) = {clamp_r2,ref1,0,pDef};*/
/*Point(4) = {clamp_r2-11e-3,ref1,0,pDef};*/
Point(5) = {clamp_r2-11e-3,ref1-2e-3,0,pDef};
Point(6) = {clamp_r1,ref1-2e-3,0,pDef};

Line(1) = {7,8};
Line(2) = {8,9};
Line(3) = {9,10};
Line(4) = {10,13};
Line(5) = {13,14};
Line(6) = {14,11};
Line(7) = {11,12};
Line(8) = {12,7};
Line(9) = {13,16};
Line(10) = {16,15};
Line(11) = {15,14};
Line(12) = {15,22};
Line(13) = {22,21};
Line(14) = {21,20};
Line(15) = {20,19};
Line(16) = {19,18};
Line(17) = {18,16};
Line(18) = {21,2};
Line(19) = {2,1};
Line(20) = {1,6};
Line(21) = {6,5};
Line(22) = {5,20};
Line Loop(23) = {6,7,8,1,2,3,4,5};
Plane Surface(24) = {23};
Line Loop(25) = {-11,-10,-9,5};
Plane Surface(26) = {25};
Line Loop(27) = {12,13,14,15,16,17,10};
Plane Surface(28) = {27};
Line Loop(29) = {-22,-21,-20,-19,-18,14};
Plane Surface(30) = {29};

Point(100)={0,-yInt,0,pInt};
Point(101)={0,-yExt,0,pExt};
Point(102)={0,yInt,0,pInt};
Point(103)={0,yExt,0,pExt};

Circle(200)={101,0,103}; /*infinity*/
Circle(201)={100,0,102}; /*internal border*/


Physical Surface(2000) = {24}; /* upper loop */
Physical Surface(2001) = {26}; /* magnet */
Physical Surface(2002) = {28}; /* slot */
Physical Surface(2003) = {30}; /* clamp */


Line(2004) = {103,102};
Line(2005) = {102,7};
Line(2006) = {8,100};
Line(2007) = {100,101};
Line Loop(2008) = {-201,2007,200,2004};
Plane Surface(2009) = {2008};
Line Loop(2010) = {9,-17,-16,-15,-22,-21,-20,-19,-18,-13,-12,11,6,7,8,-2005,-201,-2006,2,3,4};
Plane Surface(2011) = {2010};

Physical Surface(2100) = {2011}; /* air */
Physical Surface(2101) = {2009}; /* infinity air */

Physical Line(1100) = {2004,2005,1,2006,2007}; /* SurfaceGh0 */
Physical Line(1102) = {200}; /* SurfaceGInf */


--------------FBA1D314A2CC8D9EDD2009FA--