In this project a water-simulation was implemented based on the approach suggested in the presentation "Fast Water Simulation for Games Using Height Fields" by Matthias Müller-Fischer. Some effort was invested to make the result visually pleasing.
I included the [document](https://github.com/hpatjens/Watersimulation/tree/master/documentation/Fast Water Simulation for Games Using Height Fields - Matthias Müller-Fischer.pdf) in this repository as well in case of the document being offline in the future.
This is my first attempt to simulate and render water. This is not a tutorial or even an attempt to do things right. Just keep this in mind when you browse the code.
The water simulation is based on the method that Fischer proposes for 2D heightfields and implemented using the compute shader. Everything is rendered with OpenGL 4.3 at real-time frame rates applying a simple ad-hoc shading - this is not a PBR renderer. Multiple rendering passes are used to determine the depth of the water and calculating fake refractions of the ground. Light is attenuated by the water based on the distance the light has to travel inside the medium. To give the impression of subsurface-scattering and shift in color a look-up texture is used. Caustics are achieved by projecting a texture onto the ground. Nothing real here either. To give the impression of small ripples on the water a normals are read from a normal-map.
There is a compile.sh script you can execute. The needed packages are listed in the following table.
| Package | Installation on Ubuntu |
|---|---|
| GLEW | sudo apt-get install libglew-dev |
| GLFW | sudo apt-get install libglfw3-dev |
Execute sudo apt-get update before installing the libraries.
There is a visual studio project in the visual studio folder.
| Key | Function |
|---|---|
| ESC | Quit |
| F1 | Wireframe |
| F2 | Show Normals |
| F3 | Free-fly |
| F4 | - |
| F5 | Render all |
| F6 | Render ground |
| F7 | Render watermap |
| F8 | Render water |
The simulation directly operates on the vertices of the water mesh. This ties the simulation and rendering accuracy directly together which can be a drawback when it comes to dynamically scaling the rendering workload to ensure constant framerates. (Which is not happening in this application.) A huge benefit is the possibility to let the vertex and compute shader operate on the same buffer object. While rendering the buffer is interpreted via the GL_ARRAY_BUFFER target and for the simulation it is provided to the compute shader via the GL_STORAGE_BUFFER target.
These lines implement the above mentioned simulation method:
vec4 position = positionsPrev[index];
// f = c^2*(u[i+1,j]+u[i-1,j]+u[i,j+1]+u[i,j-1] – 4u[i,j])/h^2
float f = c * c * (
positionsPrev[clampedIndex(gl_GlobalInvocationID.xy + uvec2(-1, 0))].y +
positionsPrev[clampedIndex(gl_GlobalInvocationID.xy + uvec2( 1, 0))].y +
positionsPrev[clampedIndex(gl_GlobalInvocationID.xy + uvec2( 0, -1))].y +
positionsPrev[clampedIndex(gl_GlobalInvocationID.xy + uvec2( 0, 1))].y -
4.0 * position.y) / (h * h);
// v[i,j] = v[i,j] + f*∆t
position.w = position.w + f * DeltaTime;
// unew[i,j] = u[i,j] + v[i]*∆t
position.y = position.y + position.w * DeltaTime;The vertices of the water surface are altered by the compute shader using a double buffering. This ensures results that are independent of the execution order of the compute shader threads. positionsPrev refers to the vertex positions from the previous frame. position holds the position of the vertex, that is currently computed. The w component of the position stores the velocity of that water column. The last line alters the height of the current vertex, which is the y-coordinate.
To make the simulation interesting the water is pulled up every few seconds by two gaussian bells.
There are 3 framebuffers on which the rendering operates.
The first framebuffer contains the water surface which is rendered from the position of the light source using an orthographic projection. A depth buffer and a color buffer is rendered. The color buffer is made up of 3 channels that store the normals of the water surface. This buffer has a resolution of 1024x1024.
| Orthographic projection | Normal buffer |
|---|---|
 |
 |
The ground (without the water on top) is rendered from the viewers position. This framebuffer stores the color of the pixels and the depth from the viewer. This buffer as well as the third is rendered on full resolution.
| Perspective projection | Color buffer |
|---|---|
 |
 |
The water surface is rendered from the viewers position.
| Perspective projection | Color buffer |
|---|---|
 |
 |
Finally the second and third buffer are combined and rendered onto the default framebuffer.
The light gets attenuated by the water on the way to the ground and on the way to the camera. From the watermap we can get the distance the light travels on the way to the ground by projecting the currently processed fragment onto the watermap and calculating the difference to the stored value (w). The attentuation on the way to the viewer can be obtained by calulating the difference between the water fragment and the depth value that is stored in the background buffer (b).
When the ground is rendered without the water, you can see how the depth of the water affects the lighting. On the right bottom in this image you can see that individual waves are included in the depth-calculation. High waves result in a stronger attenuation and therefore a darker ground. Additionally you can see some distortion in the caustics happen.
This is of course a very approximated model that assumes homogenous density distribution and no scattering at all.
| Perspective projection | Color buffer (depth buffer is not shown) |
|---|---|
 |
|
To immitate spectral effects that change the color of the water a lookup table is used to control the change of colors in different depths. The table is implemented as a simple 1D texture.
struct Pixel { GLubyte r, g, b; };
const int size = 3;
Pixel data[3];
data[0] = Pixel{ 2, 204, 147 };
data[1] = Pixel{ 2, 127, 199 };
data[2] = Pixel{ 1, 9, 100 };
...
glTexImage1D(GL_TEXTURE_1D, 0, GL_RGBA, size, 0, GL_RGB, GL_UNSIGNED_BYTE, data);Because the simulation constantly changes the water surface, the normals have to be recalculated every frame. The geometry shader is used to verify the calculations.
