0

meshmaker error, increments problem

Hello everyone. 

i am a masters student trying to do simulation with tough+hydrate. i am trying to create a Mesh to use the ELEME NAME in my source and sink block but  i keep getting this error. The samples files run normally but when i use the increments i created from excel, the file refuse to run. 

the reservoir is suppose to be 200m wide in x dimention discretized into 90 grid blocks nonuniformly with 0.1 as the smallest size close to the well, the y dimension is just 1 and z direction have 326 block of 393.9m with the hydrate and gas zones discretized finely 0.1-0.5m and overbudden and underbudden with coarser blocks, top and bottom boundary of 180 cells with just 0.001m. the two wells are located at the first and last blocks in x dimension and from top of hydrate region to the bottom of methane gas region in z dimension. 

i made this increments but it refuse to run. someone pls help, find attached the meshmaker input 

6 replies

null
    • Staff Scientist
    • Christine_Doughty
    • 1 yr ago
    • Reported - view

    There are a bunch of lines at the start of your input file that TOUGH does not understand.  That is what the message "input conversion error" is telling you.  For MESHMAKER, the input file should start with one line that is the title line and contains whatever you want.  Then the next line has MESHM in columns 1-5.  This is a keyword that invokes MESHMAKER.  The next line has XYZ in columns 1-3, and the remainder of your input file should be okay.  EXCEPT, that you have tabs separating entries instead of spaces in some places.  TOUGH cannot read tabs.  I modified your input file to make it legal, and it runs for me.  See if it works for you.  I did not check that it is actually creating the grid you want - just that it runs.  Christine

      • Mahmud_muhammad_mahmud
      • 1 yr ago
      • Reported - view

      Christine Doughty  Thank you so much. I really appreciate it, let me try it. I am sure I will have other questions about source/sink or other subsequent blocks in the near future because I only have a month or two to finish up my thesis. It will be very helpful indeed if u could answer those as well. Thanks in advance 

    • Mahmud_muhammad_mahmud
    • 1 yr ago
    • Reported - view

    Christine Doughty

     

    The Meshmaker only assign well1 and skipped well2.

    Here, the first well perforation is at A2R00 (gridblock 1 in X at layer 96) and every other layer to  A6Q00 (gridblock 1 at layer 254).

    The reservoir domain consist of 90 (200m) grid block in the x dimension and 326 (393.9) in y dimension. The grid block 1 (0.1m) and 90 (0.1m) in x are injector and producer wells respectively, and the perfortaion starts from y dimension layer 96 (top of hydrate layer) to layer 254 (bottom of gas zone) as highlighted in red.

    A2R89 to A6Q89 were not assigned for well2

    In the block GENER, should i add these manually?

    i think in GENER i will have two ElName5C (A2R00 and A2R89) containing my sink/source, am not sure how to show where the perforation ends (A6Q00 and A6Q89).

    There is something about NSEQ and NADD in the manual, increment between the code numbers of two successive elements, pls check the attachment if i represent the well correctly

    • Staff Scientist
    • Christine_Doughty
    • 1 yr ago
    • Reported - view

    You are on the right track with your GENER block, but there are a couple of things to check.

    First, if you want to inject heat (units Watts, like with an underground heater), you use the HEAT keyword in GENER, but if you want to inject heated water, you use the keyword COM1, which means you are injecting water.  To make the water hot, you either specify the enthalpy (J/kg) in the next 10 columns after the mass injection rate (kg/s).  The alternative is to make the specific heat of the injection elements huge (so their temperature will never change), and set the desired injection temperature as an initial condition in INCON.  You do this by making the injection elements a distinct material in ROCKS, with a specific heat of something like 1.E30.

    Now, to use the NSEQ and NADD option, you need to see how your grid is ordered.  With the built-in MESHMAKER in TOUGH3, the elements are ordered by Z, so your NADD would be 1.  But, I think you must be using the MESHMAKER in TOUGHio-dash, which might order the grid differently, so check that adding 89 to the element order gets you to the next well element.  Also, I think NSEQ should be 158=254-96?  I actually haven't used the NSEQ/NADD option for a long time, so you should check the GENER output and make sure the code is doing what you want.

    Also, with the built-in MESHMAKER, you would have to specify your wells by hand.   If the TOUGHio-dash MESHMAKER assigns wells for you, you should check with @Keurfon_Luu if it is not working for your second well, or you could just assign it by hand, the way you have in your GENER block.  Usually a production well does not need to have a distinct material in ROCKS.

    I hope this helps.  For a problem like this it is usually a good idea to just run a few timesteps at first, and check the output carefully to make sure that the code is doing what you want.  For example, I would expect pressure to be increasing around your injection well and decreasing around your production well.  Another good thing to look at is the flow rates between the well elements and their neighbors, to make sure flow is going in the correct direction.  To include flow rates in the output file, you need to be sure KDATA>1 in the PARAM block.

      • Mahmud_muhammad_mahmud
      • 1 yr ago
      • Reported - view

      Christine Doughty Thank you very much.

      Actually i am using tough+ for the simulation. i was thinking the mesh creation process should be similar. i am using the meshmaker.f95 and it seems that when i use 1500m instead of 200m for the x domain and a smaller size well (0.05) instead of 0.1 the second well is assigned so i will stick to this for now.

      I appreciate the information on the GENER, i will make the correction and take note of the two different alternatives to heat the well. The manual did not specify COM1 and HEAT so i just assumed.

      Thanks

    • Staff Scientist
    • Christine_Doughty
    • 1 yr ago
    • Reported - view

    Oh, that's right - you mentioned TOUGH+ in your first post.  Since we were using the TOUGH3 forum category I forgot that.  I do not have much TOUGH+ experience.  It is supposed to be back-compatible with TOUGH3, but you should definitely be consulting the TOUGH+ documentation. I do know that the TOUGH+ MESHMAKER is more sophisticated, and I'm glad to hear its working for you now.   In the future, you will get more informed help if you post to the TOUGH+ forum category. 

Content aside

  • 1 yr agoLast active
  • 6Replies
  • 72Views
  • 2 Following