스트림 함수 -Vorticity 방법의 경계 조건

MATLAB 유한 차분 코드를 사용하여 점성의 비 압축 형 Navier Stokes 방정식의 Stream 함수 -Vorticity 공식을 기반으로 뚜껑 구동 공동 흐름을 해결했습니다. 방법에 대한 자세한 내용은 여기를 참조하십시오 .

흐름이 왼쪽에서 균일하고 흐름이 위쪽과 아래쪽의 수평 벽에 의해 제한되는 사각형 도메인의 사각형 상자 주위의 흐름을 시뮬레이션하도록 코드를 변경하고 싶습니다. 그러나 흐름이 나가는 도메인의 오른쪽에는 경계 조건을 표현하는 방법을 모릅니다.

오른쪽에 경계 조건을 어떻게 표현할 수 있습니까? 이 방법으로도 가능합니까?

(다음에서는 양의 x 방향으로 주요 흐름으로 2 차원이라고 가정하면 속도 벡터는 )$\vec{u} = (u \quad v)^T$

콘센트에 적용 할 수있는 경계 조건 :

1. $\frac{\partial u}{\partial x}$$v = 0$
2. $\frac{\vec{u}}{t} + \vec{c} \frac{\vec{u}}{\vec{x}} = 0$$\vec{c}$

이 두 조건은 압력을 무시합니다. 물리적 관점에서 분석법은 흐름을 유도하는 압력 구배를 생성해야합니다. 나는 공식이 그것을 허용한다고 생각하지만 나는 결코 그것을 사용해 본 적이 없다. 명심해야 할 또 다른 문제는 압력 구배가 장애물 뒤에서 발전 할 수 있도록 도메인이 충분히 길어야한다는 것입니다. 즉, 전체 웨이크 (또는 대부분)를 시뮬레이션해야합니다. 깨우기를 통과하면 솔버가 수렴하지 않을 수도 있습니다.

$\omega = 0$