psychrolib icon indicating copy to clipboard operation
psychrolib copied to clipboard

Calculations of wet bulb when dry bulb above 100 C converge to unreasonable value

Open dmey opened this issue 5 years ago • 3 comments

The calculations of wet bulb when dry bulb when above 100 C are not handled properly. The convergence tests implemented (test_GetTDewPointFromVapPres_convergence) converge to unreasonable values when hen above 100 C. For example GetTWetBulbFromRelHum(TDryBulb, RelHum, Pressure) with TDryBulb = 25, RelHum = 1 and Pressure = 101325. As expected the function returns 25 (when the air is saturated, wet bulb = dew point = dry bulb). With TDryBulb = 99 it returns 99.. When called with TDryBulb = 101 it returns 43.

dmey avatar Mar 21 '19 09:03 dmey

@didierthevenard we may want to revisit this at the next release -- do you have any thoughts?

dmey avatar Apr 08 '20 09:04 dmey

Hi, I'm sorry this is a bit off topic but maybe it's useful for review:

I'm testing two libraries in javascript: coolprop for thermodynamics and mathjs for units and other math stuff. When testing HAprops('DewPoint', { 'Tdb': 25 degC, 'P': 1 atm, 'RelHum': 100 % }) to degC The results is: 25.000000000064 degC

You can benchmark here if you would like: https://github.com/dvd101x/CoolPropJavascriptDemo

dvd101x avatar Jan 09 '22 03:01 dvd101x

I apologise if this is completely obvious, but it seems that GetTWetBulbFromRelHum should be evaluating Psat * RelHum at the very beginning, and verifying that it is not larger than than the actual value of Pressure, which would indicate supersaturation, and therefore "unphysical" in the world of psychrolib. If not obvious, did I misunderstand something?

PHCEAC avatar May 05 '23 19:05 PHCEAC