[alex2293]

Hello everyone,

This winter I wrote the following UDF to calculate the saturation temperature of each cell in the domain. I am trying to simulate the water vapor condensation of humid air in heat exchanger.

The interpreted UDF is working on my computer. I want to run it on a Linux HPC.

However when I initialize the solution I get Saturation_temperature <= 0 on node 0 if running in serial and node 0 and 1 if running on 2 nodes.

Do you have any clue for me?

#include "udf.h"
DEFINE_PROPERTY(saturation_temperature,c,thread)
{
real p_operating;
real p_pressure;
real abs_pressure;
real p_press;
real sat_temp;
real w;
#if !RP_HOST
Thread *t_humid_air = THREAD_SUB_THREAD(thread,1);
Thread *t_wl = THREAD_SUB_THREAD(thread,0);
w = C_YI(c,t_humid_air,0);
p_pressure = C_P(c,t_humid_air);
p_operating = RP_Get_Real("operating-pressure");
abs_pressure = p_operating + p_pressure;
p_press = ((w/18)/((1-w)/27))*abs_pressure;
sat_temp = 13.57*log(p_press)+187,67;

if (sat_temp<=0) {
sat_temp =273;
return sat_temp;
} else {
return sat_temp;
}
#endif
}