Hi Neks,
I'm trying to run a simulation with ~211.000 elements. The case runs fine on every mesh that I try as long as the number of elements is lower than 150.000. I'm using Nek5000 v19 with KTH Framework.
I am aware of the following two problems which occur if I use a mesh with more than 150.000 elements:
1. The BC for the passive scalar all seem to be zero now and therefore the PS is zero everywhere, even though it shouldn't be.
2. The more elements I have, the more of them don't seem to be part of the simulation anymore. The initial conditions for the velocity show weird gaps exactly along element lines, see attached screenshot.
I'm diagnosing all of this using ParaView, it might also be a ParaView problem, but Visit also has its problems with the files so I assume that it is more than a visualisation problem. Specifically, Visit says that it can't open the file because there is no mesh contained in the first snapshot. What makes me suspicious however, is that Nek doesn't throw a single error.
Here is what I assume could be the problem: 150.000 elements is coincidentally also the default maximum for the Nek tools. I changed this to 250.000 in the makefile and makefile.inc (in Nek5000/tools) but the problem is still there. Could it be that the maximum of 150.000 is also set somewhere else in the code where I have to change it? Also the mesh was generated with gmsh, so it might be some problem with gmsh2nek and genmap, but I executed both tools (compiled for 250.000 elements) on two different machines and the results were always the same (faulty).
What I tried to fix the problem:
1. Multiple attempts of cleaning and recompiling Nek and all Nek tools (especially gmsh2nek and genmap).
2. Compiling without HYPRE
I'm thankful for any help or suggestions and please let me know if you need any case files or outputs to make sense of the problem.
Best regards
Nils