Diffusion de vélocité

Le deuxième terme dans les équations de Navier-Stokes détermine la diffusion dans le fluide. La première étape de la simulation est la diffusion des valeurs de vélocité. L'étape de diffusion cherche une évolution de fluide qui vérifie le terme de diffusion dans les équations. Un solveur itératif de Gauss-Seidel dans la fonction lin_solve cherche une solution. Cette fonction est appliquée sur chaque composante de la vélocité.

Le solveur génère des solutions qui ne respectent pas la conservation de masse dans le fluide. Une autre fonction project utilisant aussi un solveur de Gauss-Seidel est appliquée pour transformer le champ de vélocité afin qu'il respecte la conservation de masse.