[论文解读] A quasinonlocal coupling method for nonlocal and local diffusion models

Qiang Du, Xingjie Helen Li+2 more|arXiv (Cornell University)|Apr 2, 2017

Numerical methods in engineering参考文献 35被引用 4

一句话总结

本文提出了一种拟非局部耦合方法，通过几何重构在一维中无缝融合非局部与局部扩散模型，确保通量平衡、能量守恒及最大值原理。该方法消除了界面不一致性与边界伪影，首次为非局部到局部耦合建立了严格的适定性、稳定性及建模误差估计。

ABSTRACT

In this paper, we extend the idea of "geometric reconstruction" to couple a nonlocal diffusion model directly with the classical local diffusion in one dimensional space. This new coupling framework removes interfacial inconsistency, ensures the flux balance, and satisfies energy conservation as well as the maximum principle, whereas none of existing coupling methods for nonlocal-to-local coupling satisfies all of these properties. We establish the well-posedness and provide the stability analysis of the coupling method. We investigate the difference to the local limiting problem in terms of the nonlocal interaction range. Furthermore, we propose a first order finite difference numerical discretization and perform several numerical tests to confirm the theoretical findings. In particular, we show that the resulting numerical result is free of artifacts near the boundary of the domain where a classical local boundary condition is used, together with a coupled fully nonlocal model in the interior of the domain.

研究动机与目标

• 解决非局部到局部扩散模型耦合中的界面不一致性与边界层伪影问题。
• 开发一种保持关键物理性质（通量平衡、能量守恒及最大值原理）的耦合框架。
• 为耦合的非局部-局部问题建立适定性与建模误差估计。
• 提供一种数值离散化方法，避免在使用经典狄利克雷边界条件时在边界附近产生虚假伪影。
• 将原子尺度到连续介质方法中的几何重构思想拓展至非局部扩散模型。

提出的方法

• 该方法采用几何重构，将非局部扩散算子与经典局部PDE在一维中耦合，确保界面处的一致性。
• 通过一种拟非局部算子实现耦合，该算子保持线性一致性并满足最大值原理。
• 利用拟非局部版本的庞加莱不等式，证明了适定性，确保解的存在性与唯一性。
• 通过最大值原理估计建模误差，表明当Horizon参数 δ → 0 时解收敛至局部极限。
• 提出一种一阶有限差分格式用于数值离散化，保持连续模型的物理性质。
• 该方法用经典狄利克雷条件替代体积边界条件，消除了全非局部模型中因体积约束引起的非物理边界层。

实验结果

研究问题

• RQ1 一种非局部到局部耦合方法能否同时保持通量平衡、能量守恒及最大值原理？
• RQ2 如何消除非局部-局部耦合中扩散模型的界面不一致性和鬼魅力？
• RQ3 当 δ → 0 时，耦合的非局部-局部解与局部极限之间的建模误差是多少？
• RQ4 几何重构框架能否被适配至非局部扩散模型以确保物理一致性？
• RQ5 当使用经典狄利克雷边界条件时，所提出的耦合方法在数值上表现如何？

主要发现

• 所提出的拟非局部耦合方法确保了通量平衡、能量守恒及最大值原理——这些性质在现有耦合方法中均未满足。
• 通过拟非局部庞加莱不等式，严格证明了耦合问题的适定性。
• 利用最大值原理估计了耦合解与局部极限之间的建模误差，表明当 δ → 0 时解收敛。
• 数值测试确认了在 L∞ 范数下解与梯度均呈现一阶收敛，收敛阶为 0.99–1.00。
• 该方法成功消除了在使用经典狄利克雷边界条件时，由全非局部模型中体积约束引起的虚假边界层。
• 局部-非局部-局部耦合公式能准确捕捉源项中的奇异性行为，同时在非局部区域两侧与局部解保持一致。