|
304 | 304 | " dx**2 * dy**2 / (2 * (dx**2 + dy**2)) * \n", |
305 | 305 | " b[1:-1,1:-1])\n", |
306 | 306 | "\n", |
307 | | - " p[:, -1] = p[:, -2] ##dp/dy = 0 at x = 2\n", |
308 | | - " p[0, :] = p[1, :] ##dp/dy = 0 at y = 0\n", |
309 | | - " p[:, 0] = p[:, 1] ##dp/dx = 0 at x = 0\n", |
310 | | - " p[-1, :] = 0 ##p = 0 at y = 2\n", |
| 307 | + " p[:, -1] = p[:, -2] # dp/dx = 0 at x = 2\n", |
| 308 | + " p[0, :] = p[1, :] # dp/dy = 0 at y = 0\n", |
| 309 | + " p[:, 0] = p[:, 1] # dp/dx = 0 at x = 0\n", |
| 310 | + " p[-1, :] = 0 # p = 0 at y = 2\n", |
311 | 311 | " \n", |
312 | 312 | " return p" |
313 | 313 | ] |
|
359 | 359 | " dt / dy**2 *\n", |
360 | 360 | " (vn[2:, 1:-1] - 2 * vn[1:-1, 1:-1] + vn[0:-2, 1:-1])))\n", |
361 | 361 | "\n", |
362 | | - " u[0, :] = 0\n", |
363 | | - " u[:, 0] = 0\n", |
| 362 | + " u[0, :] = 0\n", |
| 363 | + " u[:, 0] = 0\n", |
364 | 364 | " u[:, -1] = 0\n", |
365 | | - " u[-1, :] = 1 #set velocity on cavity lid equal to 1\n", |
366 | | - " v[0, :] = 0\n", |
367 | | - " v[-1, :]=0\n", |
368 | | - " v[:, 0] = 0\n", |
| 365 | + " u[-1, :] = 1 # set velocity on cavity lid equal to 1\n", |
| 366 | + " v[0, :] = 0\n", |
| 367 | + " v[-1, :] = 0\n", |
| 368 | + " v[:, 0] = 0\n", |
369 | 369 | " v[:, -1] = 0\n", |
370 | 370 | " \n", |
371 | 371 | " \n", |
|
660 | 660 | "name": "python", |
661 | 661 | "nbconvert_exporter": "python", |
662 | 662 | "pygments_lexer": "ipython3", |
663 | | - "version": "3.6.1" |
| 663 | + "version": "3.6.5" |
664 | 664 | } |
665 | 665 | }, |
666 | 666 | "nbformat": 4, |
|
0 commit comments