太阳辐射下喷雾蒸发模型,温度场问题

模型情况:小房子的模型(半透明结构,2.5×4.5×4m),模拟在太阳辐射条件下内部水喷雾情况下房子内部的温度场,有两个开口(实际对应门和排气窗)。

模型设置:两个开口均设为pressure outlet(0 Pa),其余均为wall(semi-transparent,对太阳辐射具有一定吸收率,四周和屋顶采用对流和辐射混合换热边界条件,地面采用绝热条件)。采用DO模型solar tracing加载太阳辐射,采用DPM模型设置两个喷头(injection, group, 2组,unsteady particle tracking),湍流模型为RNG k-epsilon,开启组分传输。
DPM中屋顶设置为reflect,地面和墙壁均设为trap,开口设置为escape.
其它详细设置如以下图片所示。


网格情况:结构化网格(质量>0.9),网格总量235万+,网格尺寸最大0.06m左右,内壁均设置边界层(0.0001m,1.1,6层),喷头处加密,尺寸均有过渡。如下图所示


求解用simpleC,二阶迎风格式。
松弛因子设置如下图

计算过程中continuity没有收敛的迹象。同时不断出现
temperature limited to 1.000000e+00 in 4 cells on zone 33 in domain 1
temperature limited to 5.000000e+03 in 2 cells on zone 33 in domain 1
等提示。
6000步时暂停计算,查看温度场发现壁面个别区域(特别是地板)、以及喷水位置以下区域温度很低,低于-20℃。如下图

同时,对流流场如下图(流速较小)

看了下DPM summary,结果如下:

我感觉是蒸发模型出了问题,但是不知道如何着手去查这个问题。
求大佬们指点,感谢!