[论文解读] On some computational aspects of Hermite wavelets on a class of SBVPs arising in exothermic reactions

Amit K. Verma, Diksha Tiwari|arXiv (Cornell University)|Nov 3, 2019

Fractional Differential Equations Solutions参考文献 27被引用 5

一句话总结

本文提出了一类新型奇异非线性两点边值问题（SBVPs），用于建模放热反应，并引入了四种基于埃尔米特小波与拟线性化及牛顿-拉夫森迭代相结合的计算稳定数值方法。这些方法在少量迭代内即可实现高精度，并展现出优异的数值稳定性，相较于基于哈尔小波的方法，在求解兰恩-艾米迪亚类型方程及热反应动力学中的实际SBVPs时表现更优。

ABSTRACT

We propose a new class of SBVPs which deals with exothermic reactions. We also propose four computationally stable methods to solve singular nonlinear BVPs by using Hermite wavelet collocation which are coupled with Newton's quasilinearization and Newton-Raphson method. We compare the results obtained with Hermite Wavelets with Haar wavelet collocation. The efficiency of these methods are verified by applying these four methods on Lane-Emden equations. Convergence analysis is also presented.

研究动机与目标

• 通过广义阿伦尼乌斯定律的推广形式，建立放热反应中出现的新一类奇异非线性两点边值问题（SBVPs）的数学模型。
• 针对具有奇异性与强非线性的奇异非线性SBVPs，解决基于小波的数值格式在计算上的挑战。
• 开发并比较四种新型数值方法——埃尔米特小波牛顿法（HeWNA）、埃尔米特小波拟线性化法（HeWQA）、哈尔小波牛顿法（HWNA）与哈尔小波拟线性化法（HWQA），以实现高精度求解。
• 对HeWNA方法建立收敛性分析，并在基准问题与实际SBVPs上验证所提方法的有效性。
• 证明所提方法在初始向量发生微小扰动时仍保持数值稳定性和鲁棒性。

提出的方法

• 通过多分辨率分析（MRA）使用埃尔米特小波构造谱配置的基函数，实现对奇异非线性SBVPs解的高效逼近。
• 应用拟线性化方法对微分方程中的非线性项进行线性化，将原问题转化为一系列线性子问题的序列。
• 利用牛顿-拉夫森迭代求解由小波配置产生的非线性代数系统，确保快速收敛。
• 实现埃尔米特小波的积分与微分运算矩阵，将微分方程转化为代数方程。
• 通过利用埃尔米特小波的紧支集与正交性，有效处理t=0处的奇异性，实现仅用少量网格点即可精确解析。
• 与基于哈尔小波的方法（HWNA与HWQA）进行对比，验证基于埃尔米特小波方法在精度与稳定性方面的优越性。

实验结果

研究问题

• RQ1 埃尔米特小波能否有效用于求解一类新型奇异非线性SBVPs（用于建模放热反应），并相较于现有小波方法实现更高的精度？
• RQ2 所提出的基于埃尔米特小波的方法（HeWNA、HeWQA、HWNA、HWQA）在精度与收敛速度方面相较于已有的哈尔小波方法表现如何？
• RQ3 所提方法在初始猜测向量发生微小扰动时，其数值稳定性保持程度如何？
• RQ4 HeWNA方法的收敛性分析是否证实：随着分辨率水平J的增加，误差范数||E_{k,M}||趋于零？
• RQ5 所提方法能否准确求解实际SBVPs，如兰恩-艾米迪亚方程与人类头部热分布问题（其精确解未知）？

主要发现

• 所提出的基于埃尔米特小波的方法（HeWNA与HeWQA）在精度上显著优于基于哈尔小波的方法（HWNA与HWQA），如表11与图9所示，对应示例5.4。
• 在示例5.4（非线性兰恩-艾米迪亚型SBVP）中，HeWNA与HeWQA方法计算得到的解与HWNA与HWQA方法结果在至少8位小数内一致，表明具有高度一致性和稳定性。
• 在示例5.5（人类头部热分布）中，HeWNA与HeWQA方法的解比HWNA与HWQA方法更精确，差异出现在第4至第6位小数，如表12所示。
• 收敛性分析表明，误差范数||E_{k,M}||随M→∞而趋于零，证实随着分辨率水平J的增加，解的精度不断提高。
• 初始向量的微小扰动（如[0.9,0.9,…,0.9]或[0.8,0.8,…,0.8]）不会显著改变解，表明所提方法具有良好的数值稳定性。
• 所提方法在无需精确解的前提下，成功处理了SBVPs中的奇异性与强非线性，已在五个实际应用案例中得到验证，包括查姆布、中村与达根-古德曼的问题。