question about TOUGH2 output
Hi, dear all
I'm working on geothermal model by using TOUGH2 with PetraSim.
Now I want to conduct with model calibration by using PEST.
To do so, I must output temperature from TOUGH2 model,
So I generated input data from PetraSim and modified as attached file (123.dat), and executed at TOUGH2 with output data as attached file (123.out).
Here is my question,because I would like to conduct with model calibration as I mention before.
But the output data (123.out) didn't show temperature results like TOUGH2 example as attached file (figure1),
I also modified with input data related with output control like PARAM(MOP..) and TIMES in input data(123.dat). But it didn't seem to work.
I'm wondering if someone can give me guidance can help me solve this problem.
Many Thanks
Yane
7 replies
-
Yane,
check out the following lines in your output file:
NUMBER OF ELEMENTS SPECIFIED IN DATA BLOCK "ELEME" EXCEEDS ALLOWABLE MAXIMUM OF 12000 INCREASE PARAMETER *MNEL* IN MAIN PROGRAM, AND RECOMPILE ----------------- SKIP FLOW SIMULATION ----------------- MESH HAS 0 ELEMENTS ( 0 ACTIVE) AND 129346 CONNECTIONS (INTERFACES) BETWEEN THEM GENER HAS 3846 SINKS/SOURCES
So you did not run a simulation there. Hope that helps.
Julia
-
Dear Yane Ye,
Please read the output file 123.out carefully, and you will find the error message that states that the number of elements you have is greater than the maximum number available in PetraSim.
I encourage you to go through all sample problems and develop small test problems on your own before attempting a model with 44,000 elements. Also note that iTOUGH2 would provide inverse modeling capabilities for TOUGH2 without the need to specify template and instructions files (of course, using PEST is a perfectly valid alternative).
Regards,
Stefan
-
Hi, JuDi and Stefan,
Thanks for replied.
I noticed this message in output file indeed,
By you mean, my TOUGH2 model didn't work successful.
That is because the number of elements in my model is greater than
maximum in PetraSim. But I do executed and read the temperature and others parameter successful at PetraSim .
So, dose there have another way to run a simulation on TOUGH2 successful under no change with number of elements?
Thanks replied again,
Yane
-
Yane,
the number of elements is actually limited in the EOS1.exe (or other EOS?.exe) that comes along with your Petrasim installation - unless you replaced that with your own compiled EOS?.exe.
Just guessing here: When you are reading temperature and other parameters in PetraSim make sure it is not the initial conditions you are setting. For help on Petrasim usage I also recommend the PetraSim forum: https://www.rockware.com/forum/index.php?showforum=5
Regards,
Julia
-
Yane,
The answer is: no, you cannot run the code without having it properly dimensioned.
Are you using PetraSim for pre- and postprocessing only, but run the simulation using regular TOUGH2?
Please contact PetraSim if you use PetraSim for the entire workflow.
If you have the TOUGH2 source code, simply redimension the major arrays in file T2.
Stefan
-
Hi, JuDi
I didn't compile the EOS1.exe file, so I have no ideal why the TOUGH2 simulation results is different with PetraSim.
And also I verified the PetraSim results again, it seem that there have no unusual at PetraSim results because the temperature distribution is different with initial conditions.
Hi, Stefan
At initial work, I used PetraSim for all the work including pre-, postprocessing and run simulation. The results looks normal and elements number also had no problem with simulation.
Because now I have to calibrate the parameter in my model and PetraSim have no calibration function as far as I knew.
So I chose the PEST as the calibration tool, because it can be coupled with TOUGH2 input file which can generate from my initial PetraSim file.
I actually have TOUGH2 source code, but I have no experience
to edit the source code. So, can you give me some guidance for edit source code if it's a simple work?
Many Thanks
Yane
-
Yane,
As indicated in my previous e-mail, simply edit file T2 and set MNEL and MNCON to the maximum number of elements and connections you need. Then recompile following the instructions in the read.me file.
Good luck!
Stefan