(수학이야기)
나비에 스토크스 방정식
du/dt+(u•del)u=1/rhodel(p)+nu*del^2u
비압축성조건 div u=0
은 수치해석으로 풀수 있다
단순히 유한차분법을 사용하면 (du/dt=u(n+1)-u(n)/dt)
u(n+1)값이 비압축성 조건을 만족하지 않을수도 있고,압력을 구할 방법이 없다
그래서 중간에 u*이라는 중간속도를 넣어서 비압축성 조건을 만족하게 만드는데, 중간 과정속에서는 poisson방정식의 풀이가 필요하다
poisson방정식은 유한차분법 한후 행렬의 역행렬을 구하거나, 그 행렬을 반복법으로 구한다
편미분방정식엔 경계조건이 필요한데, dirichlet조건이나, neumann조건이 있다
댓글 0