myzek 4 days ago

So once upon a time I stumbled upon simulating fluids in gamedev and I really wanted to learn how it works. Fast forward 2 months and I decided to write down everything I learned to hopefully make it easier for others in the future!

  • magicalhippo 4 days ago

    Fluid sims are just so darn fun! Nice writeup, very accessible.

    Couple of notes, I think you forgot to apply timestep when adding rocket exhaust velocity, pretty sure it should be

      u[idx] += backward.x * flame_velocity_amount * falloff * delta
      v[idx] += backward.y * flame_velocity_amount * falloff * delta
    

    You need to compensate by scaling up flame_velocity_amount, I used 85, seemed about the same.

    • myzek 4 days ago

      Yeah it seems I missed that, adding the rocket at the end was a cherry on the top and I was so exhausted already I'm surprised it even works lol

      • magicalhippo 3 days ago

        It made it very fun to play with though :D

  • nick__m 1 day ago

    excellent article, great vulgarisation and human written !

    Thank you

  • brandonpelfrey 1 day ago

    Great job. I have spent a lot of time working on fluid simulations (I still am). Glad to see more people still mesmerized. If you’re interested, this rabbit hole goes very deep.

markstock 1 day ago

Before you go adding vorticity confinement, consider performing a higher-order backward advection scheme (Runge-Kutta 2nd or similar), and using a higher-order interpolation method (triangle-shaped cloud instead of bilinear).

In my implementations I use 4th order for both and vortices stick around a lot longer.

  • magicalhippo 13 hours ago

    If you implement this on the GPU, it's my understanding you can get the 4th-order interpolation quite cheaply exploiting the bilinear texture sampling hardware[1].

    So instead of reading 16 grid values and combining them to get the interpolated sample value, you can fetch 4 bilinearly filtered samples and combine those. And thanks to the hardware filtering, those bilinear samples cost basically the same as reading an unfiltered value.

    [1]: https://developer.nvidia.com/gpugems/gpugems2/part-iii-high-...

Stevvo 1 day ago

For anyone wanting to dive further, Fluid Simulation for Computer Graphics by Robert Bridson is the definitive textbook.

jedimastert 23 hours ago

I appreciate the AI disclosure, I hope that becomes normalized

  • titanomachy 9 hours ago

    I wonder how it'll read in 20 years. Like all the covid references in blog posts from 2020, or vague references to "the current political situation" from any point in the last 10 years or so.

amarant 1 day ago

This is great! When I have some leftover time I want to try copying this implementation for 3D. I reckon I could get away with minimal modifications to support the third axis...I think...

That'll perform even worse though, hopefully my CPU can handle it or I'm gonna need a lot of leftover time to make a shader

  • markstock 1 day ago

    You are correct: Stable Fluids extends to 3d relatively easily.

frankdlc222 4 days ago

This is really cool. I love how much detail you went into explaining the setup and walking through each piece of the simulation. Definitely bookmarking this to play around with later!

amelius 1 day ago

Did they test if it satisfies the relevant conservation laws?

  • brandonpelfrey 1 day ago

    It won’t satisfy those laws. It’s also not a goal though of this post.

    • amelius 1 day ago

      I mean if you're writing a ray tracer and the reflected light has more intensity than the light sources, then that's not desired. You can have the same sort of thing going on with a fluid simulation.

      • markstock 1 day ago

        One of the nice aspects of Stable Fluids is that you don't need to iterate the pressure correction terms to convergence. Just run a fixed number of Jacobi or Gauss-Seidel sweeps and keep performance consistent. The only drawback of this is some mass loss in areas, which for the present purposes is acceptable.

      • dahart 20 hours ago

        Sure, but conservation in ray tracing is also a goal that not everyone has and isn’t required for teaching or games or making pretty & even plausible images. There are plenty of situations in both ray tracing and fluid simulation where conservation is not desirable.

CamperBob2 16 hours ago

Very cool. For your next project, how about a DIY 2D EM field solver? You've done a lot of the work already.

  • __patchbit__ 5 hours ago

    AI model latent space needs to be instrumented to a DAW and 3D helical visual effects on triangle slices to show the mathsemantic field gears grinding at the prompt: What's on your mind?

nryoo 1 day ago

Nice writeup. One thing worth adding to the limitations: without vorticity confinement, the Gauss-Seidel projection step quietly dissipates the small-scale curl that makes smoke look like smoke.

The 2001 Fedkiw/Stam/Jensen "Visual Simulation of Smoke" paper added it back as a correction force for exactly this reason. At N=16 it doesn't matter much because the grid itself can't represent fine vortices, but the moment you crank N up the missing confinement becomes visible.

analog8374 1 day ago

Oh don't let us pinch zoom. That would be a disaster.

  • zamadatix 1 day ago

    Assuming you mean the site - pinch zoom & pan works for me on Windows and iOS?

    • analog8374 1 day ago

      well android chrome doesn't

      • OsrsNeedsf2P 19 hours ago

        Web exports from Godot (and other game engines) played on mobile is a hard place between compatibility, performance, and many other factors. It's getting better, but very slowly. Try on PC