Table des matières

Prêt à explorer toutes les possibilités de GHEtool Cloud ?

Vous pouvez essayer GHEtool 14 jours gratuitement,
aucune carte de crédit n'est requise.

MuoviELLIPSE

Ce chapitre présente la modélisation du MuoviELLIPSE de MuoviTech et propose une introduction à la méthode des éléments de frontière pour le calcul de la résistance thermique effective d'un forage pour des sondes de forme irrégulière.

MuoviELLIPSE

Le MuoviELLIPSE est, comme son nom l'indique, un échangeur de chaleur elliptique mis au point par MuoviTech. À l'instar de son homologue, le TurboCollector, dont il a été question dans le chapitre précédent, le MuoviELLIPSE comporte de nombreuses petites ailettes sur sa surface interne. Ces ailettes sont orientées en alternance dans le sens horaire et antihoraire sur toute la longueur du tuyau. Agissant comme des turbulateurs passifs, elles sont conçues pour induire un écoulement turbulent à de faibles débits, améliorant ainsi le transfert thermique. Pour cette conception particulière de sonde, la transition vers la turbulence commence à environ Re = 1 850, au lieu de Re = 2 300 pour les tuyaux lisses.

Image du MuoviELLIPSE.
Image du MuoviELLIPSE.

Comme on peut clairement le constater, le MuoviELLIPSE se distingue des tuyaux ronds lisses traditionnels de deux façons :

  1. Le comportement du fluide à l'intérieur de la sonde est influencé par les ailettes internes et sa forme elliptique.
  2. Le transfert de chaleur à l'extérieur de la sonde, mais à l'intérieur du trou de forage, est influencé par la géométrie irrégulière de celui-ci.

Pour modéliser correctement le MuoviELLIPSE, ces deux aspects doivent être pris en compte. Le premier aspect peut à nouveau être traité à l'aide de la simulation numérique directe (DNS), comme cela a également été fait pour le TurboCollector dans le chapitre précédent. La méthodologie étant la même, seuls les résultats sont présentés ci-dessous.

L'autre différence majeure réside toutefois dans la géométrie de la sonde, qui modifie le transfert thermique au sein même du forage. Pour tenir compte de cet effet, on utilise la méthode des éléments de frontière (BEM). Ces deux aspects sont abordés ci-dessous.

Pour plus d'informations sur le MuoviELLIPSE lui-même, rendez-vous sur le site Site web de MuoviTech.

Développement de modèles

Afin de modéliser le MuoviELLIPSE, les corrélations entre le coefficient de frottement et le nombre de Nusselt sont d'abord examinées, sur la base des simulations DNS, puis le concept de la méthode des éléments de frontière (BEM) est présenté pour calculer le transfert thermique à l'intérieur du forage.

Corrélation pour le coefficient de frottement

Dans le graphique ci-dessous, le coefficient de frottement du modèle MuoviELLIPSE (désigné par ‘ DNS à ellipses alternées ’) est représenté en fonction des corrélations analytiques du coefficient de frottement, tant pour le régime laminaire que pour le régime turbulent.

Pour calculer le nombre de Reynolds, il faut disposer d'un diamètre caractéristique. Dans le cas d'une sonde circulaire, cela est simple, mais pour les géométries non circulaires, le diamètre hydraulique est utilisé à la place. Il s'agit du diamètre d'un tuyau circulaire fictif présentant le même comportement hydraulique que la géométrie non circulaire. Le diamètre hydraulique est défini comme suit :$$D_h=\frac{4A}{P}$$où $D_h$ est le diamètre hydraulique en (m), $A$ est l’aire de la section transversale en (m²), et $P$ est le périmètre mouillé en (m) de cette aire. Étant donné que le rapport $A/P$ est légèrement inférieur pour une sonde elliptique par rapport à une sonde circulaire, le diamètre hydraulique est également plus petit. Par conséquent, le nombre de Reynolds est légèrement plus élevé pour une sonde elliptique que pour un tuyau circulaire de même section et à débit égal.
L'aire d'une ellipse peut être calculée en multipliant le demi-grand axe $a$ par le demi-petit axe $b$ comme suit : $$A=\pi ab$$En revanche, il n'existe aucune formule analytique pour sa circonférence. Dans GHEtool, on utilise la deuxième approximation de Ramanujan pour calculer la circonférence $C$ comme suit :$$C \approx \pi(a+b) \left[ 1+\frac{3h}{10+\sqrt{4-3h}} \right]$$ où $$h=\frac{(a-b)^2}{(a+b)^2}$$
Corrélation des coefficients de frottement pour le calcul DNS du MuoviELLIPSE. (Source : (Hidman et al., 2026))
Corrélation des coefficients de frottement pour le calcul DNS du MuoviELLIPSE. (Source : (Hidman et al., 2026))

Le graphique ci-dessus présente le coefficient de frottement d'une sonde elliptique lisse (référence) ainsi que celui de la MuoviELLIPSE. Tout comme pour le TurboCollector, on constate clairement que le passage à la turbulence intervient plus tôt dans ce dernier cas, vers Re = 1 850. Tant en régime entièrement laminaire qu'en régime entièrement turbulent, le coefficient de frottement du MuoviELLIPSE est très similaire à celui de la sonde elliptique lisse.

Les corrélations exactes correspondant à la figure ci-dessus sont données par : $$w = \left(1 + exp\left[ -5 \left( \frac{Re-1850}{2300-1850} -0,5\right) \right] \right)^{-1}$$$$f=(1-w)\frac{65}{Re}+w\left[ -1,8log_{10}\left( \frac{6,9}{Re}\right) \right]^{-2}$$Pour plus de détails mathématiques, l’utilisateur est invité à se reporter à Hidman N. (2026).

Corrélation pour le nombre de Nusselt

Sur la figure ci-dessous, le nombre de Nusselt est représenté en fonction du nombre de Reynolds pour différents nombres de Prandtl.

Comme indiqué dans Partie 6.1, le nombre de Prandtl représente le rapport entre la diffusivité cinétique et la diffusivité thermique ou, pour simplifier, englobe à la fois les aspects thermiques et hydrauliques. Ce nombre peut varier entre 25 (carbonate de potassium, 30%) et 75 (MPG, 33%), l’éthanol (29%) présentant une valeur de 67 et le MEG (30%) une valeur de 36. Étant donné que le nombre de Prandtl influe également sur le nombre de Nusselt, les simulations thermohydrauliques sont réalisées pour plusieurs valeurs du nombre de Prandtl.
Corrélation du nombre de Nusselt pour le calcul DNS du MuoviELLIPSE. (Source : (Hidman et al., 2026))
Corrélation du nombre de Nusselt pour le calcul DNS du MuoviELLIPSE. (Source : (Hidman et al., 2026))

Sur le graphique ci-dessus, on observe la même transition vers le régime turbulent aux alentours de Re = 1 850. À des nombres de Reynolds plus élevés, le nombre de Nusselt converge vers la solution correspondant à une sonde elliptique lisse. Les différentes couleurs représentent différents nombres de Prandtl : le vert correspond à 20, le bleu à 40 et le rouge à 75.

Les corrélations exactes correspondant à la figure ci-dessus sont données par : $$Nu_{lam}^{reg}=\sqrt{Nu^2_{smooth, lam}+\left[ -0,321Re^{0,2}Pr^{0,21} \right]^2} pour Re < 1850$$$$Nu_{turb}^{reg}=\sqrt{Nu^2_{smooth, lam}+\left[ 1,96(Re-1849,9)^{0,295}Pr^{0,29} \right]^2} pour Re \leq 1850$$
Pour les valeurs de Re > 4 000, on utilise la corrélation de Gnielinski pour le nombre de Nusselt, avec un décalage constant $\delta$ afin de tenir compte des ailettes aux nombres de Reynolds élevés (tout comme avec le TurboCollector). Ce décalage est défini comme suit :$$\delta = Nu_{Gnielinksi}(4000) - Nu_{MuoviELLIPSE}(4 000)$$Pour plus de détails mathématiques, l’utilisateur est invité à se reporter à Niklas et al. (2026).

Méthode des éléments de frontière

Les deux corrélations précédentes décrivent le comportement thermohydraulique à l'intérieur du tuyau. Cependant, le deuxième défi consiste à modéliser l'interaction entre le tuyau et le forage, tout en tenant compte de la forme elliptique de ce dernier. Pour les tuyaux circulaires lisses classiques, les équations de transfert thermique interne peuvent être résolues de manière analytique. Or, cela n'est plus possible pour les géométries non circulaires.

Pour pallier cette limite, on a recours à une approche numérique sous la forme de la méthode des éléments de frontière (BEM).

L'application de la méthode des éléments de frontière aux champs de forages géothermiques peu profonds a été inspirée par le professeur Massimo Cimmino et développée en collaboration avec lui.

La méthode BEM est une technique numérique utilisée pour résoudre des équations aux dérivées partielles (EDP) linéaires, telles que celles régissant le transfert thermique. Dans le cas présent, elle transforme le problème bidimensionnel d'origine en un problème unidimensionnel équivalent défini le long des parois des tuyaux et de la paroi du forage. En termes simples, plutôt que de calculer l'ensemble du champ de température à l'intérieur du forage, il suffit de résoudre un problème de transfert thermique équivalent uniquement le long des surfaces des tuyaux et de la paroi du forage. Ce concept est illustré graphiquement ci-dessous.

On peut également envisager la question sous un autre angle. Supposons que l’on veuille calculer la propagation des ondes sonores à travers un lac. Une option consiste à discrétiser l’ensemble du volume en petits éléments et à les résoudre un par un. Cette méthode est extrêmement chronophage et correspond essentiellement à ce qui est fait dans les simulations DNS décrites ci-dessus. Une autre approche consiste toutefois à n’étudier que la surface de la source sonore. À l’aide des fonctions de Green, la description mathématique de l’ensemble du champ d’ondes peut être réécrite de manière à ne dépendre que de la surface de la source sonore. Cela rend le problème beaucoup plus facile à résoudre. C’est exactement ce que fait la méthode BEM.
Représentation graphique de la méthode des éléments de frontière (avec l'aimable autorisation de M. Cimmino).
Représentation graphique de la méthode des éléments de frontière (avec l'aimable autorisation de M. Cimmino).

Sur la figure ci-dessus, les différents points représentent les nœuds auxquels les équations de transfert thermique sont résolues. Les flèches indiquent les composantes tangentielles et normales du transfert thermique. En discrétisant ainsi la géométrie du tuyau, il devient possible de prendre en compte avec précision la forme elliptique réelle de la sonde.

Le principal inconvénient de la méthode des éléments de bord (BEM) réside dans le fait qu’elle est très gourmande en ressources de calcul et donc trop lente pour être utilisée directement dans GHEtool. Afin de rendre le modèle suffisamment rapide pour des simulations pratiques, un réseau neuronal artificiel (RNA) est entraîné à partir des résultats des simulations BEM. Cette approche combine le meilleur des deux mondes : un modèle précis et géométriquement représentatif pour calculer le transfert thermique au sein du forage, et un ANN qui permet d'effectuer ces calculs efficacement au sein de GHEtool.

Un Réseau neuronal artificiel (RNA)est une sous-catégorie du vaste domaine de l'IA. Le principe d’un réseau neuronal artificiel (RNA) est d’imiter le fonctionnement du cerveau humain : lorsque nous recevons une information sensorielle, qu’il s’agisse d’une odeur, d’un toucher ou d’un son, celle-ci est transmise aux neurones de notre cerveau, où le signal se propage de neurone en neurone jusqu’à ce que nous parvenions à une pensée, une action, une sensation particulière, etc. Ce comportement, qui consiste à partir d’un ensemble de données d’entrée et à passer par une série de neurones pour aboutir à une conclusion particulière, est exactement ce que nous essayons de modéliser avec un réseau neuronal artificiel. La figure ci-dessous présente une représentation schématique d’un réseau neuronal artificiel.

Représentation schématique d'un réseau neuronal artificiel. (Source : https://blog.roboflow.com/what-is-a-neural-network/)
Représentation schématique d'un réseau neuronal artificiel. (Source : https://blog.roboflow.com/what-is-a-neural-network/)

Selon l'architecture du réseau neuronal artificiel (RNA), le nombre de couches cachées et le nombre de neurones dans chaque couche peuvent varier. Dans ce cas, au niveau de chaque nœud (ou neurone), les données sont pondérées par la valeur associée à ce neurone, puis transmises aux neurones de la couche suivante. Ce processus est répété jusqu'à ce que la sortie soit atteinte.

Tout comme un bébé qui vient au monde et qui doit encore tout apprendre, le réseau neuronal ne peut rien faire dès sa mise en service. C'est pourquoi un tel modèle doit être entraîné, afin que tous les coefficients de pondération des différents nœuds soient correctement calibrés et puissent convertir les données d'entrée en une sortie correcte. Dans notre cas, cela consiste à prédire la résistance thermique d’un forage pour un diamètre de forage, une taille de MuoviELLIPSE, un espacement entre les tuyaux et une conductivité thermique du coulis donnés.

Sur la base des modèles évoqués ci-dessus, la résistance thermique effective du forage et la perte de charge du MuoviELLIPSE sont examinées ci-dessous.

Comportement du MuoviELLIPSE

Compte tenu des deux corrélations établies ci-dessus, la résistance thermique effective du forage et la perte de charge du MuoviELLIPSE sont abordées dans les sections suivantes.

Toutes les simulations ci-dessous sont réalisées à partir d'un forage d'une profondeur de 100 m et d'une profondeur d'enfouissement de 70 cm. Le diamètre du forage est de 140 mm, la conductivité thermique du coulis est de 1,5 W/(m·K) et celle du sol est de 2 W/(m·K). Les tuyaux sont placés exactement à mi-chemin entre le centre du forage et la paroi de celui-ci (c'est-à-dire à une distance de 35 mm du centre du forage pour un diamètre de forage de 140 mm). Le fluide utilisé est du MPG à 25 v/v% et à 5 °C. Tous les tuyaux ont une classe de pression PN16 (SDR11) et une conductivité thermique de 0,4 W/(mK). Sauf indication contraire, le diamètre des tuyaux est de 32 mm. Tout écart par rapport aux hypothèses ci-dessus est explicitement mentionné ci-dessous.

Résistance thermique effective du trou de forage

Le graphique ci-dessous présente la résistance thermique effective du forage pour un tube en U lisse simple et double de DN32, ainsi que pour le MuoviELLIPSE DN32.

Résistance thermique effective des forages pour les tubes en U simples et doubles, qu’ils soient lisses ou de type MuoviELLIPSE.
Résistance thermique effective des forages pour les tubes en U simples et doubles, qu’ils soient lisses ou de type MuoviELLIPSE.

Comme vous pouvez le constater, le début du régime de transition survient plus tôt pour le MuoviELLIPSE que pour les tuyaux lisses équivalents. Cela signifie que, lorsque l'on compare le MuoviELLIPSE DN32 au double DN32 lisse, la plage dans laquelle le premier offre de meilleures performances passe de 0,28 à 0,45 l/s (pour la sonde circulaire lisse) à 0,18 à 0,45 l/s. Cela implique qu’il est possible d’obtenir une résistance thermique du forage plus faible avec un seul tube en U (MuoviELLIPSE) à un débit inférieur à celui requis avec un tube lisse standard DN32.

En régime turbulent, le tuyau rond traditionnel présente une résistance thermique au niveau du forage légèrement inférieure à celle du MuoviELLIPSE. Cela s'explique par le fait que, de par sa forme, la distance entre la paroi du forage et le tuyau est plus faible pour une sonde circulaire que pour une sonde elliptique, à distance égale entre le centre du tuyau et celui du forage. Lorsque l'écoulement est turbulent, la résistance du coulis prend davantage d'importance, ce qui explique cet effet.

Sur l'image ci-dessous, cette comparaison est clairement visible : on constate que la sonde elliptique est plus éloignée de la paroi du forage que la sonde circulaire.

Comparaison de la section transversale d'un forage avec une sonde circulaire DN32 (en haut) et une sonde MuoviELLIPSE (en bas).
Comparaison de la section transversale d'un forage avec une sonde circulaire DN32 (en haut) et une sonde MuoviELLIPSE (en bas).

Une autre façon de tirer parti de ce passage plus précoce à la turbulence consiste à utiliser un tuyau de diamètre légèrement supérieur (DN40), comme le montre la figure ci-dessous.

Résistance thermique effective des forages pour des tubes en U simples et doubles de DN32, ainsi que pour un tube lisse simple et un tube MuoviELLIPSE de DN40.
Résistance thermique effective des forages pour des tubes en U simples et doubles de DN32, ainsi que pour un tube lisse simple et un tube MuoviELLIPSE de DN40.

Sur la figure ci-dessus, le MuoviELLIPSE élargit également la plage dans laquelle le tube en U simple surpasse le tube en U double DN32 équivalent, en doublant cette plage, qui passe de 0,3 à 0,45 l/s à 0,2 à 0,45 l/s.

Jusqu’à présent, les performances du MuoviELLIPSE étaient assez similaires à celles du TurboCollector, mais le principal avantage de la conception elliptique réside dans la possibilité d’utiliser un forage de plus petit diamètre, ce qui réduit la résistance du coulis. Grâce à sa forme, l’espace disponible dans le forage est légèrement plus important avec une sonde elliptique, ce qui permet de réduire le diamètre du forage tout en garantissant la pose du tuyau. Ci-dessous, nous reprenons le même graphique que ci-dessus, mais cette fois avec le MuoviELLIPSE DN40 dans un forage d’un diamètre de 100 mm au lieu de 140 mm.

Résistance thermique effective des forages pour des tubes en U simples et doubles de DN32, ainsi que pour un tube lisse simple et un tube MuoviELLIPSE de DN40 (diamètres respectifs de 140 mm et 100 mm).
Résistance thermique effective des forages pour des tubes en U simples et doubles de DN32, ainsi que pour un tube lisse simple et un tube MuoviELLIPSE de DN40 (diamètres respectifs de 140 mm et 100 mm).

Sur le graphique ci-dessus, la courbe du MuoviELLIPSE DN40 est décalée vers le bas en raison du diamètre réduit du forage et de la résistance thermique correspondante plus faible. Grâce à son diamètre de forage plus petit (100 mm contre 140 mm), elle surpasse même la configuration à double DN32 pour tous les débits supérieurs à 0,2 l/s. Cela souligne l'importance d'utiliser des diamètres de forage plus petits.

À titre d'exhaustivité, le jeu nécessaire lors de l'installation du MuoviELLIPSE DN40 dans un forage d'un diamètre de 100 mm est indiqué ci-dessous.

Coupe transversale d'un forage de 100 mm de diamètre équipé d'un MuoviELLIPSE DN40.
Coupe transversale d'un forage de 100 mm de diamètre équipé d'un MuoviELLIPSE DN40.

Perte de charge

Le graphique ci-dessous présente la perte de charge pour les configurations à tube en U simple et double de DN32, ainsi que pour le MuoviELLIPSE DN32.

Perte de charge pour une sonde simple et double DN32 et MuoviELLIPSE.
Perte de charge pour les modèles simple et double DN32 et MuoviELLIPSE.

Il apparaît clairement que la perte de charge commence à augmenter plus tôt pour le MuoviELLIPSE en raison de la turbulence accrue à de faibles débits, ce qui se traduit par une perte de charge plus élevée. De plus, la perte de charge du MuoviELLIPSE est toujours la plus élevée dans ce cas. Cela s'explique par le fait que, comme indiqué plus haut, le diamètre hydraulique du MuoviELLIPSE DN32 PN16 n'est que de 24 mm en raison de sa forme elliptique, alors qu'il est de 26 mm pour les sondes rondes traditionnelles. Cela signifie qu’à débit égal, le MuoviELLIPSE présente une vitesse d’écoulement légèrement supérieure, ce qui entraîne une perte de charge plus élevée.

Le graphique ci-dessous présente une comparaison avec un seul tuyau DN40 (à la fois de type circulaire lisse et MuoviELLIPSE).

Perte de charge pour les modèles simples et doubles DN32, ainsi que pour les modèles simples et MuoviELLIPSE DN40.
Perte de charge pour les modèles simples et doubles DN32, ainsi que pour les modèles simples et MuoviELLIPSE DN40.

Ici, la perte de charge a considérablement diminué, mais elle reste, pour les raisons évoquées plus haut, légèrement supérieure à celle de la sonde simple lisse DN40. En régime laminaire, cependant, ses performances sont identiques à celles du tube en U double DN32.

Conclusion

Dans ce chapitre, nous avons présenté le MuoviELLIPSE de MuoviTech. Il s’agit d’un échangeur de chaleur elliptique doté de la même structure interne à ailettes que le TurboCollector. Par conséquent, le comportement du fluide a également été modélisé à l’aide d’une simulation numérique directe afin de dériver des corrélations pour le coefficient de frottement et le nombre de Nusselt. Pour tenir compte de sa forme irrégulière, la méthode des éléments de frontière a été utilisée pour résoudre numériquement le transfert de chaleur à l’intérieur du forage. Afin d’accélérer ces calculs, un réseau neuronal artificiel a été entraîné à partir des résultats de ces simulations précises.

En ce qui concerne la résistance thermique effective du forage, le MuoviELLIPSE élargit la plage dans laquelle l'écoulement reste turbulent, ou du moins de type transitionnel, ce qui signifie qu’il étend la plage dans laquelle un tube en U simple surpasse un tube en U double. Lorsque l'on utilise un diamètre de forage plus petit, les performances s'améliorent encore davantage, surpassant le double tube en U sur la quasi-totalité de la plage de débits. L'inconvénient de cette turbulence accrue est une augmentation de la perte de charge, due à la fois aux ailettes internes et au diamètre hydraulique plus petit résultant de sa forme elliptique.

Références

    • Katsikadelis, J. T. (2016). La méthode des éléments de frontière pour les ingénieurs et les scientifiques. Academic Press, ISBN : 978-0-12-804493-3
    • Hidman, N. (2026). Évaluation des performances thermohydrauliques des tuyaux collecteurs géothermiques elliptiques à ailettes internes. Disponible en ligne.

Prêt à explorer toutes les possibilités de GHEtool Cloud ?

Vous pouvez essayer GHEtool gratuitement pendant 14 jours, sans carte de crédit.