How to implement anisotropic permeability tensors and unstructured mesh in TOUGH3?
Hi everyone,
I'm working on a model involving complex geological structures with anisotropic flow behavior, and I'm planning to use TOUGH3 because of its support for full permeability tensors and unstructured grids.
" I used to work with TOUGH2 the last 6 years."
Could someone please clarify:
How to define a non-diagonal permeability tensor in TOUGH3 input files?
What is the recommended mesh format for importing unstructured meshes (e.g., from Gmsh or Cubit)? i'm using TOUGHIO is there another way (tools voroni free)
Are there any specific guidelines or examples available for calculating fluxes across arbitrarily oriented faces?
Does the orientation of the permeability tensor need to be expressed in local or global coordinates?
Any example files, tutorials, or documentation references would be highly appreciated.
Thanks in advance!
1 reply
-
1. you may take a look of the function named permeability_tensor in utils.f90 if you can access the source codes. This function can be evoked by having MOP2(19)=1. It will project the input permeability to the direction of connection. This function has not been included in the user manual and I believe it has not been fully tested.
2. Any format is good if it can be transformed to TOUGH own formatted.
3. TOUGH will calculate flux of any connections between two elements. If you design your mesh coordinately, you can get fluxes across any interface.
4. I guess global coordinates.