[anandanandan86]

Hello, I am using OpenLB 1.4

[anandanandan86]

Hello,

Thank you for your tips. I checked these points and can confirm as below:

– does the problem depend on the number of cuboids/ cores?
No, for the problematic case (15m x 15m x 5m), I tried running with both 1 cuboid/core as well as 4 cuboids/cores, but observe the same error

– plot and check the material numbers in paraview, before and after each rename
I did this and I can confirm that the issue arises at the very first rename:

superGeometry.rename(0, 2, indicator);

However, as the indicator layer is created from the grid generated by the STLReader, this indicates that the problem is occurring from the STLReader itself. I cannot diagnose this further as the STLReader is a black-box to me.

– increasing the extent of the CuboidGeometry can avoid round-off errors at that stage
I tried this by increasing the extension distance in the IndicatorLayer3D<T> extendedDomain declaration five-fold, but it does not fix the issue

If possible I would like to diagnose what is going on within the STLReader further, but I am not sure how to do this, so I would like to request any possible support on further resolving this issue.

[anandanandan86]

Hello, this is quite a late response but yes, we are running in a Docker container. Is there any way to suppress the output in that case?

[anandanandan86]

Hello, Sorry for my late reply and thank you for the explanation. I was able to visualize it.

[anandanandan86]

Hello, I checked and this error seems to appear during every iteration. In this calculation setBoundaryValues is being called outside the loop (only at the start of the calculation as boundary conditions do not change during the calculation). I also commented out getResults temporarily, but even so, the error is output repeatedly during each iteration (with only collideAndStream being called for each iteration).

[anandanandan86]

I think the correct compiler settings should have been used as there is no error during compilation, and the calculation is exporting contours and result files as expected (at least below a certain core count). Is there any file that I should share to shed more light on the issue?

The program is a modification of the cylinder3d example case that has been changed to model flow within a 3d room with a few obstructions inside. The Smagorinsky turbulence model that comes with OpenLB is also implemented in this case (in contrast to the original cylinder3d case, which was laminar).

[anandanandan86]

Note: the #s I included in the error message stand in for alphanumeric characters that occasionally change depending on the line

[Author]

Posts