assigning the right enthalpy
I'm running a thermal simulation using TMVOC where I'm also including a mass rate at the second layer beneath the atmospheric boundary layer. The mass rate represents annual precipitation. I would like the temperature of this water infiltration to be at 25 deg C. So, I used a single element case and initialized it with only water at 25 deg C and then ran it. By looking at the output file I obtained the enthalpy of water at that temperature and assigned this value in the main simulation. However, when I look at the results, it seems the temperature keeps falling below the target 25 deg C. Any idea why this is happening? What is the best way to determine an enthalpy corresponding to a temperature in TOUGH2?
In principle, the process you described would be the best way to determine the enthalpy of your injectate.
However, in order to address this issue correctly, I'll need to know more details about the "main simulation" that causes your temperature to drop below 25 degree C. For example, what is the temperature of the "atmospheric" boundary layer (which I assume is a single large element or "dummy" element)? If it is below 25 C, then it seems the heat flux coming from it would dominate the heat introduced from the injection source term.
Could you provide more details about the simulation you're trying to set up?Reply
Assigning enthalpy to an injectate to achieve a desired injection temperature is not generally possible because the process of injection can change P-T conditions which then change effective heat capacity making the chosen enthalpy inconsistent with the desired injection temperature. The TOUGH trick to injecting at constant T is to inject the fluid into an extra grid block with a rock type with infinite Cp that is initialized at whatever injection T you want. This block should then be connected to your actual domain where the injection is occurring. This grid block will never change T (because of its infinite Cp) no matter what enthalpy you assign to fluid injected into it.
I cannot know if this issue is affecting your case and follow up to Mikey's question is needed. But please give the above approach a try.Reply