[Getdp] calcul de courant

[Patrick Dular]

Cher Ali,

Voici une archive avec un ensemble d'exemples de problèmes électrocinétiques et en électrostatiques à partir desquel tu pourras t'inspirer pour ton problème.

Le traitement des grandeurs globales (courant, tension, charge) y est décrit. Il est fait usage de la technique présentée dans:

Coupling of local and global quantities in various finite element formulations
and its application to electrostatics, magnetostatics and magnetodynamics 
Dular, P.; Legros, W.; Nicolet, A.; 
Magnetics, IEEE Transactions on , Volume: 34 Issue: 5 , Sep 1998 
Page(s): 3078 -3081


Au sujet de ta définition du courant pour un calcul en post-processing, i.e.

      { Name I1; 
        Value { 
          Integral { [ -sigma[] * {Grad v} ]; In Front_Ele1; 
                              Integration Int_1; Jacobian Surf; }
        }
      }

le problème est que la partie du grad v dont tu as besoin est le grandient normal à Front_Ele1. L'interpolation de grad v sur Front_Ele1 ne rend que sa composante tangentielle, c'est-à-dire qu'il est fait abstraction de la présence d'éléments volumiques et donc que le grandient n'est calculé qu'à partir du potentiel v sur les éléments de surface (ce qui donne zéro si v y est constant).

Une solution pour ce calcul en post-processing, en accord avec ce qui est proposé dans la méthode du papier ci-dessus, consiste à réaliser une integration en volume telle que appelée à la fin du fichier Plate.pro. Il faut cependant (avec la version actuelle de getdp) définir a priori tous les outils nécessaires, puisqu'il faut avoir accès à des fonctions de bases bien précises.

Bon travail et bons tests,

Patrick


> Subject: [Getdp] calcul de courant
> Date: Wed, 22 Jan 2003 17:47:34 +0100
> From: Ali ABAKAR <abakar at edf.fr>
> To: getdp at geuz.org
> 
> Bonjour,
> Nouvel utilisateur de GetDp, je voudrais calculer le courant global I
> dans un conducteur
> en électrocinetique. Après plusieurs efforts et malgres le manuel
> utilisateur et les contributions
> du forum, je n'y arrive toujours pas!
> La résolution se fait correctement et j'ai les valeurs de la densité de
> courant dans les differentes
> régions. Mais l'integrale qui devrait donner le courant global ne donne
> rien :((
> D'autre part est il possible de fixer un Dirichlet en complexe (pour
> prendre en compte un dephasage)
> si oui, comment fait on ?
> Ci joint les fichiers concernés.
> Merci d'avance et à bientôt
> Ali
> 
> --
> Office:
>  Ali Abakar                       Ph:  33 1 47 65 34 67
> EDF Division R&D                        Fax: 33 1 47 65 32 51
> 1 av. General de Gaulle                 e-mail: abakar at edf.fr
> 92141 Clamart
> 
> _______________________________________________
> GetDP mailing list
> GetDP at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

-- 
Patrick Dular, Dr. Ir., Chercheur qualifié F.N.R.S. (Research associate)
Département d'Electricité, Electronique et Informatique
Unité d'Electricité Appliquée
Université de Liège - Institut Montefiore - Bât. B28 - Parking 32
B-4000 Liège - Belgique - Tel. +32-4 3663710 - Fax +32-4 3662910
E-mail: Patrick.Dular at ulg.ac.be
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Ele_examples.tar.gz
Type: application/x-gzip
Size: 11703 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20030124/a5e28b5e/attachment.gz>