[{"data":1,"prerenderedAt":175},["ShallowReactive",2],{"$f1bdyhm7z4789t":3},{"slug":4,"publishedAt":5,"category":6,"title":7,"summary":8,"image":9,"imageAlt":14,"readingMinutes":15,"serviceId":16,"keywords":17,"imageCaption":23,"sections":24,"sources":165},"data-center-liquid-cooling-cfd","2026-10-19","methods","CFD du refroidissement liquide des centres de données : l’essentiel","Un modèle CFD de refroidissement liquide couvre la plaque froide, la boucle de baie et la salle. Ce que chacun doit saisir, avec un bilan chiffré et les pièges courants.",{"src":10,"thumbnail":11,"width":12,"height":13},"\u002Fassets\u002Fcfd-cases\u002Fplate\u002Fimages\u002F05.webp","\u002Fbrand\u002Fstudies\u002Fplate.svg",1672,941,"Illustration DEEB-FLOW du champ de température dans un échangeur compact empilé liquide et huile",15,"thermal",[18,19,20,21,22],"CFD refroidissement liquide centre de données","CFD plaque froide","simulation refroidissement direct sur puce","résistance thermique plaque froide","distribution du débit dans une baie","Illustration conceptuelle issue de la bibliothèque de cas DEEB-FLOW. Elle ne constitue pas un résultat de cet article.",[25,32,46,62,68,81,93,100,106,118,124,130,137,150],{"id":26,"title":27,"paragraphs":28},"hottest-gets-least","Le serveur le plus chaud reçoit le moins de liquide",[29,30,31],"Voici un système de refroidissement qui envoie le plus de liquide aux serveurs qui en ont le moins besoin. Personne ne l’a conçu ainsi.","Dans une boucle diphasique pompée, une plaque froide très chargée vaporise davantage de fluide frigorigène, sa perte de charge augmente, et le fluide se reporte vers les plaques qui travaillent moins. Dans un modèle de baie à 34 tiroirs réalisé par Ozguc et ses collègues, sans limiteurs de débit, le tiroir le plus défavorisé atteignait un titre de vapeur prédit de 176 % en sortie, ce qui est la façon polie du modèle de dire que la plaque s’est asséchée.","Les boucles monophasiques se comportent mieux, mais moins bien que ne le laissent croire leurs schémas. Un modèle CFD de refroidissement liquide de centre de données doit montrer de combien, et il doit couvrir trois échelles : l’écoulement et la conduction dans une plaque froide, la part de liquide que reçoit chaque plaque dans la baie, et l’air qui emporte encore une partie de la chaleur dans la salle. Cet article les parcourt une à une, avec les chiffres que l’industrie compare réellement (résistance thermique, débit par kilowatt, perte de charge), un bilan énergétique chiffré et les raccourcis de modélisation qui font paraître un modèle meilleur que le matériel.",{"id":33,"title":34,"paragraphs":35},"numbers-everyone-quotes","Quels chiffres tout le monde cite-t-il ?",[36,37,38,39,42,43,44,45],"Avant tout maillage, apprenez les trois chiffres sur lesquels votre modèle sera jugé, et la façon exacte dont chacun est défini. Deux documents de l’Open Compute Project sont les sources publiques les plus claires.","**Débit par kilowatt.** Les lignes directrices de l’OCP pour les systèmes d’accélérateurs refroidis par liquide recommandent un échauffement du liquide de 7,5 à 12 °C, ce qui correspond selon elles à 1,25 à 2,0 litres par minute et par kilowatt pour un liquide à 25 % de propylène glycol, appelé PG25. Leur objectif de conception habituel est un échauffement de 10 °C, soit 1,5 L\u002Fmin par kW. Un échauffement plus fort améliore l’efficacité du site et réduit la marge de refroidissement ; un échauffement plus faible fait l’inverse. Les mêmes lignes directrices prévoient pour chaque plaque froide d’accélérateur un débit de 1 à 2 L\u002Fmin.","**Résistance thermique.** La spécification d’essai des plaques froides de l’OCP, en vigueur depuis le 3 septembre 2026, la définit ainsi :",{"formula":40},[41],"R = (T_s − T_l) \u002F P","Elle se mesure avec la plaque froide montée sur un processeur ou un véhicule d’essai thermique, avec le matériau d’interface thermique recommandé. T_s est la température du boîtier, T_l la température du liquide en entrée et P la puissance appliquée. Le liquide d’essai est en général à 45 ± 2 °C, et l’essai est répété à plusieurs débits pour obtenir une courbe.","Les lignes directrices pour accélérateurs utilisent la même forme boîtier-entrée et donnent pour des plaques froides seules de 0,014 à 0,019 °C\u002FW pour 1,5 à 3,5 L\u002Fmin par plaque. Ces valeurs ont été obtenues en supposant une distribution uniforme du débit, et le document ajoute qu’un liquide PG25 peut coûter jusqu’à 15 % de performance par rapport à l’eau déminéralisée.","**Perte de charge.** La spécification la mesure dans une boucle d’essai comprenant tubes et raccords rapides, et soustrait un essai de tare sans plaque froide. Au moins six débits entre 25 % et 150 % du débit nominal, dont 100 %, sont présentés sous forme d’un ajustement à deux coefficients. Un essai de débit distinct monte jusqu’à un peu moins de 1,5 m\u002Fs pour vérifier l’érosion et les fuites.","Remarquez le détail de ces définitions : quelle température, mesurée où, à quelle température d’entrée, avec quel liquide, sur quelle plage de débit. Votre modèle CFD doit donner la même grandeur de la même façon. Sinon, la comparaison porte sur deux nombres qui partagent seulement un symbole.",{"id":47,"title":48,"paragraphs":49},"energy-balance","Un bilan énergétique chiffré",[50,51,56,57,58,61],"**Exemple chiffré.** Partons du couple de valeurs donné par le document lui-même : 1,5 L\u002Fmin par kW donne un échauffement de 10 K. Ce couple fixe la capacité thermique volumique du liquide, et avec elle la chaleur que chaque litre par minute peut emporter.",{"formula":52},[53,54,55],"V̇ = 1.5 L\u002Fmin = 2.5 × 10⁻⁵ m³\u002Fs","ρc_p = P \u002F (V̇ ΔT) = 1000 \u002F (2.5 × 10⁻⁵ × 10) = 4.0 × 10⁶ J\u002F(m³·K)","ρc_p V̇ = 100 W\u002FK","Chaque litre par minute emporte donc environ 67 W par kelvin d’échauffement. Servez-vous-en maintenant pour contrôler votre modèle. Si votre calcul CFD donne une sortie 12 K au-dessus de l’entrée à 1 kW et 1,5 L\u002Fmin, le refroidissement va bien et c’est le bilan thermique qui est faux. Quelque part de la chaleur apparaît, ou le débit n’est pas celui que vous croyez. Corrigez cela avant de regarder le moindre champ.","Avec une entrée à 40 °C et une résistance de 0,017 °C\u002FW, dans la plage publiée, la température du boîtier vaut :",{"formula":59},[60],"T_case = 40 + 0.017 × 1000 = 57 °C","Utilisez la fiche technique du liquide que vous modélisez réellement. La concentration en glycol et la température modifient ces chiffres, et tout ce qui suit évolue avec eux.",{"id":63,"title":64,"paragraphs":65},"what-r-includes","Ce que la résistance thermique contient discrètement",[66,67],"Une partie de ces 0,017 °C\u002FW revient au liquide. Comme R est rapportée à la température d’entrée, elle inclut l’échauffement du liquide pendant la traversée de la plaque. Dans l’exemple ci-dessus, le liquide gagne environ 10 K. Si la zone chaude de la puce se trouve à peu près à mi-parcours de l’écoulement, environ 5 K des 17 K d’écart entre boîtier et entrée viennent de l’échauffement propre du liquide, près d’un tiers. C’est une estimation grossière avec une hypothèse simple, et elle explique ce qui en intrigue beaucoup : à plus faible débit, une partie de la hausse de R vient simplement de l’échauffement plus fort du liquide, la convection restant la même.","Une définition fondée sur l’écart moyen logarithmique de température sépare les deux effets. Quelle que soit la définition que vous préférez, présentez celle qu’utilise l’essai.",{"id":69,"title":70,"paragraphs":71},"flow-is-laminar","Dans une plaque froide : l’écoulement est laminaire",[72,73,79,80],"**Exemple chiffré, avec une géométrie illustrative.** Prenons 1,5 L\u002Fmin répartis sur 100 canaux parallèles de 0,3 mm de large et 2 mm de haut, et les propriétés de l’eau à 40 °C tirées des tables usuelles, environ 992 kg\u002Fm³ et 0,65 mPa·s.",{"formula":74},[75,76,77,78],"A = 100 × 0.3 mm × 2 mm = 60 mm²","V = 2.5 × 10⁻⁵ \u002F (60 × 10⁻⁶) ≈ 0.42 m\u002Fs","D_h = 4A \u002F P_wet = 4 × 0.6 \u002F 4.6 ≈ 0.52 mm","Re = ρ V D_h \u002F μ = 992 × 0.42 × 0.00052 \u002F 0.00065 ≈ 330","Un mélange glycolé est plus visqueux que l’eau, le nombre de Reynolds réel est donc encore plus faible. Quelques centaines, c’est laminaire selon tous les critères, et beaucoup de conceptions à microcanaux et à mésocanaux tombent dans le même régime.","Désactivez donc le modèle de turbulence dans les canaux, ou sachez au moins ce qu’il y fait. Un modèle de turbulence alimenté par une intensité turbulente d’entrée par défaut peut ajouter une viscosité tourbillonnaire que l’écoulement n’a pas, ce qui relève le coefficient d’échange et rend la plaque meilleure qu’elle n’est. Vérifiez le rapport de viscosité turbulente dans les canaux. S’il n’est pas proche de zéro, vous avez votre réponse.",{"id":82,"title":83,"paragraphs":84},"flow-is-developing","Dans une plaque froide : l’écoulement est encore en développement",[85,86,91,92],"Les longueurs d’établissement laminaires valent environ 0,05 Re D_h pour le profil de vitesse et 0,05 Re Pr D_h pour le profil de température. Le produit Re Pr est le nombre de Péclet, V D_h ρc_p \u002F k, qui ne dépend pas de la viscosité ; il ne demande donc que la capacité thermique volumique du bilan énergétique et une conductivité d’environ 0,6 W\u002F(m·K).",{"formula":87},[88,89,90],"L_h ≈ 0.05 × Re × D_h ≈ 0.05 × 330 × 0.52 mm ≈ 9 mm","Pe = V D_h ρc_p \u002F k = 0.42 × 0.00052 × 4.0 × 10⁶ \u002F 0.6 ≈ 1460","L_t ≈ 0.05 × Pe × D_h ≈ 0.05 × 1460 × 0.52 mm ≈ 38 mm","Le profil de vitesse s’établit en moins d’un centimètre, et plus vite encore avec le glycol, plus visqueux. Le profil de température est encore en développement sur une longueur comparable à toute la plaque froide, et un mélange glycolé, qui conduit moins bien la chaleur que l’eau, l’allonge. Un nombre de Nusselt laminaire établi sous-estime donc l’échange sur la plus grande partie du canal. Un calcul à la main et un modèle CFD seront en désaccord, et c’est le CFD qu’il faut croire si le maillage résout la couche limite thermique. Les lignes directrices de l’OCP font la même remarque pratique : les concepteurs utilisent la simulation pour étudier le profil de température le long des ailettes et empêcher l’écoulement de s’établir complètement dans les canaux.","Maillez en conséquence : plusieurs mailles à travers chaque canal et chaque ailette, et une étude de maillage sur les deux grandeurs qui comptent, R et la perte de charge. La méthode de l’indice de convergence de grille de notre article sur l’indépendance au maillage s’applique ici comme partout.",{"id":94,"title":95,"paragraphs":96},"heat-source","La source de chaleur n’est pas un flux plat et uniforme",[97,98,99],"C’est le raccourci que je rencontre le plus souvent dans les modèles de transfert thermique conjugué en électronique, et il flatte la plaque froide à chaque fois.","Un flux uniforme sur la base de la plaque froide donne un champ de température lisse et symétrique et un pic agréablement bas. Un vrai accélérateur a des zones chaudes, et la chaleur traverse un empilement avant d’atteindre la plaque. Les lignes directrices de l’OCP énumèrent les éléments entre la jonction et le liquide : l’empilement de puces, la première interface thermique et le capot, la seconde interface thermique, la plaque froide elle-même et les échanges latéraux entre voisins. Elles ajoutent qu’une grande partie de l’échauffement peut naître à la surface du boîtier et en son sein, et que dans bien des cas la limite de refroidissement dépend davantage du boîtier que de la solution de refroidissement.","Modélisez donc l’empilement : le capot et les couches d’interface comme solides ou comme résistances thermiques minces avec des valeurs issues de données, une carte de puissance du fournisseur de la puce si vous pouvez l’obtenir, et la géométrie du véhicule d’essai thermique si vous comparez avec un essai. Si vous n’avez pas les résistances d’interface ou de contact, dites qu’il s’agit d’hypothèses et faites un calcul de sensibilité. Ce sont les deux données les plus susceptibles de déplacer la réponse.",{"id":101,"title":102,"paragraphs":103},"across-the-rack","À travers la baie : qui reçoit le liquide ?",[104,105],"Des plaques froides en parallèle se partagent le débit selon leur résistance hydraulique, et chaque collecteur, chaque flexible et chaque raccord rapide a son mot à dire. Si la répartition est inégale, la plaque qui reçoit moins de débit chauffe davantage, et aucun modèle à une seule plaque ne vous le dira.","La même physique explique pourquoi les échangeurs de chaleur n’atteignent pas leur puissance, ce que nous avons traité dans l’article sur la mauvaise distribution de l’écoulement dans les échangeurs. Les baies refroidies par liquide sont, du point de vue thermique, de petits échangeurs aux tubes très coûteux.",{"id":107,"title":108,"paragraphs":109},"series-or-parallel","En série ou en parallèle ?",[110,111,112,117],"Les lignes directrices de l’OCP pour accélérateurs comparent deux dispositions pour huit plaques froides : quatre branches parallèles de deux plaques en série (4 × 2), et les huit en parallèle (8 × 1). D’après un modèle de réseau calibré sur des essais, à 10 L\u002Fmin de propylène glycol à 25 %, la boucle 4 × 2 demandait environ 10 psi et la boucle 8 × 1 environ 7 psi. La disposition 8 × 1 maintenait toutes les puces à moins de 1 °C les unes des autres, contre 7 °C pour la 4 × 2.","**Exemple chiffré.** Pourquoi 7 °C ? Dans la disposition 4 × 2 à 10 L\u002Fmin, chaque branche reçoit 2,5 L\u002Fmin. Supposons que chaque plaque dissipe 1 kW :",{"formula":113},[114,115,116],"V̇ = 10 \u002F 4 = 2.5 L\u002Fmin = 4.17 × 10⁻⁵ m³\u002Fs","ρc_p V̇ = 4.0 × 10⁶ × 4.17 × 10⁻⁵ ≈ 167 W\u002FK","ΔT = 1000 \u002F 167 ≈ 6.0 K","La puce aval reçoit un liquide environ 6 K plus chaud que la puce amont avant tout autre effet, et l’inégalité de débit entre les branches fait le reste. La disposition entièrement parallèle supprime ce préchauffage et demande au collecteur de répartir le débit de façon égale entre huit plaques, chacune à plus faible débit, et c’est précisément là que les plaques froides doivent encore bien fonctionner.",{"id":119,"title":120,"paragraphs":121},"flow-network","Modéliser la boucle comme un réseau",[122,123],"L’étude de l’OCP montre la méthode que je recommanderais pour la plupart des boucles de baie. Calculez en CFD une plaque froide sur toute sa plage de débit pour obtenir les courbes de perte de charge et de résistance thermique en fonction du débit. Injectez ces courbes dans un réseau hydraulique unidimensionnel de toute la boucle : collecteurs, flexibles, raccords et plaques. Calibrez le réseau sur un essai de la boucle, puis servez-vous-en pour explorer dispositions et points de fonctionnement.","On pourrait mailler toute une boucle de baie en résolvant chaque ailette. C’est possible, c’est impressionnant, et ce n’est presque jamais la bonne première étape. Le réseau répond à la question de la distribution en quelques secondes, et l’effort CFD va là où la physique est difficile.",{"id":125,"title":126,"paragraphs":127},"two-phase","Le piège du diphasique",[128,129],"Revenons au début. Dans les boucles diphasiques, la perte de charge d’une plaque dépend de sa charge thermique, si bien que la distribution change avec la charge de calcul. Ozguc et ses collègues ont utilisé un modèle des collecteurs avec conservation de la masse, de la quantité de mouvement et de l’énergie, une hypothèse d’écoulement diphasique homogène et des courbes de perte de charge des tiroirs ajustées sur 104 points d’essai, la plupart prédits à 25 % près. Leur remède est un orifice en amont de chaque plaque froide, dimensionné pour que la perte de charge du limiteur domine. Bien dimensionnés, les orifices maintenaient chaque tiroir à un titre de sortie de 85 % ou moins. Sous-dimensionnés, ils ne corrigeaient pas la mauvaise distribution ; surdimensionnés, ils gaspillaient de la pression.","La leçon pour les modélisateurs est qu’une baie diphasique est un problème de réseau hydraulique avec des données empiriques bien avant d’être un problème de CFD. L’ébullition dans les canaux des plaques froides reste un sujet de recherche pour la CFD, et c’est le réseau, alimenté par des mesures, qui a permis de résoudre ce cas.",{"id":131,"title":132,"paragraphs":133},"the-room","Dans la salle : l’air compte-t-il encore ?",[134,135,136],"Oui. Le refroidissement direct sur puce couvre en général les composants les plus chauds, le reste du serveur reste refroidi par air, et une partie de la chaleur de chaque baie arrive donc toujours dans la salle. Supposons, pour l’exemple, qu’un cinquième d’une baie de 100 kW reste refroidi par air. Cela fait 20 kW de charge à air par baie, soit ce que portait à elle seule une baie dense refroidie par air.","Les lignes directrices de l’OCP notent aussi que les températures de liquide les plus élevées, de 40 à 45 °C au départ secondaire, peuvent être desservies par des équipements hybrides liquide-air comme des échangeurs en porte arrière ou latéraux. Ceux-ci renvoient volontairement la chaleur dans l’air de la salle.","Le modèle de salle survit donc. Les baies deviennent des conditions aux limites avec un débit d’air et une charge thermique, les portes arrière et échangeurs latéraux deviennent des échangeurs, et la question reste la question familière de la recirculation et des températures d’entrée. Si vous avez déjà modélisé l’air d’une salle, la méthode est la même, avec des charges à air par baie plus faibles que ne l’indique la plaque signalétique et moins de marge qu’on ne le voudrait.",{"id":138,"title":139,"paragraphs":140},"common-mistakes","Erreurs fréquentes, telles qu’elles apparaissent à l’écran",[141],{"list":142},[143,144,145,146,147,148,149],"**Une température de base lisse et symétrique** sous un flux uniforme. Les vraies puces ont des zones chaudes ; le pic du modèle est optimiste.","**Une viscosité turbulente non nulle dans des canaux laminaires.** La plaque paraît meilleure qu’elle n’est.","**Toutes les puces à la même température** dans un modèle de baie. Quelqu’un a supposé un débit uniforme. Même les bons rapports d’essai énoncent cette hypothèse, alors lisez-la.","**Une résistance de fiche technique utilisée avec un autre liquide ou un autre débit.** L’eau déminéralisée et le PG25 diffèrent jusqu’à 15 % dans les données de l’OCP, et R dépend du débit.","**Une température de sortie incompatible avec P \u002F (ρc_p V̇).** Un problème de bilan thermique. Arrêtez-vous et trouvez-le.","**Une résistance CFD calculée à partir de la température moyenne de la base,** comparée à un essai qui utilise un thermocouple au centre du boîtier. Deux nombres différents pour un seul symbole.","**Pas de modèle de salle,** parce que la baie est « refroidie par liquide ».",{"id":151,"title":152,"paragraphs":153},"what-i-would-do","Ce que je ferais",[154,155,162,163,164],"Pour un système monophasique de refroidissement direct sur puce :",{"list":156,"ordered":161},[157,158,159,160],"**Une plaque froide en détail.** Transfert thermique conjugué, écoulement laminaire dans les canaux, empilement du boîtier et carte de puissance réaliste, calculés à cinq ou six débits sur la plage de l’essai OCP, avec une étude de maillage sur R et la perte de charge.","**Une comparaison avec un essai de plaque** s’il en existe un, avec les mêmes définitions, le même liquide et la même température d’entrée.","**Un réseau hydraulique de la boucle,** fondé sur les courbes CFD et calibré sur un essai de boucle quand c’est possible.","**Un modèle de salle** où les baies sont des conditions aux limites de chaleur et de débit d’air, avec la charge résiduelle à air et les éventuelles portes arrière.",true,"Cela change pour les systèmes diphasiques, où je partirais du réseau hydraulique et de données mesurées sur les plaques, en gardant la CFD pour les questions locales. Cela change encore pour le refroidissement par immersion, où la convection naturelle dans une cuve pose le problème et où le modèle est assez différent pour mériter son propre article.","Commencez dans tous les cas par la serviette en papier : puissance, débit, échauffement et chaîne de résistances. Si les chiffres ne tombent pas juste sur la serviette, aucun maillage ne les fera tomber juste à l’écran.","Si vous concevez ou vérifiez un système refroidi par liquide et souhaitez que la plaque froide ou la boucle soit modélisée avec les contrôles ci-dessus, la page CFD des échangeurs de chaleur décrit notre approche de ce type de travail en transfert thermique conjugué, et les modes d’intervention précisent ce que comprend une étude.",[166,169,172],{"href":167,"label":168},"https:\u002F\u002Fwww.opencompute.org\u002Fwp-content\u002Fuploads\u002F2026\u002F10\u002FOCP-Cold-Plate-Test-Specification-V1.0.0.pdf","Open Compute Project (2026): Single-Phase Direct-to-Chip Cold Plate Test Specification 1.0.0",{"href":170,"label":171},"https:\u002F\u002Fwww.opencompute.org\u002Fdocuments\u002Foai-system-liquid-cooling-guidelines-in-ocp-template-mar-3-2023-update-pdf","Open Compute Project (2023 update): OAI system liquid cooling guidelines",{"href":173,"label":174},"https:\u002F\u002Fwww.electronics-cooling.com\u002F2025\u002F02\u002Fmitigating-flow-maldistribution-in-data-center-two-phase-cooling-systems-with-flow-restrictors\u002F","Ozguc, Wang, Narayanan and Bonner (2025), Electronics Cooling: mitigating flow maldistribution in data center two-phase cooling systems with flow restrictors",1792438042561]