Compte tenu du nombre et de l'importance des applications faisant appel à des échangeurs de chaleur dans les sociétés modernes, que ce soit pour la production d'énergie, la chimie, le secteur pétrolier, la réfrigération, le conditionnement de l'air, un besoin croissant s'est manifesté au cours des dernières décennies de connaître de manière précise les valeurs des coefficients d'échange thermique et de pertes de charge, en particulier lors de l'ébullition ou de la condensation des fluides.

Des travaux de recherche ont ainsi été entrepris dans de nombreux pays sur ce sujet, qui ont donné et continuent de donner lieu à de très nombreuses publications sur les corrélations permettant de les déterminer.

Une des grandes difficultés aujourd'hui est de pouvoir choisir les corrélations adaptées à un problème donné parmi l'ensemble de celles qui sont publiées. Certains articles de synthèse sont disponibles, mais ils ne sont pas toujours fiables. Des erreurs existent fréquemment dans les formules qui sont données, les définitions des variables prises en compte peuvent différer d'une publication à l'autre...

Dans cette page, nous passons en revue les corrélations qui ont été ajoutées en 2021 aux classes externes de Thermoptim, en renvoyant le lecteur intéressé aux publications dont elles sont issues, ainsi qu'aux classes externes dans lesquelles elles sont implémentées.

Les versions de Thermoptim permettant d'utiliser ces nouvelles fonctionnalités sont 2.72, 2.82 et Demo 2.82. Les versions précédentes ne sont pas compatibles.

Ces corrélations portent sur les cas suivants :

• écoulements simples

• ébullition saturée

• ébullition nucléée

• condensation

• extérieur des échangeurs à ailettes

Précisons aussi que ces corrélations sont généralement implémentées dans les classes externes de type FlowConfig, mais qu'elles le sont quelquefois dans des TechnoHX.

Cette page vous fournit quelques précisions sur la manière dont ces classes sont structurées, en complément des explications détaillées fournies dans les tomes 3 et 4 du manuel de référence de Thermoptim.

Si les implémentations de ces corrélations ne répondent pas à vos besoins, vous pouvez les compléter par d'autres classes bâties de manière analogue.

Il importe en effet de bien avoir à l'esprit que de très nombreuses corrélations ont été publiées au cours des dernières années.

Elles dépendent de nombreux paramètres, comme le fluide mis en jeu, le diamètre hydraulique, l'orientation et la géométrie des canalisations...

Nous en avons codé un certain nombre, mais il est impossible de toutes les implémenter et vous serez peut-être amenés à en introduire d'autres correspondant mieux à vos besoins.

Vous pourrez le faire en adaptant les classes existantes, ce qui ne représente généralement pas un travail très important.

Précisons que les classes actuelles ne sont pas optimisées afin de pouvoir disposer des résultats des différentes corrélations dans le fichier de résultat output.txt.

Si vous le souhaitez, vous pouvez facilement les personnaliser en supprimant les calculs qui vous sont inutiles.

Les explications sur le développement des classes externes sont données dans cette page du portail.

Les corrélations donnant le nombre de Nusselt pour les écoulements simples sont classiques : MacAdams à l'intérieur des tubes, dans la classe IntTubeConfig.java, et Colburn dans la classe ExtTubeConfig.java.

La classe TechnoEvaporator utilisait la corrélation de Gungor-Winterton calculée de manière moyennée, ce qui était criticable. Elle est maintenant implémentée en discrétisant la zone d'évaporation en 100 éléments, ce qui est beaucoup plus précis.

Cinq nouvelles corrélations d'estimation des coefficients d'échange pendant l'ébullition y ont de plus été implémentées, celles de (Shah, 1982), (Borishanskiy, 1971), (Kim et Mudawar, 2013), (Kandlikar, 2017) et (Saitoh, 2007).

Pour les pertes de charge, la corrélation de (Lockhart-Martinelli, 1949) est l'une des plus utilisée, même si elle n'est pas parmi les plus précises selon les auteurs : nombreuses même sont les publications qui s'accordent pour considérer qu'elle les surestime significativement. Il semblerait que celles de (Müller-Steinhagen & Heck, 1986) et de (Sun & Mishima, 2009) ou celle de (Friedel,1979) soient beaucoup plus précises. Les trois ont été implémentées dans la classe TechnoSteamGenerator.

Le calcul de l'ébullition nucléée est beaucoup plus délicat que celui de l'ébullition saturée. C'est un domaine d'investigation relativement récent et de nombreuses publications existent, sans que des synthèses vraiment pertinentes soient disponibles.

Lorsque la température de paroi Tw croît, trois seuils apparaissent :

• lorsque Tw atteint Tsat, des bulles peuvent commencer à apparaître,

• mais elles ne le font qu'avec un certain retard, lorsque Tw = Tw,onb (Onset of Nucleate Boiling). Un flux complémentaire du flux purement convectif q_fc apparaît alors. Il est égal à q_nb pour nucleate boiling.

• l'ébullition devient complète lorsque Tw > Tw, fdb (fully developped boiling), ce qui correspond à un titre x supérieur ou égal à 0.

Dans la classe TechnoSteamGenerator, l'ONB est détecté par une relation type (Thom, 1965), et le FNB par une relation type (Bowring, 1962).

On considère que le coefficient d'échange h_fnb est égal à celui calculé pour l'ébullition saturée, en fonction de la corrélation choisie.

Pour la transition entre l'ONB et le FDB, on fait une rampe entre le hfl local et hlvf atteint lors du FDB. Des méthodes plus sophistiquées sont proposées dans la littérature, mais elles sont complexes à implémenter et leur robustesse n'est pas garantie.

L'échangeur est découpé en 100 parties et les différents équilibres thermiques sont calculés. La température de paroi est déterminée pour chaque intervalle, avec l'hypothèse que le coefficient d'échange côté chaud reste constant.

Différentes corrélations de détection de l'ONB ont été implémentées, et celle qui conduit à la valeur la plus élevée est retenue. Il suffit de la remplacer par une autre si on le souhaite.

De la même manière, différentes corrélations de détection de l'OFB (ou FDB ou NVG) ont été implémentées, et celle qui conduit à la valeur la plus élevée est retenue. h_fdb est estimé égal à celui de l'ébullition saturée.

Le nombre de Nusselt pour la condensation laminaire à l'extérieur d'un tube horizontal peut être estimée par la théorie de Nusselt avec la modification de gravité de Dhir et Lienhard (1971). C'est cette corrélation qui est implémentée dans la classe ExtHorCondConfigDhir.java.

Une variante pour les tubes verticaux est implémentée dans la classe ExtVertCondConfig.java.

La corrélation présentée par (Lévy, 1990) pour la condensation à l'extérieur des tubes horizontaux utilisée dans les condenseurs de centrale à vapeur est implémentée dans la classe ExtHorCondConfig.java.

Pour l'intérieur des tubes horizontaux, la corrélation de (Shah, 1979) continue d'être souvent utilisée, même si l'auteur en a proposé une plus complexe par la suite. Nous avons implémenté une variante modifiée par (Bivens, 1994) dans la classe CondConfig.java.

A la différence des corrélations précédentes qui fournissent généralement le nombre de Nüsselt Nu, celles relatives au côté air des échangeurs à ailettes sont généralement rapportées au facteur de Chilton–Colburn j.

Deux corrélations peuvent être utilisées si on dispose de données expérimentales, celles de (Morisot, 2002) et celle de (Martin et Ginielski, 2000).

Trois corrélations prédictives supposent une connaissance plus ou moins détaillée de la géométrie de l'échangeur à ailettes. Il s'agit de celles de (Kim et Jacobi, 2000) de (Manglik et Bergles, 1995), et de (Wang, Chi et Chang, 2000). Elles sont implémentées dans les classes FinnedAirCoilConfigKJ.java, FinnedAirCoilConfigMB.java et FinnedAirCoilConfigWCC.java.

Elles utilisent l'écran de paramétrage de corrélations pour définir les grandeurs géométriques dont elles ont besoin, relatives à :

• l'espacement des ailettes Fp,

• le diamètre extérieur des tubes,

• le nombre N de nappes,

• la distance transversale Pt entre elles,

• leur écartement Pl.

Selon les cas, certains de ces paramètres ne sont pas pris en compte.

Un exemple de paramétrage expliqué est donné, en complément de l'échangeur simple utilisé pour illustrer le dimensionnement technologique et le non-nominal.

W.W. Akers, H. Deans, O. Crosser, Condensing heat transfer within horizontal tubes, Chem. Eng. Prog. Symp. 54 (1955) 89–90.

Bivens, D. B. and Yokozeki, A., "Heat Transfer Coefficients and Transport Properties for Alternative Refrigerants" (1994). InternationalRefrigeration and Air Conditioning Conference. Paper 263.

Borishanskiy, B.M., Andreevskij, A.A., Fromzel, V.N., Fokin, B.S., Cistgakov, V.A., Danilowa, G.N. and Bikov, G.S., (1971). Heat transfer during two-phase flows (in Russian). Teploenergetika 11, pp. 68-69.

Bowring R.W. , Physical model based on bubble detachment and calculation of steam voidage in subcooled region of a heated channel, Institute for Atomenergi, Halden, Norway, Report No. HPR-10, 1962.

Dhir V. K., Lienhard J. H., Laminar Film Condensation on Plane and Axisymmetric Bodies in Nonuniform Gravity, Journal of Heat Transfer-transactions of The Asme, 39, 97, 1971

Friedel, L., (1979). Improved friction pressure drop correlation for horizontal and vertical twophase pipe flow. European Two-Phase Flow Group Meeting, paper E2, Ispra, Italy.

Ghione A., Assessment and improvements of thermal-hydraulic correlations and methods for the analysis of the Jules Horowitz Reactor, Thesis for the Degree Of Doctor Of Philosophy, Chalmers University of Technology, Goteborg, Sweden, 2017

Gungor, K. E.; Winterton, R. H. S.: A general correlation for flow boiling in tubes and in annuli. Int. J. Heat Mass Transfer 29 (1986) 351–358

Kandlikar, Ch. 8. Boiling, Multiphase Flow Handbook, Second Edition by Crowe, Clayton T. Michaelides, Efstathios Schwarzkopf, John D., 2017

Kim G. J., Jacobi A. M. Condensate Accumulation Effects on the Air-Side Thermal Performance of Slit-Fin Surfaces, ACRC CR-26, University of Illinois at Urbana-Champaign I, 2000.

Kim, S.M. and Mudawar, I., 2013, “Universal Approach to Predicting Saturated Flow Boiling Heat Transfer in Mini/Micro-Channels Part II. Two-Phase Heat Transfer Coefficient,”International Journal of Heat and Mass Transfer, Vol. 64, pp. 1239-1256

Lévy W., Condenseurs par surface dans les centrales thermiques, Article B 1 540, Techniques de l'Ingénieur, 1990

Lockhart, R.W. and Martinelli, R.C.,Proposed correlation of data for isothermal two-phase, two-component flow in pipes. Chem. Eng. Progr. 45 (1949), pp. 39-48.

Manglik, R. M., Gergles, A.E. Heat transfer and pressure drop correlations for the rectangular offset strip fin compac heat exchanger, Experimental Thermal and Fluid Science 1995; 10:171-180, Elsevier

Martin H., Gnielinski V., Calculation of heat transfer from pressure drop in tube bundles and, 3rd European Thermal Sciences Conference 2000, E.W.P. Hahne, W. Heidemann and K. Spindler (Editors),

Morisot, O., Marchio, D. (2002). Simplified Model for the Operation of Chilled Water Cooling Coils Under Nonnominal Conditions. HVAC&R Research 8(2): 135-158.

Müller-Steinhagen, H. and Heck, K., (1986). A simple friction pressure drop correlation for twophase flow in pipes. Chem. Eng. Process. 20, pp. 297-308.

Ohrby, Fredrik, Numerical modeling of subcooled nucleate engine cooling systems, Master's thesis in the Applied Mechanics programme Chalmers University of Technology, Goteborg, Sweden 2014

R. Rabiee, M. Désilets, M. Proulx, M. Ariana, M. Julien, Determination of condensation heat transfer inside a horizontal smooth tube, International Journal of Heat and Mass Transfer 124 (2018) (816-828)

Saitoh, S., Daiguji, H. and Hihara, E., (2007). Correlation for boiling heat transfer of R134a in horizontal tubes including effect of tube diameter. Int. J. Heat Mass Transfer 50, pp. 5215-5225.

Shah, M.M. 1982. Chart correlation for saturated boiling heat transfer: equations and further study. ASHRAE Trans. 88(1):185–196.

Shah, M. A General Correlation for Heat Transfer During Flow Condensation Inside Pipes, Journal of Heat and Mass Transfer 22, 1979, 547-56

Sun, L. and Mishima, K., (2009). Evaluation analysis of prediction methods for two-phase flow pressure drop in mini-channels. Int. J. Multiphase Flow 35, pp. 47-54

Thom J. R. S. , W. M. Walker, T. A. Fallon, and G. F. S. Reising, Boiling in subcooled water during annuli. In Proceedings of Symposium on Boiling Heat Transfer in Steam Generation Units and Heat Exchangers, 1965

Thome J. R and Cioncolini, A. Encyclopedia of Two-Phase Heat Transfer and Flow I, Fundamentals and Methods, Volume 3: Flow Boiling in Macro and Microchannels, 2015

Wang, Chi-Chuan, Chi Kuan-Yu, Chang Chun-Jung Heat transfer and friction characteristics of plain fin-and-tube heat exchangers, part II: Correlation August 2000, International Journal of Heat and Mass Transfer 43(15):2693-2700 DOI: 10.1016/S0017-9310(99)00333-6

Wang, J, Y. Huang, Y. Wang, Visualized study on specific points on demand curves and flow patterns in a single-side heated narrow rectangular channel, Int. J. Heat Fluid Flow 32 (2011) 982–992.