Test Case 13 Species transport in fractrure matrix system with 4 straight fractures - GeoChemFoam/GeoChemFoam GitHub Wiki

In this example, we investigate species transport in a fracture-matrix system with 4 straight fractures. The geometry is a 5 cm x 5 cm solid grid with four fractures of aperture 2mm. The image is a 8-bit raw images where the fracture is represented by voxel of value 255, which corresponds to eps=1.0, and the matrix voxels of value 102, which corresponds to eps=0.4. The domain is meshed using a 2D uniform grid with resolution of 200 microns. A constant pressure drop is applied between the left and right boundary and transport a species is simulated

1. Navigate to the correct tutorial: cd $GCFOAM_TUTORIALS/basic/scalarTransportDBSFoam/4StraightFractures/ and copy it into 'runs' cp -r ../4StraightFractures/ $GCFOAM_RUNS/..

2. If you are using the Native version of GeoChemFoam, it is recommended that you run the test case from the runs directory cd $GCFOAM_RUNS/4StraightFractures/. If you are using Docker, it is recommended that you run the test case directly in its tutorial directory (for performance). cd $GCFOAM_TUTORIALS/basic/scalarTransportDBSFoam/4StraightFractures/ . You will then need to copy the case and the result to the 'runs' directory after each step for visualisation.

3. Run the ./createMesh.sh script to define the domain and mesh it. If you are using Docker, you need to copy the results in the runs folder to visualise them. cp -r ../4StraightFractures/* $GCFOAM_RUNS/4StraightFractures/

Fig. 1 Mesh for species transport in a fracture matrix system with 4 straight fractures.

4. Run the ./initCaseFlow.sh script to initialise the flow solver with pressure drop, fluid properties and kozeny carmen constant for the permeability of the matrix. If you are using Docker, you need to copy the results to the runs folder to visualise them. cp -r ../4StraightFractures/* $GCFOAM_RUNS/4StraightFractures/

5. Run the ./runCaseFlow.sh to run the flow solver to create the velocity field. If you are using Docker, you need to copy the results to the runs folder to visualise them. cp -r ../4StraightFractures/* $GCFOAM_RUNS/4StraightFractures/

Fig. 2 Velocity fields.

6. Run the ./processFlow.sh script to move the solution of simpleFoam to the '0' folder (ready for the heat transfer simulation) and calculate flow properties.

cat poroPerm.csv 
time poro perm Lpore Re UD
0 0.473187 2.86759e-09 0.000220185 0.00131745 2.83126e-06

7. Run the ./initCaseTransport.sh script to initialise the transport solver with diffusion coefficient and decompose if using multiple processors. If you are using Docker, you need to copy the results to the runs folder to visualise them. cp -r ../4StraightFractures/* $GCFOAM_RUNS/4StraightFractures/

8. Run the ./runCaseTransport.sh script to run the transport solver. If you are using Docker, you need to copy the results to the runs folder to visualise them. cp -r ../4StraightFractures/* $GCFOAM_RUNS/4StraightFractures/

Fig. 3 Transport of species

9. Run the ./processTransport.sh script to calculate concentration in fracture and matrix as well as the transfer coefficient

cat ST.csv 
time Tf Tm coeff
0 0 0 -3.84698e-14
500 0.0633775 0.081641 -2.68051e-07
1000 0.123305 0.136383 -5.92456e-06
1500 0.183223 0.183843 -0.000199237
2000 0.243084 0.228887 8.02512e-06
2500 0.301519 0.271617 3.3506e-06
3000 0.355653 0.309528 2.44031e-06
3500 0.405108 0.344432 2.16768e-06
4000 0.450638 0.37873 1.98659e-06
4500 0.492444 0.413048 1.82025e-06
5000 0.530193 0.44661 1.67008e-06
5500 0.563628 0.478331 1.54162e-06
6000 0.592891 0.507549 1.43551e-06
6500 0.618409 0.534076 1.34921e-06
7000 0.640683 0.557977 1.27905e-06
7500 0.660172 0.579413 1.22136e-06
8000 0.677277 0.598597 1.17308e-06
8500 0.692361 0.615775 1.1319e-06
9000 0.705746 0.631212 1.09619e-06
9500 0.717719 0.645166 1.06485e-06
10000 0.728524 0.657875 1.03709e-06
10500 0.738363 0.669543 1.01237e-06
11000 0.747399 0.680336 9.90247e-07
11500 0.755763 0.690391 9.70389e-07
12000 0.763557 0.699818 9.52517e-07
12500 0.770864 0.708704 9.36398e-07
13000 0.777751 0.71712 9.21834e-07
13500 0.78427 0.725123 9.08658e-07
14000 0.790465 0.732759 8.96725e-07
14500 0.796372 0.740067 8.8591e-07
15000 0.802021 0.747079 8.76106e-07
15500 0.807435 0.753821 8.67217e-07
16000 0.812636 0.760317 8.59162e-07
16500 0.817643 0.766585 8.51869e-07
17000 0.822469 0.772644 8.45272e-07
17500 0.82713 0.778506 8.39315e-07
18000 0.831636 0.784186 8.33948e-07
18500 0.835997 0.789695 8.29125e-07
19000 0.840224 0.795042 8.24804e-07
19500 0.844323 0.800237 8.20949e-07
20000 0.848302 0.805288 8.17525e-07
20500 0.852168 0.810202 8.14502e-07
21000 0.855926 0.814985 8.11852e-07
21500 0.859582 0.819643 8.09548e-07
22000 0.86314 0.824182 8.07567e-07
22500 0.866605 0.828607 8.05887e-07
23000 0.86998 0.832921 8.04488e-07
23500 0.873269 0.83713 8.0335e-07
24000 0.876476 0.841236 8.02457e-07
24500 0.879603 0.845243 8.01791e-07
25000 0.882653 0.849154 8.01339e-07
25500 0.885628 0.852972 8.01085e-07
26000 0.888532 0.8567 8.01017e-07
26500 0.891365 0.86034 8.01122e-07
27000 0.89413 0.863894 8.01389e-07
27500 0.896829 0.867365 8.01807e-07
28000 0.899464 0.870754 8.02366e-07
28500 0.902036 0.874064 8.03057e-07
29000 0.904547 0.877296 8.03871e-07
29500 0.906998 0.880452 8.048e-07
30000 0.90939 0.883534 8.05836e-07
30500 0.911726 0.886543 8.06972e-07
31000 0.914005 0.88948 8.08202e-07
31500 0.91623 0.892348 8.09518e-07
32000 0.918402 0.895148 8.10914e-07
32500 0.920521 0.89788 8.12387e-07
33000 0.922589 0.900546 8.13929e-07
33500 0.924607 0.903148 8.15536e-07
34000 0.926576 0.905687 8.17203e-07
34500 0.928496 0.908164 8.18927e-07
35000 0.93037 0.910581 8.20703e-07
35500 0.932197 0.912938 8.22527e-07
36000 0.93398 0.915236 8.24395e-07
36500 0.935718 0.917478 8.26304e-07
37000 0.937412 0.919664 8.28251e-07
37500 0.939064 0.921795 8.30233e-07
38000 0.940675 0.923872 8.32247e-07
38500 0.942245 0.925897 8.3429e-07
39000 0.943776 0.927871 8.3636e-07
39500 0.945267 0.929794 8.38454e-07
40000 0.946721 0.931668 8.40571e-07
40500 0.948137 0.933494 8.42707e-07
41000 0.949517 0.935273 8.44861e-07
41500 0.950861 0.937006 8.47031e-07
42000 0.952171 0.938694 8.49216e-07
42500 0.953447 0.940338 8.51413e-07
43000 0.954689 0.941939 8.53621e-07
43500 0.955899 0.943498 8.55838e-07
44000 0.957078 0.945016 8.58064e-07
44500 0.958225 0.946494 8.60296e-07
45000 0.959343 0.947933 8.62533e-07
45500 0.960431 0.949334 8.64775e-07
46000 0.96149 0.950698 8.6702e-07
46500 0.962521 0.952025 8.69266e-07
47000 0.963525 0.953317 8.71514e-07
47500 0.964502 0.954575 8.73762e-07
48000 0.965453 0.955798 8.76008e-07
48500 0.966378 0.956989 8.78253e-07
49000 0.967279 0.958148 8.80495e-07
49500 0.968156 0.959276 8.82734e-07
50000 0.969009 0.960373 8.84968e-07