[surfer]

Dear all,

I have a problem in using waves2Foam coupled with OpenFOAM 2.3.
So far I have successfully and intensely used waveFoam coupled with OpenFOAM 2.2.2 on a wide range of cases without any problem.
Now I would like to use the 2.3 version.
But using the same test cases as with 2.2.2 version, I get always the following error when running "setWaveField" utility:

Create time

Create mesh for time = 0


Reading waveProperties

Reading g

Reading field alpha

Reading field U

Reading field p

Setting the wave field ...

Face volume ratio : minimum: 0.214279 average: 0.973621
Face volume ratio check OK.
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::localCell::localizeCell(Foam::fvMesh const&, int const&) in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#4 Foam::localCell::localCell(Foam::fvMesh const&, int const&) in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#5 Foam::convexPolyhedral::dividePolyhedral(int const&, Foam::Vector<double> const&, Foam::Vector<double> const&) in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#6 Foam::setWaveField::correct() in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2FoamProcessing.so"
#7
in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/bin/setWaveField"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/surfer/OpenFOAM/surfer-2.3.0/platforms/linux64GccDPOpt/bin/setWaveField"
Segmentation fault (core dumped)

Does anybody have an idea about the reason of this error?

This happens with any case I try to run under OpenFOAM 2.3.

I thank you a lot in advance for your help.

Kind regards,

Surfer

[ngj]

Hi Surfer,

Something happens in a part of waves2Foam that I would be very unhappy to revisit (read: somewhat scared). It is therefore very important to figure out, whether it is related to waves2Foam or a inherited bug/feature from OF2.3 itself.

Therefore, please answer these questions first:

1. Does it happen for all your cases in 2.3?
2. Does the exact same case run without a problem in 2.2.x?
3. What type of cells does your mesh contain? I.e. hex's, tets, prisms, etc.
4. What happens, if you try to initialise with setWaveFields but a different source, e.g. take a random Stokes wave or shift the seaLevel a fraction of a cell.
5. Try to initialise with setFields and then run waves2Foam. Does it still happen? The relaxationZone should visit this part of the code as well, so I suppose the same error should pop up.

My experience with 2.3 is very limited, but I know that the waveFlume tutorial works, as I have tested in recently.

Kind regards,

Niels

[surfer]

Dear Niels,

I apologyse for bothering you again. Here the answers to your questions:

1. Does it happen for all your cases in 2.3?
Yes, it happens for all my cases, and for the waveFlume tutorial too;

2. Does the exact same case run without a problem in 2.2.x?
Yes, the same cases run for the 2.2.2 version. waveFlume tutorial runs perfectly for the 2.2.2 version.

3. What type of cells does your mesh contain? I.e. hex's, tets, prisms, etc.
In all my cases hexahedral cells are used, as well as in waveFlume tutorial.

4. What happens, if you try to initialise with setWaveFields but a different source, e.g. take a random Stokes wave or shift the seaLevel a fraction of a cell.
Nothing changes. I tried with different kinds of wave source as well as shifting the seaLevel by 0.00001 mm in the waveFlume case, and nothing has changed.

5. Try to initialise with setFields and then run waves2Foam. Does it still happen? The relaxationZone should visit this part of the code as well, so I suppose the same error should pop up.
I did what you suggested, and the run starts well, without any problems. But if I rerun setWaveField the same error appears as described in my previous post.

I do not have any idea about the problem.

I thank you a lot in advance for your help.

Kind regards,

surfer

[ngj]

Hi Surfer,

Thank you for the effort. I have been trying on my 2.3.0 version of OF, and I cannot reproduce the error. Not even with setWaveFields on the waveFlume tutorial.

Therefore, I suspect that the error is due to a more recent version of 2.3. Correct?

I have searched for the string "Face volume ratio" in the source, and it is found in polyMeshCheck.C, but this check is not invoked in my 2.3.0 version of OF, so I suspect that a check of the polyMesh object (primitive mesh object) created in localCell::localizeCell ruins things.

The only thing that I can suggest is that you try debugging in localCell::localiseCell and see, where it goes wrong. A very primitive approach could be to insert something like this:

Info << "A1" << endl;

etc. down through the code (and recompile waves2Foam). It will quickly give you an answer to where it goes wrong. I suppose that it is related to a call to one of the basic primitive mesh classes in OF.

Kind regards and good luck,

Niels