第四章 非局域性与相容性的崩溃
4.1 Conforming预设了局域纠缠
Conforming的"合法性"定义——共享面,不重叠——隐含假设:纠缠必须局域。两个element要交换信息必须几何上相邻。对局域PDE合理(微分算子只涉及无穷小邻域),对非局域物理错误。
4.2 非局域物理打破相容性
分数阶Laplacian (-Δ)^s u(x)每一个点依赖所有其他点。在conforming mesh上离散产生满矩阵N²。当前应对全是补丁:截断horizon、H-matrix压缩、特殊正交基——都是把非局域硬塞回局域框架。
杜强的非局域计算数学[3](渐近相容格式、非局域向量微积分)做到了顶尖,但底层mesh没有动。网格不动意味着D2不动,仍然是秩序态。这样很难进入真正的生命模拟——因为细胞的纠缠网络在随时动态改变,用一套固定的conforming网格已经将E的纠缠预设了太多,缺乏开放性。
4.3 网格是E的信息节点网络
网格不是空间的离散化——网格是E的离散化。每个节点是信息载体,承载E2三层:符号层(节点"是什么")、逻辑层(节点间的信息关系)、数学层(节点上的场)。
网格是信息网络,每个节点是三层信息处理单元,节点间的连接是信息通道。跟神经网络同构。
4.4 非局域网格:重叠的伞饼层级
非局域mesh不是conforming partition,是重叠的饼层级体系:大饼覆盖大范围(长程纠缠),小饼覆盖小范围(短程纠缠),空间上重叠——因为非局域纠缠多尺度同时发生。
传统mesh的"不重叠"公理就是局域纠缠的几何表达。SDE mesh的原则:必须重叠。不重叠不是"合法性",是不完整性——只有一个尺度的纠缠信息,丢失了所有其他尺度。
Conforming是一种最简单的特征纠缠系统——纠缠范围=零,纠缠通道数=有限,纠缠层数=一(不重叠),纠缠强度=均匀。四重退化。SDE mesh是四个参数全部打开后的完整纠缠网络。
4.5 固定网格预设了纠缠的封闭性
比conforming更深的问题:任何固定网格都预设了E。一旦固定,纠缠的拓扑、维度、分辨率全被锁定。能发现的物理不超出框架能表达的范围。
生命没有固定的感知框架。细胞的受体在移动、信号通路在切换、基因调控在跳转。细胞的"mesh"随时在重建。
SDE mesh的核心不是"生成一个好的mesh",而是永远不固定mesh。特征律:mesh永远在变(ΔSIO ≠ 0),变化的模式稳定了(d/dt ≈ 0)。
4.6 E2的信息系统如何锁死D2
不是"我们选择了秩序算法"。是conforming的E2决定了只能产生秩序算法。因果方向:E → S → D。
符号层锁死了选择空间——预定义的离散标签排除了"同时是三角形和四边形"或"同时属于两层"的可能。逻辑层锁死了决策规则——Boolean判断(refine或不refine)排除了粘性逻辑。数学层锁死了解空间——局部基函数继承了conforming的局域性。
三层加在一起:封闭系统。D2只能是秩序态。要打开D2三态,必须从E开始改。
4.7 SDE算法是必然的
非局域物理 → 非局域纠缠 → nonconforming mesh → 巨大配置空间 → 秩序态不可达 → 必须三态轮换。每步必然。
纯秩序态 → 困在局部最优。纯混沌态 → 永远不收敛。纯介生态 → 耗尽初始信息后固化。必须轮换:秩序维护 → 混沌注入 → 介生涌现 → 秩序接管 → 再次停滞 → 再次混沌 → 循环。
4.8 案例:peridynamics——非局域力学的数学框架
【A类核心证据】 Silling[4]在2000年提出的peridynamics(近场动力学)是非局域连续介质力学的代表性框架。它用积分方程替代微分方程来描述材料的变形和破坏。
传统连续介质力学: 力平衡方程 ρü = ∇·σ + b,其中σ是应力张量,∇·是散度算子。散度算子是局域的——它只依赖于无穷小邻域内的应力分布。这意味着:(1)方程在裂纹尖端(应力奇异)处失效;(2)不能自然描述裂纹的萌生和扩展——因为裂纹意味着连续性假设的破坏,而散度算子预设了连续性。
Peridynamics: 力平衡方程 ρü = ∫_{B_δ(x)} f(u(x')-u(x), x'-x) dV' + b,其中积分遍历点x的δ-邻域B_δ(x)内所有点x'。f是"bond力"——x和x'之间的非局域相互作用力。δ就是SDE伞模型中的半径r——纠缠范围。
SDE翻译: peridynamics中的每个物质点x就是一个伞中心,δ就是伞半径r。Bond f就是两个伞之间的纠缠耦合力。当δ→0时,peridynamics退化为经典连续介质力学——这存在显著的结构对应伞参数r→0的退化极限。
Peridynamics的优势在于它可以自然描述裂纹:当两个点x和x'之间的bond被拉伸超过临界值时,bond断裂——对应纠缠断开。裂纹就是一系列断裂的bond的集合。这不需要额外的断裂力学理论——断裂从非局域力的基本方程中自然涌现。
杜强(Qiang Du)的贡献: 杜强及合作者发展了peridynamics的渐近相容格式——证明了当δ→0时peridynamics解收敛到经典PDE解。这在SDE语言中是:证明了伞参数r→0时伞场方程退化为传统场方程。渐近相容性是连接非局域数学和局域数学的桥梁——它保证了撑开参数后的新理论在退化极限下不会跟旧理论矛盾。
但杜强的工作中网格仍然是conforming的——非局域的是PDE,不是mesh。SDE数学的下一步是让mesh也非局域——伞状重叠网格上的peridynamics。
4.9 分数阶微积分的伞模型解读
分数阶微积分将整数阶导数d^n/dx^n推广为非整数阶导数d^α/dx^α(α可以是0.5、1.7等任何正实数)。
Riemann-Liouville分数阶导数: D^α f(x) = (1/Γ(n-α)) · d^n/dx^n ∫_a^x (x-t)^{n-α-1} f(t) dt,其中n = ⌈α⌉。
这个定义的关键特征:它是非局域的——x处的分数阶导数依赖于从a到x的整个历史f(t)。这不是"在x点看无穷小邻域"——这是"在x点看从起点到x的整条路径"。
伞模型翻译: 分数阶导数中的积分核(x-t)^{n-α-1}就是伞的权重函数——距离x越近的点权重越大(如果n-α-1>0),距离越远权重越小但不为零。伞半径r = x-a(覆盖从起点到当前点的整个区间)。α就是伞的"非局域度"——α=1时退化为普通导数(纯局域),α→0时变成纯积分(完全非局域)。
物理意义: 分数阶导数描述的物理系统具有"记忆"——当前状态不仅取决于当前时刻的输入,还取决于过去所有时刻的输入的加权积分。这在生命系统中极其常见——免疫记忆、神经可塑性、表观遗传——都是"过去影响现在"的非局域时间纠缠。
分数阶微积分在SDE框架中的位置:它是伞微积分的一维时间特例——伞只在时间方向上展开,不在空间方向上展开。完整的伞微积分在时间和空间方向上同时展开——这比分数阶微积分更一般。
4.10 案例:地震波的非局域传播——伞模型的地球物理验证
地震波的传播是伞模型在地球物理中的天然验证。
体波(P波和S波): 在均匀介质中,体波是局域传播的——扰动从震源出发,以有限速度v向外球面扩散,影响范围随时间线性增长。这是r=0(局域)、ω=1(不重叠)、v有限的伞模型。
面波(Rayleigh波和Love波): 面波沿地表传播,振幅随深度指数衰减。面波的一个关键特征是频散——不同频率的分量以不同速度传播。低频分量穿透更深(感受更深层的地质结构),高频分量只在浅层传播。
从伞模型看频散:不同频率分量的面波就是不同半径r的"伞"——低频伞的r大(深穿透),高频伞的r小(浅穿透)。不同伞在地表重叠(ω>1——同一个地表点同时被多个频率伞覆盖),不同伞的传播速度v不同(频散——大伞慢,小伞快或反之,取决于地质结构)。
面波层析成像: 地震学家利用面波的频散特性来反演地球内部的速度结构——本质上就是从伞的表面行为(频散曲线)反推伞的内部参数(不同深度的波速)。这跟伞场的"逆问题"完全对应——给定伞场的可观测量,反推伞的参数分布。
地震波的非局域效应: 在不均匀介质(如含有断层、低速区、高速体的真实地壳)中,地震波出现散射、多次反射、波导效应——远处的地质结构可以通过波的多次反射"非局域地"影响近处的地震记录。这是r>0的伞效应——纠缠通过波的传播产生了非零的有效纠缠半径。