[논문 리뷰] Regularized Reduced Order Models for a Stochastic Burgers Equation

Traian Iliescu, Honghu Liu외 1명|arXiv (Cornell University)|2017. 01. 04.

Model Reduction and Neural Networks참고 문헌 46인용 수 19

한 줄 요약

이 논문은 비선형 곱셈 노이즈가 있는 비선형 대류 지배 영역에서 스토케스틱 버거스 방정식에 대해 갈레르킨 감소 차수 모델(G-ROM)을 안정화하기 위해 레일리 정규화된 감소 차수 모델(L-ROM)을 제안한다. 대류 항에 공간 미분 필터를 적용함으로써 L-ROM은 부정확한 수치 진동을 크게 감소시키고 노이즈 강도 변화에 대한 강건성을 향상시켜, 다양한 노이즈 수준에서 표준 G-ROM보다 정확도와 안정성 면에서 뛰어나다.

ABSTRACT

In this paper, we study the numerical stability of reduced order models for convection-dominated stochastic systems in a relatively simple setting: a stochastic Burgers equation with linear multiplicative noise. Our preliminary results suggest that, in a convection-dominated regime, standard reduced order models yield inaccurate results in the form of spurious numerical oscillations. To alleviate these oscillations, we use the Leray reduced order model, which increases the numerical stability of the standard model by smoothing (regularizing) the convective term with an explicit spatial filter. The Leray reduced order model yields significantly better results than the standard reduced order model and is more robust with respect to changes in the strength of the noise.

연구 동기 및 목표

• 대류 지배 영역에서 스토케스틱 유체 흐름에 대한 표준 갈레르킨 감소 차수 모델(G-ROM)의 수치적 불안정성을 조사하기 위해.
• 스토케스틱 PDE에서 높은 대류 지배성과 노이즈로 인해 발생하는 G-ROM의 부정확한 진동을 해결하기 위해.
• 스토케스틱 시스템에 대한 안정화 방법으로서 레일리 정규화된 감소 차수 모델(L-ROM)의 성능을 평가하기 위해.
• 다양한 노이즈 강도에서 G-ROM과 L-ROM의 정확도와 강건성을 비교하기 위해.
• 두 모델이 노이즈 강도와 감소 기저 차원의 변화에 얼마나 민감한지 평가하기 위해.

제안 방법

• 스토케스틱 버거스 방정식의 스냅샷 데이터로부터 적절한 직교 분해(POD)를 사용하여 갈레르킨 감소 차수 모델(G-ROM)을 유도한다.
• 비선형 항을 명시적으로 스무딩하기 위해 대류 항에 미분 필터를 도입함으로써 레일리 정규화된 감소 차수 모델(L-ROM)을 구현한다.
• 선형 곱셈 노이즈가 있는 스트라토노비치 적분 방정식을 사용하여 난류 변동을 모델링한다.
• 모든 시뮬레이션에서 관측된 노이즈 수준 간 유사성에 기반해, POD 기저를 결정론적 경우(σ = 0)로 고정한다.
• L-ROM에서 대류 항을 정규화하기 위해 고정된 파rameter δ = 0.12를 갖는 미분 필터를 적용한다.
• 13개의 노이즈 강도(σ ∈ [0, 0.6])에 대해 3000개의 샘플을 사용한 앙상블 시뮬레이션을 수행하여 통계적 오차 지표를 계산한다.

실험 결과

연구 질문

• RQ1 표준 갈레르킨 ROM은 대류 지배 영역에서 스토케스틱 버거스 방정식에 대해 수치적 진동을 보이는가?
• RQ2 노이즈 존재 하에서 L-ROM은 표준 G-ROM에 비해 수치적 진동을 어떻게 완화하는가?
• RQ3 노이즈 강도가 증가함에 따라 G-ROM과 L-ROM의 강건성은 어떻게 변화하는가?
• RQ4 고정된 POD 기저(σ = 0에서 유도)를 사용할 경우, 다양한 노이즈 수준에서 ROM의 정확도에 어떤 영향을 미치는가?
• RQ5 스토케스틱 힘 작용 하에서 낮은 차원의 감소 공간에서도 L-ROM은 정확성과 안정성을 유지할 수 있는가?

주요 결과

• 표준 갈레르킨 ROM(G-ROM)은 대류 지배 영역에서 심각한 부정확한 수치적 진동을 유발하며, 노이즈 강도가 증가할수록 악화된다.
• L-ROM은 수치적 진동을 크게 감소시켜, 특히 낮은 감소 모델 차원(r = 10)에서 G-ROM보다 더 정확한 해를 제공한다.
• 모든 노이즈 수준에서 L-ROM은 상대 L² 오차의 평균과 분산이 훨씬 낮아, 노이즈 변화에 대한 우수한 강건성을 보인다.
• 그림 7의 오차 막대는 G-ROM 오차의 표준편차가 노이즈 강도 증가에 따라 급격히 증가하는 반면, L-ROM 오차는 안정적이고 좁은 범위 내에 유지됨을 보여준다.
• 3000개의 앙상블 샘플을 사용한 통계적 수렴 분석을 통해 오차 지표의 신뢰성이 확인되었으며, 보고된 오차 통계치의 신뢰성이 보장된다.
• 결정론적 경우(σ = 0)에서 유도한 고정된 POD 기저를 사용할 경우, 노이즈 경로별로 기저를 갱신하는 것과 비교해 오차 표준편차에 미미한 변화만을 보이며 유사한 결과를 얻는다.