I am attempting to simulate a 1-D Stefan problem (with phase change) using Fluent. The issue at hand is condensation. The domain is meshed using a single cell layer in the y-direction. Other simulation parameters are illustrated in Figure (a). However, the obtained pressure distribution, depicted in Figure (b), does not match the theoretical calculation. The pressure distribution in the vapor region (black line) differs from the expected result (red line). Is there a possible solution to this discrepancy? Thank you very much!