lauchmelder/NumericalMethods
1
0
Fork 0
Code Issues Pull requests Actions Packages Projects Releases Wiki Activity

fixed projection bug

Browse source
This commit is contained in:
Lauchmelder 2021-12-16 16:38:07 +01:00
parent b947b0b263
commit 8bfed7b189
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/EulerFluid/FluidField.cpp
View file

@ -222,7 +222,7 @@ void FluidField::Project()
{ {
for (int j = 1; j <= N; j++) for (int j = 1; j <= N; j++)
{ {
velocity[1].vertical[IDX(i, j, size)] = -0.5 * h * (velocity.Current().vertical[IDX(i + 1, j, size)] - velocity.Current().vertical[IDX(i - 1, j, size)] + velocity.Current().vertical[IDX(i, j + 1, size)] - velocity.Current().vertical[IDX(i, j - 1, size)]); velocity[1].vertical[IDX(i, j, size)] = -0.5 * h * (velocity.Current().horizontal[IDX(i + 1, j, size)] - velocity.Current().horizontal[IDX(i - 1, j, size)] + velocity.Current().vertical[IDX(i, j + 1, size)] - velocity.Current().vertical[IDX(i, j - 1, size)]);
velocity[1].horizontal[IDX(i, j, size)] = 0; velocity[1].horizontal[IDX(i, j, size)] = 0;
} }
} }