[论文解读] Uniform exponential stability of Galerkin approximations for damped wave systems

Herbert Egger, Thomas Kugler|arXiv (Cornell University)|Nov 26, 2015

Advanced Numerical Methods in Computational Mathematics被引用 5

一句话总结

本文通过证明半离散和全离散格式在最小正则性假设下保持物理能量的指数衰减，建立了阻尼波动系统伽辽金逼近的统一指数稳定性。分析依赖于修正的能量估计和混合有限元空间的一般相容性条件，确保在空间和时间离散化参数下均保持统一稳定性。

ABSTRACT

We consider the numerical approximation of linear damped wave systems by Galerkin approximations in space and appropriate time-stepping schemes. Based on a dissipation estimate for a modified energy, we prove exponential decay of the physical energy on the continuous level provided that the damping is effective everywhere in the domain. The methods of proof allow us to analyze also a class of Galerkin approximations based on a mixed variational formulation of the problem. Uniform exponential stability can be guaranteed for these approximations under a general compatibility condition on the discretization spaces. As a particular example, we discuss the discretization by mixed finite element methods for which we obtain convergence and uniform error estimates under minimal regularity assumptions. We also prove unconditional and uniform exponential stability for the time discretization by certain one-step methods. The validity of the theoretical results as well as the necessity of some of the conditions required for our analysis are demonstrated in numerical tests.

研究动机与目标

• 在阻尼和解的正则性假设最小的条件下，建立阻尼波动系统伽辽金逼近的统一指数稳定性。
• 分析半离散和全离散格式，包括空间和时间离散化，同时保持能量衰减特性。
• 为应用于阻尼波动系统的混合有限元方法推导统一的误差估计和稳定性界。
• 证明稳定性与收敛性结果在网格尺寸和时间步长上保持一致。
• 通过数值实验验证理论结果，确认关键分析条件的必要性。

提出的方法

• 使用修正的能量泛函推导连续问题的耗散估计，证明所有能量 $ E^k(t) $（$ k \geq 0 $）的指数衰减。
• 在伽辽金逼近空间上应用一般相容性条件，以确保适定性及半离散能量的统一指数衰减。
• 采用混合有限元方法进行空间离散化，在最小正则性假设下推导收敛性和统一误差估计。
• 使用单步时间积分格式（如Crank-Nicolson、后向欧拉）并证明全离散系统的无条件、统一指数稳定性。
• 依赖离散能量估计和离散版本的杨不等式，以控制时间导数并限制离散能量演化中的增长。
• 建立离散能量范数之间的等价性，并通过归纳法证明全离散格式下修正能量的指数衰减。

实验结果

研究问题

• RQ1 在阻尼参数正则性假设最小的条件下，能否保证阻尼波动系统伽辽金逼近的统一指数稳定性？
• RQ2 速度与压力有限元空间之间的何种相容性条件可确保离散能量的统一指数衰减？
• RQ3 采用单步时间积分器的全离散格式是否能以时间步长和网格尺寸无关的方式保持指数稳定性？
• RQ4 理论误差估计与稳定性界如何依赖于阻尼参数的下界与上界？
• RQ5 通过数值实验验证，稳定性与收敛性的理论条件是否必要？

主要发现

• 当阻尼 $ a(x) $ 一致为正时，连续阻尼波动系统的所有能量 $ E^k(t) $（$ k \geq 0 $）表现出统一指数衰减，包括物理能量 $ E^0(t) $。
• 伽辽金逼近空间的一般相容性条件可确保适定性，并使所有 $ k \geq 0 $ 的半离散能量实现统一指数衰减。
• 混合有限元方法在最小正则性假设下，可获得统一误差估计，并实现全离散格式的无条件、统一指数稳定性。
• 某些单步方法（如Crank-Nicolson、后向欧拉）的时间离散化可独立于时间步长大小保持统一指数稳定性。
• 构造了修正的能量泛函 $ E_{h,\varepsilon}^{1,n} $，以证明离散能量衰减速率 $ \varepsilon_0 > 0 $，从而得到 $ E_{h,\varepsilon}^{1,n} \leq e^{-\varepsilon\tau/3} E_{h,\varepsilon}^{1,n-1} $，确保指数衰减。
• 数值实验确认了理论条件的必要性，并验证了分析所预测的统一稳定性和收敛速率。