Non isothermal simulations using ECO2N module

I am trying to run some simulations using TOUGH2-MP and compare the results with TOUGH2.

The simulation runs successfully by TOUGH2 but when I tried to run it with TOUGh2-MP, I get convergence issues and simulation fails to run.

### no convergence in NSIN at * 1890   * for (T, XS, XLCO2) = (102.5102, 0.200000E+00, 0.100000E-05) P00=0.113029E+06
     latest parameters are   PS=-.778998E+08   X2S=0.170500E-01
 +++++++++   CANNOT FIND PARAMETERS AT ELEMENT * 1890   *          XX(M) = 0.267755E+08 0.200000E+00 0.990000E-06 0.102510E+03
 +++++++++ REDUCE TIME STEP AT (  149, 1) ++++++++ ++++++++   NEW DELT =0.273642E-04

But the same file works if I make it a isothermal case. Are there any issues related to non isothermal simulations using ECO2N module in TOUGH2-MP

2replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • Please provide full outputs from both runs. Otherwise, it is not possible to diagnose the problem.  For TOUGH2-MP, you may need to increase MOP(16) to a larger number, like 6 instead of 4. 

    Reply Like
  • Thank You George Pau

    By changing MOP(16) to 6 the test problem and uniform grid worked for non isothermal cases, up to temperatures 100°C.

    For non-uniform grid the simulations does not seem to converge even for temperatures less than 100°C and with MOP(16)=6.

    I have attached Non-uniform grid, isothermal and non isothermal input files and 

    RCC3 test problem at temperature greater than 100°C and for non-isothermal simulations.

    Could you please verify the attached non-isothermal cases. 

    Reply Like
Like Follow
  • 2 yrs agoLast active
  • 2Replies
  • 740Views
  • 2 Following