# Simulation des propriétés thermodynamiques du fluide R134a

Djamila Diaf et Ali Malek

Division Energie Solaire Photovoltaïque
Centre de Développement des Energies Renouvelables
B.P. 62, Route de l'Observatoire, Bouzaréah, Alger

(reçu le 10 Juin 2006 – accepté le 25 Mars 2007)

**Résumé** - Nous présentons, dans cet article, des relations simples de forme polynomiale pour déterminer les propriétés thermodynamiques du fluide frigorigène R134a par simulation sans avoir recours aux abaques. Les relations proposées ont été déterminées par des méthodes numériques et sont applicables dans l'intervalle de température de saturation de – 40°C à 90°C pour l'état de saturation et pour une température réelle de + 120°C dans le domaine de la phase surchauffée. Cette étude nous est indispensable pour connaître les performances du système frigorifique à compression d'un conservateur médical alimenté par énergie solaire photovoltaïque.

**Abstract** - We present in this article, the simples relations of polynomial form to determine the thermodynamic properties of the refrigerant fluid R134a by simulation without having recourse to the abacuses. The equations suggested were determined by numerical methods and are applicable in the interval of saturation temperature of -40°C to 90°C for the saturation state and a real temperature of + 120°C for the field of the overheated phase. This study is necessary for us to know the performances of a refrigerating system with compression (medical conservator) supplied with photovoltaic solar energy.

**Mots clés:** Fluides frigorigènes - Réfrigération - Conservateur médical - Simulation.

## 1. INTRODUCTION

Pour la plupart des régions du monde, la mise à disposition de sources de froid est une base considérable de progrès. Le besoin du froid pour la conservation des fruits et légumes, et pour la conservation des produits pharmaceutiques, surtout des vaccins, est une nécessité de première importance.

Pour subvenir à la demande en production de froid, de nombreux systèmes frigorifiques ont été mis en pratique, selon divers cycles thermodynamiques. Parmi ces cycles, le cycle à compression de vapeur, le cycle à absorption et le cycle à adsorption sont les plus utilisés dans les systèmes de réfrigération et de la climatisation.

Les systèmes de réfrigération les plus usités sont des systèmes à compression de vapeur et fonctionnant avec des fluides frigorigènes, R12, R22, R502 etc. Ils nécessitent de l'énergie électrique pour leur fonctionnement. Dans les régions où l'énergie électrique est disponible, ce système de réfrigération est habituellement adéquat pour satisfaire la plupart des exigences en réfrigération. Cependant, plusieurs régions ne disposent pas de l'énergie électrique, principalement des villages ruraux et des sites isolés dans les pays en voie de développement.

La réfrigération dans ces régions, par conséquent, exige une solution alternative à ces systèmes de réfrigération. L'utilisation de l'énergie solaire comme source d'énergie électrique pour alimenter des systèmes frigorifiques est intéressante, puisqu'elle est inépuisable et est abondante dans les régions, où le besoin en froid est énorme.

Dans ce contexte, plusieurs systèmes ont été développés, particulièrement les systèmes frigorifiques photovoltaïques restent les plus appropriés.

Pour contribuer au développement de systèmes frigorifiques photovoltaïques performants, un nouveau type de réfrigérateur destiné à la conservation des vaccins a été proposé, étudié et
développé. Et c'est dans le but de simuler les performances de ce réfrigérateur solaire, qu'une étude des propriétés thermodynamiques du fluide frigorigène utilisé s'est avérée nécessaire [1].

Le fluide le plus utilisé dans la réfrigération à compression de vapeur était le R12. Or du fait de l'interdiction des fluides chlorofluorocarbones 'CFC', dû à leur 'ODP -Ozone Depletion Potential' et 'GWP -Global Warming Potential' élevés [2] et de la grande difficulté à trouver des substituts qui aient les mêmes propriétés thermodynamiques que les CFC, la plupart des industriels se sont convertis à utiliser ce fluide, qui est le R134a, dans les installations de moyenne et haute température (réfrigération domestique et industrielle).

Afin de simuler et de prédire les performances du réfrigérateur médical proposé fonctionnant à base du fluide frigorigène R134a, une étude synthétique sur les propriétés thermo physiques de ce fluide est indispensable.

Dans la pratique, les propriétés thermodynamiques et physico-chimiques des fluides se présentent sous forme d'abaques ou de tableaux, où sont donnés l'ensemble des paramètres indicatifs, comme la pression P, la température T, le volume spécifique v, l'enthalpie H et l'entropie S. L'expérience montre que l'utilisation de ces tableaux est très peu pratique. Pour rendre exploitables ces données, il est proposé dans cet article, de les approcher par des formes polynomiales simples établies par le biais de méthodes numériques.

Les états étudiés de ce fluide frigorigène sont les suivants: état du fluide sous forme liquide saturé et vapeur saturée, état du fluide sous forme de vapeur surchauffée.

Dans les états saturés, les grandeurs physico-chimiques sont exprimées en fonction de la température de saturation $T_s$. Par contre, pour les états hors saturation (vapeurs surchauffées, par exemple), les grandeurs s'expriment en fonction de la température de saturation $T_s$ et la température réelle T.

L'intervalle étudié et reproduit par des équations polynomiales est représenté sur le diagramme enthalpique de Mollier (Fig. 1) par les deux droites isobares correspondantes respectivement à $T_s = -40°C$ et $T_s =90°C$ et par les deux courbes correspondantes respectivement à la courbe de saturation à l'état liquide et la courbe de la vapeur surchauffée à 120°C.

<figure><img src="image_1.png" /><figcaption>Fig. 1: Diagramme de Mollier du fluide frigorigène R134a</figcaption></figure>
## 2. PROPRIETES GENERALES DU FLUIDE FRIGORIGENE R134a

Le fluide frigorigène R134a est un hydrofluorocarbone 'HFC' de la famille de l'éthane. Il a pour formule organique CH₂F-CF₃. Ses principales propriétés sont résumées dans le tableau suivant:

<table><caption>Tableau 1: Caractéristiques physiques du R134a [3]</caption><thead><tr><td>Propriétés</td><td>Valeur</td></tr></thead><tbody><tr><td>Température d'ébullition</td><td>- 26.07 °C</td></tr><tr><td>Température de fusion</td><td>- 103.3 °C</td></tr><tr><td>Température critique</td><td>+ 101.06 °C</td></tr><tr><td>Pression critique</td><td>40.59 Bar abs.</td></tr><tr><td>Masse molaire</td><td>0.102 kg/mol</td></tr><tr><td>Enthalpie de vaporisation à T<sub>s</sub> = -10°C</td><td>206 kJ/kg.K</td></tr><tr><td>Masse volumique du liquide à T<sub>s</sub> = -10°C</td><td>1328 kg/m<sup>3</sup></td></tr><tr><td>Masse volumique de la vapeur à T<sub>s</sub> = -10°C</td><td>10.04 kg/m<sup>3</sup></td></tr></tbody></table>

## 3. SIMULATION DES PROPRIETES THERMODYNAMIQUES DU FLUIDE R134a A L'EQUILIBRE LIQUIDE - VAPEUR

Les principales propriétés thermodynamiques du fluide frigorigène R134a, représentées par des équations polynomiales à l'état de saturation sont: la pression de saturation ou la tension de vapeur P, l'enthalpie H des phases liquide et vapeur, le volume massique V et la masse volumique ρ des deux phases, l'entropie S des phases liquide et vapeur.

Les équations de ces grandeurs, données ci-après, et décrivant au mieux les points expérimentaux [3], ont été développées grâce à des méthodes numériques.

* **Pression de saturation [bar]**

$$P = p_1 T_s^4 + p_2 T_s^3 + p_3 T_s^2 + p_4 T_s + p_5 \quad (1)$$

$$p_1 = 2.12766 \times 10^{-8}; p_2 = 8.9522 \times 10^{-8}; p_3 = 0.00147189; p_4 = 0.106194; p_5 = 2.93554$$

* **Enthalpie de la phase liquide [kJ/kg]**

$$H_{liq} = h_{l1} T_s^4 + h_{l2} T_s^3 + h_{l3} T_s^2 + h_{l4} T_s + h_{l5} \quad (2)$$

$$h_{l1} = 1.43312 \times 10^{-7}; h_{l2} = -1.86962 \times 10^{-6}; h_{l3} = 1.9844 \times 10^{-3}; h_{l4} = 1.30623; h_{l5} = 200.041$$

* **Enthalpie de la phase vapeur [kJ/kg]**

$$H_{vap} = h_{v1} T_s^4 + h_{v2} T_s^3 + h_{v3} T_s^2 + h_{v4} T_s + h_{v5} \quad (3)$$

$$h_{v1} = -1.70029 \times 10^{-5}; h_{v2} = -0.844368 \times 10^{-3}; h_{v3} = 0.588189; h_{v4} = 396.476$$

* **Masse volumique de la phase liquide [kg/l]**

$$\rho_{liq} = \rho_{l1} T_s^4 + \rho_{l2} T_s^3 + \rho_{l3} T_s^2 + \rho_{l4} T_s + \rho_{l5} \quad (4)$$

$$\rho_{l1} = -1.12536 \times 10^{-9}; \rho_{l2} = -2.5417 \times 10^{-8}; \rho_{l3} = -7 \times 10^{-6}; \rho_{l4} = -3.28683 \times 10^{-3}; \rho_{l5} = 1.296547$$

* **Masse volumique de la phase vapeur [kg/m³]**

$$\rho_{vap} = \rho_{v1} T_s^4 + \rho_{v2} T_s^3 + \rho_{v3} T_s^2 + \rho_{v4} T_s + \rho_{v5} \quad (5)$$

$$\rho_{v1} = 1.0439 \times 10^{-6}; \rho_{v2} = 3.4835 \times 10^{-5}; \rho_{v3} = 5.66156 \times 10^{-3}; \rho_{v4} = 0.526964; \rho_{v5} = 14.9159$$
### • Volume massique de la phase liquide [l/kg]

$$V_{liq} = v_{l1} T_s^4 + v_{l2} T_s^3 + v_{l3} T_s^2 + v_{l4} T_s + v_{l5} \quad (6)$$

$$v_{l1} = 2.14945 \ 10^{-9}; v_{l2} = 2.3399 \ 10^{-8}; v_{l3} = 6.98645 \ 10^{-6}; v_{l4} = 1.99174 \ 10^{-3}; v_{l5} = 0.772442$$

### • Volume massique de la phase vapeur [l/kg]

$$V_{vap} = v_{v1} T_s^5 + v_{v2} T_s^4 + v_{v3} T_s^3 + v_{v4} T_s^2 + v_{v5} T_s + v_{v6} \quad (7)$$

$$v_{v1} = -7.40927 \ 10^{-8}; v_{v2} = 1.41186 \ 10^{-5}; v_{v3} = -1.05562 \ 10^{-3}; v_{v4} = 5.08282 \ 10^{-2}; v_{v5} = -2.28546; v_{v6} = 68.3095$$

### • Entropie de la phase liquide [kJ/kg.K]

$$S_{liq} = s_{l1} T_s^3 + s_{l2} T_s^2 + s_{l3} T_s + s_{l4} \quad (8)$$

$$s_{l1} = 2.58087 \ 10^{-8}; s_{l2} = -1.31995 \ 10^{-6}; s_{l3} = 4.71104 \ 10^{-3}; s_{l4} = 1$$

### • Entropie de la phase vapeur [kJ/kg.K]

$$S_{vap} = s_{v1} T_s^3 + s_{v2} T_s^2 + s_{v3} T_s + s_{v4} \quad (9)$$

$$s_{v1} = -7.17036 \ 10^{-8}; s_{v2} = 7.23032 \ 10^{-6}; s_{v3} = -5.36793 \ 10^{-4}; s_{v4} = 1.71932$$

## 4. SIMULATION DES PROPRIETES THERMODYNAMIQUES DE LA VAPEUR SURCHAUFFEE DU R134a

Les deux grandeurs thermodynamiques du fluide R134a, du domaine de la phase vapeur surchauffée, qui nous intéressent dans cette étude sont l'enthalpie et le volume spécifique.

### 4.1 Détermination de la relation de l'enthalpie de la phase vapeur surchauffée

La forme polynomiale proposée, relative à l'enthalpie de la vapeur surchauffée du fluide frigorigène R134a, est la suivante:

$$H = AT^2 + BT + C \quad (10)$$

L'intervalle d'étude est compris entre -40 et +90 °C pour la température de saturation et la température réelle maximale utilisée est de 120 °C.

Sur le tableau 2, sont présentées les constantes calculées A, B et C de l'expression (10) de l'enthalpie de la phase vapeur surchauffée. Ces coefficients sont aussi fonction de la température de saturation $T_s$.

<table><caption>Tableau 2: Coefficients de l'expression de l'enthalpie de la phase vapeur surchauffée</caption><thead><tr><th>T<sub>s</sub> (°C)</th><th>P<sub>v</sub> (Bar)</th><th>A</th><th>B</th><th>C</th></tr></thead><tbody><tr><td>- 40</td><td>0.521</td><td>0.00084872</td><td>0.821593</td><td>403.887</td></tr><tr><td>- 30</td><td>0.852</td><td>0.00079678</td><td>0.833264</td><td>402.918</td></tr><tr><td>- 20</td><td>1.335</td><td>0.00072816</td><td>0.849649</td><td>401.497</td></tr><tr><td>- 10</td><td>2.012</td><td>0.00063839</td><td>0.872131</td><td>399.489</td></tr><tr><td>0</td><td>2.933</td><td>0.00052241</td><td>0.902526</td><td>396.699</td></tr><tr><td>+ 10</td><td>4.151</td><td>0.00036886</td><td>0.943631</td><td>392.898</td></tr><tr><td>+ 20</td><td>5.722</td><td>0.00016687</td><td>0.999072</td><td>387.755</td></tr><tr><td>+ 30</td><td>7.707</td><td>-0.00010680</td><td>1.07481</td><td>380.804</td></tr><tr><td>+ 40</td><td>10.17</td><td>-0.00048777</td><td>1.18075</td><td>371.287</td></tr><tr><td>+ 50</td><td>13.18</td><td>-0.00105325</td><td>1.33639</td><td>357.850</td></tr><tr><td>+ 60</td><td>16.813</td><td>-0.00178298</td><td>1.55111</td><td>339.083</td></tr><tr><td>+ 70</td><td>21.156</td><td>-0.00141492</td><td>1.58392</td><td>325.349</td></tr><tr><td>+ 80</td><td>26.312</td><td>-0.00213763</td><td>1.86366</td><td>295.681</td></tr><tr><td>+ 90</td><td>32.417</td><td>-0.00398695</td><td>2.4881</td><td>238.528</td></tr></tbody></table>
La variation des constantes A, B et C de l'expression (10) en fonction de la température de saturation T_s est donnée dans les relations (11-13) suivantes:

Variation de A en fonction de T_s

$$A = A_1 T_s^3 + A_2 T_s^2 + A_3 T_s + A_4 \quad (11)$$
$$A_1 = -2.25 \times 10^{-9}; \quad A_2 = -1.49 \times 10^{-7}; \quad A_3 = -1.3 \times 10^{-5}; \quad A_4 = 0.477 \times 10^{-3}$$

Variation de B en fonction de T_s

$$B = B_1 T_s^3 + B_2 T_s^2 + B_3 T_s + B_4 \quad (12)$$
$$B_1 = 1.3 \times 10^{-6}; \quad B_2 = 3.5 \times 10^{-5}; \quad B_3 = 2.44 \times 10^{-3}; \quad B_4 = 0.917412$$

Variation de C en fonction de T_s

$$C = C_1 T_s^3 + C_2 T_s^2 + C_3 T_s + C_4 \quad (13)$$
$$C_1 = -0.000147; \quad C_2 = -0.003073; \quad C_3 = -0.17755; \quad C_4 = 395.399$$

## 4.2 Détermination de l'équation du volume spécifique de la phase vapeur suchauffée

Pour la variation du volume spécifique en fonction de la température réelle, la forme polynomiale suivante est proposée :

$$V = \frac{1}{\rho} = DT^3 + ET^2 + FT + G \quad (14)$$

Les coefficients D, E, F et G calculés de l'expression (14) du volume spécifique pour la phase vapeur suchauffée sont présentés dans le tableau 3. Ces coefficients sont également fonction de la température de saturation T_s.

<table><caption>Tableau 3: Coefficients calculés de l'équation du volume spécifique de la vapeur suchauffée</caption><thead><tr><th>T<sub>s</sub> (°C)</th><th>Pv (Bar)</th><th>D</th><th>E</th><th>F</th><th>G</th></tr></thead><tbody><tr><td>- 40</td><td>0.521</td><td>0</td><td>-3.69437.10<sup>-7</sup></td><td>0.00164474</td><td>0.420233</td></tr><tr><td>- 30</td><td>0.852</td><td>2.24843.10<sup>-9</sup></td><td>-6.48696.10<sup>-7</sup></td><td>0.00103962</td><td>0.254432</td></tr><tr><td>- 20</td><td>1.335</td><td>1.88925.10<sup>-9</sup></td><td>-5.99542.10<sup>-7</sup></td><td>0.00069337</td><td>0.15981</td></tr><tr><td>- 10</td><td>2.012</td><td>1.81159.10<sup>-9</sup></td><td>-5.94648.10<sup>-7</sup></td><td>0.000489749</td><td>0.103706</td></tr><tr><td>0</td><td>2.933</td><td>1.90983.10<sup>-9</sup></td><td>-6.18908.10<sup>-7</sup></td><td>0.000364568</td><td>0.0687191</td></tr><tr><td>+ 10</td><td>4.151</td><td>1.71887.10<sup>-9</sup></td><td>-6.04293.10<sup>-7</sup></td><td>0.000286045</td><td>0.0462147</td></tr><tr><td>+ 20</td><td>5.722</td><td>1.68273.10<sup>-9</sup></td><td>-6.06575.10<sup>-7</sup></td><td>0.000234871</td><td>0.0312576</td></tr><tr><td>+ 30</td><td>7.707</td><td>1.7575.10<sup>-9</sup></td><td>-6.67.10<sup>-7</sup></td><td>0.000207923</td><td>0.0207409</td></tr><tr><td>+ 40</td><td>10.17</td><td>2.47677.10<sup>-9</sup></td><td>-8.92155.10<sup>-7</sup></td><td>0.000206107</td><td>0.0128494</td></tr><tr><td>+ 50</td><td>13.18</td><td>2.49708.10<sup>-9</sup></td><td>-9.59317.10<sup>-7</sup></td><td>0.00020316</td><td>0.00692612</td></tr><tr><td>+ 60</td><td>16.813</td><td>4.662.10<sup>-9</sup></td><td>-1.62637.10<sup>-6</sup></td><td>0.000260283</td><td>0.000647457</td></tr><tr><td>+ 70</td><td>21.156</td><td>9.79022.10<sup>-9</sup></td><td>-3.31236.10<sup>-6</sup></td><td>0.000435921</td><td>-0.00901819</td></tr></tbody></table>

La variation de ces coefficients avec T_s est présentée par les expressions suivantes (15-18).

Variation de D en fonction de T_s

$$D = D_1 T_s^7 + D_2 T_s^6 + D_3 T_s^5 + D_4 T_s^4 + D_5 T_s^3 + D_6 T_s^2 + D_7 T_s + D_8 \quad (15)$$
$$D_1 = 7.35046 \times 10^{-21}; \quad D_2 = -8.01992 \times 10^{-19}; \quad D_3 = 1.63598 \times 10^{-17}; \quad D_4 = 1.02313 \times 10^{-15};$$
$$D_5 = -2.99745 \times 10^{-14}; \quad D_6 = -2.20572 \times 10^{-14}; \quad D_7 = 2.59011 \times 10^{-12}; \quad D_8 = 1.77957 \times 10^{-9}$$

Variation de E en fonction de T_s

$$E = E_1 T_s^6 + E_2 T_s^5 + E_3 T_s^4 + E_4 T_s^3 + E_5 T_s^2 + E_6 T_s + E_7 \quad (16)$$
$E_1 = -3.5515 \ 10^{-17}$; $E_2 = -1.01016 \ 10^{-15}$; $E_3 = 2.16742 \ 10^{-13}$; $E_4 = -2.46083 \ 10^{-12}$; $E_5 = -2.67134 \ 10^{-10}$; $E_6 = -1.24046 \ 10^{-9}$; $E_7 = -5.71967 \ 10^{-7}$

Variation de F en fonction de T<sub>s</sub>

$$F = F_1 T_s^4 + F_2 T_s^3 + F_3 T_s^2 + F_4 T_s + F_5 \quad (17)$$

$$F_1 = 5.96895 \ 10^{-11}; F_2 = -5.76126 \ 10^{-9}; F_3 = 2.52697 \ 10^{-7}; F_4 = -8.84529 \ 10^{-6}; F_5 = 3.58414 \ 10^{-4}$$

Variation du coefficient G en fonction de T<sub>s</sub>

$$G = G_1 T_s^4 + G_2 T_s^3 + G_3 T_s^2 + G_4 T_s + G_5 \quad (18)$$

$$G_1 = 1.05921 \ 10^{-8}; G_2 = -1.52981 \ 10^{-6}; G_3 = 7.747 \ 10^{-5}; G_4 = -2.5723 \ 10^{-3}; G_5 = 0.0645791$$

## 5. COMPARAISON DES RESULTATS ENTRE LES VALEURS EXPERIMENTALES ET CELLES CALCULEES PAR LES MODELES

Sur les figures 2-4, on constate que les relations établies représentent correctement les propriétés thermodynamiques du fluide frigorigène R134a sur la courbe de saturation dans l'intervalle de température compris entre -40 à + 90°C. Il en est de même, pour les relations représentant l'enthalpie et l'entropie de la vapeur surchauffée du fluide R134a (Fig. 5, 6).

<figure><img src="image_2.png" /><figcaption>Fig. 2: Courbe de la pression de saturation du fluide R134a en fonction de T<sub>s</sub></figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 3: Courbe de l'enthalpie du fluide R134a en fonction de T<sub>s</sub></figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Fig. 4: Courbe de l'entropie du R134a en fonction de la température de saturation T<sub>s</sub></figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Fig. 5: Enthalpie de la vapeur surchauffée pour $T_s = -40°C$, $0°C$ et $40°C$</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 6: Volume spécifique de la vapeur surchauffée pour différentes températures $T_s$</figcaption></figure>

Pour une comparaison des résultats obtenus par le biais des expressions précédentes des principales grandeurs thermodynamiques du fluide frigorigène R134a et les résultats expérimentaux obtenus par le logiciel ‘SOLKANE’ [3], nous déterminons, par calcul, les erreurs relatives ‘ER’, ainsi que les écarts types moyens ‘ETM’ pour chacune des expressions polynomiales. La formulation des ‘ER’ et des ‘ETM’ est représentée par les relations suivantes:

$$ER = \frac{1}{N} \sum_{i=1}^{N} \frac{(Y_i^C - Y_i^{Exp})}{Y_i^{Exp}} \times 100 \quad (19)$$

$$ETM = \left[ \frac{1}{N} \sum_{i=1}^{N} \left( \frac{Y_i^C - Y_i^{Exp}}{Y_i^{Exp}} \right)^2 \right]^{1/2} \times 100 \quad (20)$$

où $Y_i^C$ et $Y_i^{Exp}$ sont respectivement les valeurs calculées à l'aide des équations et les valeurs expérimentales, N étant le nombre de points des valeurs expérimentales.

Les résultats obtenus sont résumés dans le tableau suivant:

<table><caption>**Tableau 4:** Comparaison entre valeurs expérimentales et valeurs calculées par les modèles établis à l'état de saturation</caption><thead><tr><td>Grandeurs thermodynamiques</td><td>Érreur Relative (%)</td><td>Écart Type (%)</td></tr></thead><tbody><tr><td>Pression de saturation P</td><td>0.09</td><td>0.157</td></tr><tr><td>Enthalpie phase liquide H<sub>1</sub></td><td>0.02</td><td>0.023</td></tr><tr><td>Enthalpie phase vapeur H<sub>v</sub></td><td>0.04</td><td>0.04</td></tr><tr><td>Masse volumique phase liquide ρ<sub>1</sub></td><td>0.011</td><td>0.20</td></tr><tr><td>Masse volumique phase vapeur ρ<sub>v</sub></td><td>2.40</td><td>4.40</td></tr><tr><td>Volume massique phase liquide V<sub>1</sub></td><td>0.087</td><td>0.01</td></tr><tr><td>Volume massique phase vapeur V<sub>v</sub></td><td>2.70</td><td>4.18</td></tr><tr><td>Entropie phase liquide S<sub>1</sub></td><td>0.035</td><td>0.043</td></tr><tr><td>Entropie phase vapeur S<sub>v</sub></td><td>0.02</td><td>0.025</td></tr></tbody></table>
Les résultats obtenus par les expressions proposées précédemment, en comparaison avec les résultats expérimentaux sont satisfaisants et ont montré une très bonne approche avec une erreur relative maximale de 2.7 % pour le modèle de calcul du volume massique de la vapeur.

Tandis que pour le cas de la vapeur surchauffée du R134a, on remarque d'après le tableau 5 que l'erreur maximale obtenue, lors du calcul de l'enthalpie, n'est que de 0.5 %, ce qui confère au modèle établi une bonne fiabilité dans l'intervalle de température étudié (-40 à 80°C).

<table><caption>**Tableau 5:** Erreurs relatives à l'enthalpie et du volume spécifique de la vapeur surchauffée</caption><thead><tr><th>T<sub>s</sub></th><th>Érreur Relative H (%)</th><th>Érreur Relative V (%)</th></tr></thead><tbody><tr><td>- 40</td><td>0.52</td><td>0.82</td></tr><tr><td>- 30</td><td>0.16</td><td>2.74</td></tr><tr><td>- 20</td><td>0.38</td><td>0.57</td></tr><tr><td>- 10</td><td>0.074</td><td>0.97</td></tr><tr><td>0</td><td>0.14</td><td>4.86</td></tr><tr><td>+ 10</td><td>0.058</td><td>0.96</td></tr><tr><td>+ 20</td><td>0.21</td><td>6.67</td></tr><tr><td>+ 30</td><td>0.27</td><td>11.8</td></tr><tr><td>+ 40</td><td>0.22</td><td>5.70</td></tr><tr><td>+ 50</td><td>0.086</td><td>14.14</td></tr><tr><td>+ 60</td><td>0.10</td><td>30.21</td></tr><tr><td>+ 70</td><td>0.25</td><td>25</td></tr></tbody></table>

Par contre, pour le modèle du volume spécifique, on constate que plus la température de saturation $T_s$ augmente et s'approche du point critique, plus l'erreur relative 'ER' devient importante. Néanmoins, on peut dire que le modèle établi peut être utilisé dans l'intervalle de température de saturation de -40 à 40 avec une erreur relative maximale de 11.7 % à $T_s = 40°C$, domaine utilisé surtout pour les applications du froid où généralement le volume spécifique est déterminé à l'entrée de la compression, donc à des températures basses.

## 6. CONCLUSION

Les expressions établies dans cet article, vérifient correctement les fonctions thermodynamiques fondamentales avec une bonne précision, puisque les erreurs relatives et les écarts types sont respectivement de 2.7 % et 4.4 % pour l'état de saturation et de 0.5 % et 11 % pour le cas de la vapeur surchauffée. Ces équations peuvent être facilement introduites dans les programmes de calcul pour la génération des diagrammes et des cycles thermodynamiques, indispensables dans le cadre de la prédiction des performances d'un réfrigérateur médical, telle est notre but.

## REFERENCES

[1] D. Diaf, 'Etude de Simulation d'un Conservateur Médical Solaire', Mémoire de Magister, Université M'Hamed Bouguerra, Boumersès, Novembre 2003.

[2] Montréal Protocol 1991. 'Report of the Refrigeration, Air Conditioning and Heat Pumps Technical Options Committee', Décembre 1991, PNUE.

[3] Logiciel Solkane, 'Solvay Fluor und Derivate GMBH', Technical service refrigerants, Germany.

[4] J. Raznjevic, 'Handbook of Thermodynamic Tables and Charts', McGraw Hill, 1982.

[5] Logiciel REFPROP, 'Thermodynamic and Transport Properties of Refrigerants and Refrigerant Mixture', Version 6.01, NIST.