stable_fluid 改写boundary condition

image

原来是四周是solid wall,想再加上两条杠,目的是让气体从这两条杠之间喷出。我的问题是我不知道该怎么改写linear equation Ap=-d 。

根据第11课中的公式,如果是杠边沿的grid,打个比方,下图中的蓝色grid,设置为p_i,j, 那他的p_{i,j+1}是有值的,左边p_{i-1, j} 是没有值的,用velocity 来代替。那最后对应这个grid,p_i,j 对应的是3, 上右下对应的是-1 吗? 对于-d的设置,对应的v_{i-1/2, j} 和 v_{i, j+1/2} 改成solid wall 的velocity。 如果是橙色的grid,那就是 p_i,j 对应4, 上下左右对应-1 吗?只需要改正v_{i, j+1/2} = solid wall velocity 就可以了吗?

image