v7.31.111 WTNV111 – Flux thermique sur un milieu poreux saturé#
Résumé :
On considère un problème tridimensionnel de thermique dans un milieu poreux saturé.
Ce test consiste à étudier l’effet de l’application d’un flux de température sur la face supérieure du modèle sur la distribution de la température dans l’élément. On se limite à l’étude du premier pas de temps. On bloque la pression et les déplacements.
Les modèles étudiés sont \(\mathrm{2D}\) plans (DPQ8 et DPTR6) et \(\mathrm{3D}\) volumique (HEXA20) avec un comportement linéaire.
La solution de référence est alors la solution d’un calcul en thermique pure (THER_LINEAIRE) avec Code_Aster.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
La solution de référence est la solution d’un calcul en thermique pure (THER_LINEAIRE) avec le code_Aster .
Grandeur de référence#
\(\mathrm{TEMP}\) : température
Grandeur et résultat de référence#
Point |
Type de valeur |
Instant \((s)\) |
Référence \((°C)\) |
\(C,D\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(2.9412\times {10}^{-1}\) |
\(A,B\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(\approx 0.\) |
Modélisation A#
Caractéristiques de la modélisation A#
Modélisation volumique : 3D_THM
1 maille HEXA20 de la modélisation 3D_THM : THM_HEXA20
Résultat de la modélisation A#
Discrétisation en temps : un seul pas de temps grand : \({10}^{13}s\) .
Le schéma en temps est implicite \((\vartheta =1)\) .
Nœud |
Type de valeur |
Instant \((s)\) |
Référence \((°C)\) |
Tolérance |
\(\mathrm{NO1}\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(2.9412\times {10}^{-1}\) |
\(1.\text{\%}\) |
\(\mathrm{NO20}\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(\approx 0.\) |
\({10}^{-6}\) |
Modélisation B#
Caractéristiques de la modélisation B#
Modélisation plane: D_PLAN_THM
1 maille DPQ8 de la modélisation D_PLAN_THM : THM_DPQ8
Résultat de la modélisation B#
Discrétisation en temps : un seul pas de temps grand : \({10}^{13}s\) .
Le schéma en temps est implicite \((\vartheta =1)\) .
Nœud / Point |
Type de Valeur |
Instant \((s)\) |
Référence \((°C)\) |
Tolérance |
\(\mathrm{N3}\text{/}\text{C}\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(2.9412\times {10}^{-1}\) |
\(1.\text{\%}\) |
\(\mathrm{N1}\text{}/\text{}A\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(\approx 0.\) |
\({10}^{-6}\) |
Modélisation C#
Caractéristiques de la modélisation C#
Modélisation plane: D_PLAN_THM
9
2 mailles DPTR6 de la modélisation D_PLAN_THM : THM_DPTR6
Résultat de la modélisation C#
Discrétisation en temps : un seul pas de temps grand : \({10}^{13}s\) .
Le schéma en temps est implicite \((\vartheta =1)\) .
Nœud / point |
Type de Valeur |
Instant \((s)\) |
Référence \((°C)\) |
Tolérance |
\(\mathrm{N3}\text{/}\text{C}\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(2.9412\times {10}^{-1}\) |
\(1.\text{\%}\) |
\(\mathrm{N1}\text{}/\text{}A\) |
\(\mathrm{TEMP}\) |
\({10}^{13}\) |
\(\approx 0.\) |
\({10}^{-6}\) |
Synthèse des résultats#
Les écarts observés entre la solution de référence et la solution Code_Aster sont très faibles. Il y donc une bonne concordance des résultats.