[논문 리뷰] Holistic discretisation ensures fidelity to dynamics in two spatial dimensions

Tony MacKenzie, A. J. Roberts|ArXiv.org|2009. 04. 06.

Advanced Mathematical Modeling in Engineering참고 문헌 28인용 수 5

한 줄 요약

이 논문은 두 차원 반응-확산 PDE에 대한 통합 이산화 방법을 제안하며, 겹치는 유한 요소와 국소 외 상호작용 조건을 통해 미세 구조 과정을 체계적으로 해소함으로써 원래의 역학에 대한 충실도를 보장한다. 이 방법은 중심다양체 이론과 일致성 분석을 조합하여 渐近적으로 정확하고 고차수의 이산 모델을 유도하며, 해석적 방법이 실패하는 곳에서 수치적 미세 구조 해상도를 통해 고차수 정확도를 달성한다.

ABSTRACT

Developments in dynamical systems theory provides new support for the discretisation of \pde{}s and other microscale systems. By systematically resolving subgrid microscale dynamics the new approach constructs asymptotically accurate, macroscale closures of discrete models of the \pde. Here we explore reaction-diffusion problems in two spatial dimensions. Centre manifold theory ensures that slow manifold, holistic, discretisations exists, are quickly attractive, and are systematically approximated. Special coupling of the finite elements ensures that the resultant discretisations are consistent with the \pde to as high an order as desired. Computer algebra handles the enormous algebraic details as seen in the specific application to the Ginzburg--Landau equation. However, higher order models in 2D appear to require a mixed numerical and algebraic approach that is also developed. Being driven by the residuals of the equations, the modelling here may be straightforwardly adapted to a wide class of reaction-diffusion differential and lattice equations in multiple space dimensions.

연구 동기 및 목표

• 원래 PDE의 역학을 유지하는 두 차원 반응-확산 PDE에 대한 체계적이고 渐近적으로 정확한 매크로스케일 이산화 모델을 개발하기 위해.
• 일차원에서의 통합 이산화 기법을 다차원 공간으로 확장하는 데 도전하는 데서, 고차수 모델의 대수적 구성이 비가능해지는 문제를 해결하기 위해.
• 고차수 정확도를 유지하는 상호작용 조건을 유도함으로써 이산 모델이 마이크로스케일 PDE 역학과 일致하도록 보장하기 위해.
• 해석적 해가 얻어지지 않을 경우 고차수 미세 구조 필드를 구성하기 위한 수치-대수적 융합 접근법을 탐색하기 위해.
• 기준으로 사용된 2D 실 Ginzburg–Landau 방정식을 통해 방법을 검증하여 예측 정확도와 오차 제어 능력을 입증하기 위해.

제안 방법

• 미세 구조 역학의 다중스케일 해상도를 위해 2D 영역을 겹치는 정사각형 유한 요소로 나누기.
• γ로 매개변수화된 비국소적 요소 간 상호작용 조건 (4)–(5)를 적용하여 PDE와의 일치성과 중심다양체 이론의 적용 가능성을 확보하기.
• 중심다양체 이론을 활용하여 정확히 닫혀 있고 渐近적으로 정확한 이산 모델을 제공하는 느린 다류가 존재함을 정당화하기.
• 컴퓨터 대수를 활용하여 이산 모델의 2차 및 3차 渐近적 근사식을 유도하며, 특히 수반 연산자를 사용해 추가로 한 단계의 정확도를 추출하기.
• 고차수 모델을 위해 대수적 표현식을 유지하는 수치적 미세 구조 필드 계산 방법을 도입하여 효율적인 모델 구축 가능하게 하기.
• 일치성 정리에 의해 검증된 혼합 수치-대수적 접근법을 통해 오차 척도 O(h^{2p−2})를 갖는 p차수 정확도의 이산 모델을 구성하기.

실험 결과

연구 질문

• RQ1 일차원 PDE에서 사용된 통합 이산화 프레임워크를 두 차원 공간으로 일반화할 수 있는가? 이때 원래의 역학에 대한 충실도는 유지되는가?
• RQ2 겹치는 유한 요소 간의 어떤 상호작용 조건이 2D에서 PDE와의 일치성과 중심다양체 이론의 적용 가능성을 동시에 보장하는가?
• RQ3 순수한 대수적 방법이 왜 두 차원에서 고차수 통합 모델을 구성하는 데 실패하는가? 이 제약은 어떻게 극복할 수 있는가?
• RQ4 저해상도의 미세 구조 필드가 2D 반응-확산 시스템에서 매크로스케일 역학을 얼마나 정확하게 예측할 수 있는가?
• RQ5 수치적 미세 구조 해상도와 대수적 계수 표현식의 조합이 고차수이고 일치하는 이산 모델을 어떻게 가능하게 하는가?

주요 결과

• Corollary 5에 의해 증명된 linely, 상호작용을 O(γ^p)로 잘라내었을 때, 결과적으로 매크로스케일 이산 모델이 원래 PDE와 오차 O(h^{2p−2})로 일치함을 보장한다.
• 중심다양체 이론은 지수적으로 빠르게 유도되는 느린 다류의 존재를 보장하여 이산 모델의 정확도에 대한 이론적 기반을 제공한다.
• 결합된 요소에서 확산 연산자의 수반 연산자를 활용하여 추가로 한 단계의 정확도를 추출할 수 있었으며, 이는 Ginzburg–Landau 방정식에 대해 3차 이산 모델 (17)을 도출하게 하였다.
• 미세 구조 구조와 요소 간 상호작용이 충실하게 해상화되므로, 조밀한 미세 구조 해상도가 아니더라도 매크로스케일 역학을 정확히 예측할 수 있다.
• 순수한 대수적 방법은 2차원에서 고차수 모델을 구성하는 데에는 실패하므로, 고차수 모델을 구성하기 위해 수치-대수적 융합 접근법이 필수적이다.
• 컴퓨터 대수 시스템 Reduce는 모델 구축 속도에서 Mathematica를 한 단계 이상 뛰어넘는 등 계산적으로 효율적이다.