[논문 리뷰] Postprocessing Galerkin method applied to a data assimilation algorithm: a uniform in time error estimate

Cecilia F. Mondaini, Edriss S. Titi|arXiv (Cornell University)|2016. 12. 21.

Advanced Mathematical Modeling in Engineering참고 문헌 34인용 수 8

한 줄 요약

이 논문은 2차원 나비에-스토크스 방정식에 대한 연속적 데이터 아울러션 알고리즘을 향상시키기 위해 후처리 갈레르킨 방법을 도입하며, 수치적 근사와 진짜 해 사이의 시간에 관계없이 일관된 오차 추정을 달성한다. 노칭 파ram터, 메쉬 해상도, 수치적 방법에 적합한 조건 하에서, 이 방법은 모든 시간에 걸쳐 참조 해로의 지수적 수렴을 보장하며, 푸리에 모드와 유한체적 평균을 포함한 광범위한 종류의 보간 연산자에 대해 유효하다.

ABSTRACT

We apply the Postprocessing Galerkin method to a recently introduced continuous data assimilation (downscaling) algorithm for obtaining a numerical approximation of the solution of the two-dimensional Navier-Stokes equations corresponding to given measurements from a coarse spatial mesh. Under suitable conditions on the relaxation (nudging) parameter, the resolution of the coarse spatial mesh and the resolution of the numerical scheme, we obtain uniform in time estimates for the error between the numerical approximation given by the Postprocessing Galerkin method and the reference solution corresponding to the measurements. Our results are valid for a large class of interpolant operators, including low Fourier modes and local averages over finite volume elements. Notably, we use here the 2D Navier-Stokes equations as a paradigm, but our results apply equally to other evolution equations, such as the Boussinesq system of Benard convection and other oceanic and atmospheric circulation models.

연구 동기 및 목표

• 2차원 나비에-스토크스 방정식에 대한 연속적 데이터 아울러션 알고리즘의 엄밀한 수치 분석을 제공하기 위해 후처리 갈레르킨 방법을 사용한다.
• 초기 조건이 알려지지 않은 경우에도 수치 근사와 진짜 해 사이의 시간에 관계없는 오차 경계를 확립한다.
• 데이터 아울러션 프레임워크의 적용 가능성을, 부스인스크 방정식 및 대기/해양 모델을 포함한 일반적인 소산성 진화 방정식으로 확장한다.
• 저주파 푸리에 모드와 국소 공간 평균을 포함한 다양한 보간 연산자에 대해 방법의 강건성을 검증한다.
• 초기 자료에 관계없이 오차 추정이 유한하고 시간에 따라 균일하게 감소함을 보장한다.

제안 방법

• 이 방법은 피드백 제어(노칭) 기반의 연속적 데이터 아울러션 알고리즘에 후처리 갈레르킨 기법을 적용하여, 군더미 스케일 측정값으로부터 세밀한 스케일 해를 재구성한다.
• 알고리즘은 보간 측정값 $ I_h(\textbf{v}) $ 와 $ I_h(\textbf{u}) $ 의 차이에 비례하는 피드백 항을 사용하며, 이는 이완 파ram터 $ \beta $ 를 통해 조절된다. 이를 통해 근사 해 $ \textbf{v} $ 는 참값 해 $ \textbf{u} $ 쪽으로 이동한다.
• 공간 이산화 방법은 유한차원 부분공간 위에 갈레르킨 사영을 수행한 후, 표준 갈레르킨 해를 초월해 정확도를 향상시키기 위한 후처리를 수행한다.
• 분석은 나비에-스토크스 방정식의 함수형 표현을 기반으로 하며, $ A $ (라플라시안) 및 $ B $ (비선형 대류 항) 연산자를 사용하고, $ L^2 $ 및 소볼레프 노름에서 에너지 추정을 포함한다.
• 보간 연산자 $ I_h $ 는 저주파 푸리에 모드와 유한체적 요소 위의 국소 평균을 포함하여 일반화되어, 광범위한 적용 가능성을 확보한다.
• 사전 추정과 에너지 방법을 사용하여 오차 $ \| \textbf{v} - \textbf{u} \|_{L^2} $ 의 경계를 유도하며, $ \beta $, $ h $, 수치 해상도에 적절한 조건 하에서 시간에 관계없이 오차가 감소하는 것을 보여준다.

실험 결과

연구 질문

• RQ1 후처리 갈레르킨 방법은 2차원 나비에-스토크스 방정식에 적용된 데이터 아울러션 알고리즘에 대해 시간에 관계없이 일관된 오차 추정을 제공할 수 있는가?
• RQ2 공간 후처리를 통한 노칭 기반 데이터 아울러션 기법을 사용할 경우, 수치 근사와 진짜 해 사이의 오차는 시간에 따라 어떻게 변화하는가?
• RQ3 진짜 해로의 지수적 수렴을 보장하기 위해 이완 파aram터 $ \beta $, 메쉬 해상도 $ h $, 수치적 방법에 어떤 조건이 필요한가?
• RQ4 오차 추정은 푸리에 모드 또는 유한체적 평균을 포함한 보간 연산자 $ I_h $ 의 선택에 얼마나 의존하는가?
• RQ5 제안된 방법은 2차원 나비에-스토크스 방정식을 초월한 다른 소산성 진화 방정식으로 확장할 수 있는가?

주요 결과

• 후처리 갈레르킨 방법은 수치 근사 $ \textbf{v} $ 와 진짜 해 $ \textbf{u} $ 사이에 시간에 관계없이 일관된 오차 추정을 제공하며, 오차는 시간에 따라 지수적으로 감소한다.
• 오차 경계는 초기 자료에 독립적이므로, $ \textbf{v} $ 의 초기 추측에 관계없이 강건함을 보장한다.
• 이 방법은 저주파 푸리에 모드와 유한체적 요소 위의 국소 공간 평균을 포함한 광범위한 종류의 보간 연산자에 대해 수렴성을 확보한다.
• 오차 추정은 군더미 메쉬의 해상도 $ h $, 이완 파aram터 $ \beta $, 스펙트럴 간격 $ \lambda_{N+1} $ 에 따라 달라지며, 추정치는 $ \mathcal{O}(L^{1/2}/h^2 \cdot \lambda_{N+1}^{-1/4}) $ 의 비율로 척도화된다.
• 분석은 피드백 제어 항 $ -\beta(I_h(\textbf{v}) - I_h(\textbf{u})) $ 가 수치 해를 효과적으로 안정화시키고 참값 해로 이동시킴을 확인한다.
• 결과는 부스인스크 방정식 및 대규모 대기 및 해양 모델과 같은 다른 시스템으로도 확장 가능하며, 넓은 적용 가능성을 보여준다.