Zone gridpoint fix velocity-x range position-x 350000.The development and application of roadway group layout methods in coal mines have become more common and the mutual disturbance of blasting and driving of roadway groups has also become more prominent at depth.

Zone gridpoint fix velocity-x range position-x -350000.0 Zone gridpoint fix velocity-z range position-z -3000.1 Zone gridpoint initialize temperature 10.0 gradient 0 0 -0.025 Zone gridpoint initialize pore-pressure 1.e5 gradient 0 0 -9810.0 Zone initialize stress-zz 0.0 gradient 0 0 22170.6 Zone initialize stress-yy 0.0 gradient 0 0 22170.6 Zone initialize stress-xx 0.0 gradient 0 0 22170.6

Zone gridpoint initialize fluid-modulus 1.0 You can install toughio in FLAC3D by running this command in a Windows console: cmd /C ''C:/Program Files/Itasca/FLAC3D700/exe64/python36/Scripts/pip.exe'' install toughio -user Reading TOUGH input and output files can be done using toughio. # Interpolate my_pp and my_temp from SAVE file To update a zone pore pressure and temperature, I don't remember if there is a zone array for these properties, but you can loop over each zone and update the parameters as follow: import itasca as it For instance, you can call an external program in Python as follow: import subprocess Loop back to step 1 until end of simulation.įLAC3D has a Python interface which should make this a bit easier to implement compared to using FISH. Import pressure and temperature from SAVE file and update zone pore pressure and temperature in FLAC3D,ĥ. Solve fluid and heat flow in TOUGH to time step Dt,Ĥ. Calculate stress-induced changes in porosity and/or permeability and update TOUGH input files (MESH and/or INCON) with these new values,ģ. Solve the mechanical state in FLAC3D to a time step Dt,Ģ. However, the TOUGHFLAC implementation is more optimized as it does not require to read and write TOUGH input/output files at each time step.Īssuming that you have already setup your problem for TOUGH and FLAC3D independently, a quick algorithm would look like as below:ġ. It's even more complicated with TOUGH3 since you should now also consider that arrays are distributed in parallel. I am suggesting it since the current implementation of TOUGHFLAC (TOUGH calling FLAC) alters TOUGH source code a lot which requires a really good knowledge of the Fortran language. Floriane Youzan I didn't try the FLACTOUGH (FLAC calling TOUGH) implementation myself.

Groups flac3d