[논문 리뷰] Uniform in Time Error Estimates for a Finite Element Method Applied to a Downscaling Data Assimilation Algorithm for the Navier--Stokes Equations

García-Archilla, Julia Novo외 1명|arXiv (Cornell University)|2018. 07. 23.

Advanced Mathematical Modeling in Engineering참고 문헌 33인용 수 1

한 줄 요약

이 논문은 2차원 및 3차원 나비에-스토크스 방정식에 대해 연속적인 다운스케일링 데이터 아울러션 알고리즘에 적용된 유한요소법에 대해 시간에 관계없이 일관된 오차 추정을 제시한다. 혼합 유한요소를 사용하거나 그라드-디브 안정화를 적용함으로써, 이 방법은 점도의 역수에 의존하지 않는 최적 수렴 속도와 경계를 확보하여, 굵은 스케일 측정치로부터 참조 해의 안정적인 장기적 근사가 가능하다.

ABSTRACT

In this paper we analyze a finite element method applied to a continuous downscaling data assimilation algorithm for the numerical approximation of the two- and three-dimensional Navier--Stokes equations corresponding to given measurements on a coarse spatial scale. For representing the coarse mesh measurements we consider different types of interpolation operators including a Lagrange interpolant. We obtain uniform-in-time estimates for the error between a finite element approximation and the reference solution corresponding to the coarse mesh measurements. We consider both the case of a plain Galerkin method and a Galerkin method with grad-div stabilization. For the stabilized method we prove error bounds in which the constants do not depend on inverse powers of the viscosity. Some numerical experiments illustrate the theoretical results.

연구 동기 및 목표

• 나비에-스토크스 방정식에 대한 다운스케일링 데이터 아울러션 알고리즘의 유한요소 근사에 대해 시간에 관계없이 일관된 오차 경계를 확립한다.
• 실제 유한요소 가정 하에 그라드-디브 안정화 유무에 따라 갈레르킨 방법의 수렴 행동을 분석한다.
• 오차 경계가 점도 매개변수 ν에 어떻게 의존하는지, 특히 저점도(고레이놀즈 수) 영역에서 어떻게 되는지를 조사한다.
• 저점도 측정치에 대해 라그랑주 보간을 사용할 때 수렴이 유지되는 조건을 규명한다.
• 이전 연구에서 가정된 뉴징 매개변수 β에 대한 제한적인 상한을 제거하여 분석의 강건성을 향상시킨다.

제안 방법

• IH(v) − IH(u)에 대한 뉴징 항 −β(IH(v) − IH(u))을 나비에-스토크스 방정식에 추가하여 연속적인 다운스케일링 데이터 아울러션 알고리즘을 수립한다. 여기서 IH는 저점도 스케일 보간 연산자이다.
• 속도와 압력 근사에 대해 인프-스터블한 혼합 유한요소를 사용하는 공간에 대한 반연속 유한요소 방법을 적용한다.
• 두 가지 변형을 고려한다: 표준 갈레르킨 방법과 그라드-디브 안정화를 적용한 갈레르킨 방법으로, 질량 보존성과 안정성을 향상시킨다.
• 저점도 관측치를 표현하기 위해 보간 연산자(라그랑주 보간 포함)를 사용하며, H는 저점도 메쉬 크기를 나타낸다.
• 에너지 추정과 그로나발 유형 부등식을 적용하여, 유한요소 해와 참조 해 사이의 시간에 관계없이 일관된 오차 경계를 도출한다.
• 이론적 수렴 속도를 검증하고 β, ν, H/h, 안정화의 영향을 평가하기 위해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 나비에-스토크스 방정식에 대한 다운스케일링 데이터 아울러션 알고리즘의 유한요소 근사에 대해 시간에 관계없이 일관된 오차 추정을 유도할 수 있는가?
• RQ2 그라드-디브 안정화가 오차 경계가 점도 ν에 어떻게 의존하는지, 특히 저점도 영역에서 어떻게 되는가?
• RQ3 라그랑주 보간을 사용할 때 수렴을 보장하기 위해 저점도-세밀점도 메쉬 비율 H/h이 수행해야 할 역할은 무엇인가?
• RQ4 이전 연구에서 가정한 바와 같이 뉴징 매개변수 β에 상한을 두지 않더라도 분석이 여전히 유효한가?
• RQ5 ν, β, 메쉬 정밀도 변화 조건 하에서 이론적 수렴 속도가 수치 실험으로 확인되는가?

주요 결과

• 그라드-디브 안정화 없이 갈레르킨 방법을 사용할 경우, 매끄러운 해에 대해 오차 경계는 최적 수렴 속도 O(h^3)를 확보하며, 이는 유한요소 공간의 최선 근사 오차와 일치한다.
• 그라드-디브 안정화를 적용할 경우, 오차 경계는 점도 ν의 역수에 의존하지 않으며, 이는 고레이놀즈 수 유동에서 매우 중요하다.
• 수치 실험 결과, 작은 ν 값에서 안정화된 방법은 ν에 관계없이 O(h^2) 수렴을 달성하는 반면, 안정화되지 않은 방법은 테스트한 h 값에서 수렴하지 못함을 확인하였다.
• 라그랑주 보간을 사용할 경우 H/h가 유한하게 유지되어야 수렴이 보장되며, H를 고정한 채로 h→0으로 갈 경우 수렴하지 않는 행동을 보인다.
• 더 큰 뉴징 매개변수 β는 초기 오차의 감쇠를 가속화하며, 이 행동은 시간에 따라 변화하는 ∥IH(eh)∥₀ / ∥eh∥₀의 하한 경계를 통해 분석으로 설명된다.
• 이전 연구에서 가정한 β에 대한 제한적인 상한을 제거하여 알고리즘 설계의 유연성을 높였다.