0

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

6replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • 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----*----8

     

    To 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

    Like
  • 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

    Like
  • ...don't forget to add gravity if you want to see different values in each layer...
     

    Like
  • 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

    Like
  • Gregori,

    All I wanted to remind you of is to set variable GF in block PARAM.2 to 9.81.

    Stefan

    Like
  • Hello Stefan,

    Yes, of course, thank you!!

    Cheers,

    Gregori

    Like
Like Follow
  • 7 yrs agoLast active
  • 6Replies
  • 1018Views
  • 2 Following