[论文解读] Singular limits of reaction diffusion equations and geometric flows with discontinuous velocity

Cecilia De Zan, Pierpaolo Soravia|arXiv (Cornell University)|May 23, 2019

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 20被引用 5

一句话总结

本文研究具有不连续速度场的双稳态反应-扩散方程的奇异极限，证明解在时间全局收敛至移动界面外的稳定平衡态。该界面以不连续法向速度演化，其演化由几何哈密顿-雅可比方程的粘性解控制，通过粘性理论将水平集方法推广至非光滑速度函数的情形。

ABSTRACT

We consider the singular limit of a bistable reaction diffusion equation in the case when the velocity of the traveling wave solution depends on the space variable and converges to a discontinuous function. We show that the family of solutions converges to the stable equilibria off a front propagating with a discontinuous velocity. The convergence is global in time by applying the weak geometric flow uniquely defined through the theory of viscosity solutions and the level-set equation.

研究动机与目标

• 分析具有依赖于 ε 的双稳非线性项及不连续波速的反应-扩散方程在时间全局范围内的渐近行为。
• 将几何流框架扩展至极限速度函数在超曲面上不连续的情形。
• 在速度梯度在 ε → 0 时趋于无穷大的情况下，建立解在移动界面外收敛至稳定平衡态的结果。
• 通过粘性解定义并表征极限前驱面的运动，其运动由具有不连续哈密顿量的哈密顿-雅可比方程控制。
• 利用比较原理与逼近技术，将水平集方法与弱几何流理论推广至不连续速度场的情形。

提出的方法

• 使用粘性解理论定义并确保即使在速度函数不连续时，几何流仍具有唯一性。
• 应用水平集方法，将演化界面表示为哈密顿-雅可比方程粘性解的零水平集。
• 通过行波解形貌与障碍函数构造上下解，以控制前驱面的传播。
• 通过光滑函数序列逼近不连续速度函数，以应用文献中已知的比较原理。
• 采用 Barles-Souganidis 的奇异极限框架，通过与光滑演化局部比较，将其适配于处理不连续哈密顿量的情形。
• 使用截断函数与正则化距离函数，将局部上下解扩展至整个时空区域的全局解。

实验结果

研究问题

• RQ1 能否在时间全局范围内表征具有不连续速度的反应-扩散方程奇异极限？
• RQ2 当速度函数不连续时，界面如何演化？其运动能否由几何PDE描述？
• RQ3 极限速度的上半连续包与下半连续包在定义初始前驱面时起什么作用？
• RQ4 当速度不连续时，粘性解是否仍能提供唯一且稳定的前驱面运动描述？
• RQ5 在何种条件下，反应-扩散方程的解会收敛至界面外的稳定平衡态？

主要发现

• 解族 $ u^\nu $ 局部一致收敛至 $ 1 $，当且仅当几何方程的粘性解 $ u $ 为正；在 $ u $ 为负的区域则收敛至 $ -1 $。
• 极限前驱面 $ \Gamma_t $ 以法向速度 $ -\alpha(x) $ 移动，其中 $ \alpha $ 为不连续的极限速度，其由哈密顿-雅可比方程 $ u_t + \alpha(x)|Du| = 0 $ 的粘性解定义。
• 初始前驱面 $ \Gamma_0 $ 由 $ \{ x : \alpha_*(x)/2 \leq g(x) \leq \alpha^*(x)/2 \} $ 给出，其中 $ \alpha_* $ 与 $ \alpha^* $ 分别为 $ \alpha $ 的下半连续包与上半连续包。
• 即使当 $ \|Dc^\varepsilon\|_\infty $ 与 $ \|D^2c^\varepsilon\|_\infty $ 在 $ \varepsilon \to 0 $ 时趋于无穷大，解在前驱面外仍收敛至稳定平衡态。
• 几何流通过粘性解唯一定义，且水平集演化在解的递增变换下保持不变，从而确认其几何本质。
• 证明依赖于构造逼近前驱面运动的上下解，方法包括在全时空区域内使用正则化行波形貌与截断技术。