Hello guys . Does anybody know how to simulate the deformation of the soil caused by fluidflow in both TOUGH and FLAC3D? Thanks

Floriane Youzan You shouldn't save the data to FLAC3D file, you should import pressure and temperature data arrays from TOUGH in FLAC3D's Python, and then use itasca Python's module built in FLAC3D to import these arrays directly into FLAC3D's zones.

Assuming that your grid is loaded, the following code will import pressure and temperature calculated by TOUGH to FLAC3D's zones:

import toughio
import itasca as it

# Load SAVE file
save = toughio.read_output("SAVE")

# Loop over FLAC3D's zones and set pressure and temperature
for i, zone in enumerate(it.zone.list()):
    zone.set_pp(save.data["X1"][i])
    zone.set_temp(save.data["X6"][i])

You must change the primary variable names depending on the EOS you are using. From your SAVE file, I guessed that you are using EOS7R so pressure corresponds to variable X1 and temperature to variable X6.

Floriane Youzan 

import itasca as it

it.command("zone fluid zone-based-pp on")

If you are using module thermal, also add:

it.command("zone thermal zone-based-temp on")

Floriane Youzan To put it briefly, you can use 3 different languages in FLAC3D: FLAC3D commands, FISH and Python. Obviously, Python only accepts Python instructions.

Since

model configure fluid

is a FLAC3D command, you cannot use it as-is in a Python script. However, FLAC3D's commands can be used within a Python script using the function "it.command":

import itasca as it

it.command("model configure fluid")

Floriane Youzan Can you share your script?

Floriane Youzan What do you get with?

it.zone.count()

and

save.data["X1"].size

If your TOUGH MESH file has been created with toughio from your FLAC3D grid, you should have the same number of zones in both.

Have you imported your FLAC3D grid before running the script?

Floriane Youzan Looks like you imported your grid in FLAC3D thrice (18648 = 3 x 6216).

Run the FLAC3D command "model new" in FLAC3D's console (or "it.command("model new")" in FLAC3D's iPython's console), then reimport your grid once, and run your script again.

Note that at this current stage, you have merely imported the steady state results from TOUGH into FLAC3D. You should now loop over FLAC3D time steps (see my early answer).

Floriane Youzan What do you mean? The sample code so far only imports pore pressure calculated by TOUGH in FLAC3D. As I said earlier, you should now loop over time steps and sequentially couples TOUGH and FLAC3D:

t = 0.0
dt = 3600.0  # Time step, let's say one hour time step
tmax = 24.0 * 3600.0  # Simulation time, let's say one day
while t < tmax:
    # Code to solve fluid flow with TOUGH (e.g., using subprocess to call TOUGH)
    subprocess.call(...)

    # Import pore pressure in FLAC3D
    save = toughio.read_output("SAVE")
    for i, zone in enumerate(it.zone.list()):
        zone.set_pp(save.data["X1"][i])

    # Solve mechanical state in FLAC3D
    it.command(f"model solve mechanical time {dt}")

    # Prepare next time step for TOUGH
    # Modify input file using toughio
    # Rename SAVE file to INCON and modify porosity and/or permeability for each zone
    # (see Rutqvist et al., 2002 for equations)
    ...

    # Increment simulation time
    t += dt

Floriane Youzan No response, no problem. It means that the code ran successfully.

Now, FLAC3D can use pore pressure from TOUGH to update the mechanical state.