Home » How do I create 2D water with dynamic waves?

How do I create 2D water with dynamic waves?


I tried it.

Splashes (springs)

As that tutorial mentions, the surface of water is like a wire: If you pull on some point of the wire, the points next to that point will be pulled down too. All points are also attracted back to a baseline.

It’s basically lots of vertical springs next to each other that pull on each other also.

I sketched that in Lua using LÖVE and got this:

animation of a splash

Looks plausible. Oh Hooke, you handsome genius.

If you want to play with it, here is a JavaScript port courtesy of Phil! My code is at the end of this answer.

Background waves (stacked sines)

Natural background waves look to me like a bunch of sine waves (with different amplitudes, phases and wavelengths) all summed together. Here’s what that looked like when I wrote it:

background waves produced by sine interference

The interference patterns look pretty plausible.

All together now

So then it’s a pretty simple matter to sum together the splash waves and the background waves:

background waves, with splashes

When splashes happen, you can see small grey circles showing where the original background wave would be.

It looks a lot like that video you linked, so I’d consider this a successful experiment.

Here’s my main.lua (the only file). I think it’s quite readable.

-- Resolution of simulation
-- Width of simulation
WIDTH = 400
-- Spring constant for forces applied by adjacent points
-- Sprint constant for force applied to baseline
-- Vertical draw offset of simulation
Y_OFFSET = 300
-- Damping to apply to speed changes
DAMPING = 0.98
-- Number of iterations of point-influences-point to do on wave per step
-- (this makes the waves animate faster)

-- Make points to go on the wave
function makeWavePoints(numPoints)
    local t = {}
    for n = 1,numPoints do
        -- This represents a point on the wave
        local newPoint = {
            x    = n / numPoints * WIDTH,
            y    = Y_OFFSET,
            spd = {y=0}, -- speed with vertical component zero
            mass = 1
        t[n] = newPoint
    return t

-- A phase difference to apply to each sine
offset = 0

-- Amounts by which a particular sine is offset
sineOffsets = {}
-- Amounts by which a particular sine is amplified
sineAmplitudes = {}
-- Amounts by which a particular sine is stretched
sineStretches = {}
-- Amounts by which a particular sine's offset is multiplied
offsetStretches = {}
-- Set each sine's values to a reasonable random value
    table.insert(sineOffsets, -1 + 2*math.random())
    table.insert(sineAmplitudes, math.random()*BACKGROUND_WAVE_MAX_HEIGHT)
    table.insert(sineStretches, math.random()*BACKGROUND_WAVE_COMPRESSION)
    table.insert(offsetStretches, math.random()*BACKGROUND_WAVE_COMPRESSION)
-- This function sums together the sines generated above,
-- given an input value x
function overlapSines(x)
    local result = 0
        result = result
            + sineOffsets[i]
            + sineAmplitudes[i] * math.sin(
                x * sineStretches[i] + offset * offsetStretches[i])
    return result

wavePoints = makeWavePoints(NUM_POINTS)

-- Update the positions of each wave point
function updateWavePoints(points, dt)
    for i=1,ITERATIONS do
    for n,p in ipairs(points) do
        -- force to apply to this point
        local force = 0

        -- forces caused by the point immediately to the left or the right
        local forceFromLeft, forceFromRight

        if n == 1 then -- wrap to left-to-right
            local dy = points[# points].y - p.y
            forceFromLeft = SPRING_CONSTANT * dy
        else -- normally
            local dy = points[n-1].y - p.y
            forceFromLeft = SPRING_CONSTANT * dy
        if n == # points then -- wrap to right-to-left
            local dy = points[1].y - p.y
            forceFromRight = SPRING_CONSTANT * dy
        else -- normally
            local dy = points[n+1].y - p.y
            forceFromRight = SPRING_CONSTANT * dy

        -- Also apply force toward the baseline
        local dy = Y_OFFSET - p.y
        forceToBaseline = SPRING_CONSTANT_BASELINE * dy

        -- Sum up forces
        force = force + forceFromLeft
        force = force + forceFromRight
        force = force + forceToBaseline

        -- Calculate acceleration
        local acceleration = force / p.mass

        -- Apply acceleration (with damping)
        p.spd.y = DAMPING * p.spd.y + acceleration

        -- Apply speed
        p.y = p.y + p.spd.y

-- Callback when updating
function love.update(dt)
    if love.keyboard.isDown"k" then
        offset = offset + 1

    -- On click: Pick nearest point to mouse position
    if love.mouse.isDown("l") then
        local mouseX, mouseY = love.mouse.getPosition()
        local closestPoint = nil
        local closestDistance = nil
        for _,p in ipairs(wavePoints) do
            local distance = math.abs(mouseX-p.x)
            if closestDistance == nil then
                closestPoint = p
                closestDistance = distance
                if distance <= closestDistance then
                    closestPoint = p
                    closestDistance = distance

        closestPoint.y = love.mouse.getY()

    -- Update positions of points
    updateWavePoints(wavePoints, dt)

local circle = love.graphics.circle
local line   = love.graphics.line
local color  = love.graphics.setColor

-- Callback for drawing
function love.draw(dt)

    -- Draw baseline
    line(0, Y_OFFSET, WIDTH, Y_OFFSET)

    -- Draw "drop line" from cursor

    local mouseX, mouseY = love.mouse.getPosition()
    line(mouseX, 0, mouseX, Y_OFFSET)
    -- Draw click indicator
    if love.mouse.isDown"l" then
        love.graphics.circle("line", mouseX, mouseY, 20)

    -- Draw overlap wave animation indicator
    if love.keyboard.isDown "k" then
        love.graphics.print("Overlap waves PLAY", 10, Y_OFFSET+50)
        love.graphics.print("Overlap waves PAUSED", 10, Y_OFFSET+50)

    -- Draw points and line
    for n,p in ipairs(wavePoints) do
        -- Draw little grey circles for overlap waves
        circle("line", p.x, Y_OFFSET + overlapSines(p.x), 2)
        -- Draw blue circles for final wave
        circle("line", p.x, p.y + overlapSines(p.x), 4)
        -- Draw lines between circles
        if n == 1 then
            local leftPoint = wavePoints[n-1]
            line(leftPoint.x, leftPoint.y + overlapSines(leftPoint.x), p.x, p.y + overlapSines(p.x))

For the solution (mathematically speaking you can solve the problem with the solving of differential equations, but im sure they don’t do it that way) of creating waves you have 3 possibilities(depending on how detailed it should get):

  1. Calculate the waves with the trigonometric functions (most simple and the fastest)
  2. Do it like Anko has proposed
  3. Solve the differential equations
  4. Use texture lookups

Solution 1

Really simple, for each wave we calculate the (absolute) distance from each point of the surface to the source and we calculate the ‘hight’ with the formula

1.0f/(dist*dist) * sin(dist*FactorA + Phase)


  • dist is our distance
  • FactorA is a value which means how fast/dense the waves should be
  • Phase is the Phase of the wave, we need to increment it with time to get an animated wave

Note that we can add as many terms together as we like (superposition principle).


  • Its really fast to calculate
  • Is easy to implement


  • For (simple) reflections on a 1d Surface we need to create “ghost” wave sources to simulate reflections, this is more complicated at 2d surfaces and it is one of the limitations of this simple approach

Solution 2


  • Its simple too
  • It allows to calculate reflections easily
  • It can be extended to 2d or 3d space relativly easily


  • Can get numerically instable if the dumping value is too high
  • needs more calculation power than Solution 1 (but not so much like Solution 3)

Solution 3

Now i hit a hard wall, this is the most complicated solution.

I didn’t implement this one but it is possible to solve these monsters.

Here you can find a presentation about the mathematics of it, its not simple and there exists also differential equations for different kinds of waves.

Here is a not complete list with some differential Equations to solve more special cases (Solitons, Peakons, …)


  • Realistic waves


  • For most games not worth the effort
  • Needs the most calculation time

Solution 4

A bit more complicated than solution 1 but not so complicated a solution 3.

We use precalculated textures and blend them together, after that we use displacement mapping (actually a method for 2d waves but the principle can also work for 1d waves)

The game sturmovik has used this approach but i don’t find the link to the article about it.


  • it is more simple than 3
  • it gets good looking results (for 2d)
  • it can look realistic if the artists good a great job


  • difficult to animate
  • repeated patterns could get visible on the horizon

To add constant waves add a couple of sine-waves after you have calculated dynamics. For simplicity I would make this displacement a graphical effect only and not let it affect the dynamics themselves but you could try both alternatives and see which works out the best.

To make the “splashhole” smaller I would suggest altering the method Splash(int index, float speed) so that it directly affects not only index but also some of the close vertices, so as to spread out the effect but still have the same “energy”. The number of vertices affected could depend on how wide your object is. You’ll probably need to tweak the effect a lot before you have a perfect result.

To texture the deeper parts of the water you could either do as described in the article and just make the deeper part “more blue” or you could interpolate between two textures depending on the depth of the water.

Related Solutions

Extract file from docker image?

You can extract files from an image with the following commands: docker create $image # returns container ID docker cp $container_id:$source_path $destination_path docker rm $container_id According to the docker create documentation, this doesn't run the...

Transfer files using scp: permission denied

Your commands are trying to put the new Document to the root (/) of your machine. What you want to do is to transfer them to your home directory (since you have no permissions to write to /). If path to your home is something like /home/erez try the following:...

What’s the purpose of DH Parameters?

What exactly is the purpose of these DH Parameters? These parameters define how OpenSSL performs the Diffie-Hellman (DH) key-exchange. As you stated correctly they include a field prime p and a generator g. The purpose of the availability to customize these...

How to rsync multiple source folders

You can pass multiple source arguments. rsync -a /etc/fstab /home/user/download bkp This creates bkp/fstab and bkp/download, like the separate commands you gave. It may be desirable to preserve the source structure instead. To do this, use / as the source and...

Benefits of Structured Logging vs basic logging

There are two fundamental advances with the structured approach that can't be emulated using text logs without (sometimes extreme levels of) additional effort. Event Types When you write two events with log4net like: log.Debug("Disk quota {0} exceeded by user...

Interfaces vs Types in TypeScript

2019 Update The current answers and the official documentation are outdated. And for those new to TypeScript, the terminology used isn't clear without examples. Below is a list of up-to-date differences. 1. Objects / Functions Both can be used to describe the...

Get total as you type with added column (append) using jQuery

One issue if that the newly-added column id's are missing the id number. If you look at the id, it only shows "price-", when it should probably be "price-2-1", since the original ones are "price-1", and the original ones should probably be something like...

Determining if a file is a hard link or symbolic link?

Jim's answer explains how to test for a symlink: by using test's -L test. But testing for a "hard link" is, well, strictly speaking not what you want. Hard links work because of how Unix handles files: each file is represented by a single inode. Then a single...

How to restrict a Google search to results of a specific language?

You can do that using the advanced search options: http://www.googleguide.com/sharpening_queries.html I also found this, which might work for you: http://www.searchenginejournal.com/how-to-see-google-search-results-for-other-locations/25203/ Just wanted to add...

Random map generation

Among the many other related questions on the site, there's an often linked article for map generation: Polygonal Map Generation for Games you can glean some good strategies from that article, but it can't really be used as is. While not a tutorial, there's an...

How to prettyprint a JSON file?

The json module already implements some basic pretty printing in the dump and dumps functions, with the indent parameter that specifies how many spaces to indent by: >>> import json >>> >>> your_json = '["foo", {"bar":["baz", null,...

How can I avoid the battery charging when connected via USB?

I have an Android 4.0.3 phone without root access so can't test any of this but let me point you to /sys/class/power_supply/battery/ which gives some info/control over charging issues. In particular there is charging_enabled which gives the current state (0 not...

How to transform given dataset in python? [closed]

From your expected result, it appears that each "group" is based on contiguous id values. For this, you can use the compare-cumsum-groupby pattern, and then use agg to get the min and max values. # Sample data. df = pd.DataFrame( {'id': [1, 2, 2, 2, 2, 2, 1, 1,...

Output of the following C++ Program [closed]

It works exactly like this non-recursive translation: int func_0() { return 2; } int func_1() { return 3; } int func_2() { return func_1() + func_0(); } // Returns 3 + 2 = 5 int func_3() { return func_2() + func_1(); } // Returns 5 + 3 = 8 int func_4() { return...

Making a circle out of . (periods) [closed]

Here's the maths and even an example program in C: http://pixwiki.bafsoft.com/mags/5/articles/circle/sincos.htm (link no longer exists). And position: absolute, left and top will let you draw: http://www.w3.org/TR/CSS2/visuren.html#choose-position Any further...

Should I use a code converter (Python to C++)?

Generally it's an awful way to write code, and does not guarantee that it will be any faster. Things which are simple and fast in one language can be complex and slow in another. You're better off either learning how to write fast Python code or learning C++...

tkinter: cannot concatenate ‘str’ and ‘float’ objects

This one line is more than enough to cause the problem: text="რეგულარი >> "+2.23+ 'GEL' 2.23 is a floating-point value; 'GEL' is a string. What does it mean to add an arithmetic value and a string of letters? If you want the string label 'რეგულარი...