[논문 리뷰] A discrete data assimilation scheme for the solutions of the 2D Navier-Stokes equations and their statistics

Ciprian Foiaş, Cecilia F. Mondaini외 1명|arXiv (Cornell University)|2016. 02. 18.

Meteorological Phenomena and Simulations참고 문헌 27인용 수 17

한 줄 요약

이 논문은 노이즈가 섞인 시간 이산 측정값을 나우딩 기반 피드백 제어를 통해 통합하는 2차원 나비에-스토크스 방정식을 위한 이산 데이터 통합 기법을 제안한다. 나우딩 파aram터, 메쉬 크기, 시간 간격에 대한 조건 하에서 측정 오차에 비례하는 항을 제외하고는 근사 해가 참값으로 지수 수렴함을 입증한다.

ABSTRACT

We adapt a previously introduced continuous in time data assimilation (downscaling) algorithm for the 2D Navier-Stokes equations to the more realistic case when the measurements are obtained discretely in time and may be contaminated by systematic errors. Our algorithm is designed to work with a general class of observables, such as low Fourier modes and local spatial averages over finite volume elements. Under suitable conditions on the relaxation (nudging) parameter, the spatial mesh resolution and the time step between successive measurements, we obtain an asymptotic in time estimate of the difference between the approximating solution and the unknown reference solution corresponding to the measurements, in an appropriate norm, which shows exponential convergence up to a term which depends on the size of the errors. A stationary statistical analysis of our discrete data assimilation algorithm is also provided.

연구 동기 및 목표

• 시간 이산적이고 오차가 포함된 측정값을 사용하여 2D 나비에-스토크스 방정식에 실용적인 데이터 통합 기법을 개발하기 위해.
• 연속 시간 나우딩 알고리즘을 현실적인 시간 이산적 환경과 관측 오차를 고려하여 확장하기 위해.
• 일반적인 보간 연산자 하에서 근사 해가 참값으로 수렴하는 것을 분석하기 위해.
• 정적 통계 분석을 포함하여 이산 데이터 통합 알고리즘의 정적 성질을 분석하고, 불변 측도 수렴을 제공하기 위해.
• 해의 근사 오차를 측정 오차, 메쉬 해상도, 시간 간격의 함수로 정량화하기 위해.

제안 방법

• 근사 해의 거시적 스케일이 기준 해의 것과 일치하도록 유도하는 피드백 제어(나우딩) 항을 나비에-스토크스 방정식에 도입한다.
• 유한 체적 요소 위에서 국소 공간 평균을 통해 정의된 보간 연산자 $ I_h $ 를 통해 이산 측정값을 통합한다.
• 측정 오차를 유계 벡터 $ \varepsilon_{n,j} $ 로 모델링하며, $ |\varepsilon_{n,j}| \leq \varepsilon $ 를 만족하고, 이들이 해 오차에 미치는 영향을 분석한다.
• 기저 함수 $ \psi_j $ 를 사용한 속도 장의 유한 차원 사영을 적용하며, $ I_h(\varphi) = \sum \overline{\varphi}_j \psi_j $ 를 만족한다. 여기서 $ \overline{\varphi}_j $ 는 국소 평균이다.
• 에너지 추정과 보간 성질을 활용하여 참 해와 근사 해의 차이에 대한 $ L^2 $ 및 $ H^1 $ 노름에서의 오차 추정을 유도한다.
• 이완 파aram터 $ \beta $, 메쉬 크기 $ h $, 시간 간격 $ \Delta t $ 에 대한 조건 하에서 수렴 결과를 확립하며, 이는 $ \varepsilon / h $ 에 비례하는 항을 제외하고 지수 수렴을 보장한다.

실험 결과

연구 질문

• RQ1 노이즈가 섞인 측정값을 사용하는 이산 시간 데이터 통합 기법이 2D 나비에-스토크스 방정식의 참값으로 지수 수렴을 달성할 수 있는가?
• RQ2 측정 오차와 시간 이산화가 나우딩 기반 데이터 통합 프레임워크 내에서 근사 해의 정확도에 어떤 영향을 미치는가?
• RQ3 나우딩 파aram터, 메쉬 크기, 시간 간격에 어떤 조건이 성립해야 기법의 안정성과 수렴성이 보장되는가?
• RQ4 이산 데이터 통합 기법이 참 시스템의 통계적 성질, 예를 들어 불변 측도를 유지하는가?
• RQ5 해의 근사 오차를 측정 오차와 공간 해상도의 함수로 유계로 제시할 수 있는가?

주요 결과

• 근사 해는 $ H^1 $ 노름에서 참값으로 지수 수렴하며, 이는 $ \varepsilon / h $ 에 비례하는 항을 제외하고는 성립한다. 여기서 $ \varepsilon $ 는 최대 측정 오차이고 $ h $ 는 메쉬 크기이다.
• $ L^2 $ 노름에서의 오차 bound는 $ |\eta_n|_{L^2} \leq \varepsilon |\Omega|^{1/2} $ 이고, $ H^1 $ 노름에서는 $ \|\eta_n\|_{H^1} \leq c \frac{\varepsilon}{h} |\Omega|^{1/2} $ 를 만족하며, $ c = C_0(2C_1)^{1/2} $ 이다.
• 에너지 추정에서 유도된 특정 조건을 만족할 경우, 이완 파aram터 $ \beta $, 메쉬 크기 $ h $, 시간 간격 $ \Delta t $ 가 지수 수렴을 보장한다.
• 정적 통계 분석 결과, 적절한 조건 하에서 이산 통합 기법의 불변 측도가 참 시스템의 불변 측도로 수렴함을 보였다.
• 측정 오차에 대한 체계적 오차에 대해 기법이 강건하며, 저주파 수식 모드와 국소 공간 평균을 포함한 일반적인 보간 연산자에 적용 가능하다.
• 지역 평균과 분할 정수 기저 함수를 기반으로 한 일관성 있는 보간 연산자 $ I_h $ 를 구성함으로써 이론적 오차 유계가 검증되었다.