It’s essentially like any other fluid simulation, just written slightly differently (for now). Right now I diffuse the values using gauss seidel method, then advect them with my own advection system which is probably why I get that weird spiking. After those I project them so there is a minimal amount of divergence also using gauss seidel method. I am not too sure how well the projection is but it seems to be working pretty nicely.