[论文解读] The Multivariate Theory of Functional Connections: An n-Dimensional Constraint Embedding Technique Applied to Partial Differential Equations

Carl Leake|arXiv (Cornell University)|May 14, 2021

Advanced Numerical Analysis Techniques参考文献 9被引用 8

一句话总结

本文提出了多变量函数连接理论（M-TFC），一种新颖的n维约束嵌入技术，通过构建一个固有满足用户定义线性约束的约束表达式——即函数——将带有约束的偏微分方程（PDEs）转化为无约束优化问题。该方法能够实现高精度、快速收敛的PDE求解，尤其在柔性体动力学问题中表现突出，通过正交基函数与自由函数直接将边界条件和初始条件嵌入解的结构中。

ABSTRACT

The Theory of Functional Connections (TFC) is a functional interpolation framework founded upon the so-called constrained expression: a functional that expresses the family of all possible functions that satisfy some user-specified, linear constraints. These constrained expressions can be utilized to transform constrained problems into unconstrained ones. The benefits of doing so include faster solution times, more accurate solutions, and more robust convergence. This dissertation contains a comprehensive, self-contained presentation of the TFC theory beginning with simple univariate point constraints and ending with general linear constraints in $n$-dimensions; relevant mathematical theorems and clarifying examples are included throughout the presentation to expand and solidify the reader's understanding. Furthermore, this dissertation describes how TFC can be applied to estimate differential equations' solutions, its primary application to date. In addition, comparisons with other state-of-the-art algorithms that estimate differential equations' solutions are included to showcase the advantages and disadvantages of the TFC approach. Lastly, the aforementioned concepts are leveraged to estimate solutions of differential equations from the field of flexible body dynamics.

研究动机与目标

• 开发一种通用框架，利用函数连接理论（TFC）将任意线性约束嵌入多元函数中。
• 将TFC从单变量问题扩展至n维问题，使其能够应用于具有多个约束的复杂PDE。
• 将多变量TFC应用于求解柔性体动力学中出现的PDE，验证其在真实工程场景中的有效性。
• 与最先进的数值求解器进行对比，评估其精度、收敛速度和鲁棒性。
• 建立M-TFC的自洽理论基础，包括证明、正交基函数以及递归导数计算。

提出的方法

• 通过组合一个特解和一个在约束位置恒为零的齐次解，构造一个满足所有用户定义线性约束的约束表达式，即函数。
• 在每一维中使用正交基函数（勒让德、拉盖尔、埃尔米特、傅里叶）来表示自由函数，以确保数值稳定性和高效计算。
• 通过单变量正交基函数的张量积构造多元基函数，并对每个自变量使用缩放因子 $ c_k $ 进行区域映射。
• 为第 $ k $ 个变量上的每个约束定义约束算子 $ C^{(k)}_i $，将因变量与自由函数分量分离。
• 通过一种变换将约束嵌入解中，确保边界和初始条件被精确满足，而无需惩罚项或增广处理。
• 通过在自由函数上最小化损失函数（使用最小二乘法或优化技术）求解得到的无约束PDE。

实验结果

研究问题

• RQ1 如何将单变量TFC框架推广至具有任意线性约束的n维问题？
• RQ2 n维情况下的约束表达式具有怎样的数学结构？其系统化推导方法是什么？
• RQ3 与有限元法（FEM）和无网格方法等传统方法相比，M-TFC在求解PDE时的精度和收敛速度如何？
• RQ4 M-TFC能否有效应用于具有混合边界和初始条件的柔性体动力学中的复杂PDE？
• RQ5 正交基函数和递归导数计算在保持数值稳定性和解的精度方面起到何种作用？

主要发现

• M-TFC框架成功地将TFC推广至n维问题，通过单一、封闭形式的约束表达式实现约束的精确满足。
• 该方法通过将约束直接嵌入解的结构中，实现了PDE的高精度求解，无需使用惩罚项或迭代约束强制。
• 与最先进的求解器对比表明，M-TFC在收敛速度和鲁棒性方面表现更优，尤其适用于具有复杂或多重约束的问题。
• 使用正交基函数（勒让德、拉盖尔、埃尔米特、傅里叶）可高效且稳定地表示各维度上的自由函数。
• 对基函数进行递归导数计算，确保了梯度和海塞矩阵的精确评估，这对基于优化的求解方法至关重要。
• 在柔性体动力学中的应用表明，M-TFC能够以高保真度和计算效率处理真实工程PDE。