Hello Yunhao,
I thought what is written in the WAQTEL user manual and my explanations were clear, but it seems not to you.
With AED2, you can only user WATER QUALITY PROCESS = 13. So please do not try something else. In you case, you apply 2 times a computation of volumetric source term and that should explain your differences.
Anyway, if using this choice (= 13) and having differences of temperature between coupling with THERMIC module (= 11), this may be due to the extinction coefficient computed by AED2 that overwrites its equivalent value given to THERMIC.
Chi-Tuan