-
Mesh Creation
Hi everyone. I recently started learning ToughReact. I’m trying to learn Mesh creation from scratch and how to create a 1D or 2D mesh. how to I go about this? What resources are available?
-
TOUGH+HYDRATES convergence issue in phase change
I want to calculate the depth interval at which gas hydrates forms, given P and T initial conditions, fixed on top. The idea is to run the model until it reaches the equilibrium,…
-
A problem about the direction of the anisotropic permeability
Hello everyone, I am using the EOS3 module to simulate the air and water flow in a hypothetical dome structure. I developed a 2-D radial model, and the PER in the ROCKS were set as 1e-12 (PER1),…
-
Error while running a T2LBM model
Does anyone know what is this error? forrtl: severe (157): Program Exception - access violation
-
MESHMAKER CORNERS
Good Morning All, I believe that a version of MESHMAKER (Meshmaker 2?) exists that writes a CORNERS attribute defining the nodes of an unstructured grid.…
-
run time error M6201:math
Hello everyone, I don't know where this error occurs and how to solve it, hopefully you can help me，thanks!
-
How to visualize TOUGH MESH file in pyvista?
Hello TOUGH3 Forum, I got the above files, and want to visualize it with some python packages, e.g., pyvista or pytecplot. Now I can only read "MESH" file with the following code through toughio:…
-
SMCO2 is not increased
Hello I am using Toughreact v3 for simulation for CO2 injection and our simulation runs well, however, CO2 is not mineralized after the simulation has finished.…
-
Empty SAVE file
Dear All, I have run the executable tough2 with the input is a DAT file (got it from Petrasim). When I run the simple model, it worked. But when I made a field model in Petrasim and got the DAT file,…
-
EOS3-sink/source rates are piecewise constant, how to set up?
As mentioned in the title, I want to set the piecewise constant sink/source rate as shown in the attached figure, but I can't achieve the desired result after setting it many times in petrasim,…
-
Last bock in CO2TAB file
Hello everyone: I have difficulties to understand the last block in CO2TAB file and search for help from here. I know the previous blocks are pressure, temperature, viscosity,…
-
TOUGHREACT/ECO2N:CANNOT FIND PARAMETERS AT ELEMENT * 1*
Hi Everyone, I'm new at using ECO2N of TOUGHREACT. And the following problem is encountered when I used PetraSim to simulate. It's a simple model. And I model a strata from -2260m to -2900m below…
-
Tough2-MP post-processing and visualization
Hi everybody, can anyone, please, advise on software to extract and visualize results from TOUGH2-MP with ECO2N? I realise that EXT is not an option to extract results to tecplot,…
-
CONVP block in solute.inp
Hello everyone: I have a question about implementing CONVP (chemical convergence bounds for time stepping) in solute.inp. I firstly make MOPR(10)=2 to make simulator read CONVP block.…
-
Why MESHMAKER not writing MATERIAL to MESH file: 2D Radial heterogenous Mesh: what is wrong?
hello moderator/everyone i adapted TSN_Radial example from MESHMAKER V1.5 example: but when i run it, it makes mesh file but it it it does not write material to mesh file.…
-
call TOUGH in jupyter notebook and PyTOUGH
Hi All, I am using PyTOUGH to automate TOUGH simulation. However, when I try to call TOUGH in Jupyter notebook with this command subprocess.call(["path_to_my_tough3_executable(tough3-eos1)",…
-
Run ends early
Hi, I was working on a larger block mesh using a doublet configuration for an aquifer, over 10 years. When I made the mesh finer, the simulation timed out and quit after 16 hours.…
-
Pre-Processing using Matlab
Hai everyone, I am new to using tough2 simulation. Do anyone has experience in using Matlab as pre-processing to run the tough2? I tried to follow the code from the internet which is…
-
zero pore pressure boundary condition
Can TOUGH3 define pore pressure Dirichlet boundary condition as zero values? In order to validate the coupled THM code, I have to compare the numerical solution to analytical solution.…
-
Ice in porous media T+H
Is it possible to model gas hydrate simulations in ice-bearing layers? Although there are different ice phases in TOUGH+HYDRATE, when we add ice into the system we got convergence failure problem.
-
TOUGH3-FLAC3D time step setting
Hi, I am trying to couple TOUGH3 and FLAC3Dv7.0 through python scripting engine built in FLAC3D software. I am confused about the time step control after reading some published papers (e.g.…
-
How to realize an air-injection well with constant bottom-hole pressure?
Dear all, I am using the EOS3 module to simulate the injection of air into a confined aquifer via a fully penetrating well. In order to prevent potential damage to the rock formation,…
-
MAX.RES. is small but DELTEX keeps decreasing
Hi All, I am running TOUGH3 EOS7C. The simulation seems stuck at extremely small time step. Please see below for the failure message. You can see that the residual is already very small but the time…
-
Is TOUGHREACT only suitable for NaCl salt?
Hello everyone, I am conducting an experiment using different salts and distilled water to prepare brine. I was wondering if I can simulate the experiment using TOUGHREACT as it appears Toughreact…
-
Inserting Permeability
Dear all, I am using Petrasim E0S5. I would like to know how to input different permeability and porosity for each grid block. My model is 3600 elements.…