RELP as a function of a VOC mole fraction

Hi everyone,

I would like to make the relative permeability of the gas phase a function of the mole fraction of the the first VOC. Actually, the first VOC, in this case, is identical to water and is just here to model a concentration in the aqueous phase.

I do not know wich is the very first VOC mole fraction : CMOLEFW(1) or CMOLEFW(2) as the AIR mole fraction might be CMOLEFW(1)...

Thank you in advance for your help!

Cheers,
Quentin

4replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • Quentin,

    The VOCs are in slots CMOLEFW(1) to CMOLEFW(NHC); the NCGs are in slots CMOLEFW(NHC+1) to CMOLEFW(NHC+NCG). This is indicated in the TMVOC manual in the paragraph following Table E.1 (on page 161).

    Stefan

    Reply Like
      • Quentin Giraud
      • PhD Student - Hydrogeology - Contaminated Sites - Modelling
      • Quentin_Giraud
      • 7 mths ago
      • Reported - view

      Stefan Finsterle 
      Thank you very much Stefan for your answer!
      It is all clear now!

      Cheers,

      Quentin

      Reply Like
      • Quentin Giraud
      • PhD Student - Hydrogeology - Contaminated Sites - Modelling
      • Quentin_Giraud
      • 6 mths ago
      • Reported - view

      Stefan Finsterle 

      Hi Stefan,
      When I use CMOLEFW(1) in the modified RELP subroutine, it is a constant through the simulation run.
      I was wondering if you had any hints on how to make sure that during the "CALLs" of the subroutine RELP, the CMOLEFW(1) takes into account the iterative process.

      I tried by replacing CMOLEFW(1) with PAR(NKL2W+NBK+1+NCG+1) but it does not compile...

      Thank you very much in advance for your help,

      Cheers,

      Quentin

      Reply Like
  • Quentin,

    What compilation error did you get?

    Yes, use PAR(NKL2W+NBK+1+NCG+1) (even though it seems this is a hardwired fix; make it general by using appropriate variable indices to indicate the right VOC or NCG). I also think it should be NB instead of NBK. Finally, make sure you USE the right modules to have all the needed variables available in RELP. looking at your compilation error should help you identify the issue rather quickly.

    Stefan

    Reply Like
Like Follow
  • 6 mths agoLast active
  • 4Replies
  • 17Views
  • 2 Following