[论文解读] Holistic discretisation ensures fidelity to dynamics in two spatial dimensions

Tony MacKenzie, A. J. Roberts|ArXiv.org|Apr 6, 2009

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动力学的一致性。
• 探索在解析解不可得时，构建更高阶亚网格场的混合数值-代数方法。
• 以二维实Ginzburg–Landau方程为基准，验证该方法，展示其预测精度与误差控制能力。

提出的方法

• 将二维区域划分为重叠的正方形有限元，以实现对亚网格动力学的多尺度解析。
• 应用由γ参数化的非局部单元间耦合条件（4）–（5），以确保与PDE的一致性，并支持中心流形理论。
• 利用中心流形理论证明存在一个慢流形，从而提供一个精确封闭的、渐近精确的离散模型。
• 使用计算机代数推导离散模型的二阶与三阶渐近近似，特别是利用伴随算子提取额外一阶精度。
• 对于更高阶模型，引入一种保持系数代数表达式的数值亚网格场计算方法，从而实现高效模型构建。
• 采用混合数值与代数方法，构建误差缩放为O(h^{2p−2})的离散模型，其中p阶精度，通过一致性定理加以验证。

实验结果

研究问题

• RQ1 能否将一维PDE中使用的综合离散化框架推广至二维空间，同时保持对底层动力学的保真性？
• RQ2 在二维中，何种重叠有限元之间的耦合条件既能确保与PDE的一致性，又能保证中心流形理论的适用性？
• RQ3 为何纯代数方法在二维中无法构建更高阶综合模型？该局限性如何克服？
• RQ4 低分辨率亚网格场在多大程度上能准确预测二维反应-扩散系统中的宏观尺度动力学？
• RQ5 数值亚网格解析与代数系数表达式的结合，如何实现高阶、一致的离散模型？

主要发现

• 如引理5所证明，当相互作用截断至O(γ^p)时，综合离散化方法可确保离散模型与原始PDE的一致性，误差为O(h^{2p−2})。
• 中心流形理论保证了慢流形的存在，且其吸引速度呈指数级，为离散模型的精度提供了理论基础。
• 通过在耦合单元上对扩散算子的伴随算子进行操作，可提取额外一阶精度，从而为Ginzburg–Landau方程得到三阶离散模型（17）。
• 即使亚网格分辨率较粗，仍足以准确预测宏观尺度动力学，原因在于亚网格结构与单元间相互作用的忠实解析。
• 在二维中构建高阶模型时，混合数值-代数方法至关重要，因为纯代数方法在低阶以上即失效。
• 该方法计算高效，计算机代数系统Reduce在模型构建速度上比Mathematica快一个数量级以上。