Please help with 2D version of 1D radial flow ECO2N example
Hello Everyone,
I am very new to TOUGHREACT, and in all honesty, I know I am in over my head, but I refuse to give up. I apologize in advance, although I find my problem very difficult for me it might be plainly obvious to you.
I am working with the ECO2N module, specifically the example 1D Radial flow problem for CO2 injection into a saline aquifer. I was able to run the example and plot some data. I was able to go through the manual and figure out what most of the values in the flow input file represent.
My difficulty, is that I want to make a 2D version of this problem. But I can't understand the user manual enough to see how to even setup the Meshmaker for a 2-D radial. I just don't know what to add to the flow.inp file. I don't even know if it's possible.
I have tried looking on the internet for examples, but I can't find anything that shows an example 2-D radial setup. If someone can point me to an example, or even better share how to change the 1D radial example into a 2D radial, then I can learn that way and get past this hurdle.
Thank you all for your time and patience.
Cheers,
Gregori
6 replies
-
Gregori,
(In case it makes you feel any better: TOUGHREACT is way above my head as well - but I may be able to help with MESHMAKER).
Assuming your current 1D radial model is horizontal (as in sample problem rcc3) and you want to make it 2D by having multiple layers, all you need to do is specify variable NLAY to the number of layers you want, and then specify the thicknesses of these layers. The following is the MESHMAKER block from rcc3, extended to include 4 instead of just one horizontal layer:
MESHMAKER1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
RZ2D
RADII
1
0.
EQUID
1 .3
LOGAR
200 1.E3
LOGAR
100 3.E3
LOGAR
100 1.E4
LOGAR
34 1.E5
LAYER----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
4
5. 15. 35. 35.
ENDFI----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8To be honest, making the mesh 2D is the simple step - the tricky part is to specify appropriate initial and boundary conditions, as you now need to account for gravity/buoyancy/capillary equilibrium effects.
Hope this helps,
Stefan
-
Hello Stefan,
Thank you for the quick response, I do feel better thank you. I thought it might have something to do with the layers, but I didn't make the connection. No one on my advisor's team has any experience with the software, so they all voted that the new student (me) become the "expert" on TR. This is really helpful, since what I am doing (for now) is proof of concept perhaps I can simply change the number of layers and see if I can actually see values for the different layers in the output file. I am sure I will have more questions as I progress. Thank you once again!
Cheers,
Gregori
-
...don't forget to add gravity if you want to see different values in each layer...
-
Hello Stefan,
Thanks for the reminder. I assume I would be adding this gravity data as a new layer type (e.g. shale) under the ROCKS section of the flow input file?
Thanks again,
Gregori
-
Gregori,
All I wanted to remind you of is to set variable GF in block PARAM.2 to 9.81.
Stefan
-
Hello Stefan,
Yes, of course, thank you!!
Cheers,
Gregori