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 structurerun 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:
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.
Tip
Run many simulations for different fault widths and permeabilities
Explore how vent temperature changes
…and when the plume is captured by the fault!