toujours des problemes...

trophime christophe trophime at labs.polycnrs-gre.fr
Mon Mar 5 11:52:56 CET 2001


> 
> Print[ b, OnPlane {{1.e-3,0,0}{250.e-3, 0, 0}{1.e-3,250.e-3,0}} {50,50}
> , File "CoreSta_b2.pos"] ;
> 
> Pour ameliorer les resultats en Axi, je te conseille egalement
> d'utiliser, non pas VolAxi et VolAxiSphShell, mais VolAxiSqu et
> VolAxiSquSphShell (cf. F. Henrotte et al., Finite element modelling with
> transformation techniques. IEEE Transactions on Magnetics,
> 35(3):1434-1437, May 1999).
> 

Effectivement ca me permet de faire le calcul. Mais je  continue a pas
comprendre
ce qui se passe. En fait dans ce cas simple on peut expliciter les
solutions
en termes de potentiel et de champ magnetique a l'aide d'integrales
elliptiques
completes... Ce que j'essaie de faire c'est retrouver la solution
analytique bien evidemment avec
une precision qui depend du maillage. Pour cela je fais tracer a_{theta}
et b_z
suivant l'axe Or:

    Print[ b, OnLine {{1.e-06,0,0}{250.e-3, 0, 0}} {50}, File
"CoreSta_b3.cut" , Format Gnuplot ] ;
    Print[ a_theta, OnLine {{0,0,0}{250.e-3, 0, 0}} {50}, File
"CoreSta_a3.cut" , Format Gnuplot ] ;

avec comme definition de a_theta et b :

        { Name a_theta ; Value { Local { [ CompZ[{a}] ]   ; In
Domain_Mag ; } } }
        { Name b  ; Value { Local { [ {d a} ] ; In Domain_Mag ; Jacobian
VolAxi ; } } }

Dans le cas d'une distribution uniforme je ne retrouve pas du
tout a_{theta} mais par contre b_z est "juste". 
Dans le cas d'une distribution de bitter (i.e. en 1/r) rien ne va plus.
J'arrive pas a savoir s'il s'agit d'un probleme de definition chez moi
ou d'un Bug?

Pour le maillage que je t'ai communique dans mon precedent mail
la solution analytique (au signe pres) est donnee pour l'axe Or
dans le fichier que je te mets en attachment. La valeur du courant a ete
fixee
a 1.e5. La densite de courant est "recalculee" a chaque fois pour le cas
"Bitter"
et "Uniforme".

> >
> > Par ailleurs j'ai une remarque. Quand vous definissez les fonctions de
> > Green
> > je pense qu'il y a une precision a apporter. La fonction de Green 2D
> > n'est uniquement
> > valable qu'en 2D plan. Dans le cas de l'axi elle est differente. En
> > particulier on a :
> >
> > r = \sqrt{r^2-2 a r \cos\theta + (z-z_a)^2},
> >
> > et la fonction de Green est alors de type 1/r a un facteur 2 \pi pres
> > (bien sur il
> > faut encore intégrer sur \theta).
> >
> > Avez-vous essayer des couplages avec des methodes integrales de
> > frontieres pour le
> > cas axisymetrique?
> >
> 
> Oui, c'est exact, mais je n'ai pas encore introduit le noyau pour le cas
> axisymetrique. Si mes souvenirs sont bons, il est assez complexe
> integrer analytiquement... Si tu possedes une experience avec ce type de
> noyaux (et si tu es interesse, bien sur !), on pourrait songer a une
> collaboration.

Je ne suis pas vraiment un specialiste de la question mais j'ai
travaille
sur le calcul analytique des inductances pour des geometries axi simples
qui se trouvent "ressembler" aux integrales a introduire dans le modele
numerique.
Aussi je pense que ca pourrait etre interressant de songer a
quelquechose... 
mais je souhaiterais dans un premier temps
arriver a faire tourner getdp pour nos types de problemes pour ensuite
convaincre
mes chefs plus facilement.

Encore merci.
Ch. Trophime
Tel 33 (0) 4 76 88 90 02