[gsg001]

Hi ngj,

After reading H Jasak's document[1] where they have taken waves to foam and implicitly coupled it to navalFoam, I began experimenting if this was possible using waveFoam.

The main difference I found was that they[2] add the following line to UEqn.h
" ==
- alpha1*fvm::Sp(beach.damping(), U)"
while you add relaxation.correct() and interface.correct() to waveFoam.C.

I tried modifying the UEqn.h in wave foam as follows
" ==
- alpha1*fvm::Sp(relaxationZone.targetVelocityField( ), U)"

It spits up "error: expected primary-expression before '.' token." I am guessing the way it is written relaxationZone.H is not callable from within UEqn.H

The reason for trying to modify the solver to implicit is to allow the simulation to run at a Courant No of 25 -100. waves2Foam fails to make the correct wave at higher courant Numbers despite switching courantCorrection ON.

Thanks in advance

[1] http://www.turing-gateway.cam.ac.uk/documents/Jasak.pdf
[2] navalFoamVariant2 https://github.com/Unofficial-Extend...alFoamVariant2

[ngj]

Good morning,

Note that the version of navalFoam that you are looking at is not with the implicit relaxation zone technique. The approach by Hrv is very different from the simple addition of a numerical beach.

The numerical beach is available in waves2Foam, but it has never been populated with actual functionality. Therefore, it is likely that I will retire it in the future.

You are looking at the right place in the code to obtain the implicit relaxation, but you are not using the right equations. I am not aware, whether these equations have been published by Hrv. I do have a working version of the implicit relaxation, but it will not be coming out any time in the near future.

Furthermore, note that the wave propagation in itself is very diffusive, so it could be that you only gain dissipation of the wave energy by using larger Courant numbers.

Kind regards,

Niels