******************* This email originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list https://spam.ic.ac.uk/SpamConsole/Senders.aspx to disable email stamping for this address. ******************* Hi Everyone What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far. I appreciate your kind help. Kind regards syavash
Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient. You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well. Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ________________________________ From: nektar-users-bounces@imperial.ac.uk <nektar-users-bounces@imperial.ac.uk> on behalf of Ehsan Asgari <eh.asgari@gmail.com> Sent: Wednesday, March 22, 2023 3:54 PM To: nektar-users <nektar-users@imperial.ac.uk> Subject: [Nektar-users] High-order mesh This email from eh.asgari@gmail.com originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list<https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address. Hi Everyone What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far. I appreciate your kind help. Kind regards syavash
Hi Kaloyan Thank you for your kind reply. I guess just sent my last email before receiving your email. I will try your solution. But I think I will have to recompile at least part of the code, right? Kind regards syavash On Wed, Mar 22, 2023, 19:31 Kirilov, Kaloyan S < kaloyan.kirilov19@imperial.ac.uk> wrote:
Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient.
You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well.
Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ------------------------------ *From:* nektar-users-bounces@imperial.ac.uk < nektar-users-bounces@imperial.ac.uk> on behalf of Ehsan Asgari < eh.asgari@gmail.com> *Sent:* Wednesday, March 22, 2023 3:54 PM *To:* nektar-users <nektar-users@imperial.ac.uk> *Subject:* [Nektar-users] High-order mesh
This email from eh.asgari@gmail.com originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list <https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address.
Hi Everyone
What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far.
I appreciate your kind help.
Kind regards syavash
Hi Syavash, No worries at all. Yes, you will need to recompile. ( it should take just seconds because there are no dependencies on the module ) Kind Regards, Kaloyan ________________________________ From: Ehsan Asgari <eh.asgari@gmail.com> Sent: Wednesday, March 22, 2023 4:05 PM To: Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk>; nektar-users <nektar-users@imperial.ac.uk> Subject: Re: [Nektar-users] High-order mesh Hi Kaloyan Thank you for your kind reply. I guess just sent my last email before receiving your email. I will try your solution. But I think I will have to recompile at least part of the code, right? Kind regards syavash On Wed, Mar 22, 2023, 19:31 Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk<mailto:kaloyan.kirilov19@imperial.ac.uk>> wrote: Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient. You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well. Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ________________________________ From: nektar-users-bounces@imperial.ac.uk<mailto:nektar-users-bounces@imperial.ac.uk> <nektar-users-bounces@imperial.ac.uk<mailto:nektar-users-bounces@imperial.ac.uk>> on behalf of Ehsan Asgari <eh.asgari@gmail.com<mailto:eh.asgari@gmail.com>> Sent: Wednesday, March 22, 2023 3:54 PM To: nektar-users <nektar-users@imperial.ac.uk<mailto:nektar-users@imperial.ac.uk>> Subject: [Nektar-users] High-order mesh This email from eh.asgari@gmail.com<mailto:eh.asgari@gmail.com> originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list<https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address. Hi Everyone What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far. I appreciate your kind help. Kind regards syavash
Hi Kaloyan I modified the code and recompiled it. Now I only get a bunch of warnings like: [ProcessProjectCAD] WARNING: surface 0.000537058 unknown (tolerance: 0.0005) When comparing the initial mesh with the optimized one, I noticed no cell was added. Only some minor changes in some regions (Highlighted in the attached image): [image: Untitled.jpg] I used order=4 for this geometry. How can I decrease the tolerance and improve the performance? Kind regards syavash On Wed, Mar 22, 2023 at 7:40 PM Kirilov, Kaloyan S < kaloyan.kirilov19@imperial.ac.uk> wrote:
Hi Syavash, No worries at all. Yes, you will need to recompile. ( it should take just seconds because there are no dependencies on the module ) Kind Regards, Kaloyan ------------------------------ *From:* Ehsan Asgari <eh.asgari@gmail.com> *Sent:* Wednesday, March 22, 2023 4:05 PM *To:* Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk>; nektar-users <nektar-users@imperial.ac.uk> *Subject:* Re: [Nektar-users] High-order mesh
Hi Kaloyan
Thank you for your kind reply. I guess just sent my last email before receiving your email. I will try your solution. But I think I will have to recompile at least part of the code, right?
Kind regards syavash
On Wed, Mar 22, 2023, 19:31 Kirilov, Kaloyan S < kaloyan.kirilov19@imperial.ac.uk> wrote:
Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient.
You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well.
Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ------------------------------ *From:* nektar-users-bounces@imperial.ac.uk < nektar-users-bounces@imperial.ac.uk> on behalf of Ehsan Asgari < eh.asgari@gmail.com> *Sent:* Wednesday, March 22, 2023 3:54 PM *To:* nektar-users <nektar-users@imperial.ac.uk> *Subject:* [Nektar-users] High-order mesh
This email from eh.asgari@gmail.com originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list <https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address.
Hi Everyone
What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far.
I appreciate your kind help.
Kind regards syavash
Dear All, I could manage to run "projectcad" on a 3-D hex mesh without significant issues. The resulting mesh is a smooth one, however, the link below has made me believe that additional cells should be added to the optimized mesh. In my case, I only get a smooth mesh without additional cells representing a high-order grid (Based on the mesh examination in tecplot). https://www.nektar.info/mesh-conversion-from-gridgen-pointwise-to-nektar/ Noteworthy, I tried to run the examples provided in the link above, but I am getting some errors. I am confused here. How can I have a high-order mesh with additional cells? I played around with the "order" parameter but the result is the same. I would appreciate your kind help. Kind regards syavash On Thu, Mar 23, 2023 at 11:31 AM Ehsan Asgari <eh.asgari@gmail.com> wrote:
Hi Kaloyan
I modified the code and recompiled it. Now I only get a bunch of warnings like:
[ProcessProjectCAD] WARNING: surface 0.000537058 unknown (tolerance:
0.0005)
When comparing the initial mesh with the optimized one, I noticed no cell was added. Only some minor changes in some regions (Highlighted in the attached image):
[image: Untitled.jpg]
I used order=4 for this geometry. How can I decrease the tolerance and improve the performance?
Kind regards syavash
On Wed, Mar 22, 2023 at 7:40 PM Kirilov, Kaloyan S < kaloyan.kirilov19@imperial.ac.uk> wrote:
Hi Syavash, No worries at all. Yes, you will need to recompile. ( it should take just seconds because there are no dependencies on the module ) Kind Regards, Kaloyan ------------------------------ *From:* Ehsan Asgari <eh.asgari@gmail.com> *Sent:* Wednesday, March 22, 2023 4:05 PM *To:* Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk>; nektar-users <nektar-users@imperial.ac.uk> *Subject:* Re: [Nektar-users] High-order mesh
Hi Kaloyan
Thank you for your kind reply. I guess just sent my last email before receiving your email. I will try your solution. But I think I will have to recompile at least part of the code, right?
Kind regards syavash
On Wed, Mar 22, 2023, 19:31 Kirilov, Kaloyan S < kaloyan.kirilov19@imperial.ac.uk> wrote:
Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient.
You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well.
Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ------------------------------ *From:* nektar-users-bounces@imperial.ac.uk < nektar-users-bounces@imperial.ac.uk> on behalf of Ehsan Asgari < eh.asgari@gmail.com> *Sent:* Wednesday, March 22, 2023 3:54 PM *To:* nektar-users <nektar-users@imperial.ac.uk> *Subject:* [Nektar-users] High-order mesh
This email from eh.asgari@gmail.com originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list <https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address.
Hi Everyone
What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far.
I appreciate your kind help.
Kind regards syavash
Hi Syavash, I am not exactly sure what you mean by additional cells. What the Master version of the ProjectCAD module does is projecting the Edge Quadrature points to the CAD. It is not creating any new elements/cells. If you see a smooth ( curved ) boundary in Teckplot, then the module has worked as intended. You can also send me the xml and I will check it later. If you want to split the element around the surface for Boundary Layer y+ control, you need to call: -m bl:surf=x,y,z:nq=5 , nq is the number of points on the edge (p_order +1). I hope that helps. Kind Regards, Kaloyan PS: Keep in mind the version in the Master works for 3D elements, not 2D. The link that you shared has the 2D version of ProjecrCAD and needs to be build separately. ________________________________ From: Ehsan Asgari <eh.asgari@gmail.com> Sent: Saturday, March 25, 2023 8:35:21 AM To: Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk>; nektar-users <nektar-users@imperial.ac.uk> Subject: Re: [Nektar-users] High-order mesh Dear All, I could manage to run "projectcad" on a 3-D hex mesh without significant issues. The resulting mesh is a smooth one, however, the link below has made me believe that additional cells should be added to the optimized mesh. In my case, I only get a smooth mesh without additional cells representing a high-order grid (Based on the mesh examination in tecplot). https://www.nektar.info/mesh-conversion-from-gridgen-pointwise-to-nektar/ Noteworthy, I tried to run the examples provided in the link above, but I am getting some errors. I am confused here. How can I have a high-order mesh with additional cells? I played around with the "order" parameter but the result is the same. I would appreciate your kind help. Kind regards syavash On Thu, Mar 23, 2023 at 11:31 AM Ehsan Asgari <eh.asgari@gmail.com<mailto:eh.asgari@gmail.com>> wrote: Hi Kaloyan I modified the code and recompiled it. Now I only get a bunch of warnings like: [ProcessProjectCAD] WARNING: surface 0.000537058 unknown (tolerance: 0.0005) When comparing the initial mesh with the optimized one, I noticed no cell was added. Only some minor changes in some regions (Highlighted in the attached image): [Untitled.jpg] I used order=4 for this geometry. How can I decrease the tolerance and improve the performance? Kind regards syavash On Wed, Mar 22, 2023 at 7:40 PM Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk<mailto:kaloyan.kirilov19@imperial.ac.uk>> wrote: Hi Syavash, No worries at all. Yes, you will need to recompile. ( it should take just seconds because there are no dependencies on the module ) Kind Regards, Kaloyan ________________________________ From: Ehsan Asgari <eh.asgari@gmail.com<mailto:eh.asgari@gmail.com>> Sent: Wednesday, March 22, 2023 4:05 PM To: Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk<mailto:kaloyan.kirilov19@imperial.ac.uk>>; nektar-users <nektar-users@imperial.ac.uk<mailto:nektar-users@imperial.ac.uk>> Subject: Re: [Nektar-users] High-order mesh Hi Kaloyan Thank you for your kind reply. I guess just sent my last email before receiving your email. I will try your solution. But I think I will have to recompile at least part of the code, right? Kind regards syavash On Wed, Mar 22, 2023, 19:31 Kirilov, Kaloyan S <kaloyan.kirilov19@imperial.ac.uk<mailto:kaloyan.kirilov19@imperial.ac.uk>> wrote: Hi Syavash, As per my last email about ProjectCAD with Hexes - we usually use it with Prisms/Tets/Pyramids. In theory, nothing stops you from projecting the Hex Edge-Nodes on the CAD. I just had a quick look into the code and the only reason why it fails is that we haven't specified Hex function for ProcessProjectCAD.cpp / IsNotValid(). This is an element validity after the linear vertices are projected on the CAD. So If you use thick Hexes around the cylinder, this will be sufficient. You can comment out Line #226 ProcessProjectCAD (m_log(FATAL)<< "Only prisms and tetrahedra supported" << endl; ) function, rebuild and ProjectCAD should work for Hexes as well. Let me know how it goes. I will mark it as functionality that we need to add soon. I hope this helps. Kind Regards, Kaloyan ________________________________ From: nektar-users-bounces@imperial.ac.uk<mailto:nektar-users-bounces@imperial.ac.uk> <nektar-users-bounces@imperial.ac.uk<mailto:nektar-users-bounces@imperial.ac.uk>> on behalf of Ehsan Asgari <eh.asgari@gmail.com<mailto:eh.asgari@gmail.com>> Sent: Wednesday, March 22, 2023 3:54 PM To: nektar-users <nektar-users@imperial.ac.uk<mailto:nektar-users@imperial.ac.uk>> Subject: [Nektar-users] High-order mesh This email from eh.asgari@gmail.com<mailto:eh.asgari@gmail.com> originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list<https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address. Hi Everyone What is the best strategy to have a high-order mesh starting from a Hex mesh around a wall-mounted cylinder? So far, I have tried spherigon, high-order cylinder, and projectcad but I couldn't get a valid mesh. When I try projectcad, it gives me an error message that it only works with tets and prisms. The initial Hex mesh is created in Nastan format and then revised in Gmsh. Subsequently, NekMesh is used to convert it to xml. It is a valid mesh but building a high-order grid has been a challenge for me so far. I appreciate your kind help. Kind regards syavash
participants (2)
-
Ehsan Asgari
-
Kirilov, Kaloyan S