退化交叉扩散下不连续入侵种群的自相似平滑及移动界面不存在性

《The ANZIAM Journal》：SELF-SIMILAR SMOOTHING OF A DISCONTINUITY BY DEGENERATE CROSS-DIFFUSION

　　

研究人员考虑一个理想化交叉扩散（cross-diffusion）模型，用于描述生物入侵，其中入侵种群的扩散受到驻地种群抑制，而驻地种群又被入侵种群降解。研究人员建立并数值求解了一个问题，该问题描述初始为分段常数的入侵种群的自相似平滑（self-similar smoothing）。在极限情形下，初始不连续面前方的高度消失，传播速度也消失，但仅以对数方式缓慢消失。该结果由匹配渐近分析（matched asymptotic analysis）予以证实。该极限的奇异性表明，该系统不允许存在表现出移动前沿或界面（moving front or interface）的紧支撑解（compactly supported solutions），这对模拟同样具有该退化扩散（degenerate diffusion）机制的更复杂模型具有重要意义。

研究背景与问题。在物理和生物扩散现象的数学建模中，退化扩散（degenerate diffusion）是技术上极具挑战性且颇受关注的问题。典型例子为多孔介质方程（porous medium equation），其扩散率v^m在v=0处消失；当m>0时，从分段阶跃初值出发可出现解为零的奇异点，即界面或自由边界（interface/free boundary），并可用弱解（weak solution）定义。向多孔介质方程加入逻辑斯蒂反应项后得到多孔Fisher方程，它是经典Fisher–Kolmogorov–Petrovsky–Piskunov（Fisher–KPP）问题的退化扩散版本；其行波（travelling wave）解具有界面。近年来，多物种退化反应–扩散系统被用于生物细胞侵袭与相互作用，其中交叉扩散（cross-diffusion）使一个物种影响另一物种的扩散率。一类两物种模型包含非运动驻地细胞种群?和扩散性入侵细胞种群v，?抑制v扩散，v降解局部?并间接提高自身扩散率，非无量纲形式如(1.3)。已有研究将此类系统用于酸介导肿瘤生长近似的Gatenby–Gawlinski模型，以及忽略驻地种群恢复或考虑对承载能力影响的细胞侵袭模型。此类系统在驻地细胞达到容量时出现退化。一个与细胞侵袭相关的初值配置是驻地种群处于容量、入侵种群初始局域于紧支撑区域，例如分段常数初值(1.4)。与多孔介质方程类比，关键问题是该系统是否存在具有传播界面的弱解，即驻地种群?处于承载能力且入侵种群v为零的点。该问题此前未有定论：一方面，有数值模拟显示形式(1.3)在初值(1.4)下可出现传播界面，并常演化为行波解；另一方面，Gallay和Mascia表明其模型的一参数族行波必须在波前具有V的代数衰减，而非有限位置界面，且波速趋于零时波前变陡；Mascia、Moschetta和Simeoni由偏微分方程（PDE，partial differential equation）模拟得到的波速估计对空间分辨率有轻微但不可忽略的依赖，引起数值结果是否收敛的疑问。为此，研究人员深入考虑形式(1.3)的理想化版本，包含退化交叉扩散与降解的耦合效应，即(1.5)：?u/?t=v，?v/?t=?/?x(u?v/?x)。其中u=1??为驻地种群亏缺，可视为时空相关扩散率，u=0表示驻地种群处于容量。研究人员特别关注扩散率初始为零且入侵种群紧支撑的解，但初始紧支撑v会带来问题；因此采用初始条件(1.6)：u(x,0)=0，v(x,0)=1（x0^*）、ε（x>x₀^*），并考察前方高度ε∈(0,1)趋于零的极限。

研究人员开展的工作与结论。研究人员建立并数值求解了上述理想化系统的自相似解，并进行小ε匹配渐近分析。主要发现是极限ε→0为奇异极限：对小ε，会出现表观界面（apparent interface），即v较小且v斜率随ε增大而很大的点；但该点的传播速度消失，剖面斜率在ε→0时变为无穷。因此，系统(1.5)不支持具有移动界面的解，这与多孔介质方程不同；该极限以对数方式缓慢趋近，具体为O([log(1/ε)]^?1/2)。这解释了为何此类模型的数值模拟可能看似收敛：数值方法在O(ε)尺度上正则化问题，而传播速度等观测量表面上可呈一阶量级。该结果对模拟具有同类退化扩散机制的更复杂模型具有重要启示，尤其强调需仔细检验数值收敛性。

关键技术方法（不超过250字）。研究人员先以相似变量η=(x?x₀^*)/t、u=tf(η)、v=g(η)将PDE化为常微分方程（ODE，ordinary differential equation）系统，并施加远场条件；在η=0奇异性处用展开式并分左右两套ODE，通过AUTO-07p软件和同伦延拓（homotopy continuation）计算高精度相似解。PDE数值模拟采用MATLAB中的有限体积法（finite volume method）离散空间、ode15s隐式时间推进，用于验证相似解。小ε分析采用匹配渐近分析，将负η外层、0<η<η^*压缩区及η^*附近内层进行匹配。本研究无实验样本队列。

研究结果。2 数值解与自相似性：通过PDE数值模拟展示ε=0.1时u与v的耦合演化，并将PDE解缩放到相似变量，结果显示精确塌缩到相似解，且与ODE相似解一致；因此后续可聚焦相似解计算。2.1 相似解系统的构造：通过相似变换得到ODE系统及远场条件；在η=0处系统奇异，展开显示f可不可微，g需跨奇异点可微，a₂^±可不同，从而提供额外自由度以满足四个边界条件；其后果是u在初始间断点x₀^*处随时间仍不可微。2.2 相似解的数值计算：通过将系统分为η=0左右两段、在奇异点附近用展开作边界条件，并应用同伦延拓先将误差项和远场值逐步调整到目标值，再扩大计算域；所得相似解与PDE塌缩剖面完全一致，确认其为原PDE在给定初值下的精确解。2.3 小ε数值结果：通过数值延拓小ε发现a₀、a₁在约ε=0.1处有转折；小ε时η=0右侧形成陡峭表观界面，η^*表示该点在PDE中的传播速度；数值表明η^*、1?a₀和a₁均趋于零但仅以[log(1/ε)]^?1/2量级缓慢趋零；还发现扩散率u为零时问题不一定满足比较原理（comparison principle），且两个不同ε可有相同的η<0剖面而在η>0不同。3 小ε渐近分析：3.1 负η：以δ=1?a₀展开，得到g₁、f₁及斜率展开系数a₁₁、a₁₂。3.2 表观界面外正η：令η=δξ并展开，得到g₀=1、g₁=?1、f₀线性等，以及g₂、g₃在ξ→ξ^*处的对数奇异项，为内层匹配准备。3.3 表观界面附近内层问题：定义内变量z、f?、?，得到自治系统，其解与(1.5)的行波型内层问题相关；积分得df?/dz=?ξ^*log(f?)，隐式解涉及对数积分（li，logarithmic integral）和指数积分（Ei，exponential integral）；匹配得到δ～√(π/2)[log(1/ε)]^?1/2，并给出η^*与a₁的两项渐近式。3.4 与数值结果比较：通过将两项渐近展开与数值计算比较，发现对数项收敛很慢；ε=10^?6时[log(1/ε)]^?1/2≈0.27仍不很小；使用数值δ和η^*拟合外层剖面，内层缩放后与渐近内层解几乎完全一致，内层误差为O(ε)。

讨论总结与研究结论翻译。讨论部分总结：本研究主要发现是，在耦合退化扩散与周围环境降解的系统中，当试图趋近可表现移动界面（即扩散率消失点）的情形时，极限具有奇异性；传播速度消失且解的梯度变为无穷，是系统不能表现此类移动界面的强证据。该理想化系统(1.5)的分析也适用于形式(1.3)的复杂模型：在早期时间，同样相似缩放给出主导阶行为，反应项和?衰减变化只引入O(t)修正；内层问题(3.8)预期适用于任何表观移动界面附近，因此速度对O(ε)正则化的对数依赖可推广到行波等情形。对数依赖解释了文献中主要来自PDE模拟的观察：数值格式中ε可代表空间网格尺寸，网格尺寸常难小于约10^?3，此时表观传播速度可为正，因此仔细检验数值收敛性至关重要。ε也可等价于离散过程如个体细胞模型中的无量纲细胞尺寸等正则化参数；不同内层细节会变，但与外层匹配仍给出相同对数依赖。该正则化效应与界面流体动力学中的接触线悖论有类比，其中接触线速度对小正则化参数呈对数依赖；也与Fisher–KPP方程反应项扰动对行波速度的对数依赖（如cut-off效应）类比。未来扩展包括含φ(v)?u/?x的交叉扩散系统，预期有相同对数依赖；若将退化性弱化为u^m且m∈(0,1)，则可能允许具有移动界面的弱解，但该结果在m→1^?时奇异。研究结论翻译：研究人员的主要发现是，一个耦合退化扩散和周围环境降解的系统在试图趋近能够表现出移动界面（即扩散率消失的点）的系统时，其极限具有奇异性质；传播速度在该极限中消失且解的梯度变为无穷，这有力表明该系统不能表现出这种移动界面。该结果还表明，表观界面传播速度仅以对数方式缓慢消失，说明紧支撑的移动前沿解不存在，并对具有该退化扩散机制的更复杂模型的数值模拟具有重要含义。