基于Yang变换混合方法的时间分数阶Richards方程数值分析

《Fractal and Fractional》：Numerical Analysis of the Caputo Fractional Richards Equation Using Yang Transform-Based Hybrid Methods

　　

Richards方程是研究非饱和土体（如多孔介质）入渗特性的著名方程。本研究的主要目标是揭示非饱和土体中水分运移问题的行为。研究人员借助Yang变换迭代法（YTIM）与同伦扰动变换法（HPTM），考察了非线性时间分数阶Richards方程（TFRE）的近似解。所用方法新颖且具有吸引力，成功将Yang变换方法、新迭代法及同伦扰动法相结合。两种技术均采用计算量更少的迭代过程，高效给出快速收敛的级数型解。为阐明所施方法的有效性，研究人员通过Caputo导数考察了TFRE的两个算例。绘制了不同分数阶下的图形以展示所得解的物理行为。并将导出的近似解与精确解进行数值比较，以证明YTIM与HPTM的效能。将所提方法获得的绝对误差与q-同伦分析变换法（q-HATM）所得者对照。非整数阶、空间与时间参数变化对所得解行为的影响亦以图形展示。使用Mathematica软件获取近似级数解并绘制若干分数阶的图形显示。结果表明，所提技术易于实施，可用于考察由时间分数阶非线性偏微分方程控制的复杂物理系统。

研究背景与问题提出

分数阶微积分（Fractional Calculus, FC）将积分与微分算子推广至非整数幂，其非局部性可刻画记忆与遗传效应，近三十年在连续介质力学、混沌理论、信号处理、复杂网络及最优控制等领域被广泛采用。经典Richards方程由连续性方程与Darcy–Buckingham定律耦合而成，是非饱和多孔介质（如土壤）水分入渗时空演化的基本模型；但因模型参数依赖于含水量、水力传导度与压力水头之间的强非线性关系，且涉及Brooks–Corey、Van Genuchten等土壤水滞留函数，精确解极难获得。传统数值或半解析法（如变分迭代、Adomian分解）常需线性化、离散化或计算复杂多项式。为此，研究人员开展本研究，旨在将Yang变换（Yang Transform, YT）分别与新迭代法、同伦扰动法结合，提出YTIM与HPTM两种低计算量、快收敛的混合算法，求解Caputo意义下的时间分数阶Richards方程（TFRE），并通过与精确解及q-HATM对照验证精度。

关键技术与方法

研究人员采用的核心方法为：①Yang变换（YT），其对Caputo导数具有直接代数映射Y{CDtαu}=?αuˉ(?)?∑k=0n?1?α?k?1u(k)(0)而免卷积；②YTIM，将YT与新迭代法融合，通过逆YT得迭代公式um+1=Y?1{?α1Y{N(um)}}；③HPTM，将YT与同伦扰动法及He多项式结合，以嵌入参数p展开u=∑pmum并逐项配平。算例参数取自Brooks–Corey型K(θ),D(θ)特殊化后的TFRE两情形（α∈(0,1]），使用Mathematica推演级数解并绘图，无外部样本队列。

研究结果

8.1 分数阶敏感性分析

绘制Case I（式49）与Case II（式64）在固定深度x下，α=0.5,0.7,0.9,1.0时的3D/2D剖面。结果显示：非饱和土含水量随α增大而上升，α→1时逼近精确解；物理上意味着孔隙介质湿度随分数阶增高逐步趋近饱和。

8.2 物理洞察

图形证实YTIM/HPTM近似解与精确解几乎重合；含水率随分数阶增加单调增大，表明记忆效应减弱、系统回归经典扩散。

8.3 精度与收敛验证

表1–表6列示α取不同值时YTIM、HPTM与精确解的数值对照，以及二者与q-HATM的绝对误差比较。YTIM与HPTM在α=1时误差趋零；HPTM在Case I误差量级10^-6，Case II近零；YTIM在Case II大参数下误差升至10^-2负向。

8.4 误差分析

表7–表8汇总α=1时CPU时间、最大绝对/相对误差：HPTM CPU略低于YTIM；HPTM最大绝对误差Case I约10^-6、Case II近0，YTIM随截断项n增多缓慢改善。

8.5 定量结果

Case I中HPTM绝对误差最小（如t=0.5,x=0.5处达10^-7级），YTIM在低x处劣于q-HATM但高x追平；Case II中HPTM全程10^-8~10^-9，YTIM在x=1处误差恶化为-10^-2。

8.6 数值稳定性

HPTM在所有α,t,x范围光滑增长无振荡；YTIM在Case II大α失稳，迭代展开无法抑制误差传播。

8.7 收敛速度

HPTM残差10^-6~10^-9仅需少数项；YTIM收敛慢，Case II中甚至偏离。

8.8 计算效率

HPTM迭代少、CPU与内存占用低；YTIM在稳定区可与q-HATM竞争，但在大参数区因无效迭代浪费资源。

8.9 方法局限

YTIM/HPTM依赖级数截断，长时程或强非线性下截断误差累积；复杂非线性项需反复逆YT与He多项式，初值偏离真实解时收敛不保；对Caputo–Fabrizio或Atangana–Baleanu–Caputo导数体系可能收敛变差。

8.10 未来方向

可换用Caputo–Fabrizio、Atangana–Baleanu等非奇异非局部算子；拓展至多维、变阶、随机外力及严格收敛性证明；将YTIM/HPTM与物理信息神经网络（PINN）耦合。

讨论与结论翻译

讨论指出，HPTM在精度、稳定性、收敛速度与CPU耗时上均优于YTIM与q-HATM，YTIM仅在部分参数窗内具竞争力。两种方法均免线性化/离散化、直接代入初边值，适用于强非线性时间分数阶偏微分方程组。

结论：研究人员针对Caputo型TFRE两算例给出YTIM与HPTM近似解，各分数阶α及t,x下结果与精确解高度吻合，α=1回归经典Richards解；2D/3D图显示分数阶调控入渗剖面；数值表证伪HPTM误差最低、YTIM条件稳定。所提框架易实现、可推广至科学工程中各类非线性分数阶系统。本文发表于《Fractal and Fractional》。