[논문 리뷰] Multiscale finite elements through advection-induced coordinates for transient advection-diffusion equations

Konrad Simon, Jörn Behrens|arXiv (Cornell University)|2018. 02. 21.

Advanced Mathematical Modeling in Engineering인용 수 3

한 줄 요약

이 논문은 일시적 확산-대류 방정식에서 코arse 그리드를 사용한 기후 시뮬레이션을 위한 업스케일링을 향상시키기 위해 대류 유도 좌표를 사용하는 새로운 다중스케일 유한요소 방법을 제안한다. 대류 효과가 강한 코어스 엘리먼트 경계 근처의 영향을 줄이기 위해 좌표를 변환함으로써, 진동하는 속도와 확산 조건을 가진 1차원 문제에서 높은 정확도를 달성하며, 최적 수렴 속도와 스케일 분리 조건 하에서의 강건성을 입증한다.

ABSTRACT

Long simulation times in climate sciences typically require coarse grids due to computational constraints. Nonetheless, unresolved subscale information significantly influences the prognostic variables and can not be neglected for reliable long term simulations. This is typically done via parametrizations but their coupling to the coarse grid variables often involves simple heuristics. We explore a novel up-scaling approach inspired by multi-scale finite element methods. These methods are well established in porous media applications, where mostly stationary or quasi stationary situations prevail. In advection-dominated problems arising in climate simulations the approach needs to be adjusted. We do so by performing coordinate transforms that make the effect of transport milder in the vicinity of coarse element boundaries. The idea of our method is quite general and we demonstrate it as a proof-of-concept on a one-dimensional passive advection-diffusion equation with oscillatory background velocity and diffusion.

연구 동기 및 목표

• 코어스 그리드를 사용한 장기 기후 시뮬레이션에서 해상도가 낮은 서브그리드 스케일 과정의 문제를 해결하기 위해.
• 직관적 파rametrization를 피하는 수학적으로 일관된 업스케일링 프레임워크를 개발하기 위해.
• 기후 모델링에서 흔한 대류 지배 문제에 다중스케일 유한요소 방법(MsFEM)을 적용하기 위해.
• 대류 유도 왜곡을 줄이기 위해 좌표를 변환함으로써 일시적 대류-확산 방정식의 안정성과 정확도를 향상시키기 위해.
• 진동하는 배경 속도와 확산 조건을 가진 1차원 문제에서 방법의 효과성을 입증하기 위해.

제안 방법

• 대류 특성에 기반한 좌표 변환을 도입하여 코어스 엘리먼트 경계 근처의 운반 효과를 완화한다.
• 세밀한 스케일 행동을 반영하면서도 안정성을 유지하는 변환된 좌표계에서 다중스케일 기저 함수를 구성한다.
• 대류 유도 좌표에서 국소 문제를 해결하여 유도된 비다항 기저 함수를 사용한 갈레르킨 형식을 적용한다.
• 진동하는 속도와 확산 계수를 가진 1차원 일시적 대류-확산 방정식에 방법을 적용한다.
• 좌표 변환을 통해 기저 함수를 주로 흐름 방향에 맞추어 질량 보존과 안정성을 확보한다.
• 수렴성과 정확도 평가를 위해 다양한 그리드 해상도와 확산 주기를 가진 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 대류 특성에 기반한 좌표 변환은 대류 지배 일시적 문제에서 다중스케일 유한요소 방법의 안정성과 정확도를 향상시키는가?
• RQ2 확산 계수가 코어스 그리드에 비해 빠르게 진동할 경우, 제안된 방법은 표준 유한요소 방법보다 어떻게 성능을 발휘하는가?
• RQ3 배경 속도에 주된 평균이 없거나 시간에 따라 변하지 않는 영점이 포함되어 있을 경우, 방법의 한계는 무엇인가?
• RQ4 이 방법은 다차원 또는 보존 형태의 대류-확산 방정식으로 일반화될 수 있는가?
• RQ5 확산 계수의 진동 주파수가 변화하고 그리드 해상도가 증가할 때, 방법의 수렴 행동은 어떠한가?

주요 결과

• L² 및 L∞ 노름에서 최적 수렴 속도를 달성하며, 그리드 해상도를 두 배로 늘일 경우 상대 오차가 한 계급 감소한다.
• k=10일 때, N=24에서 L² 오차는 2.732×10⁻⁴에서 N=384로 증가하면서 3.939×10⁻⁷로 감소하여 2차 수렴을 나타낸다.
• 고주파 확산 진동 조건에서도 잘 작동하여 스케일 분리 조건 하에서 강건성을 입증한다.
• 배경 속도에 시간에 따라 변하지 않는 영점이 존재할 경우 안정성 문제가 발생하며, 이는 특성 붕괴를 유도하는 안착점(아트랙터)을 생성하기 때문이다.
• 보존 형태의 방정식에서는 방법이 실패한다. 이는 유량 항을 대류 및 반응 항으로 분해할 필요가 있어 각각 다른 업스케일링 메커니즘이 요구되기 때문이다.
• 한계가 있음에도 불구하고, 이 방법은 관련 시나리오에서 뛰어난 성능을 보이며, 향후 반라그랑주 방법과 같은 하이브리드 접근법을 통한 확장에 기초를 마련한다.