一阶广义差分模型_向后差分公式法、广义 α 法和龙格-库塔法

解决方法

瞬态求解器提供三种不同的时间步进方法：隐式向后差分公式、广义 α 方法以及显式龙格-库塔系列方法。

“向后差分公式”(BDF) 求解器是一种使用向后差分公式的隐式求解器，其精度在一阶(也称为向后欧拉法)到五阶之间变化。向后差分公式法的使用已经有很长的历史，并以其稳定性而著称。然而，这些方法会产生严重的阻尼效应，低阶方法尤其如此。“向后欧拉法”会严重衰减任何高频。即使您期望得到具有陡峭梯度的解，但由于该方法的阻尼效应，您可能会得到在时间上非常平滑的解。在 COMSOL 中，向后差分公式的实现具有阶数可变的性质。也就是说，在可能的情况下会使用高阶，但在必要时将自动采用低阶以获得稳定性。向后差分公式法也是一种“微分代数方程”(DAE) 求解器。

“广义 α”的属性与二阶向后差分公式类似，但底层技术不同。该方法包含一个参数(文献中称为 α)，用于控制高频阻尼的程度。与向后差分公式(最高二阶)相比，“广义 α”具有更小的阻尼，因此更准确。出于相同的原因，其稳定性也较差。在 COMSOL 中，“广义 α”的实现表现为检测哪些变量在时间上是一阶，哪些是二阶，然后对这些变量应用正确的公式。即使“广义 α”常用于求解微分代数方程，但该方法并不是正式的微分代数方程求解器。“广义 α”于 1993 年被首次提出用于结构力学二阶系统，即具有二阶时间导数的系统(参考文献 1)。结构力学领域是该方法应用最为广泛的领域。人们于 2000 年提出了一个用于一阶系统(纳维-斯托克斯方程)的公式(参考文献 2)。

向后差分公式比“广义 α”更稳定、更通用，但同时会引入更多的阻尼。因此，对于结构力学、声学和瞬态电磁波等对过度耗散敏感的应用，“广义 α”是默认的瞬态求解器。在传递等其他应用中，通常需要较高的稳定性，