【差分法】SMAC法で中心差分を用いてNavier-Stokes方程式を解きました C++コード付き

SMAC法(Simplified Maker And Cell method)を用いてNavier-Stokes方程式を数値的に解いてみました。テスト問題としてよく取り上げられる2次元キャビティー流れです。水で満たされた正方形があり、上部の壁だけ右へ一定の速さで動いていきます。最初は静止していた水が、上部の壁に引きずられて動き始めます。そして、最…