[Nek5000-users] Using mshconvert With Curves

Thu Aug 10 12:10:51 CDT 2017

Dear Nicolas,

I think it will be helpful if you email your 2D msh format.

Thanks,
Marco
On Aug 10, 2017, at 1:03 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:

Dear Marco,

I have actually very little experience with the script. Some people used it here (at KTH) for 2D and 3D wings and we never encountered such an issue.
If you cannot get it to work for simple cases and if that helps, I can send you a simple 2D msh file that works. From there, you could compare the files format... I am sorry that I cannot help much more than that.

Best,
Nicolas

On Thu, Aug 10, 2017 at 5:44 PM, <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
Dear Nicolas,

thank you for your reply. We tested the python script on different meshes to get familiar with it before diving into our working mesh. We can generate a rea file with velocity and temperature boundary conditions fine. However, whenever we try to use the curved side option we get the following error message:

Traceback (most recent call last):
  File "makemesh.py", line 6, in <module>
    convert('Single-Pitch-Mesh.msh',curves={13:{'type': 'm'},14:{'type': 'm'},15:{'type': 'm'}},bcs={13:'W',14:'W',15:'W',16:'O',17:'v'})
  File "/Users/v9g/Documents/Nek5000/mshconvert.py", line 1525, in convert
    read_curved_sides(curves)
  File "/Users/v9g/Documents/Nek5000/mshconvert.py", line 751, in read_curved_sides
    cell_face_m.append((cell, face_map[cell][(nds[0], nds[1])], face_number))
KeyError: (20355, 20211)

I am not sure how to proceed next, besides digging in the source code and try to resolve this error. We are currently trying to debug on a simpler geometry with curved boundaries, but we would really appreciate any help.

Thanks,

Marco

On Aug 10, 2017, at 10:41 AM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:

Dear Tommy,

Besides the short REAMDE file, there is a detailed explanation on how to use the script at the end of the file 'mshconvert.py' (starting from line 1385).
Basically, you need to identify the IDs of the curved boundaries in the .msh file and give it a Nek type (either midpoint 'm' or curve 'C' + some info: I think that only the radius is required, the rest is optional) and pass this info to the script in the form of a Python dictionary called 'curves'.
If the problem persists, feel free to share some details about the error that you get.

Best,
Nicolas

On Wed, Aug 9, 2017 at 10:43 PM, <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>> wrote:

Greetings,

I am trying to use the mshconvert script to produce a rea file for a helical geometry. I am running into errors when I try to define the boundaries that have curved surfaces. Has anyone used it to generate a rea file with curved data?

Thanks,

Tommy

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170810/86d26d49/attachment.html>