Keurfon Luu's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: Output the spatial results for specific layer
@Refaat G Hashish Here is a script to over a line in a more general fashion. First, we create a fake mesh from the cloud points and import the output data:…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Output the spatial results for specific layer
@Refaat G Hashish My bad, it should be: import matplotlib.pyplot as plt
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Output the spatial results for specific layer
@Refaat G Hashish The function `toughio.read_output` returns a list of all time steps in your file. In Python, `[-1]` refers to the last element in the list ([0] being the first element,…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Output the spatial results for specific layer
You can use a Python module of mine (https://github.com/keurfonluu/toughio). To install, you need a Python distribution (e.g. Anaconda) and then run in Windows console:…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: How to identify the cause of an incomplete simulation run
@Refaat G Hashish You have reached the maximum number of cycles (9999). Rename the output SAVE file as INCON, remove the INCON block in your input file,…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3 on PC: CPU timing routine
I installed TOUGH3 on a fresh WSL setup (Ubuntu 20.04) and it installed smoothly without any modification (with the latest files from the repository). I just installed the following packages…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Commingled injection well modeling
@Refaat G Hashish As I said, you can explicity list all these elements instead. GENER----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8 A1 1inj 1 COM3 1.00 A2 1inj 2 COM3 1.…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Commingled injection well modeling
It's correct, depending on what you are trying to do. In your case, elements "A1 1", "A1 2", "A1 3", "A1 4" and "A1 5" will inject COM3 at a rate of 1 kg/s each.…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Try this command: zone create brick point 0 0 0 0 point 1 7403 0 0 point 2 0 5471 0 point 3 0 0 600 size 37 28 6
- Keurfon_Luu
- 4 yrs ago
- 1
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan From my understanding, the four points must form a right-handed coordinate basis, which means that a unit brick with origin (0, 0, 0) can be defined by this set of points:…
- Keurfon_Luu
- 4 yrs ago
- 1
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Your command is wrong: zone create brick size (37,28,6) point 0 (1890,2300,6) point 1 (7400,2300,6) point 2 (7400,5471,6) point 3 (1890,5471,…
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan There are 6156 elements in your grid, the CONNE block starts at line 6159 of the MESH file.
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan You are not using the latest release of FLAC3D7. According to the .f3grid file you provided, you are using FLAC3D7 R091 while the latest is R118.
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan It does not return any variable. It creates a MESH file in your current working directory.
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan I got this file with your command: zone create brick size 38,27,6 I am using a more recent version of FLAC3D, though. Can you update to the latest version (https://www.itascacg.…
Attachment
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Can you provide a screenshot of the grid in FLAC3D? The .f3grid file you provided has all the z-components of the grid points equal to 6 (you can open it with notepad).
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Your mesh is flat, is it expected? TOUGH is a integral finite-volume simulator, so you need a 3D mesh (FLAC3D too, as far as I know). Edit:…
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Can you provide the export .f3grid file for debugging purpose?
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan No, it's the .f3grid file that you can export from FLAC3D (File > Grid > Export to FLAC3D). It will export all the information from your grid as a text file (by default).
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan Yes, if you have a Python distribution installed (the most popular being Anaconda), you can simply run in Windows console the following command line to install the package:…
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: Tough2-FLAC 3D
@Floriane Youzan As stated by Christine, there is no way to convert a TOUGH mesh to a FLAC3D grid since we are lacking the cell connectivity data required to reconstruct the original cells (TOUGH…
- Keurfon_Luu
- 4 yrs ago
- Questions & Answers
- Reported - view
Re: How to interpret the data in CO2TAB file
Note that rcc3 sample problem is running isothermally. Don't forget to modify the second parameter NEQ in block MULTI to disable the isothermal option:…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Problem with running a sample case using TOUGH3/ECO2N
@Refaat G Hashish The input file for rcc3 sample problem is using AZTEC then switching to PETSc solvers which are not available on Windows platform. You have to edit your input file,…
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Problem with running a sample case using TOUGH3/ECO2N
I am not sure if it will solve your problem, but as far as I know, the command line to run a TOUGH3 simulation should be something like: tough3-eco2n INFILE OUTPUT.out
- Keurfon_Luu
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
@Christine Doughty Also, for Dirichlet boundary conditions, nodal distances should be set to a really small value. Here is a Python script to modify your input file with Christine's suggestions:…
Attachment
- Keurfon_Luu
- 4 yrs ago
- 1
- TOUGH2
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website