2. Solution de référence#

2.1. Méthode de calcul#

2.1.1. Calcul de la pression de vapeur à partir de la température#

Nous supposons la courbe de saturation linéaire. Elle s’écrit donc:

_images/Object_7.svg

éq 2.1.1-1

L’équation [éq 2.2.3.3-2] du document de référence [R7.01.11] donne alors:

_images/Object_8.svg

éq 2.1.1-2

On écrit que la masse totale d’eau et la masse totale d’air sont conservées (car il n’y a pas de flux d’eau ni de gaz au bord) et on obtient:

_images/Object_9.svg _images/Object_10.svg

éq 2.1.1-3

_images/Object_11.svg _images/Object_12.svg

éq 2.1.1-4

[R7.01.11] [éq 4.1.4-1] donne par ailleurs:

_images/Object_13.svg

éq 2.1.1-5

Le couplage des équations [éq 2.1.1-3], [éq 2.1.1-4] et [éq 2.1.1-5], auquel il faut ajouter l’équation des gaz parfaits pour la vapeur, l’air sec et l’air dissous ainsi que la loi de Henry est un système fortement non linéaire que nous résoudrons en petites perturbations, ce qui permet de le linéariser.

Tous calculs faits, on obtient:

_images/Object_14.svg

éq 2.1.1-6

2.1.2. Calcul de la température#

L’équation [éq 3.2.4.3-1] du document de référence [R7.01.11] donne:

_images/Object_15.svg

éq 2.1.2-1

(puisque les autres coefficients de dilatation sont nuls).

L’équation [éq 3.2.4.3-2] donne:

_images/Object_16.svg

éq 2.1.2-2

On obtient donc:

_images/Object_17.svg

éq 2.1.2-3

Dans ce problème,

_images/Object_18.svg

n’est rien d’autre que la chaleur apportée par unité de volume.

En appelant

_images/Object_19.svg

le volume total de la pièce et

_images/Object_20.svg

sa surface latérale et

_images/Object_21.svg

le temps d’application des flux:

_images/Object_22.svg

éq 2.1.2-4

2.1.3. Système à résoudre#

_images/Object_23.svg

éq 2.1.3-1

_images/Object_24.svg _images/Object_25.svg _images/Object_26.svg _images/Object_27.svg _images/Object_28.svg _images/Object_29.svg

(calculé)

_images/Object_30.svg

5,00E-01

-1,00E-12

3,00E+02

3,70E+03

2,50E+06

2,67E-02

1,00E+03

_images/Object_31.svg _images/Object_32.svg _images/Object_33.svg

(calculé)

_images/Object_34.svg _images/Object_35.svg

l

_images/Object_36.svg _images/Object_37.svg

(calculé)

2,20E+03

3,00E-01

2,93E+03

1,05E+03

4,18E+03

1,90E+03

2,78E+06

_images/Object_38.svg _images/Object_39.svg _images/Object_40.svg _images/Object_41.svg

1,00E+06

10

6,00E+04

1,00E+06

On obtient les résultats suivants:

Après résolution de ce système, on obtient:

_images/Object_42.svg

Ce qui donne en terme de résultat Aster(incrément) :

\(\mathrm{PRE1}\)

\(\mathrm{PRE2}\)

\(\mathrm{DT}\)

\(\mathrm{PVP}\) (\(\mathrm{V3}\) )

1.49E5

113

0.216

44

2.2. Incertitudes#

Les incertitudes sont assez grandes parce-que la solution analytique est une solution approchée du fait de la linéarisation des équations.