Solve the 2D wave equation with the finite difference method in real time. Click to add wave sources and watch interference, diffraction, and reflection patterns emerge.
The core physics is described by the 2D wave equation, which relates how the wave displacement changes in time to how it curves in space.
$$ \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) $$Here, \( u(x,y,t) \) is the wave height (displacement) at position (x,y) and time \(t\). \(c\) is the wave speed, a constant for the medium. The left side is the acceleration of the wave, and the right side sums its curvature in the x and y directions.
The simulator solves this equation numerically using the explicit Finite Difference Time Domain (FDTD) method. Stability requires satisfying the Courant窶擢riedrichs窶鏑ewy (CFL) condition.
$$ \frac{c \Delta t}{\Delta x}< \frac{1}{\sqrt{2}}$$Here, \( \Delta t \) is the time step and \( \Delta x \) is the grid spacing. This condition ensures the numerical calculation doesn't outrun the physical wave, preventing simulation blow-up. The factor \(1/\sqrt{2}\) is specific to 2D simulations.
Acoustic Room Design: Engineers use similar 2D wave simulations to model how sound waves reflect and interfere in concert halls or auditoriums. By testing different wall shapes and materials virtually, they can optimize the design for perfect acoustics before construction begins.
Ultrasound Imaging: Medical ultrasound machines rely on the physics of wave reflection and interference. Simulating how high-frequency sound waves propagate through different tissues helps in improving image resolution and developing new diagnostic techniques.
Antenna & RF Design: The propagation of radio waves from an antenna array follows the same wave principles. CAE simulations are crucial for designing antenna layouts that maximize signal strength and directionality for cell towers and satellite communications.
Seismic Analysis: Geophysicists use wave propagation models to understand how earthquake waves (seismic waves) travel through the Earth's layers. This analysis is vital for assessing earthquake risks and for oil and gas exploration.
When you start using this simulator, there are several points that are easy to misunderstand, especially for CAE beginners. The first one is the assumption that "setting the wave speed (c) to closely match real water waves will result in a more accurate simulation." While wave speed is indeed a physical parameter, it is strongly tied to the CFL condition ($c \Delta t / \Delta x \lt 1$), which determines numerical stability. For instance, if you set the wave speed c too high, you will need to make the time step $\Delta t$ extremely small to satisfy this condition, leading to heavier computations or instability. In practical work, finding a balance between capturing the essence of the phenomenon and keeping computational costs manageable is key.
The second point concerns the setting of the damping coefficient. While this tool mimics "water resistance," there are actually choices in how this damping is modeled. For example, in the analysis of vibrating structures, models like "viscous damping" (proportional to velocity) or "hysteretic damping" (proportional to displacement) are used depending on the phenomenon. If you set the simulator's damping too strong, interference fringes may become unclear, so the trick is to adjust it according to your observation goals.
Finally, there is the understanding that "boundary conditions are simply walls." This tool uses a "fixed end" (wave reflects at a wall), but in practice, diverse conditions like "free ends" (wall vibrates) or "absorbing ends" (wave dissipates at the boundary) exist. For example, in designing an anechoic chamber, "absorbing end" conditions are set on the walls to absorb sound as much as possible. When interpreting simulation results, always being aware of what boundary conditions were assumed is the first step to avoiding pitfalls in real-world applications.
Configure a 1 m × 1 m tank with wave speed = 1.2 m/s, damping = 0.05, frequency = 3 Hz. Place two coherent sources 0.4 m apart on opposite sides. At constructive interference zones, amplitude peaks reach 2× single-source amplitude; at destructive nodes, amplitude drops to near-zero. Adding a 0.1 m wide barrier in the center demonstrates diffraction, with wavelength λ = v/f = 1.2/3 = 0.4 m creating measurable diffraction spreading beyond the obstacle shadow zone.