Exercise 2

Ok, time to put our new knowledge to work! Let’s explore a more realistic case and investigate how a more permeable fault zone can deviate hydrothermal upflow. You can find background on such a setup in [Andersen et al., 2015].

Steps

• check all settings and modify model run time so that the plume arrives at the surface

• modify (or create) a setFieldsDict to create the permeability structure.

• use a function called rotatedBoxToCell to set the permeability structure

• run the case and explore how the results change for different fault widths and permeability contrasts.

Step 1

Copy the $HOME/hydrothermalfoam-master/cookbooks/2d/Regular2DBox/ into your working directory.

cd $HOME/HydrothermalFoam_runs
cp -r /home/openfoam/hydrothermalfoam-master/cookbooks/2d/Regular2DBox/ ./Fault_zone_2D

Step 2

Now modify the case and setup the permeability structure. It should look like this:

Fig. 36 Permeability structure of our detachment fault example cause.

We can use the rotatedBoxToCell function in the setFieldsDict for this. Check the openfoam documentation !

Tip

It’s all about rotating the coordinate system to figure out the origin and i,j,k vectors for rotatedBoxToCell!

regions
(
    rotatedBoxToCell
    {
        origin   ( ? ? ?);
        i        ( ? ? ?);
        j        ( ? ?  ?);
        k        ( ? ?  ?);

        fieldValues
        (
            volScalarFieldValue permeability 1e-13
        );
    }
);

Step 3

Explore the results in paraview! And investigate how the results change for different fault widths and permeability contrasts.

Fig. 37 Hydrothermal flow along a preferential pathway.

Tip

• Run many simulations for different fault widths and permeabilities

• Explore how vent temperature changes

• …and when the plume is captured by the fault!