[논문 리뷰] Time dependent simulation of the Driven Lid Cavity at High Reynolds Number

Nuno Vieira Cardoso, Pedro Bicudo|ArXiv.org|2008. 09. 18.

Computational Fluid Dynamics and Aerodynamics참고 문헌 11인용 수 3

한 줄 요약

이 연구는 1024×1024 격자와 vorticity-stream function 설정에서의 SSPRK(5,4) 스킴을 사용하여 레이놀즈 수가 30,000에 이르는 2차원 구름이 이동하는 캐비티 유동의 시간에 따라 변하는 시뮬레이션을 제시한다. 双정밀도 CUDA 가속 계산을 통해 결과는 안정된 상태 솔루션으로 수렴하며, 이는 이전 연구에서 관측된 비정상적 솔루션이 거친 격자에 기인한 수치적 오류일 뿐이며, 물리적 행동이 아님을 밝혀내어 오랫동안 지속된 논란을 해결한다.

ABSTRACT

In this work, numerical solutions of the two dimensional time dependent incompressible flow, in a driven cavity at high Reynolds number Re, are presented. At high Re, there is a controversy. Some studies predicted that the flow is steady, others found time dependent non-steady flow, either periodic or aperiodic. In this study, the driven lid cavity is successfully solved using a very fine grid mesh, for Re up to $30 000$. We discretize the Vorticity-Stream formulation of the Navier-Stokes equation with the SSPRK(5,4) scheme in a $1024 imes1024$ grid. Using this very fine grid, the results obtained converge to a stationary solution. Detailed results for Re between $5 000$ and $30 000$ are presented. The driven lid cavity problem is solved with a NVIDIA GPU using the CUDA programming environment with double precision.

연구 동기 및 목표

• 이동 캐비티 유동이 고레이놀즈 수에서 비정상적인 행동을 보이는지 여부에 대한 오랫동안 지속된 논란을 해결하기 위해.
• 매우 미세한 격자로 수행된 시간에 따라 변하는 시뮬레이션의 수렴 여부가 정상 상태 솔루션인지, 아니면 진동하거나 주기적인 행동인지 확인하기 위해.
• 이전 문헌에서 보고된 비정상적 솔루션이 부족한 공간 해상도로 인한 수치적 오류인지 조사하기 위해.
• 충분히 미세한 격자가 진정한 물리적 솔루션을 포착하고 거친 메쉬에서 발생하는 잘못된 진동을 피하기 위해 필요하다는 주장을 검증하기 위해.
• GPU 가속을 통한 고해상도 직접 수치 시뮬레이션의 효과성을 입증하기 위해.

제안 방법

• 연속성 조건을 자동으로 만족시키기 위해 비압축성 2차원 나비에-스토크스 방정식의 비선형성에 기반한 코어스-스트림 함수 설정을 사용한 수치적 해법.
• 시간에 따라 변하는 코어스 운반 방정식의 이산화에 대해 시간 정확도와 안정성을 확보하기 위해 SSPRK(5,4) 룬게-쿠타 스킴을 사용.
• 특히 작은 소용돌이가 형성되는 모서리 영역에서 고공간 해상도를 확보하기 위해 1024×1024 정렬된 직각 격자를 사용.
• 오차 최소화 및 해의 정확도 향상을 위해 이중 정밀도 산술을 사용한 NVIDIA GPU에서의 CUDA 기반 구현.
• 모든 벽에 대해 슬립 조건을 적용하고, 상단 뚜껑만 일정 속도로 이동하여 유동을 구동하는 일관된 경계 조건 적용.
• 잔차 노름과 코어스 필드의 진화를 모니터링하여 해가 정적 상태에 도달할 때까지 시간 적분 수행.

실험 결과

연구 질문

• RQ1 레이놀즈 수가 30,000에 이르는 고레이놀즈 수의 이동 캐비티 유동에서 시간에 따라 변하는 행동을 보이거나, 정상 상태로 수렴하는가?
• RQ2 이전의 직접 수치 시뮬레이션(DNS)에서 보고된 비정상적 솔루션은 물리적 불안정성 때문인지, 아니면 거친 격자 해상도로 인한 수치적 오류 때문인가?
• RQ3 매우 미세한 격자 메쉬(1024×1024)가 주요 소용돌이뿐만 아니라 작은 모서리 소용돌이(예: BL3, TL2 등)까지 모두 해상하고 해를 안정화시킬 수 있는가?
• RQ4 GPU에서 이중 정밀도 산술을 사용하면 낮은 정밀도 방법에 비해 정확도와 수렴 성능이 향상되는가?
• RQ5 정적 나비에-스토크스 방정식에서 정상 상태 솔루션이 존재한다고 해서 물리적 실재성을 결론 내릴 수 있는가, 아니면 시간에 따라 변하는 시뮬레이션을 통해 이를 확인해야 하는가?

주요 결과

• 레이놀즈 수 5,000에서 30,000 사이에서는 1024×1024 격자에서 시간에 따라 변하는 시뮬레이션이 정상 상태 솔루션으로 수렴하며, 이는 진정한 물리적 솔루션이 정상 상태임을 시사한다.
• Re = 5,000일 때는 주요 소용돌이와 일부 모서리 소용돌이(BL1, BR1, TL1, TR1)만 명확하게 관찰 가능하며, BL3 및 TL2 소용돌이는 해상도 부족으로 인해 누락된다.
• Re = 15,000일 때에는 BL3 및 TL2 소용돌이를 포함한 모든 주요 소용돌이가 해상도가 충분하여 완전히 해결되고 안정된 상태를 유지한다.
• Re = 20,000, 25,000 및 30,000일 때도 해는 정적 상태를 유지하며 고레이놀즈 수에서 정상 상태로 수렴함을 확인한다.
• 이 연구는 이전 연구에서 보고된 비정상적 또는 주기적인 흐름이 거친 격자로 인한 수치적 진동 때문이며, 물리적 불안정성 때문이 아니라고 기인한다.
• 결과는 에르투르크(Erturk)의 주장, 즉 매우 미세한 격자가 작은 척도의 소용돌이를 해상하고 정확한 정상 상태 솔루션을 확보하기 위해 필수적이라는 주장을 지지한다.