Spogis / Psychrometry

Psychrometry Template for DWSIM - Version 1.1
https://www.youtube.com/@CascaGrossaSuprema
GNU General Public License v3.0
3 stars 0 forks source link

Psychrometric calculation errors #1

Open Scorpiomom69 opened 8 months ago

Scorpiomom69 commented 8 months ago

Hi. Thanks for posting this nice tool but I wanted to point out apparent errors in the calculation. To clarify, I have a masters degree strongly focussed on humid airflow calculations, which is where my comment originates. If you go to a site like https://www.kwangu.com/work/psychrometric.htm you can calculate the psychrometric properties of most humid airstreams pretty accurately. What's particularly interesting is the humidity ratio given. It must be noted that this is not a ratio of vapor mass to total mass, this is a ratio of vapor mass to dry air mass. The psychrometric calculations use this is because the dry air mass is constant and only the vapor mass changes. Thus the vapor mass ratio should be calculated as: (vapor mass/kg da)/(1kg dry air/kg da + vapor mass/kg da) [da = dry air, vapor mass/kg da is the humidity ratio]. This seems to be incorrectly done in your calculator as I find it impossible in many instances to generate a 100% humid air stream and also the mass fraction of the vapor ['vapor mass ratio'] is currently = to the humidity ratio in the psychrometric popup, which is not possible. I do find it interesting that the density calc works, since it requires the same foreknowledge, but perhaps this is due to rounding or different formulas, idk. Please feel free to contact me about this and thanks again for the tool. image image

Spogis commented 8 months ago

Dear @ Scorpiomom69,

Your informed points are indeed very important, but I believe you might be confusing two different pieces of information.

What you're mentioning to me is what I refer to as Absolute Humidity (g water / kg of dry air). Note that this value in the figure you sent me is 80.29 g water/kg of dry air, which is practically identical to the value obtained from (https://www.kwangu.com/work/psychrometric.htm), which is 0.0803 kg water / kg of dry air. Notice there is a difference in units! Thus, if you convert the units, you'll see that the values are identical.

The Degree of Saturation, used in Psychrometrics, refers to the ratio of the actual amount of water vapor present in the air to the maximum amount of water vapor the air can hold at a given temperature, when it is saturated. Essentially, it measures how "full" the air is with moisture compared to its maximum capacity for moisture at that specific temperature. The degree of saturation is a dimensionless number and is often expressed as a percentage. When air is fully saturated (100% relative humidity), the degree of saturation is 1.0 or 100%, indicating that the air cannot hold any more moisture at that temperature.

Relative Humidity is a measure of the amount of water vapor present in the air compared to the maximum amount of water vapor the air can hold at a given temperature. At 100% relative humidity, the air is saturated, meaning it cannot hold any more water vapor at that temperature.

Thus, the values obtained by this plugin have been extensively compared and validated with the data from the 2009 ASHRAE Handbook—Fundamentals, and with the ASHRAE Psychrometric Chart App (https://www.ashrae.org/technical-resources/bookstore/psychrometrics), global references in the field.

Please contact me if you still have any questions.

Atenciosamente / Best regards,

Prof. Nicolas Spogis, Ph.D. Phone/WhatsApp: https://wa.me/5519998440460 +55 (19) 99844-0460

E-mail: @. @.>

LinkTree: https://linktr.ee/CascaGrossaSuprema

De: Scorpiomom69 @.> Enviada em: sábado, 2 de março de 2024 08:15 Para: Spogis/Psychrometry @.> Cc: Subscribed @.***> Assunto: [Spogis/Psychrometry] Psychrometric calculation errors (Issue #1)

Hi. Thanks for posting this nice tool but I wanted to point out apparent errors in the calculation. To clarify, I have a masters degree strongly focussed on humid airflow calculations, which is where my comment originates. If you go to a site like https://www.kwangu.com/work/psychrometric.htm you can calculate the psychrometric properties of most humid airstreams pretty accurately. What's particularly interesting is the humidity ratio given. It must be noted that this is not a ratio of vapor mass to total mass, this is a ratio of vapor mass to dry air mass. The psychrometric calculations use this is because the dry air mass is constant and only the vapor mass changes. Thus the vapor mass ratio should be calculated as: (vapor mass/kg da)/(1kg dry air/kg da + vapor mass/kg da) [da = dry air, vapor mass/kg da is the humidity ratio]. This seems to be incorrectly done in your calculator as I find it impossible in many instances to generate a 100% humid air stream and also the mass fraction of the vapor ['vapor mass ratio'] is currently = to the humidity ratio in the psychrometric popup, which is not possible. I do find it interesting that the density calc works, since it requires the same foreknowledge, but perhaps this is due to rounding or different formulas, idk. Please feel free to contact me about this and thanks again for the tool. image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/f0554c12-eb99-49db-b869-a9cc7562542c image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/06aa01c2-5e5b-43c0-ab55-313c4f626fc6

— Reply to this email directly, view it on GitHub https://github.com/Spogis/Psychrometry/issues/1 , or unsubscribe https://github.com/notifications/unsubscribe-auth/AIHRJYKKXQ4XBFZGSQ3FOADYWGYCXAVCNFSM6AAAAABEDAROECVHI2DSMVQWIX3LMV43ASLTON2WKOZSGE3DINZWGQYTIMY . You are receiving this because you are subscribed to this thread. https://github.com/notifications/beacon/AIHRJYJC36Q7QTZDNGX3KV3YWGYCXA5CNFSM6AAAAABEDAROECWGG33NNVSW45C7OR4XAZNFJFZXG5LFVJRW63LNMVXHIX3JMTHICB5N54.gif Message ID: @. @.> >

Spogis commented 8 months ago

(vapor mass/kg da)/(1kg dry air/kg da + vapor mass/kg da) [da = dry air, vapor mass/kg da is the humidity ratio]. -> This is not Humidity Ratio! It’s Absolute Humidity (Correct Name by ASHRAE).

Absolute Humidity is a measure of the actual amount of water vapor present in the air, regardless of the air's temperature. Typically expressed in grams of water / kg dry air. Unlike relative humidity, which indicates how saturated the air is with water vapor, absolute humidity quantifies the actual concentration of water vapor in the air.

Atenciosamente / Best regards,

Prof. Nicolas Spogis, Ph.D. Phone/WhatsApp: https://wa.me/5519998440460 +55 (19) 99844-0460

E-mail: @. @.>

LinkTree: https://linktr.ee/CascaGrossaSuprema

De: Scorpiomom69 @.> Enviada em: sábado, 2 de março de 2024 08:15 Para: Spogis/Psychrometry @.> Cc: Subscribed @.***> Assunto: [Spogis/Psychrometry] Psychrometric calculation errors (Issue #1)

Hi. Thanks for posting this nice tool but I wanted to point out apparent errors in the calculation. To clarify, I have a masters degree strongly focussed on humid airflow calculations, which is where my comment originates. If you go to a site like https://www.kwangu.com/work/psychrometric.htm you can calculate the psychrometric properties of most humid airstreams pretty accurately. What's particularly interesting is the humidity ratio given. It must be noted that this is not a ratio of vapor mass to total mass, this is a ratio of vapor mass to dry air mass. The psychrometric calculations use this is because the dry air mass is constant and only the vapor mass changes. Thus the vapor mass ratio should be calculated as: (vapor mass/kg da)/(1kg dry air/kg da + vapor mass/kg da) [da = dry air, vapor mass/kg da is the humidity ratio]. This seems to be incorrectly done in your calculator as I find it impossible in many instances to generate a 100% humid air stream and also the mass fraction of the vapor ['vapor mass ratio'] is currently = to the humidity ratio in the psychrometric popup, which is not possible. I do find it interesting that the density calc works, since it requires the same foreknowledge, but perhaps this is due to rounding or different formulas, idk. Please feel free to contact me about this and thanks again for the tool. image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/f0554c12-eb99-49db-b869-a9cc7562542c image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/06aa01c2-5e5b-43c0-ab55-313c4f626fc6

— Reply to this email directly, view it on GitHub https://github.com/Spogis/Psychrometry/issues/1 , or unsubscribe https://github.com/notifications/unsubscribe-auth/AIHRJYKKXQ4XBFZGSQ3FOADYWGYCXAVCNFSM6AAAAABEDAROECVHI2DSMVQWIX3LMV43ASLTON2WKOZSGE3DINZWGQYTIMY . You are receiving this because you are subscribed to this thread. https://github.com/notifications/beacon/AIHRJYJC36Q7QTZDNGX3KV3YWGYCXA5CNFSM6AAAAABEDAROECWGG33NNVSW45C7OR4XAZNFJFZXG5LFVJRW63LNMVXHIX3JMTHICB5N54.gif Message ID: @. @.> >

Scorpiomom69 commented 8 months ago

Good day Prof.

Thanks for the feedback. I do admit that I wish they bothered to name the absolute humidity something a little more discerning as the terminology gets lost, even on myself. Thus from here I will refer to them as such: The value indicating the absolute humidity in [kg moisture/kd dry air] will be 'apparent specific humidity' and the value indicating the fraction of the total massflow consistsing of water vapor will be 'vapor fraction'. I will also use g/kg to avoid decimals when stating the values but calculations should, of course, be done in the correct units.

-According to the Barenburg charts, (I do not have the ASHRAE charts, but I'm sure they are the same), at 35 C, 100kPa, saturated air should have an apparent specific humidity of 35.4 [g/kg da]. See Image 1. The vapor fraction is calculated as: Vapor fraction = (app. spec. hum.)/(app. spec. hum. + 1 kg da). In this case it comes to 34.2 [g vapor mass/kg total mixture mass]. -I can also calculate this myself using various equations provided in literature and I find the value for the apparent specific humidity to be 37.1 [g/kg da]. This is fairly close already. See Image 2. Subsequently the vapor fraction is 37.1/(1000+37.1) = 35.7 [g vapor mass/kg total mixture mass]. -Your applet is unable to produce a saturated air stream when I specify the 100% humidity condition. It does some calculations, unknown to me, and then reduces the humidity to approx. 96 % (sometimes less). See Image 3 & 4. From Image 4 you can also see that it predicts the vapor fraction to be 37.091345 [g vapor mass/kg total mass]. But as you can see from the two samples above, this is NOT actually the vapor fraction. In reality this is the apparent specific humidity value. The true vapor fraction has already been shown above, but to drive it home let's reverse the calculation to get the "apparent specific humidity" value from this "vapor fraction": (vapor fraction)/(1000-vapor fraction) = apparent specific humidity = 38.52 [g/kg da], which is even higher than my calculated values.

Now I don't doubt that your Psychrometric charts may have some increased accuracy over mine, but coupling the applet's inability to calculate the saturated conditions and the fact that I have fed it several points for which it produces a "vapor fraction" that exactly matches what I expected for the apparent specific humidity, I have concluded that it is not functioning as intended.

I hope to have informed you sufficiently and look forward to your response. We all learn at different times and places and perhaps you have some information I am missing. I am not too proud to be taught. Kind Regards. Ir. Austin Labuschagne

Image 1: Image 1

Image 2: Image 2

Image 3: Image 3

Image 4: Image 4

Spogis commented 8 months ago

Open again for review.

Spogis commented 8 months ago

Dear @Scorpiomom69

I upload the Template again with some changes. Could you download the new version and make some checks again?

Thank you very much for your help.

Best regards,

Nicolas

Spogis commented 8 months ago

I'll need to take a deeper look also in the best Thermodynamic models to use on these cases.

Scorpiomom69 commented 8 months ago

Good day Nicolas

Thanks for taking another look! I also checked and it seems to be working now, except for the mass-frac value in the simulation: image The value being displayed on the left is the 'absolute humdity' but you should be displaying the water mass fraction (vapor fraction). To display the right value use: (abs hum)/(abs hum+1) = vapor fraction [In this case 35.73g] image

Spogis commented 8 months ago

Dear @Scorpiomom69,

How are you getting on?

Thank you very much for your very important contribution.

I just uploaded a new version (1.1), with your last request.

Now the absolute humidity is fine (thank you very much again).

Please add me on WhatsApp (if is possible). I would like to discuss with you other important questions about this template.

Atenciosamente / Best regards,

Prof. Nicolas Spogis, Ph.D. Phone/WhatsApp: https://wa.me/5519998440460 +55 (19) 99844-0460

E-mail: @. @.>

LinkTree: https://linktr.ee/CascaGrossaSuprema

De: Scorpiomom69 @.> Enviada em: segunda-feira, 4 de março de 2024 14:40 Para: Spogis/Psychrometry @.> Cc: Nicolas Spogis @.>; State change @.> Assunto: Re: [Spogis/Psychrometry] Psychrometric calculation errors (Issue #1)

Good day Nicolas

Thanks for taking another look! I also checked and it seems to be working now, except for the display on the splash screen: image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/d4fe32e2-5e9a-4d28-a008-f23024923cb5 As we discussed earlier, you corrected the calculations but the value in the simulation is still not corrected. The value being displayed is the 'absolute humdiity' but you should be displaying the water mass fraction (vapor fraction). To display the right value use: (abs hum)/(abs hum+1) = apparent specific humidity [In this case 35.73g] image.png (view on web) https://github.com/Spogis/Psychrometry/assets/76644903/c2307725-21b2-428e-aabc-fdc28922e84d

— Reply to this email directly, view it on GitHub https://github.com/Spogis/Psychrometry/issues/1#issuecomment-1977125819 , or unsubscribe https://github.com/notifications/unsubscribe-auth/AIHRJYPPKEQFC3I5TB2W25DYWSWXHAVCNFSM6AAAAABEDAROECVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNZXGEZDKOBRHE . You are receiving this because you modified the open/close state. https://github.com/notifications/beacon/AIHRJYMZRT4G3RAZEIOKV6LYWSWXHA5CNFSM6AAAAABEDAROECWGG33NNVSW45C7OR4XAZNMJFZXG5LFINXW23LFNZ2KUY3PNVWWK3TUL5UWJTTV3CF3W.gif Message ID: @. @.> >

Scorpiomom69 commented 8 months ago

Hi Prof. I am doing well. How are you?

Indeed the calculations now work but they are going in the wrong direction. I.E. you are using x=f(y) but is should be y=f(x). See: image

I have added you on WA and sent you a message. I'm not sure if you wanted me to post this last issue there or here?

Kind Regards