设置
统计
0.00
最大振幅
0.00
总能量
使用有限差分法实时求解二维波动方程。点击添加波源，观察干涉、衍射和反射图案的形成。
本模拟器的核心是求解描述波传播的二维波动方程。它告诉我们，空间中某一点水面高度的加速度，与周围水面高度的曲率（即凹凸程度）成正比。
$$\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$ 是波速，决定了波传播的快慢。方程右边括号里的两项分别代表了在 $x$ 和 $y$ 方向上的“弯曲”程度。
为了保证数值计算的稳定性，模拟采用了显式有限差分法（FTCS）进行离散求解，并必须满足CFL稳定性条件。这个条件限制了时间步长和空间步长的关系。
$$c \cdot \frac{\Delta t}{\Delta x}< \frac{1}{\sqrt{2}}$$这里，$\Delta t$ 是计算的时间步长，$\Delta x$ 是网格的空间步长。如果不满足这个条件，模拟结果会迅速发散（爆炸），无法看到正确的波动。这个条件是所有基于显式时间积分的CAE动力学分析（如碰撞、爆炸模拟）都必须严格遵守的数学基础。
光学设计与干涉测量：模拟器中双缝产生的干涉条纹，直接对应光波的干涉现象。工程师利用这一原理设计光学薄膜（如相机镜头镀膜）来增强或减弱特定波长的反射，也用于精密测量表面形貌和距离。
声学工程与噪声控制：声波同样遵循波动方程。通过模拟声波的干涉和衍射，可以设计音乐厅的声学结构，优化音响布局，或为汽车、飞机驾驶舱设计有效的主动降噪系统。
地震波分析与地质勘探：地震波在地层中的传播可以用更复杂的波动方程描述。通过分析地震波在不同岩层界面的反射、折射和衍射信号，可以推断地下结构，用于石油勘探或评估地震灾害。
CAE软件中的核心算法基础：本模拟器使用的有限差分法是数值求解偏微分方程的基石。在商业CAE软件（如ANSYS, Abaqus）中，虽然更多采用有限元法（FEM）或有限体积法，但其在时间域上的求解思路，以及对稳定性条件的处理，都与这里展示的原理一脉相承，是理解高级仿真工具的入门钥匙。
在开始使用本模拟器时，有几个CAE初学者容易误解的要点。首先是“只要将波速(c)设置得与现实水波完全一致，就能得到更精确的模拟结果”这种观念。虽然波速确实是物理参数，但它与决定数值计算稳定性的CFL条件（$c \Delta t / \Delta x < 1$）紧密相关。例如，若将波速c设置得过大，为满足该条件就需要将时间步长$\Delta t$设置得极小，可能导致计算负荷加重或结果不稳定。在实际工程中，关键在于平衡现象本质的捕捉与计算成本的控制。
第二点关于衰减系数的设置。本工具虽模拟了“水的阻力”，但实际上这种衰减的施加方式存在多种建模选择。例如在分析振动传播结构时，会根据现象选用与速度成正比的“粘性衰减”或与位移成正比的“滞回衰减”等不同模型。若将模拟器的衰减设置过强，干涉条纹会变得模糊不清，因此根据观察目的进行调整是关键技巧。
最后是“边界条件仅仅意味着存在墙壁”这种理解。本工具采用“固定端”（在墙壁处反射），但实际工程中还存在“自由端”（墙壁可振动）、“吸收端”（在边界处消除波）等多种形式。例如在设计隔音室时，墙壁会设置为尽可能吸收声音的“吸收端”条件。在解读模拟结果时，时刻注意所假设的边界条件类型，是避免实际工程失误的第一步。
这个二维波动方程的模拟，以其原始形式或扩展形式，支撑着众多尖端工程领域的核心。其一是无损检测领域。例如向飞机复合材料或桥梁钢架发射超声波，通过测量分析内部裂纹的反射波，可确定缺陷位置和尺寸。其物理原理与本模拟器中观察墙壁反射的现象相同，对测量数据的解读离不开基于模拟的“预研”。
另一领域是医学图像处理，特别是超声诊断设备。探头送入体内的超声波脉冲在脏器边界反射或在组织中传播的现象，正是由波动方程描述的。通过将不同组织间的声阻抗（声音传播难度）差异设为边界条件，可实现更贴近现实的模拟。这有助于开发新型图像重建算法，推动诊断技术进步。
此外还应用于岩土工程与资源勘探。通过向地下发送人工地震波，分析在地层边界反射返回的波（地震波勘探），可揭示地下结构或探测石油天然气储层。这种情况下介质非均匀，各地层波速c不同，模拟更为复杂，但其基础物理原理与本工具所体验的波的传播与反射完全一致。
熟悉本工具后，可进一步探索“为何通过这种计算能模拟波传播”这一数值解法的核心。首先应理解有限差分法的思想，学习如何用相邻点的值差（差分）来近似波动方程中的二阶微分 $\frac{\partial^2 u}{\partial t^2}$ 和 $\frac{\partial^2 u}{\partial x^2}$。例如时间二阶微分可离散化为：$$ \frac{\partial^2 u}{\partial t^2} \approx \frac{u^{n+1}_{i,j} - 2u^{n}_{i,j} + u^{n-1}_{i,j}}{(\Delta t)^2} $$ 将此式代入波动方程并求解未来值 $u^{n+1}$，即可计算出下一时间步的波形。
下一步建议尝试向三维扩展。现实中的声波或电磁波均在三维空间传播，其控制方程为：$$ \frac{\partial^2 u}{\partial t^2}= c^2 \left( \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}\right) $$ 计算量将显著增加，但核心解法与二维相同。这直接关联到汽车室内噪声（三维空间声压分布）模拟等应用。
最终可将视野拓展至其他数值解法。有限差分法虽直观，却不适用于复杂形状（如弯曲隔音壁或不规则地形）。此时需要学习有限元法（FEM）或边界元法（BEM）等方法。这些方法虽均通过将计算域分割为微小单元（网格）进行近似，但为波动问题提供了更灵活强大的解决方案。通过本工具扎实掌握波动现象的物理图像，将成为理解这些高级数值方法的最佳基础。