Diagramme psychrométrique 3D

[BastienBastien]

Bonjour à tous,

J'étudie les systèmes HVAC (clim'), en ce moment, et voudrais - juste pour le fun - voir un diagramme psychrométrique en 3D.

La 3ème dimension étant la pression.

Habituellement, ce genre de diagramme est toujours donné à pression atymosphérique.

Je me demande ce que donne une telle représentation....

Merci.
-----

[BastienBastien]

Bonjour à tous,

J'ai trouvé, sur FuturaScience, le fichier Excel ci-joint. Il s'agit d'une feuille de calcul avec les Entrées/Sorties suivantes :

Entrées :

• Tambiante (ex. 20.0 °C)

• Hr (ex. 50%)

Sorties :

• 100%Hr (ex. 17.12 g/m3 d'air)

• Ha (ex. 8.56 g/m3 d'air)

• Th (ex. 13.8 °C)

Je me suis amusé à tracer les courbes suivantes (avec 't' la température d'air sec et 'p' la pression) :

• 2D, pression constante : 100%Hr = f(t);

• 3D, pression variable : 100%Hr = f(t, p);

Voici le code MAXIMA :

## Plot 2D (à pression constante) :

A : 18.3036$
B : 3816.44$
C : -46.13$
P : 1.01325$
rH(t) := ((( ((exp(A - B/(t+273.16+C))) * (P / 760 ))/P ) * 18.015 ) * ( 1000 / (((8.3145*(t+273))/(P*100000))*1000) ))$
plot2d(rH(t), [t,-15,30], [legend,"Courbe de saturation"], [xlabel,"Temperature d'air sec (Celcius)"], [ylabel,"(g d'eau et 1 kg d'air sec)"]);

_________________________________________________________________________________________________

## Vérification de valeurs :

rH(t,alpha) := ((( ((exp(A - B/(t+273.16+C))) * (P*alpha / 760 ))/P*alpha ) * 18.015 ) * ( 1000 / (((8.3145*(t+273))/(P*alpha*100000))*1000) ));
float( rH(t,1) );

_________________________________________________________________________________________________

## Plot 3D (à pression variable) :

A : 18.3036$
B : 3816.44$
C : -46.13$
P : 1.01325$
rH(t,alpha) := ((( ((exp(A - B/(t+273.16+C))) * (P*alpha / 760 ))/P*alpha ) * 18.015 ) * ( 1000 / (((8.3145*(t+273))/(P*alpha*100000))*1000) ))$
plot3d( rH(t,alpha), [t,-15,30], [alpha,0,2], [grid,40,40]);

Malheureusement, je me suis rendu compte, après coup, que la sortie 100%Hr est exprimée en g/m3 d'air et non pas en g/kg_d_air.
Mais je note que mon diagramme obtenu ressemble fortement au diagamme psychrométrique, même si j'ai toujours 2 ou 3 grammes d'erreur dans mes vérifications de valeur.

Quelqu'un peut-il m'aider à obtenir un vrai diag' psychro' ?

Merci.

[BastienBastien]

Bonjour à tous,

Voici deux plots plus clairs (plus grands et avec des labels sur les axes).

Si quelqu'un sait comment je dois modifier l'équation pour avoir un vrai diagramme psychrométrique, je suis preneur !

Le problème actuel, c'est que l'axe 'z' est exprimé en g d'eau / m3 d'air alors que ça devrait être en g d'eau / kg d'air.

L'allure est la même, mais toujours j'ai quelques g d'erreur...

Merci.

[BastienBastien]

Hello,

J'abandonne. Je suis passé à autre chose. Si quelqu'un propose la solution, je suis preneur ! Mais je ne peux pas passer plus de temps sur ce sujet.

+

[A voir en vidéo sur Futura]

[mayonaise]

Bonjour BastienBastien,

Je te promet, j'avais fait une superbe réponse, mais il est tard et j'ai cliqué sur ce @$#%£@&! de bouton répondre à la discussion en haut à gauche du cadre de réponse rapide et mon chef d’œuvre à disparu.
Je vais donc être bref. Ma source (en anglais, mais c'est très complêt pour les curieux) : http://www.natmus.dk/cons/tp/atmcalc/atmoclc1.htm

Tu as besoin de connaitre Hr (Humidité relative), P (pression ambiante (l'unité de ton choix, on va travailler avec des rapports)) et T (température ambiante en °C)

D’abord calcules la pression de saturation PS :
Ps = 610.78 *exp( T / ( T + 238.3 ) *17.2694 )

Ensuite la pression partielle de la vapeur d'eau dans ton air p :
p = Ps * Hr

Ensuite, la masse molaire de l'air sec étant 0.029kg et celle de l'eau de 0.018kg, Avogadro te dis que :
kg de vapeur d'eau/kg d'air sec =
=0.018 *p / ( 0.029 *(P - p ) )
= 0.62 *p / (P - p )

Voilà, je sais pas si ça te servira 2 ans après ton post, mais le tien m'a bien servi dans mes recherches alors je complète, si c'est pas pour toi, ce sera pour un autre.

Merci encore et bonne soirée,

Mayonaise

[BastienBastien]

Hello !

Je vais voir ça dès que je peux et faire un joli "plot" dans Scilab/Octave/Gnu_plot !

Merci beaucoup !

Bonne journée à toi et à tous les forumeurs !