[논문 리뷰] A Novel Model Order Reduction Approach for Navier-Stokes Equations at High Reynolds Number

Maciej Balajewicz, Earl H. Dowell외 1명|arXiv (Cornell University)|2012. 11. 07.

Model Reduction and Neural Networks참고 문헌 53인용 수 106

한 줄 요약

이 논문은 고레이놀즈수에서 나비에-스토크스 방정식에 대한 새로운 모델 순서 감소(MOR) 기법을 제안한다. 이 기법은 데이터 표현뿐만 아니라 에너지 소산 특성까지 최적화된 공간 기저 함수를 도입함으로써 경험적 난류 모델을 피한다. 기존의 표준 Proper Orthogonal Decomposition (POD)-Galerkin ROMs는 난류 운동에너지의 과대 예측을 겪는 반면, 새로운 기저 함수는 음의 에너지 생성률을 보이며 안정적이고 정확한 ROM을 생성하여 최소한의 기저 함수로 직접 수치 해석(DNS) 해에 수렴한다.

ABSTRACT

A new approach to model order reduction of the Navier-Stokes equations at high Reynolds number is proposed. Unlike traditional approaches, this method does not rely on empirical turbulence modeling or modification of the Navier-Stokes equations. It provides spatial basis functions different from the usual proper orthogonal decomposition basis function in that, in addition to optimally representing the training data set, the new basis functions also provide stable and accurate reduced-order models. The proposed approach is illustrated with two test cases: two-dimensional flow inside a square lid-driven cavity and a two-dimensional mixing layer.

연구 동기 및 목표

• 고레이놀즈수 난류 유동에서 표준 POD-Galerkin 감소된 순서 모델(Rom)의 불안정성과 부정확성을 해결하기 위해.
• 나비에-스토크스 방정식의 역학을 수정하는 경험적 난류 모델에 의존하지 않도록 하기 위해.
• 데이터 표현과 에너지 소산 특성을 동시에 최적화하는 공간 기저 함수를 개발하기 위해.
• 기존의 표준 POD-Galerkin 접근 방식에 비해 더 적은 모드로도 안정적이고 정확한 ROM을 가능하게 하기 위해.
• 레이놀즈수 기반의 표준 난류 유동 기준 문제인 리드드 캐비티와 혼합층에서의 방법의 효과성을 입증하기 위해.

제안 방법

• 기존의 표준 POD 모드와 다름없는 새로운 공간 기저 함수를 유도하기 위해 소규모 제약 최적화 문제를 수립한다.
• POD에서 재구성된 시간 평균 난류 운동에너지와 새로운 기저 함수 사이의 차이를 최소화하도록 기저 함수를 최적화한다.
• 물리적 일관성을 확보하기 위해 시간 평균 난류 운동에너지에 제약 조건을 부과한다.
• 비볼록 최적화 문제를 해결하기 위해 표준 MATLAB 최적화 루틴(fmincon)을 사용하며, QR 유사 정규화를 통한 정규화를 적용한다.
• 새로운 기저 함수를 사용하여 갈레르킨 투영을 적용하여 비압축성 나비에-스토크스 방정식의 ROM을 유도한다.
• 변환 행렬 X는 새로운 시간 계수의 정규직교성을 보장하여 수치적 안정성을 유지한다.

실험 결과

연구 질문

• RQ1 경험적 난류 모델을 피하면서도 고레이놀즈수 유동에서 정확성을 유지할 수 있는 모델 순서 감소 접근법을 개발할 수 있는가?
• RQ2 학습 데이터를 잘 표현할 뿐 아니라 소규모 난류 구조에서의 에너지 소산 특성을 내재적으로 포착할 수 있는 공간 기저 함수를 구성할 수 있는가?
• RQ3 새로운 기저 함수의 에너지 생성률은 기존의 표준 POD 모드와 비교해 어떻게 다른가? 이는 ROM의 안정성에 어떤 영향을 미치는가?
• RQ4 새로운 ROM은 기존의 표준 POD-Galerkin ROM에 비해 더 적은 모드로 난류 운동에너지의 정확한 예측을 달성할 수 있는가?
• RQ5 새로운 기저 함수는 기준 난류 흐름에서 직접 수치 해석(DNS) 결과로의 수렴을 어느 정도 향상시키는가?

주요 결과

• 새로운 기저 함수는 표준 POD 모드가 양의 ϵ을 가지는 것과는 달리 음의 난류 운동에너지 생성률(ϵ < 0)을 보이며, 이는 안정화 효과를 설명한다.
• Re = 20,000인 리드드 캐비티에서, 새로운 기저 함수를 사용한 ROM은 정확한 평균 난류 운동에너지로 수렴하지만, 기존의 표준 POD-Galerkin ROM은 이를 과대 예측한다.
• Reδω = 500인 혼합층에서, n = 5, 10, 20 모드를 사용한 새로운 ROM은 각각 50.21%, 66.67%, 79.14%의 시간 평균 난류 운동에너지의 비율을 포괄하며, 기존의 표준 POD-ROM을 능가한다.
• 기존의 표준 POD-Galerkin 접근 방식에 비해 훨씬 적은 모드로도 정확하고 안정적인 ROM을 달성하여 계산 비용을 감소시킨다.
• 제안된 방법은 나비에-스토크스 방정식을 수정하거나 경험적 에디 유동점성 항을 도입하지 않고도 ROM의 안정성을 성공적으로 확보한다.
• 최적화에서 유도된 변환 행렬 X는 새로운 시간 계수의 정규직교성을 보장하여 수치적 강건성을 유지한다.