Dimension of geometric curvature to include for a 3D simulation
******************* 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. ******************* Dear Nektar users, I hope you're all doing fine, i'm using the incompressible navier-stokes solver, in fully 3D case with modified basis, and I have a question about the level at which adding curved geometry information, under <CURVED>....</CURVED> xml node, should be done. *Question:* is it necessary to add surface curvature information or should one be just fine with adding curved edges information only? Regarding this, I have noticed that exporting gmsh higher order mesh produces curved edges and curved faces under <CURVED>....</CURVED> section, which typically reads:
<E ID="45953" EDGEID="26270" NUMPOINTS="3" TYPE="PolyEvenlySpaced">2.54571429e-01 -7.50000000e-03 -6.42857143e-04 2.50785714e-01 -7.50000000e-03 -6.96428571e-04 2.47000000e-01 -7.50000000e-03 -7.50000000e-04 </E> <F ID="45954" FACEID="26260" NUMPOINTS="9" TYPE="PolyEvenlySpaced">-1.46714286e-01 7.50000000e-03 3.00000000e-02 -1.50500000e-01 7.50000000e-03 3.00000000e-02 -1.54285714e-01 7.50000000e-03 3.00000000e-02 -1.46714286e-01 1.12500000e-02 3.00000000e-02 -1.50500000e-01 1.12500000e-02 3.00000000e-02 -1.54285714e-01 1.12500000e-02 3.00000000e-02 -1.46714286e-01 1.50000000e-02 3.00000000e-02 -1.50500000e-01 1.50000000e-02 3.00000000e-02 -1.54285714e-01 1.50000000e-02 3.00000000e-02 </F>
Since I'm adding curvature information using an external python script, should I be worried about explicitly adding the face curvature as produced inside an auto-generated session file through NekMesh? A first naïve thought I had was that the Edges curvature information would propagate the curvature to higher dimensions through the tensor product !? However, I'm not even close to being formal on this statement. Thank you all in advance, Kind regards, Saad
******************* 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. ******************* Dear Nektar users and developers, I would really appreciate any feedback, thoughts or clarification on this question, i.e. for a fully 3d simulation using the incompressible N-S solver, should the curved section include face curvature or would it be fine to define curvature for edges only. Thank you all, best regards Saad
Hi Saad, Certainly just using edge curvature will get you most of the way to representing the geometry and in the past we have also just implemented this first. We do observe we can get some element patterning in shear stress if we do not also include face curvature. However you probably will not be able to see many issues in the velocity field. Hope that helps, Spencer. From: nektar-users-bounces@imperial.ac.uk <nektar-users-bounces@imperial.ac.uk> on behalf of Debbahi Saad <saad.debbahi@gmail.com> Date: Sunday, 3 November 2024 at 14:49 To: nektar-users <nektar-users@imperial.ac.uk> Subject: [Nektar-users] Dimension of geometric curvature to include for a 3D simulation This email from saad.debbahi@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. Dear Nektar users, I hope you're all doing fine, i'm using the incompressible navier-stokes solver, in fully 3D case with modified basis, and I have a question about the level at which adding curved geometry information, under <CURVED>....</CURVED> xml node, should be done. Question: is it necessary to add surface curvature information or should one be just fine with adding curved edges information only? Regarding this, I have noticed that exporting gmsh higher order mesh produces curved edges and curved faces under <CURVED>....</CURVED> section, which typically reads: <E ID="45953" EDGEID="26270" NUMPOINTS="3" TYPE="PolyEvenlySpaced">2.54571429e-01 -7.50000000e-03 -6.42857143e-04 2.50785714e-01 -7.50000000e-03 -6.96428571e-04 2.47000000e-01 -7.50000000e-03 -7.50000000e-04 </E> <F ID="45954" FACEID="26260" NUMPOINTS="9" TYPE="PolyEvenlySpaced">-1.46714286e-01 7.50000000e-03 3.00000000e-02 -1.50500000e-01 7.50000000e-03 3.00000000e-02 -1.54285714e-01 7.50000000e-03 3.00000000e-02 -1.46714286e-01 1.12500000e-02 3.00000000e-02 -1.50500000e-01 1.12500000e-02 3.00000000e-02 -1.54285714e-01 1.12500000e-02 3.00000000e-02 -1.46714286e-01 1.50000000e-02 3.00000000e-02 -1.50500000e-01 1.50000000e-02 3.00000000e-02 -1.54285714e-01 1.50000000e-02 3.00000000e-02 </F> Since I'm adding curvature information using an external python script, should I be worried about explicitly adding the face curvature as produced inside an auto-generated session file through NekMesh? A first naïve thought I had was that the Edges curvature information would propagate the curvature to higher dimensions through the tensor product !? However, I'm not even close to being formal on this statement. Thank you all in advance, Kind regards, Saad
participants (2)
-
Debbahi Saad
-
Sherwin, Spencer J