Defining multiple time dependent sources in GENER block

Dear all,

            I am attempting a methane migration problem in TOUGH3 EOS7CA, and I need to define multiple time dependent sources in the GENER block.

            From the descriptions given in the manual, it was understood that we can input the number of additional elements (as similar sources) and the increment in the code name of the element through NSEQ and NADD.

            For the meshing done, the default element names for the intended sources are AC1 1, AD1 1, AE1 1...... AL1 1. Since the difference (/ increment) here happens in the EL format and not in the NE, how can I define the NADD for the elements mentioned above?

Can someone help me in this regard?

13replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • Very quick answer:

    You cannot automatically generate element/source-code names in the EL portion of the element name; it only works for the integer section (i.e., NE) of the element name. Just copy the sources, or regenerate the mesh (e.g., rotated by 90 degrees), or (depending on the problem), inject into a single dummy element that is connected to all the source element, or ...

    Reply Like 1
    • Stefan Finsterle 

      Thank you Stefan..!! You helped me once again..!

      Just out of curiosity, why am I not being able to just list multiple time-dependent sources in the GENER block?

      Something like the one below...
       

      GENER----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
      AC1 1CH4                     5     COM3   6.41E-101.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AD1 1CH4                           COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AE1 1CH4                           COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AF1 1CH4                           COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      
      
      Reply Like
    • Stefan Finsterle

      Sir,

      Kindly ignore the errors in the format of the above excerpt.

      Just wanted to know about listing the time dependent sources only.

      Like

      GENER----*----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
      AC1 1CH4
                                  Time
                                  Generation rates
      AD1 1CH4
                                  Time
                                  Generation rates
      AE1 1CH4
                                  Time
                                  Generation rates

      Reply Like
    • Parameswaran T G 

      Currently, you're labeling all of your sources "CH4."  You should given each one its own label.  Something like this...

      GENER----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
      AC1 1CH4C                    5     COM3   6.41E-101.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AD1 1CH4D                          COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AE1 1CH4E                          COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      AF1 1CH4F                          COM3           1.2520E+04
                   0       4320000       8640000      12960000      17280000
            6.41E-10      2.56E-09      7.69E-09      8.98E-09      8.34E-09
      

      I'm not entirely sure this will address your issue, but give it a shot!

      Reply Like 1
    • Mikey Hannon 

      Thank you Sir.. I'll definitely try this out. Thanks a lot..!! 

      Reply Like
      • Stefan Finsterle
      • Finsterle GeoConsulting
      • Stefan_Finsterle
      • 7 mths ago
      • 1
      • Reported - view

      Parameswaran T G 

       

      Moreover, the format of source code names should be AAAII (see manual p. 174), so use something like:

      CH4 1

      CH4 2

      ...

      instead of:

      CH4A

      CH4B

      ...

      Stefan

      Reply Like 1
    • Whoops!  Thanks Stefan Finsterle!

      Reply Like
    • Stefan Finsterle 

      Thank you, Sir... Thanks a lott..!!!

      Reply Like
  • Hi All and  Mikey Hannon

    What will happen if same source/sink name is given? Will it gets override with latest entry in GENER? Or will it consider all entries of same source/sink name?

    Thanks,

    Amol

    Reply Like
  • Stefan Finsterle

    Is it not possible in TOUGH3 EOS7CA to enter multiple elements of same constant methane generation rate with NSEQ and NADD?

    Something like

    GENER----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
    A31 1CH4     49    1               COM3  1.160E-041.2520E+04

    Reply Like
      • Stefan Finsterle
      • Finsterle GeoConsulting
      • Stefan_Finsterle
      • yesterday
      • Reported - view

      Parameswaran T G 

       

      Please check the manual and set NADS. Then run again and inform your colleagues on the Forum if it works, and in case it does not work, please document your findings so everybody can look into it.

      Thanks,

      Stefan

      Reply Like
  • Already tried that.

    Gave values for NADS, tried without specifying, kept LTAB 0/1, Tried without entering any of them, etc.

    Every time, this is what is shown in the command prompt-

    “forrt1: severe (408): fort: (2) Subscript #1 of the array F1 has value 1 which is greater than the upper bound of -1.”

    The funny part is, this multiple source entering runs smoothly when I enter time-dependent methane generation rates.

    Reply Like
  • Thanks for revealing the error message (this always helps narrowing down the issue)! I'm stumped by the "funny part". Something goes wrong with dynamic array allocation (which should not be related at all to which component is later used for injection). Please attach your input file, so somebody from the TOUGH3 support team can have a look at it.

    Stefan

    Reply Like
Like Follow
  • 19 hrs agoLast active
  • 13Replies
  • 64Views
  • 4 Following