February 2021 - Nektar-users

Roots of Jocobi or Legendre polynomial?
by Dongmyung Suh 07 Mar '21

07 Mar '21

******************* 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 All, I am a user of Nektar ++, and learning its foundation with Sherwin's book, Spectral hp element method for CFD [2ed]. I have a small question, while I study the book. When the book explains the Gauss-Lobatto-Legendre formula for integration (pg.61) and differentiation (pg 68), xi are the roots of Jocobi polynomial. However, in other references ( https://en.wikipedia.org/wiki/Gaussian_quadrature) they are the roots of Legendre polynomial. Which one is right? or both are correct? Please let me be clear about that. Thank you. Dongmyung Suh

2 2

setting-up bc's for Compressible Navier-Stokes solver
by wasim sarwar 07 Mar '21

07 Mar '21

******************* 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 All, Good Morning, and wishing all of you very good day! I am trying to solve the problem of convergent-divergent nozzle with pressure at inlet and outlet. Does anyone have the idea how to setup the pressure inlet and pressure outlet boundary conditions in compressible Navier-Stokes solver. Also, give an idea of how-to setup the symmetry bc about x-axis. The input file is also attached herein for reference. Mant thanks in advance Wasim Sarwar UPC, Barcelona

2 1

Turbulent channel flow test case
by Eli Frydman 07 Mar '21

07 Mar '21

******************* 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 all, I tried to run the turbulent channel flow test case from the user guide (session file TurbChFl_3DH1D.xml) but i didn't get any turbulent flow. I used a body force in the u direction as explained and added also a perturbation in the spanwise direction and didn't get any transition. Regards, Eli

3 2

3D case problem
by Alessandro Alberto de Lima 25 Feb '21

25 Feb '21

******************* 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. My name is Alessandro and I work with prof. Julio Meneghini at University of São Paulo in Brazil. I am trying to model a flow past a cylinder with the Quasi-3D and 3D approach. I have not succeeded yet in both ways. The 3D cases run only in serial? I am getting this message when try to run the case: Fatal : Level 0 assertion violation unable to determine face orientation What is wrong? Thanks. Regards.

2 1

Warning: Level 0 assertion violation Using Default variable expansion definition for field: u
by ztdepyahoo 25 Feb '21

25 Feb '21

2 1

Processing gmsh file with multiple element types
by Vidhi Zala 24 Feb '21

24 Feb '21

******************* 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 all, I am trying to use NekMesh to convert a 2D mesh made using gmsh to XML. I am successful in converting a msh file which has only 1 type of element (only quad or tri). However, when I try to convert a mesh that has quads and triangles together, it gives me the following error: Fatal : Level 0 assertion violation Where : nektar++/library/SpatialDomains/MeshGraph.cpp[610] Message : Unable to read composite index range: terminate called after throwing an instance of 'Nektar::ErrorUtil::NekError' what(): Level 0 assertion violation Where : nektar++/library/SpatialDomains/MeshGraph.cpp[610] Message : Unable to read composite index range: Aborted (core dumped) Attached files .geo and .msh are what I am working on. ,-~~-.___._. / | ' \ | |--------. Vidhi Zala ( ) 0 | | | School of Computing \_/-, ,----' | | | Scientific Computing and Imaging Institute ==== !_!-v---v--. University of Utah / \-'~; .--------. / __/~| ._-""|| | vidhi.zala(a)utah.edu =( _____|_|____||________| http://sci.utah.edu/people/vidhi.html I support Ruff Haven Crises Shelter for Homeless Pets <https://www.ruffhaven.org/donate>. Consider a donation today :)

1 1

Error when using compressible flow solver
by Octavio Morales 23 Feb '21

23 Feb '21

******************* 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, I am trying to use the compressible flow solver on a very coarse mesh so I can have an idea of how does it works, it is just a cube. I am attachng the xml that i am running with the solver and the screenshot of what appears once I run the command "CompressibleFlowSolver oct.xml". Besides that, I have another question. What is the maximum Mach number that the solver can handle in 3D? Is the compressible flow solver able to simulate in the supersonic regime? and, how can I get the drag and lift force or the shear stress in a hull using Nektar++? Thank you for your time! Best regards, Octavio A. Morales.

1 0

Arpack error
by Matt Duran 22 Feb '21

22 Feb '21

******************* 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. ******************* Hello, my name is Matt Duran and I am from the University of Central Florida. I have recently begun using Nektar++ 5.0.1 on Ubuntu and was going through the ‘Global Stability Analysis: Cylinder Flow’ tutorial. However, when I ran the ‘IncNavierStokesSolver’ with ‘Cylinder_Direct.xml’ I got the following error: [image: image.png] Arpack is already installed and the relevant sections of the .xml file is shown below. [image: image.png] I was wondering if I could get some assistance on why this error may be happening and how to fix it. Thanks in advance, Matt Duran

2 1

help setting up a Quasi-3D case.
by Alessandro Alberto de Lima 22 Feb '21

22 Feb '21

******************* 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. I am trying to run a simulation of a flow past a cylinder. I am using the mesh of the flow stability analysis tutorial and my results are completely two dimensional. No three dimensionalities are present. The xml is attached to this message. Can anyone check if I missed something or what I did wrong? Thanks. Alessandro

2 1

FieldConvert -m addfld (Not working in the recent version)
by Nikesh Yadav 21 Feb '21

21 Feb '21

******************* 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, I have updated my Nektar++ and FieldConvert -m addfld feature is not working. Segmentation fault (core dumped) keeps coming. Are you facing the same issue? -- *Thanks and regardNikesh, Ph.D.* *Assistant at Aerodynamics Dept, * *MEiLWarsaw University of Technology, **Warsaw, Poland* *Contact: - **+48 22 234 7523* *Email-id: **nikeshspartan(a)gmail.com <nikeshspartan(a)gmail.com> / * *nyadav(a)meil.pw.edu.pl <nyadav(a)meil.pw.edu.pl>*

2 1