mass conservation problem in interCondensingEvaporatingFoam

Dear OF Developers,

I am working with interCondensingEvaporatingFoam solver. I have encountered a possibly very important problem. I am calculating flow with a phase change in-between heated plates. Domain has one inlet, one outlet, and a body limited by two heated plates. At the inlet the liquid enters the domain, is heated, partially evaporates and leaves the domain via outlet.

The problem is that the mass is not conserved. There is approximately 6x larger mass flow of the liquid/vapor mixture at the outlet that at the inlet. Please see the attached picture, which shows some filed values and integrated 'u_x*rho' over the inlet and outlet patch (u_x in normal to the inlet/outlet, so it gives mass flow).

upper row of the figure shows the data for the inlet; the value "U_x_times_rho_inlet" in the table on the right corresponds to inlet mass flow lower row of the figure shows the data for the outlet; the value "U_x_times_rho_outlet" in the table on the right corresponds to outlet mass flow

Regards Jim