... | ... | @@ -115,8 +115,8 @@ |
|
|
courant sur le maillage d'entrée. La seconde image montre la même ligne
|
|
|
sur notre maillage interpolé.
|
|
|
|
|
|
![Image:Vit1 3d.001.png](Vit1_3d.001.png "Image:Vit1 3d.001.png")
|
|
|
![Image:Vit0 3d.001.png](Vit0_3d.001.png "Image:Vit0 3d.001.png")
|
|
|
![Vit1 3d.001.png](Vit1_3d.001.png "Vit1 3d.001.png")
|
|
|
![Vit0 3d.001.png](Vit0_3d.001.png "Vit0 3d.001.png")
|
|
|
|
|
|
### Projet particules
|
|
|
|
... | ... | @@ -160,8 +160,7 @@ |
|
|
égale à +5. (On ne peut mettre que des conditions de Dirichlet pour
|
|
|
l'instant).
|
|
|
|
|
|
![Image:Orzel2d 2trous
|
|
|
maillage.png](Orzel2d_2trous_maillage.png "Image:Orzel2d 2trous maillage.png")
|
|
|
![Orzel2d 2trous maillage.png](Orzel2d_2trous_maillage.png "Orzel2d 2trous maillage.png")
|
|
|
|
|
|
Afin de mieux comprendre ce qu'il se passe, j'ai cassé la symétrie en
|
|
|
imposant une source asymétrique dans l'épaisseur du tuyau. Pour ce
|
... | ... | @@ -169,8 +168,7 @@ |
|
|
selon l'axe y=-x, négatif en bas à droite, positif en haut à gauche. (La
|
|
|
formule est du type -x+y)
|
|
|
|
|
|
![Image:Orzel2d
|
|
|
2trous.png](Orzel2d_2trous.png "Image:Orzel2d 2trous.png")
|
|
|
![Orzel2d 2trous.png](Orzel2d_2trous.png "Orzel2d 2trous.png")
|
|
|
|
|
|
La température sur l'ensemble de l'élément est majorée par 0 et 10,
|
|
|
c'est-à-dire les températures extrêmes.
|
... | ... | @@ -182,8 +180,7 @@ |
|
|
Contrairement à l'exemple précédent, la température dépasse les valeurs
|
|
|
extrêmes, en négatif (bas-gauche) ou en positif (haut-droit).
|
|
|
|
|
|
![Image:Orzel2d
|
|
|
2trous2.png](Orzel2d_2trous2.png "Image:Orzel2d 2trous2.png")
|
|
|
![Orzel2d 2trous2.png](Orzel2d_2trous2.png "Orzel2d 2trous2.png")
|
|
|
|
|
|
Documents:
|
|
|
|
... | ... | @@ -248,8 +245,7 @@ |
|
|
intérieur, je dois donc obtenir des isothermes perpendiculaires à la
|
|
|
frontière. Ce qui est bien le cas :
|
|
|
|
|
|
![Image:Condition neumann
|
|
|
test.png](Condition_neumann_test.png "Image:Condition neumann test.png")
|
|
|
![Condition neumann test.png](Condition_neumann_test.png "Condition neumann test.png")
|
|
|
|
|
|
*Tests*
|
|
|
|
... | ... | @@ -264,8 +260,7 @@ |
|
|
cercle vaut (x,y)/R2, donc la derivée normale vaut 2\*R2). J'obtiens
|
|
|
ceci:
|
|
|
|
|
|
![Image:Condition neumann
|
|
|
test2.png](Condition_neumann_test2.png "Image:Condition neumann test2.png")
|
|
|
![Condition neumann test2.png](Condition_neumann_test2.png "Condition neumann test2.png")
|
|
|
|
|
|
Ce qui est déjà prometteur, car la température est distribuée
|
|
|
radialement entre 4\*4 et 10\*10.
|
... | ... | @@ -398,8 +393,7 @@ |
|
|
Pour traiter un exemple, j'ai choisi de modéliser la température d'une
|
|
|
pièce en L dont voici le maillage :
|
|
|
|
|
|
![Image:Exemple2d sans robin
|
|
|
maillage.png](Exemple2d_sans_robin_maillage.png "Image:Exemple2d sans robin maillage.png")
|
|
|
![Exemple2d sans robin maillage.png](Exemple2d_sans_robin_maillage.png "Exemple2d sans robin maillage.png")
|
|
|
|
|
|
L'intérieur représente l'air de la pièce, les trois petites surfaces
|
|
|
représentent des murs (oui, assez épais, disons qu'il s'agit d'un
|
... | ... | @@ -421,8 +415,7 @@ |
|
|
|
|
|
Voici ce que j'obtiens:
|
|
|
|
|
|
![Image:Exemple2d sans robin
|
|
|
result.png](Exemple2d_sans_robin_result.png "Image:Exemple2d sans robin result.png")
|
|
|
![Exemple2d sans robin result.png](Exemple2d_sans_robin_result.png "Exemple2d sans robin result.png")
|
|
|
|
|
|
On note que le gradient est bien plus grand dans les murs, qui jouent
|
|
|
évidemment le rôle d'isolant par rapport à l'extérieur. On observe bien
|
... | ... | @@ -430,6 +423,5 @@ |
|
|
|
|
|
*Exemple complet avec condition de robin*
|
|
|
|
|
|
![Image:Exemple2d avec robin
|
|
|
maillage.png](Exemple2d_avec_robin_maillage.png "Image:Exemple2d avec robin maillage.png")
|
|
|
![Exemple2d avec robin maillage.png](Exemple2d_avec_robin_maillage.png "Exemple2d avec robin maillage.png")
|
|
|
|
... | ... | @@ -435,6 +427,5 @@ |
|
|
|
|
|
![Image:Exemple2d avec robin
|
|
|
result.png](Exemple2d_avec_robin_result.png "Image:Exemple2d avec robin result.png")
|
|
|
![Exemple2d avec robin result.png](Exemple2d_avec_robin_result.png "Exemple2d avec robin result.png")
|
|
|
|
|
|
La température est alors plus homogène dans la pièce (Regarder la
|
|
|
surface où la température est entre 20 et 18 degrés)
|
... | ... | @@ -608,5 +599,5 @@ |
|
|
donc modélisé sous Mefisto un cylindre (le contenant), et un hexaèdre
|
|
|
oblong (résistance chauffante).
|
|
|
|
|
|
![Image:Exemple3d.cylindre.png](Exemple3d.cylindre.png "Image:Exemple3d.cylindre.png")
|
|
|
![Exemple3d.cylindre.png](Exemple3d.cylindre.png "Exemple3d.cylindre.png")
|
|
|
|
... | ... | @@ -612,5 +603,5 @@ |
|
|
|
|
|
![Image:Exemple3d.resistance.png](Exemple3d.resistance.png "Image:Exemple3d.resistance.png")
|
|
|
![Exemple3d.resistance.png](Exemple3d.resistance.png "Exemple3d.resistance.png")
|
|
|
|
|
|
J'ai triangulé ces deux surfaces..., j'ai 'amelioré' (option 30 du menu
|
|
|
surfaces) ces deux surfaces. J'ai ensuite créé un volume (option 8) avec
|
... | ... | @@ -626,4 +617,4 @@ |
|
|
**UPDATE** M.Perronnet m'a expliqué le problème : la surface inférieur
|
|
|
de l'hexaèdre de la résistance est confondu avec le disque inférieur du
|
|
|
cylindre, c'est ce qui pose problème à Mefisto. A la rigueur, si le
|
|
|
maillage correspondait, cela serait peut-être passé... |
|
|
\ No newline at end of file |
|
|
maillage correspondait, cela serait peut-être passé... |