[논문 리뷰] Accurate macroscale modelling of spatial dynamics in multiple dimensions

A. J. Roberts, Tony MacKenzie외 1명|arXiv (Cornell University)|2011. 03. 07.

Advanced Mathematical Modeling in Engineering참고 문헌 41인용 수 4

한 줄 요약

이 논문은 중심다중면 이론과 일반화된 갭-톱 메서드를 활용하여 두 개 이상의 공간 차원에서 반응-확산 시스템의 정확한 매크로스케일 모델링을 위한 엄밀하고 다스케일 프레임워크를 제시한다. 이는 미세스케일 시뮬레이터(예: PDE 또는 패치 기반 시뮬레이션)로부터 특별히 설계된 패치 간 결합 조건을 통해 일致하고 고차 정확도를 갖춘 이산 모델이 유도된다는 것을 입증하며, 이는 이론적 타당성과 계산 효율성을 모두 확보한다.

ABSTRACT

Developments in dynamical systems theory provides new support for the macroscale modelling of pdes and other microscale systems such as Lattice Boltzmann, Monte Carlo or Molecular Dynamics simulators. By systematically resolving subgrid microscale dynamics the dynamical systems approach constructs accurate closures of macroscale discretisations of the microscale system. Here we specifically explore reaction-diffusion problems in two spatial dimensions as a prototype of generic systems in multiple dimensions. Our approach unifies into one the modelling of systems by a type of finite elements, and the `equation free' macroscale modelling of microscale simulators efficiently executing only on small patches of the spatial domain. Centre manifold theory ensures that a closed model exist on the macroscale grid, is emergent, and is systematically approximated. Dividing space either into overlapping finite elements or into spatially separated small patches, the specially crafted inter-element/patch coupling also ensures that the constructed discretisations are consistent with the microscale system/PDE to as high an order as desired. Computer algebra handles the considerable algebraic details as seen in the specific application to the Ginzburg--Landau PDE. However, higher order models in multiple dimensions require a mixed numerical and algebraic approach that is also developed. The modelling here may be straightforwardly adapted to a wide class of reaction-diffusion PDEs and lattice equations in multiple space dimensions. When applied to patches of microscopic simulations our coupling conditions promise efficient macroscale simulation.

연구 동기 및 목표

• 다양한 차원에서 공간 동역학의 이론적으로 타당한 체계적인 매크로스케일 모델링 방법을 개발하는 것.
• 중심다중면 이론을 활용하여 유한요소 방법과 '방정식 기반' 패치 기반 모델링을 통합하는 공통 프레임워크를 수립하는 것.
• 매크로스케일 이산화와 기초가 되는 미세스케일 동역학(예: PDE 또는 시뮬레이터) 간의 고차 정확도 일致성을 확보하는 것.
• 공간적으로 분리된 패치에서의 소규모 미세시뮬레이션을 결합하여 대규모 시스템의 효율적 시뮬레이션을 가능하게 하는 것.
• 간격-톱 방법을 정확도와 안정성이 입증된 방식으로 두 개 이상의 차원으로 일반화하는 것.

제안 방법

• 중심다중면 이론을 사용하여 느린 다중면의 존재성과 그가 매크로스케일 모델로서의 닫힌 형태로 존재함을 증명한다.
• 공간적으로 분리된 패치 또는 겹치는 유한요소 간의 비국소적 결합 조건을 갖춘 일반화된 갭-톱 방법을 적용한다.
• 점진 전개와 대수적 연산을 통해 미세구조 동역학을 일치시켜 고차 정확도를 갖춘 이산화를 유도한다.
• 고차 폐쇄 모델에 필요한 복잡한 기호 계산을 처리하기 위해 컴퓨터 대수 기법을 활용한다.
• 미세구조 필드의 수치 계산과 기호 계수 표현을 조합하여 다차원에서의 고차 정확도를 달성한다.
• 매크로스케일 모델이 격자 간격 H → 0일 때 진짜 PDE로 수렴함을 보여줌으로써 일치성을 검증하며, 오차 한계는 O(δi²ᵖ + δj²ᵖ)이다.

실험 결과

연구 질문

• RQ1 중심다중면 이론은 두 개의 공간 차원에서 반응-확산 시스템의 매크로스케일 이산화를 정당화하기 위해 확장될 수 있는가?
• RQ2 기초가 되는 미세스케일 동역학과의 고차 정확도 일치를 보장하기 위해 패치 간 또는 요소 간의 결합 조건은 어떻게 설계할 수 있는가?
• RQ3 수반 시스템은 매크로스케일 모델의 정확도를 향상시키는 데 어떤 역할을 하는가?
• RQ4 순수 기호 방법이 다차원에서 실패할 경우, 대수적 기법과 수치 기법을 어떻게 조합하여 고차 모델을 구성할 수 있는가?
• RQ5 정확도와 효율성을 유지하면서 갭-톱 방법을 다차원으로 체계적으로 일반화할 수 있는가?

주요 결과

• 중심다중면 이론을 통해 유한 격자 간격 H에서 느린 다중면의 존재성이 엄밀히 입증되어 매크로스케일 모델이 유도되고 잘 정의됨을 보장한다.
• 제안된 결합 조건은 기초가 되는 PDE와의 고차 정확도 일치를 보장하며, 오차 항은 O(δi²ᵖ + δj²ᵖ)로 제한된다.
• 긴츠부르크-란다 PDE의 경우, 미세시뮬레이션이 작은 고립된 패치에 국한되어 있어도 일致한 매크로스케일 진화를 달성하며, 전체 미세스케일 시뮬레이션 대비 계산 비용은 약 16%로 감소한다.
• 수치적 기법과 기호 기법의 융합적 접근은 다차원에서의 고차 모델 구축을 가능하게 하여 순수 기호 방법의 한계를 극복한다.
• 이 방법은 다차원에서 다양한 반응-확산 PDE와 격자 방정식에 일반화 가능하다.
• 이 프레임워크는 명시적 폐쇄 방정식이 필요 없이도 복잡한 다스케일 시스템의 효율적이고 방정식 기반의 매크로스케일 시뮬레이션을 위한 이론적 기반을 제공한다.