2D波动方程求解器

此模拟器的基础是以下2D波动方程：

$$ \frac{\partial^2 u}{\partial t^2}= c^2 \left( \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2}\right) $$

其中，$u(x,y,t)$ 是位置$(x,y)$、时刻$t$的水面位移，$c$ 是介质中波的相速度（波速）。此微分方程通过计算机使用有限差分法（FDM, Finite Difference Method）求解，并应用时间中心差分格式（TCS）。

为了确保数值求解的稳定性，时间步长$\Delta t$与空间网格间隔$\Delta x, \Delta y$之间需要满足CFL条件（Courant-Friedrichs-Lewy条件）：

$$ c \frac{\Delta t}{\Delta x}\lt \frac{1}{\sqrt{2}} $$

违反此条件会导致数值计算不稳定，波会无限增大，导致计算结果完全错误。在此模拟器中，系统会自动调整确保此条件满足，但如果波速过大而网格步长太小，求解可能会变得不稳定。同时，阻尼系数（通常为0.98-0.995）用于表示在每个时间步中物理的能量损耗，来模拟真实水的阻尼。

振动与应力解析：在顾问公司工具中，水面波的反射、干涉、衍射等现象的模拟可用于优化产品设计。通过数值计算可视化波动的传播模式，工程师能在设计早期发现潜在问题。在地震波分析、水下建筑防波堤设计、电子设备散热器中的流体振动等应用中，都采用类似方法。

非线性问题与参数化：在振动与应力分析中，破坏波的反弹、非线性效应的衍射等现象通过模拟确定。通过参数化分析，可获得材料特性对波动特征影响的详细信息，为材料选择和工艺优化提供指导。

声波传播：超声波传感、核磁共振、超声检测等技术中，声波在介质中的传播模式（包括散射、衍射、多重反射等）在数值模拟中被精确计算。通过改变物理参数，可评估在不同条件下的超声波信号强度。

深海波浪与海岸工程：在海洋工程中，使用波动方程模拟波的传播速度，计算不同深度下的波高变化，以及在防波堤上的反射与衍射。这些计算结果直接指导海岸防护工程设计。

使用这个模拟器时，有几个关键的CAE误区值得注意。首先是"波速(c)的准确设置至关重要"——如果使用了错误的物理参数，即使数值计算完全正确，模拟结果也会完全错误。确保波速是物理上真实的参数，但CFL条件——即$c \Delta t / \Delta x < 1$——是硬性约束。例如，如果波速过大而你想使用较大的网格步长，就会出现时间步长Δt必须减小的情况，导致计算变得缓慢或不稳定，出现不同的失败形式。

第二个误区是"阻尼系数的设置"。在这个工具中，我们通过在每个时间步乘以阻尼系数来实现能量散逸，但在真实建模中选择阻尼模式有多种方式。例如，可以让阻尼与频率成正比（"比例阻尼"），或让阻尼与位移成正比（"刚度阻尼"）等。不同的阻尼模型会导致完全不同的物理行为。此模拟器采用了简单的统一阻尼以保持计算成本低廉，但这只是对现实的近似表示。

最后，"边界条件是看不见但至关重要的"。在此模拟器中应用的是固定边界条件（即波在边界处被完全反射），但在真实工程中，自由边界（波通过）或透射边界（波部分通过、部分反射）等多种边界条件都可能存在。理解模拟结果时，必须明确地注意采用了哪种边界条件——不同的边界条件会产生截然不同的结果，这对理解数值计算的有效性至关重要。