use ENDFI to create mesh before initializing simulation , but how??? tried several combinations.....Mesh file remains empty
hello everyone, I am using TOUGH+Hydrate v 1.5
I bought it from LBNL(executable only)
I think, that for mesh creation, ENDFI must be put in place of ENDCY
now running this, creates mesh file, others as well INCON ALLOC etc and skips the simulation .
but all of these are empty??😥😥😥😥😥😥😥😥
I tried several options but none of them produces MESH
I have attached the files, both (One containing mesh, other without mesh)
Please advise and have a look
So I don't run TOUGH+ simulations, but it looks like you have an ENDCY keyword entered before ENDFI in your 'Test_2Qp_without mesh' file. I'm pretty sure that TOUGH codes will ignore anything in your input files after that.
In your other file (Test_2Qp) you have the ELEME and CONNE blocks in there, too. I'm not sure what TOUGH's reaction will be to having those blocks active in addition to an active MESHM block. Just to be safe, I would deactivate ELEME and CONNE. The easiest way to do that is to make the first letter of each lower case (i.e., eLEME and cONNE). That'll get TOUGH to ignore them. Give those a shot and let me know if they work. If this ends up being a TOUGH+ specific problem, I won't be able to help you much.
For TOUGH+HYDRATE, you need to generate the meshes (ELEME and CONNE blocks) independently, using the separate MeshMaker code, and include them in the input file before the ENDCY tag (or as an external MESH file) . The ENDFI tag (and internal MeshMaker) has been deprecated.
For any hydrate case beyond the simplest horizontal 1D cases, you would need to perform separate equilibration runs to generate an INCON block that is in correct P-T-x equilibrium before any production simulation is run, anyway.
Matthew Reagan said:
perform separate equilibration runs to generate an INCON block that is in correct P-T-x equilibrium
How to do this ??
I would be grateful, if you can attach a sample for these, it would also benefit the community
This is covered in our multi-day training courses. The location of the bottom of the hydrate stability zone is a function of the (hydrostatic) P, T, and hydrate saturation profiles, and the T and S profiles are interrelated due to the thermal conductivity of the hydrate. Typically one performs a series of 1D column simulations until a stable steady-state condition is determined. This is beyond the level of what can be written in a forum post!