[论文解读] Treating Smoothness and Balance during Data Exchange in Explicit Simulator Coupling or Cosimulation

Dirk Scharff, Thilo Moshagen+1 more|arXiv (Cornell University)|Mar 16, 2017

Numerical methods for differential equations参考文献 3被引用 5

一句话总结

本文提出了一种新方法，通过将交换信号的外推与专门设计的光滑帽函数结合使用，以增强显式共仿真中的稳定性和精度。该方法确保了交换信号的连续导数，减少了非物理动力学，同时保持了系统范围内的守恒定律，从而提高了耦合物理系统仿真结果的可靠性。

ABSTRACT

Cosimulation methods allow combination of simulation tools of physical systems running in parallel to act as a single simulation environment for a big system. As data is passed across subsystem boundaries instead of solving the system as one single equation system, it is not ensured that systemwide balances are fulfilled. If the exchanged data is a flow of a conserved quantity, approximation errors can accumulate and make simulation results inaccurate. The problem of approximation errors is typically addressed with extrapolation of exchanged data. Nevertheless balance errors occur as extrapolation is approximation. This problem can be handled with balance correction methods which compensate these errors by adding corrections for the balances to the signal in next coupling time step. This work aims at combining extrapolation of exchanged data and balance correction in a way that the exchanged signal not only remains smooth, meaning the existence of continuous derivatives, but even in a way reducing the derivatives, in order to avoid unphysical dynamics caused by the coupling. To this end, suitable switch and hat functions are constructed and applied to the problem.

研究动机与目标

• 解决显式共仿真中因数据交换导致的能量或质量等守恒量累积近似误差与平衡误差的问题。
• 开发一种方法，确保交换信号的光滑性（连续导数），同时纠正平衡违规问题。
• 构建合适的切换函数与帽函数，实现校正信号的平滑重注入，而不会引入高频动力学。
• 在不损害方法显式并行结构的前提下，提高共仿真的稳定性和物理一致性。
• 通过最小化由耦合伪影引起的非物理瞬态，实现复杂耦合系统的可靠仿真。

提出的方法

• 作者设计了光滑、紧支撑的帽函数，其来源于两个对称帽函数的积分，以确保C1连续性及单位积分。
• 这些帽函数用于在多个时间步内分布平衡校正项，确保校正信号以平滑方式重新注入子系统。
• 该方法结合了信号外推与平衡校正，其中校正通过加权、光滑的贡献方式施加，以最小化导数突变。
• 通过将基函数ψ(t)进行变换，构造出对称、C1连续且多项式次数最小的帽函数。
• 该方法确保在时间区间上贡献之和保持为单位一，从而防止引入人工转折点或振荡。
• 该方法通过理论分析与数值示例得到验证，显示出更高的平滑性与更少的非物理动力学。

实验结果

研究问题

• RQ1 如何结合平衡校正与信号外推，以在不引入高频动力学的前提下保持共仿真中的光滑性？
• RQ2 哪类光滑、紧支撑的函数可用于在时间上分布校正，同时保持守恒定律？
• RQ3 能否构造出多项式次数最小的C1连续帽函数，使其满足单位积分与对称性要求？
• RQ4 对校正重注入过程进行平滑处理，在多大程度上可减少显式共仿真中的非物理瞬态？
• RQ5 所提出的方法是否在不损害共仿真显式并行结构的前提下，提升了仿真稳定性和精度？

主要发现

• 所提出的帽函数实现了C1连续性与单位积分，确保了校正信号的平滑且守恒的重注入。
• 该方法通过最小化校正信号的导数，减少了因导数突变引发的不稳定性，从而降低了非物理动力学。
• 通过“两帽函数积分”构造方法，作者推导出满足所有要求特性的最小次数多项式帽函数。
• 仅靠平衡校正无法完全消除误差，但当与平滑处理结合时，可显著延长仿真产生有意义结果的时间。
• 该方法将类似欧拉前向的共仿真转化为更具鲁棒性的方法，适用于更广泛的耦合系统类别。
• 该方法不保证对非线性 DAE 的收敛性，但显著提升了实际稳定性与物理一致性。