[lxwd]

Dear Jacobsen,

Thank you for providing such a useful toolbox! I have successfully compiled waves2Foam on foam-extend-4.0, but I cannot run tutorials in parallel. I have tested several cases and found that the parallel can only run without the use of waveGaugesNProbes, so I think there might be something wrong in sampling. I attached the log file for your information.

I also encounter a problem that the solitary wave cannot keep its height outside the inlet relaxation zone (wave height should be 0.15m) and there even exist a small amplitude wave after the pass of the solitary wave at the middle of wave flume. I am sure my mesh is fine enough (2cm in x direction, 1cm in y direction) and there is no structure within the domain. I also attached wave elevation time series for one wave gauge right after inlet relaxation zone and one wave gauge at the middle of wave flume. Could you please give me some suggestions on this problem?

Thanks in advance and best regards,

Lee

[ngj]

Hi Lee,

Thank you for reporting the wave gauge issue. Could I ask you to test, whether you can execute in parallel, if all wave gauges are within one processor?

With respect to the solitary wave, then I assume that it is caused by the low accuracy (first order) of the implemented wave theory. I have seen myself that it creates extra waves and lose height.

Kind regards

Niels

[lxwd]

Hi Niels,

Thank you for your reply. I will do the test this morning and report back.

Have you ever tested the applicability of using a stream function wave to represent solitary wave or Do you have any plan to implement a higher order solitary wave, e.g. 3rd order, in waves2Foam?

Thanks again for your help.

Best regards,

Lee

[ngj]

No, I have never really tried anything with solitary waves, as it has never been a topic in my research.

Stream function waves could be a possibility, but you will have a negative velocity in the trough, so I do not know, whether that would be a problem for your case?

If you make an implementation of a higher order wave, I will gladly include it in the repository.

Kind regards

Niels

[lxwd]

It still doesn't work even all wave gauges are within one processor. This time it seems to get stuck at the first time step and I have to terminate the job.

I think maybe a higher oder solitary wave theory is more suitable for me. I will try to make one and let you know if it can work.

Best Regards,

Lee

[ngj]

Hi Lee,

I have just tried running the waveFlume case in parallel (2 and 4 processors) and it runs smoothly, so I cannot reproduce your error. Could I please ask you to start from the waveFlume case and check, whether the problem also exists there on your compilation?

If you do not already do it, try to give this environmental property, before executing the test:

export FOAM_ABORT=1

This normally gives a bit more information on the source of the error. As it is now, I cannot help you from this end.


Kind regards,

Niels

[lxwd]

Hi Niels,

I have just run the waveFlume case without any change in parallel (6 processors). I also gave that environmental property before I running the case. Below is the link to the test case in dropbox.

https://www.dropbox.com/sh/mly1hy5i0...VsAAQg33a?dl=0

I am now working on implementation of a higher order solitary wave. I have already checked the solitaryFirst.C file, but I am confused on how you obtain the expression for pressure term. I found several references about higher order solitary wave and they only gave the expressions for wave profile and water particle velocity. Could you please let me know where you find this equation or how you derive it, if it is derived by yourself.

Many thanks and best regards,

Lee

[ngj]

Hi Lee,

I do the same thing as you, but I cannot reproduce the problem. I can only suggest that you try to replace the wave gauges with the native sampling utility. The surface elevation tool is (long ago) derived from line sampling. If it still fail as well, the problem is unrelated to surfaceelevation.

Concerning pressure field, then I believe that I simply inserted the velocity field and the surface elevation in the Bernoulli equation to obtain the pressure. It is a long time ago, so I am not certain.

Kind regards

Niels