[AbdoMusaad]

Hi everyone! i am new to OpenFoam (first Thread).

I am trying to simulate a 2D Numerical wave tank with cyclic boundary conditions for the inlet and the outlet. i started from the stokes 1 example of OpenFoam v2306. in the example the waves are generated at the inlet.


So, my problem now is how can i create the waves without the waveAlpha condition at the inlet? should i create a nonuniform alpha.water , p_rgh , U files? and how can i do that in a practical way? i tried running the stokes 1 example until it reaches a developed stage, then change the boundary conditions in the fields files. but i could not make it to work and it is not ideal for my intended use. ideally i want to impose the theoretical stokes 1 solution and see how much the numerical solution will deviate form the theoretical solution due to spatial and temporal discretisation excluding the wave generation and absorption errors.


sorry if it is a basic question.



it would be very helpful if you can provide me with examples. Thanks in advance.

[AbdoMusaad]

I am posting my solution to this problem for reviewing, suggestions and hopefully the benefit of others. Note that i am not 100% sure if this is the best and most correct way to do it. but it looks fine from what i saw.



what i did basically is initializing the VOF (alpha) and the velocity at time zero using codeStream and let openfoam deal with calculating the pressure. i relied on "water wave mechanics for engineering and scientists" as a reference for the theoretical equations of the velocity and water surface shape.


Cheers