Calcul numérique de capacité
Calcul numérique de capacité
Bonjour,
Je dois déterminer numériquement la capacité lineique d'un câble coaxial de forme cylindrique. Ainsi, j'ai considéré une tranche de ce câble avec un potentiel de 1V au niveau du conducteur central et 0V sur le deuxieme conducteur. Avec la methode des differences finies appliquée à l'équation de Laplace, j'ai déterminé une carte des potentiels sur la section.
Cependant, je dois maintenant en déduire la capacité lineique mais je ne vois pas comment faire. J'ai essayé d'utiliser le théorème de gauss, mais je ne vois pas comment calculer le flux du champ E à partir des potentiels.
Si quelqu'un peut m'expliquer, je suis preneur.
Merci.
Je dois déterminer numériquement la capacité lineique d'un câble coaxial de forme cylindrique. Ainsi, j'ai considéré une tranche de ce câble avec un potentiel de 1V au niveau du conducteur central et 0V sur le deuxieme conducteur. Avec la methode des differences finies appliquée à l'équation de Laplace, j'ai déterminé une carte des potentiels sur la section.
Cependant, je dois maintenant en déduire la capacité lineique mais je ne vois pas comment faire. J'ai essayé d'utiliser le théorème de gauss, mais je ne vois pas comment calculer le flux du champ E à partir des potentiels.
Si quelqu'un peut m'expliquer, je suis preneur.
Merci.
Taupe (201)5
Re: Calcul numérique de capacité
bonjour,
le calcul est faisable analytiquement (pour un cylindre infini en tout cas)
mais numériquement il "suffit" de calculer le gradient de V pour en déduire le champ électrique...
après en fonction du maillage que tu as pris il y a moyen que ça devienne pénible ...
sinon tu peux regarder du côté de femm qui doit permettre de faire cela
le calcul est faisable analytiquement (pour un cylindre infini en tout cas)
mais numériquement il "suffit" de calculer le gradient de V pour en déduire le champ électrique...
après en fonction du maillage que tu as pris il y a moyen que ça devienne pénible ...
sinon tu peux regarder du côté de femm qui doit permettre de faire cela
Sciences Physiques,MP*-ex PSI* Corneille Rouen
Re: Calcul numérique de capacité
Oui, j'ai la formule analytique, mais je voulais tester une méthode numérique.
J'ai calculé le gradient et j'obtient des valeurs de C moyennement précises.
La methode des éléments finis m'a l'air un peu trop compliquée pour moi actuellement.
Merci de votre aide.
J'ai calculé le gradient et j'obtient des valeurs de C moyennement précises.
La methode des éléments finis m'a l'air un peu trop compliquée pour moi actuellement.
Merci de votre aide.
Taupe (201)5
Re: Calcul numérique de capacité
essaie de voir ce qui se passe en resserrant le maillage
en pratique tu pars d'un maillage assez lache et tu pars de la solution trouvée pour reserrer le maillage (ça devrait converger plus vite)
bon le mieux ce serait un maillage multi échelle mais là c'est une autre paire de manche ...
en pratique tu pars d'un maillage assez lache et tu pars de la solution trouvée pour reserrer le maillage (ça devrait converger plus vite)
bon le mieux ce serait un maillage multi échelle mais là c'est une autre paire de manche ...
Sciences Physiques,MP*-ex PSI* Corneille Rouen
Re: Calcul numérique de capacité
Je vais essayer de resserer le maillage. Le problème est qu'en agrandissant ma matrice, les calculs commencent à prendre un bon moment à tourner.
Taupe (201)5
Re: Calcul numérique de capacité
Tu utilises quoi comme algorithme/code/logiciel ?
Sinon oui, E = -grad V.
Ta valeur de C change-t-elle beaucoup lorsque tu changes le maillage ? Si oui c'est que celui-ci n'est pas assez petit. Idem si tu changes le critère de convergence (tu dois bien avoir un truc qui te dit quand stopper le code ?).
Sinon oui, E = -grad V.
Ta valeur de C change-t-elle beaucoup lorsque tu changes le maillage ? Si oui c'est que celui-ci n'est pas assez petit. Idem si tu changes le critère de convergence (tu dois bien avoir un truc qui te dit quand stopper le code ?).
Re: Calcul numérique de capacité
Jusque là, j'ai utilisé des tableaux sur python et j'applique par récursivité: V(i, j) (n) =1/4(V(i+i,j)(n)+V(i-1,j)(n)+...).
Mais j'ai pas pas fixé de condition d'arrêt de l'algorithme, je lui dit juste combien d'itérations il doit faire.
Si je lui fixe une condition j'ai peur qu'il ne s'arrête jamais car les valeurs de C que je trouve au bord du deuxième conducteur me semblent assez fausses..
Mais j'ai pas pas fixé de condition d'arrêt de l'algorithme, je lui dit juste combien d'itérations il doit faire.
Si je lui fixe une condition j'ai peur qu'il ne s'arrête jamais car les valeurs de C que je trouve au bord du deuxième conducteur me semblent assez fausses..
Taupe (201)5
Re: Calcul numérique de capacité
En général pour ce genre d'algorithme on définit $ e = sum_{i,j}|V_\text{iteration~k}[i,j]-V_\text{iteration~k-1}[i,j]| $, et on indique à l'algorithme de s'arreter lorsque $ e $ est inférieur à une certaine valeur epsilon. 10^-4 par exemple.
Il faut au minimum faire une étude de la convergence des résultats, c-a-d est-ce qu'ils dépendent de epsilon. Pour cela calcule les valeurs de C pour epsilon = 10^-5, 10^-4, 10^-3, et voit si ça varie beaucoup ou pas.
Il faut faire pareil avec la taille de la grille.
Après pour accélérer la convergence il y a des techniques, mais déjà commence par ça.
Il faut au minimum faire une étude de la convergence des résultats, c-a-d est-ce qu'ils dépendent de epsilon. Pour cela calcule les valeurs de C pour epsilon = 10^-5, 10^-4, 10^-3, et voit si ça varie beaucoup ou pas.
Il faut faire pareil avec la taille de la grille.
Après pour accélérer la convergence il y a des techniques, mais déjà commence par ça.