Laurent Berquez

Olivier Eichwald
Vincent Bley
Hubert Caquineau
Yann Cressault
COMSOL V5.5
1. Présentation du problème .................................................................................................... 3
2. Présentation de COMSOL ..................................................................................................... 3
3. Assistant de modélisation ..................................................................................................... 4
4. Définition globale ................................................................................................................. 6
5. Construction de la géométrie................................................................................................ 6
6. Définition des matériaux ...................................................................................................... 7
7. Définition de la physique sur les domaines ........................................................................... 8
8. Définition des conditions aux limites sur les frontières ......................................................... 9
9. Maillage ............................................................................................................................... 9
10. Résolution....................................................................................................................... 10
11. Post traitement ............................................................................................................... 10
12. Modification du maillage ................................................................................................ 12
12.1. Raffinement du maillage ........................................................................................................................................ 12
12.2. Raffinement du maillage dans les zones proches des aretes.................................................................. 13
13. Modification des fonctions d’interpolation ..................................................................... 14
14. Modification de la géométrie .......................................................................................... 14
15. Application à un autre problème physique ...................................................................... 15

2
1. PRESENTATION DU PROBLEME
On se propose de calculer le champ électrique entre 2 électrodes cylindriques constituant un
condensateur plan.

d
e

2. PRESENTATION DE COMSOL
Le bureau Comsol est organisé autour de multiples fenêtres (Constructeur du modèle, Réglages,
Graphique, Messages, Progressions …) qui peuvent être affichées ou non en fonction du besoin.

Un clic droit sur un nœud de l’arbre du modèle, permet d’afficher un menu contextuel présentant
toutes les fonctionnalités disponibles que l’on retrouve dans les menus et barres d’outils.

3
Si vous choisissez une fonction qui requiert des informations, une fenêtre des réglages permet de
spécifier les paramètres de cette fonction.
Chaque action permettant la modélisation est enregistrée dans l’arbre du modèle et peut être
modifiée à tout moment.

3. ASSISTANT DE MODELISATION
1. Dans le dossier Applications, aller dans le
dossier COMSOL55 puis dans le dossier
Multiphysics et double cliquer sur Comsol
Multiphycics (Classkit License) puis cliquer
sur Assistant de Modélisation.

4
2. Dans l’assistant de modélisation sélectionner la
dimension de l’espace 2D axisymétrique.
Justifier les choix de cette dimension d’espace.

3. Dans la fénêtre Sélectionner la physique,

développer AC/DC, développer champs et
courants électriques puis selectionner
Electrostatique et cliquer sur Ajouter

4. Cliquez sur le bouton Etude

5. Dans la fénêtre Sélectionner l’étude,

sélectionner Stationnaire puis cliquer sur
le bouton Terminé .
Les études prédéfinies ont un solveur et des
réglages adaptés à la physique sélectionnée.

5
4. DEFINITION GLOBALE
La branche « Définitions globales » de l’arbre du modèle permet de définir des paramètres, des
variables et des fonctions. Si on veut par la suite faire une étude paramétrique sur la géométrie, il
faut définir la géométrie à partir de paramètres. On peut
également définir des paramètres permettant de contrôler le
maillage ou certaines propriétés physiques du matériau.

1. Cliquer sur Paramètres.

2. Dans le tableau de paramètres, cliquer

sur la première ligne sous « nom » et
entrez d. Sur la première ligne sous
« Expression » entrer la valeur de d :
1[mm].
On peut entrer n’importe quelle unité
entre les crochets.

3. Entrer les autres paramètres : e, Re, c, U

Afin de ne pas oublier la signification des variables, il est en effet pertinent d’entrer une description
des variables.

5. CONSTRUCTION DE LA GEOMETRIE
1. Sur Géométrie 1 vous pouvez choisir l’unité de longueur m
ou mm
2. Dans l’arbre du modèle, sous « Composant 1 » faire un clic
droit sur « géométrie » et sélectionner Rectangle

3. Dans la zone Taille et forme de la fenêtre de réglage, entrer

Re pour la largeur et e pour la hauteur

4. Dans la zone position, entrer 0 pour r et –e-d/2 pour z

5. Cliquer sur Construire jusqu’à la sélection (premier bouton
en haut à droite)

Vous venez de créer le premier rectangle nommé rectangle 1

6. Répéter la même opération (refaire un clic sur géométrie,

rectangle etc..), pour créer un rectangle nommé « rectangle
2 » de largeur Re et de hauteur e. Il est positionné à d/2 en z
et 0 en r

6
7. Répéter la même opération, pour créer un rectangle nommé « retangle 3 » de largeur 3*Re et de
hauteur 10*(2*e+d). Il est positionner à -5*(2*e+d) en z et 0 en r

On obtient alors le tracé de 3 rectangles dans la fenêtre graphique. Les 2 premiers réprésentent les
électrodes et le troisième les limites du domaine de calcul.

Le calcul du potentiel électrique va se faire sur un domaine correspondant au troisième rectangle

sans les 2 premiers. On va donc faire une opération booléene pour définir ce domaine de calcul.

8. Dans l’arbre du modèle, sous « modèle 1 » faire

un clic droit sur « géométrie » et sélectionner
Booléens et partitions > Différences

9. Sélectionner le rectangle 3 (r3) dans la fenêtre

graphique pour l’ajouter dans la zone objet à
ajouter

10. Cliquer sur Actif près de la zone objets à

soustraire et sélectionner les rectangles 1 et 2 (r1
et r2) dans la fenêtre graphique pour les ajouter
dans la zone objet à soustraire

11. Cocher Conserver les objets d’origine

12. Terminer la fonction en cliquant sur Construire

jusqu’à la sélection

On a donc maintenant construit 3 domaines, l’air et

les 2 électrodes, pensez à sauvegarder votre
modèle.

6. DEFINITION DES MATERIAUX

1. Dans l’arbre du modèle, sous « Composant 1 »
faire un clic droit sur « Matériaux » et
sélectionner Ajouter un matériau

7
2. Dans la fenêtre Ajouter un matériau,
rechercher Air, faire un clic droit et
sélectionner Ajouter à Composant 1

3. Faire la même chose avec le matériau Copper

Il faut maintenant affecter le bon matériau au bon domaine

4. Dans l’arbre du modèle, cliquer sur « Composant 1 > Matériaux > Air »

5. Dans la fenêtre réglage, sélectionner le domaine 2

(l’électrode inférieure) et cliquer sur pour enlever
l’affectation air au domaine 2

6. Faire la même chose pour le domaine 3.

Il faut maintenant affecter le matériau copper aux

domaines 2 et 3

7. Dans l’arbre du modèle, cliquer sur « modèle 1 > matériaux > Copper »

8. Dans la fenêtre réglage, sélectionner le domaine 2 (l’électrode inférieure) et cliquer sur pour
affecter copper au domaine 2

9. Faire la même chose pour le domaine 3.

7. DEFINITION DE LA PHYSIQUE SUR LES DOMAINES

Les équations régissant les phénomènes d’électrostatique n’interviennent que dans le domaine 1
(air). Elles ne s’appliquent par sur les domaines 2 et 3 (Electrodes) qui sont conductrices.

1. Dans l’arbre du modèle, cliquer sur

« Composant 1 > Electrostatique»

2. Dans la fenêtre réglage, sélectionner le

domaine 2 (l’électrode inférieure) et cliquer sur
pour supprimer la simulation du potentiel
sur l’électrode

3. Faire la même chose pour le domaine 3

8
8. DEFINITION DES CONDITIONS AUX LIMITES SUR LES FRONTIERES
Dans Comsol, les conditions aux limites sont regroupées par type de conditions (axysymétrie, valeur
donnée d’une grandeur sur la frontière etc..)
Par défaut, il y a une axisymétrie sur les frontières de l’axe de symétrie et une charge nulle sur les
autres frontières.
On va donc garder une charge nulle sur les limites de la zone de calcul et ajouter un potentiel nul
(masse) sur l’électrode inférieure et une tension U sur l’électrode supérieure.

1. Dans l’arbre du modèle, sous « Composant 1 »

faire un clic droit sur « Electrostatique» et
sélectionner Potentiel électrique

2. Dans la fenêtre réglage, sélectionner les

contours de l’électrode supérieure (8, 10 et 13)
pour affecter une condition au limite de
potentiel électrique sur ces arêtes.

3. Saisir la valeur U pour le potentiel électrique V

4. Faire de même affecter une condition aux

limites de masse sur les frontières de
l’électrode inférieure

9. MAILLAGE
Dans un premier temps, nous allons utiliser les options de maillage proposées par défaut. Nous
ferons par la suite des modifications.

1. Dans l’arbre du modèle, sous « Composant

1 » cliquer sur « Maillage 1»
2. Garder l’option maillage contrôlé par la
physique et l’option fin pour la taille de
l’élément
3. Générer le maillage en cliquant sur Construire
toute la séquence

9
10. RESOLUTION
1. Dans l’arbre du modèle, sous « Composant 1 » cliquer sur « Etude 1 »
2. Cliquer sur Calculer

11. POST TRAITEMENT

Après calcul, des graphiques par défaut sont générés : dans notre cas, une représentation du
potentiel électrique dans le plan méridien.

Il faut maintenant avoir un regard critique sur la solution ! Pour cela, il faut se poser quelques
questions.

10
 - Le potentiel est il nul sur les contours de l’électrode inférieure ?
- Le potentiel est il égal à la tension appliquée sur l’électrode supérieure ?
- Le potentiel est-il calculé uniquement dans l’air ?
- Le potentiel calculé est-il « continu » ?

Nous allons maintenant afficher le champ électrique

1. Dans l’arbre du modèle, sous « Composant 1 »

cliquer avec le bouton droit sur « Résultat »
2. Sélectionner Groupe de graphique 2D
3. Dans l’arbre du modèle, sous « Composant 1 »
cliquer avec le bouton droit sur « Résultat >
Groupe de graphiques 2D 2 »
4. Sélectionner Renommer
5. Saisir Champ électrique (es) dans le boite de
dialogue
6. Dans l’arbre du modèle, sous « Composant 1 » cliquer avec le bouton droit sur « Résultat >
Champ électrique (es) »
7. Sélectionner Surface
8. Dans la fenêtre de réglage, cliquez sur remplacer l’expression
9. Choisir Composant 1 > Electrostatique > Electrique > Champ électrique > es.EZ - Champ
électrique, composante z
10. Cliquer sur afficher

On obtient alors une représentation du

champ électrique qu’il faut critiquer !

- Le champ est-il constant entre les

électrodes ?
- Le champ est-il nul en dehors de
la zone inter électrode ?
- Y a-t-il des zones de variation de
champ ?

11
Si on zoome maintenant sur la zone à l’extrémité des
électrodes, on voit manifestement un souci de calcul !

Cela s’explique par 2 raisons :

- Les coins des électrodes génèrent des zones de
variation rapide du champ, il faut donc affiner le
maillage dans les zones proches de ces arrêtes.
Dans la réalité, ces coins existent ils ? (faire un
affichage en fil de fer pour comprendre)
- Le champ électrique est dérivé à partir du
potentiel, a-t-on utilisé des éléments qui
garantissent la continuité du gradient de potentiel
entre les éléments ?

11. Avant d’aller plus loin, sauvegarder les images des deux champs de grandeur (potentiel et champ
électrique suivant z) que vous avez calculés.

12. MODIFICATION DU MAILLAGE

On se propose de tester 2 configurations de maillage :
- Raffinement du maillage
- Raffinement du maillage dans des zones proches des arêtes

12.1. RAFFINEMENT DU MAILLAGE

1. Dans l’arbre du modèle, sous « Composant 1 » cliquer sur « Maillage»
2. Garder l’option maillage contrôlé par la physique et l’option extrêmement fin pour la taille de
l’élément
3. Générer le maillage en cliquant sur Construire toute la séquence
4. Relancer le calcul de la solution.

12
Il semble y avoir une amélioration mais pas très convaincante et on est passé de 1648 éléments à
43345 éléments. Sauvegarder l’image du champ que vous avez obtenue.

12.2. RAFFINEMENT DU MAILLAGE DANS LES ZONES PROCHES DES ARETES

1. Modifier l’option maillage contrôlé par la
physique par maillage contrôlé par l’utilisateur
2. Dans l’arbre du modèle, sous « Composant 1 »
cliquer sur « Maillage 1 > Taille » et renseigner
les paramètres de taille d’élément avec les
paramètres ci contre.

3. Dans l’arbre du modèle, sous « Composant 1 »

cliquer sur « Maillage > Triangle libre 1 » et
utiliser Delaunay en méthode de pavage.
4. Générer le maillage en cliquant sur Construire
toute la séquence

13
5. Relancer le calcul de la solution.

L’amélioration du calcul du champ est visible avec une zone de renforcement du champ sur les
coins des électrodes. On est passé à 78949 éléments.

6. Sauvegarder l’image du champ que vous avez obtenue.

7. Raffiner encore le maillage en diminuant la taille de l’élément maximum (2.5e-4) et la résolution

des régions étroites à 60. Relancer un calcul et conclure quant au compromis taille de maillage /
précision / temps de calcul.

13. MODIFICATION DES FONCTIONS D’INTERPOLATION

Il est possible pour un maillage donné de modifier les fonctions d’interpolation.

1. Dans l'arbre du modèle, cliquer sur Electrostatique, dans la fenêtre de réglage sélectionner
linéaire pour la discrétisation du potentiel électrique.
2. Relancer les calculs demandés précédemment et observer les changements.

14. MODIFICATION DE LA GEOMETRIE

On a vu que le calcul du champ électrique dans les zones proches des arêtes des électrodes posait
problème. La question est : A-t-on trop simplifié la géométrie des électrodes ?
On se propose donc d’ajouter un congé au niveau des arêtes.

1. Dans l’arbre du modèle, sous « Composant 1 » cliquer sur « Géométrie 1 > Différence 1»
14
2. Décocher Conserver les objets d’origine
3. Terminer la fonction Construire jusqu’à la sélection
4. Faire un clic droit sur « Géométrie 1» et sélectionner Congé
5. Ajouter les quatres sommets des électrodes dans la zone Points à raccorder par un congé
6. Utiliser la variable c pour le rayon
7. Utiliser les différents paramètres de maillage et comparer les solutions

15. APPLICATION A UN AUTRE PROBLEME PHYSIQUE

Un industriel souhaite à présent étudier le phénomène de claquage d’un plasma entre deux
électrodes en fonction de la géométrie utilisée. Il souhaite réaliser ce claquage pour une tension U
entre les électrodes la plus faible possible. Sachant qu’à la température ambiante, et à la pression
atmosphérique le champ de claquage dans l’air est d’environ 30 kV/cm, il souhaite faire un choix
entre les deux configurations présentées sur la figure ci – dessous.

10 mm
Electrode
supérieure
10 mm 10 mm
5 mm 5 mm
5 mm 3mm 5 mm
5 mm 5 mm
Electrode
Inférieure
V=0

Mettez en œuvre une simulation qui puisse l’aider à réaliser son choix et à quantifier par exemple
quelle tension il devra appliquer dans chaque cas pour obtenir ce claquage.

15