I'm still in the process of studying computational fluid simulation, but I took two days off to implement this algorithm that simulate the motion of water surface using a simpler method (basically using sine waves).
The algorithm still doesn't give me the desired result, especially in the normals calculation, but it's as close as I could get in this short time.
And it's in a state in which I think makes sens to share it.
To develop this I used this two resources:
GPU Gems article from NVidia -
and attached source code.
And this implementation in OpenGL:
Hope you can find this useful, thanks for your support!