Thank you for your reply. Could you please tell me your opinion about changing the Advection scheme to Upwind and the transient scheme to first-order backward Euler? I tried these options and the results are more close to the experiment. I was wondering why the more precise scheme gives me the wrong answer.