- Nektar-users - imperial.ac.uk

WSS taking too much time
by Ehsan Asgari 20 May '24

20 May '24

******************* 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 trying to extract the wall shear stress for a 2.5-D case using the wss module: FieldConvert -m wss:bnd=3 <mesh.xml> <session.xml> <time-averaged.fld> <time-averaged-wss.fld> However, it is taking too much time to execute, more than 1 hour with no output. I noticed another report on this issue in the mail archive, not sure if my case has the same problem. I appreciate any help you can provide. syavash

1 0

NekMesh: Inconsistent triangular face orientation
by Alexander Schukmann 20 May '24

20 May '24

******************* 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 everybody, I tried converting a 3D tetra/prism mesh (generated with gmsh) of a simple cylinder from msh to xml using NekMesh. The domain has periodic boundaries along the cylinder axis. However, I receive the following error: Orientation of triangular face (id = 113069) is inconsistent with face 3 of prism element (id = 55392) since Dir2 is aligned with Dir1. Mesh setup needs investigation terminate called after throwing an instance of 'Nektar::ErrorUtil::NekError' I've seen a similar post from May 12th, where it was suggested to try the peralign option, so I used: NekMesh -m peralign:surf1=700:surf2=800:dir=z:orient CylReD3900.msh CylReD3900.xml with no success. The two surface IDs 700 and 800 of course refer to the periodic boundary faces. I attached my msh file and appreciate any help. All the best Alex Sicher versendet mit [Proton Mail](https://proton.me/).

2 1

ARPACK is not set-up for parallel execution
by 🌸那慢慢🌸 15 May '24

15 May '24

******************* 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 get the error "ARPACK is not set-up for parallel execution" when running the incompressbile flow solver in parallel with MPI，but in the file build/CMakeCache.tst have a line saying:NEKTAR_USE_MPI:BOOL=ON.How do I fix this? Hongtao Cheng

4 4

Acoustic solver Source term
by dauood Tariq 13 May '24

13 May '24

******************* 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 understand the source term provided in the user guide however the equation for the source term is very confusing and I cannot find any information on it. I would appreciate it. Regards, Dauood

2 1

(no subject)
by Wang Tongkuai 12 May '24

12 May '24

******************* 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. ******************* 从 Windows 版邮件发送 Dear All, I meshed a cylinder with hexas and prisms at gmsh with no error. Here is the problem occurred when I NekMesh it: Orientation of triangular face (id = 10638) is inconsistent with face 3 of prism element (id = 3588) since Dir2 is aligned with Dir1. Mesh setup needs investigation terminate called after throwing an instance of 'Nektar::ErrorUtil::NekError' I don’t know which surface is face 3 and how does the inconsistency happen. Dir 1 and Dir 2 are x- and y-axis direction, right? A half of y-z plane is shown. Regards, Tongkuai.

2 1

PowerKernel for not available along homogeneous directions?
by Alexander Schukmann 12 May '24

12 May '24

******************* 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 everybody, I am currently running Quasi-3D CG-SVV iLES with the IncNavierStokesSolver and noticed that the SVV PowerKernel is not activated/available for Fourier discretizations along homogeneous directions. However in the paper "Industry-relevant ILES of a high-performance road car via S/HPEM" (Mengaldo et al. 2021) it says: "The current best practice in Nektar has been to use the "DG" kernel in general while employin the "power-law" kenel only along homogeneous directions in case periodic boundary conditions are present" So does this mean the PowerKernel is only available if S/HP is also used along homogeneous directions? Why not use the DGKernel then? And how does one setup different SVV kernels for different spatial directions? If it DOES refer to the quasi-3D case, however, there seems to be an issue, since no output is given regarding SVV, if the "SpectralVanishingViscosityHomo1D" flag is set to "PowerKernel" within the session file. Also I noticed that it seems to be impossible to activate SVV only for the S/HP discretizations in quasi-3D cases, i.e. without using SVV for the Homo1D directions, but vice versa it works. All the best Alex Sicher versendet mit [Proton Mail](https://proton.me/).

2 4

Global vs local polynomial de-aliasing
by Alexander Schukmann 12 May '24

12 May '24

******************* 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 everybody, I have a quick question regarding the "SPECTRALHPDEALIASING" flag for the IcnNavierStokesSolver within the "SOLVERINFO" tag of the session file. Does this activate local de-aliasing only, i.e. de-aliasing with respect to the nonlinear convection term or are other parts of the equation (mass matrix, diffusive flux etc.) also "overintegrated"? Thanks and all the best Alex Sicher versendet mit [Proton Mail](https://proton.me/).

2 1

Taylor Hood approximation causing divergence
by Alexander Schukmann 11 May '24

11 May '24

******************* 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 everybody, I am currently working on a CG-SVV iLES of a cylinder flow at ReD = 3900 with the IcnNavierStokesSolver and wanted to run some 2D tests first, in order to find an appropriate polynomial order for the wall shear stress to converge in a cheap fassion, as suggested by Guglielmo in another post. I read that for under-resolved simulations it can be beneficial to resort to an equivalent of the Taylor Hood approximation and use a polynomial order for pressure that is reduced by 1 compared to velocity. However, the simulation immediately crashes, if "SPECTRALHPDEALIASING" is activated. If however I define the number of quadrature points manually within the "EXPANSIONS" tag to be at least equal to the amount needed to achieve consistent integration of quadratic nonlinearities or higher, the simulations remain stable. I assume by setting the number of quadrature points like this, global de-aliasing is performed, i.e. all the equation terms are overintegrated, instead of only the convective term being consistently integrated. Whereas by activating "SPECTRALHPDEALIASING" local de-aliasing is performed. Has anybody observed this kind of behavior and has an explanation? The simulation remains stable with "SPECTRALHPDEALIASING" switched on, if I rely on P_pressure = P_velocity for the polynomial order. Thanks and all the best Alex Sicher versendet mit [Proton Mail](https://proton.me/).

1 0

CFD to CAA coupling issue
by dauood Tariq 11 May '24

11 May '24

******************* 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, I have been attempting to use the cwipi function to form a coupled process with the nektar++ incompressible navier stokes solver and the acoustic solver however there is a repeated error that states: Fatal : Level 0 assertion violation Need to set up pressure field definition I am unsure what this means and why it is occuring. Any help would be greatly appreciated Kind regards, Dauood Tariq

1 0

Cwipi coupling not working
by dauood Tariq 08 May '24

08 May '24

******************* 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 attempting to run a cwipi coupling of the APE acoustic and incompressible navier stokes solver however they are not running correctly. The NS solver is sending out velocity files to the APE solver and is running correctly however the APE solver is meant to be reading files but cannot find a the file I have named despite it being done correctly. Any help would be appreciated. Regards, Dauood Tariq

1 0