Simulation d'Écoulement dans un Tube en T

Simulation d'Écoulement dans un Tube en T : Mélange d'Eau Chaude et Froide avec SolidWorks Flow Simulation

 

Dans cet article, nous explorons en profondeur la simulation d'écoulement dans un tube en T, un composant courant en ingénierie des fluides, où de l'eau chaude à 120°C et de l'eau froide à 20°C se mélangent avec des débits massiques égaux de 2 kg/s. Basé sur une vidéo que j'ai créée utilisant SolidWorks Flow Simulation, nous décortiquerons le processus de modélisation, les principes physiques sous-jacents, des calculs analytiques pour la température et la vitesse, ainsi que des codes MATLAB pour reproduire ces analyses. Nous inclurons également des graphiques illustrant la répartition de la vitesse et de la température le long de deux axes du tube. Cette analyse vise à optimiser la compréhension des phénomènes de mélange en CFD (Computational Fluid Dynamics), en comblant les lacunes souvent observées dans les ressources en ligne, telles que le manque de calculs théoriques et de visualisations détaillées.

 Figure 1: Tube en T

Comprendre les Principes de l'Écoulement dans un Tube en T

Un tube en T, ou jonction en T, est une configuration de tuyauterie où deux flux entrants convergent vers une sortie unique. Dans notre cas, l'une des entrées reçoit de l'eau chaude à 120°C et l'autre de l'eau froide à 20°C, chacune avec un débit massique de 2 kg/s. Les dimensions du tube, basées sur la pièce jointe (diamètre interne de 26 mm, épaisseur de 3 mm, branches d'entrée de 100 mm de longueur et branche de sortie de 200 mm), influencent grandement le comportement de l'écoulement.

L'écoulement dans un tel système implique plusieurs phénomènes physiques. D'abord, le mélange convectif : les deux flux se heurtent au niveau de la jonction, créant des zones de turbulence qui favorisent l'échange thermique. La diffusion moléculaire joue également un rôle, bien que mineur comparé à la convection forcée. En termes de régime d'écoulement, calculons le nombre de Reynolds (Re) pour déterminer si l'écoulement est laminaire ou turbulent.

La formule du nombre de Reynolds est : $Re = \frac{\rho v D}{\mu}$ où $\rho$ est la densité de l'eau (environ 1000 kg/m³), $v$ la vitesse moyenne, $D$ le diamètre interne (0.026 m), et $\mu$ la viscosité dynamique (environ 0.001 Pa·s pour l'eau à température ambiante).

Le débit massique $\dot{m} = 2$ kg/s par entrée. La section transversale $A = \pi (D/2)^2 \approx 5.31 \times 10^{-4}$ m². Ainsi, la vitesse $v = \dot{m} / (\rho A) \approx 3.77$ m/s.

En remplaçant : $Re \approx \frac{1000 \times 3.77 \times 0.026}{0.001} \approx 98,000$

Puisque Re > 4000, l'écoulement est turbulent, ce qui implique des vortex et une meilleure homogénéisation du mélange. Sans simulation, on peut estimer la température de sortie par bilan énergétique, assumant un mélange parfait et des capacités thermiques spécifiques égales (cp ≈ 4180 J/kg·K pour l'eau) : $T_{sortie} = \frac{\dot{m}_{chaude} T_{chaude} + \dot{m}_{froide} T_{froide}}{\dot{m}_{total}} = \frac{2 \times 120 + 2 \times 20}{4} = 70^\circ C$

Cependant, en réalité, des gradients persistent en raison de la géométrie.

Configuration de la Simulation dans SolidWorks Flow Simulation

Pour reproduire cette simulation, comme dans ma vidéo, commencez par modéliser le tube en T dans SolidWorks. Utilisez des esquisses pour créer les cercles concentriques (diamètre externe 32 mm, interne 26 mm) et extrudez les branches. Fermez les extrémités avec des couvercles minces (1 mm) pour définir les frontières.

Paramètres du Projet Flow Simulation

Activez l'add-in Flow Simulation via Outils > Add-Ins. Créez un nouveau projet avec l'assistant :

• Type d'analyse : Interne (écoulement dans le tube).
• Fluide : Eau (liquide).
• Conditions initiales : Température ambiante, pression atmosphérique.
• Modèle de turbulence : k-ε pour capturer les effets turbulents.

• Activez « Heat conduction in solids » si vous voulez étudier le transfert dans les parois. 

Définissez les conditions aux limites :

• Entrée chaude : Débit massique 2 kg/s, température 120°C.
• Entrée froide : Débit massique 2 kg/s, température 20°C.
• Sortie : Ouverture à pression environnementale.

Générez un maillage automatique avec raffinement local à la jonction (niveau 5 pour précision). Lancez la simulation ; elle converge après environ 100 itérations, en fonction de la puissance de calcul.

Cette configuration permet de visualiser l'écoulement via des trajectoires de flux, des coupes planes pour les contours de température et de vitesse, et des animations montrant le mélange progressif.

Analyse des Résultats de la Simulation

Les résultats de SolidWorks Flow Simulation révèlent un mélange asymétrique initial à la jonction, avec une stratification : l'eau chaude (densité moindre) tend à monter, tandis que l'eau froide descend. La température moyenne à la sortie approche 70°C, mais avec des variations locales de ±10°C dues à la turbulence.

Pour la vitesse, le profil est parabolique dans les branches d'entrée (vitesse max au centre ~4.5 m/s), mais devient irrégulier à la jonction en raison des vortex. La perte de charge totale est estimée à environ 50 kPa, calculée via la différence de pression statique.

Pour illustrer, considérons la répartition le long de deux axes :

• Axe horizontal (branche de sortie, de 0 à 0.2 m) : La température décroît exponentiellement de 100°C (près de la jonction) à 70°C, reflétant l'homogénéisation. La vitesse montre un profil parabolique modulé par des oscillations turbulentes.

 

 Figure 2: La température le long du grand axe 

 

• Axe vertical (à travers la jonction, de -0.1 m à 0.1 m) : La température transite de 20°C (froide) à 120°C (chaude) via une fonction tanh-like, indiquant une zone de mélange abrupte. La vitesse oscille avec des pics aux entrées.

 

 Figure 3: La température le long du petit axe

 

Ces distributions sont basées sur des données simulées ; en pratique, exportez-les depuis SolidWorks pour un post-traitement.

Calculs Théoriques et Codes MATLAB pour l'Analyse

Pour valider les résultats, effectuons des calculs analytiques. Par exemple, la vitesse moyenne à la sortie : avec $\dot{m}_{total} = 4$ kg/s et même section, $v_{sortie} \approx 7.54$ m/s.

Utilisons MATLAB pour modéliser ces distributions. Voici un code simple pour calculer et plotter la température le long de l'axe horizontal (modèle exponentiel simplifié) :

Visualisation du Champ de Température

% Code MATLAB pour température le long axe horizontal
x = linspace(0, 0.2, 100); % Position en m
T0 = 100; % Temp initiale près jonction
Tf = 70; % Temp finale
k = 20; % Constante de décroissance
T = Tf + (T0 - Tf) * exp(-k * x);
figure;
plot(x, T);
title('Répartition de la Température le Long de l''Axe Horizontal');
xlabel('Position (m)');
ylabel('Température (°C)');
grid on;

 

Explication : Ce code génère un vecteur de positions x, calcule la température via une décroissance exponentielle (modélisant la diffusion), et trace le graphique. Pour arriver à cette solution, partez de l'équation de transport de chaleur simplifiée en 1D : dT/dx = -k (T - Tf), résolue analytiquement.

De même, pour la vitesse le long de l'axe vertical (modèle oscillant) :

% Code MATLAB pour vitesse le long axe vertical
y = linspace(-0.1, 0.1, 100); % Position en m
v_max = 4; % Vitesse max
v = v_max * abs(cos(5 * pi * y)); % Profil oscillant
figure;
plot(y, v);
title('Répartition de la Vitesse le Long de l''Axe Vertical');
xlabel('Position (m)');
ylabel('Vitesse (m/s)');
grid on;

 

Ce code simule des oscillations dues à la turbulence. La fonction abs(cos) crée des pics périodiques, représentant des vortex. Pour la solution, intégrez numériquement ou analytiquement à partir des équations de Navier-Stokes simplifiées.

Enfin, un code pour calculer Re et T_sortie :

% Calculs basiques
rho = 1000; % kg/m3
mu = 0.001; % Pa.s
D = 0.026; % m
mdot = 2; % kg/s par entrée
A = pi * (D/2)^2;
v = mdot / (rho * A);
Re = rho * v * D / mu;
T_chaude = 120;
T_froide = 20;
T_sortie = (T_chaude + T_froide) / 2;
disp(['Reynolds: ', num2str(Re)]);
disp(['Temp sortie: ', num2str(T_sortie), ' °C']);

 

Exécutez-le pour obtenir Re ≈ 98020 et T_sortie = 70°C.

 

Optimisation et Applications Pratiques

En intégrant ces éléments, cette simulation avec SolidWorks Flow Simulation offre des insights précieux pour des applications comme les systèmes de chauffage, les mélangeurs industriels ou les réseaux de distribution d'eau. Pour optimiser davantage, variez les débits ou ajoutez des baffles pour améliorer l'homogénéité.

En conclusion, cette analyse dépasse les tutoriels basiques en ajoutant des calculs rigoureux et des outils comme MATLAB, rendant l'article un guide complet pour les ingénieurs et étudiants en mécanique des fluides.

 

Lien vers le téléchargement du tube: Tube 

 

 Autres articles utiles:

Simulation d'Écoulement dans un Tube

Essai de traction non linéaire | SolidWorks Simulation

 

Conclusion

Cette simulation SolidWorks Flow Simulation d'un tube en T avec mélange eau chaude (120 °C) et eau froide (20 °C) à débit égal (2 kg/s) démontre la puissance de la CFD pour prédire température, vitesse et qualité du mélange. Les résultats confirment un régime turbulent intense (Re ≈ 178 000), une température de sortie proche de 70 °C et une homogénéisation progressive sur la branche aval.

Grâce aux calculs analytiques, aux visualisations et aux codes MATLAB fournis, vous pouvez reproduire et approfondir cette étude. Que ce soit pour un projet universitaire, une optimisation industrielle ou simplement pour mieux comprendre les phénomènes de mélange thermique, cet exemple concret illustre parfaitement l'intérêt de coupler modélisation 3D et simulation numérique.

N'hésitez pas à regarder la vidéo associée et à tester vous-même sur votre pièce ! Partagez vos résultats en commentaire.

FAQ – Questions Fréquentes sur la Simulation de Mélange dans un Tube en T

Quelle est la différence entre mélange parfait et mélange réel dans un tube en T ?

Dans un mélange parfait (théorique), la température est uniformément 70 °C dès la jonction. En réalité, des gradients persistent sur plusieurs diamètres à cause de la stratification et de la turbulence incomplète.

Pourquoi utiliser le modèle k-ε plutôt que k-ω SST dans SolidWorks Flow Simulation ?

Le k-ε est robuste et rapide pour les écoulements libres turbulents. Le k-ω SST est plus précis près des parois et dans les zones de recirculation, recommandé si vous analysez les pertes de charge ou les fluctuations thermiques.

Comment exporter les données de température et vitesse pour les tracer dans MATLAB ?

Utilisez « Surface Plot » ou « Flow Trajectories », puis « Export » → « Results » → format CSV. Importez ensuite avec csvread ou readtable dans MATLAB.

Peut-on simuler un débit massique différent entre les deux entrées ?

Oui ! Modifiez simplement les conditions limites (ex. 1.5 kg/s et 2.5 kg/s) → la température de sortie devient pondérée (ex. 62 °C dans cet exemple).

La température d’entrée à 120 °C est-elle réaliste pour de l’eau liquide ?

À pression atmosphérique, l’eau bout à 100 °C. Pour éviter l’ébullition, la simulation suppose une pression suffisante (typique dans les systèmes fermés) ou utilise de l’eau surpressurisée. SolidWorks gère bien cela via les tables thermodynamiques.